关于Matlab怎么计算Lambert投影m值？(请教)

肥鹅

rt.图里是想用格点书（i、j）和数值天气书上的Lambert公式循环计算出每个格点的m，将其储存一起再输出图片。循环一块矩阵方面也有很多问题搞不清…在此请教大佬，感谢。

以下是循环处的程序：

for i=1,172;
    for j=1,232;
        In=i-(m-1)/2;
        Jn=(l-j*d)/d;
        l=sqrt((In.^2+Jn.^2)*d^2);
        sinfai=(le^(2/k)-l.^(2/k))/(le^(2/k)+l.^(2/k));
        an=asin(sinfai)/pi*180;
        bn=b+atan(In/Jn)/k/pi*180;
        m(i,j)=k*l/(a*sqrt(1-sinfai^2));
    end
end

肥鹅

已经捣鼓出来解决了