matlab下标是用括号表示的
L = 4; % L划分层数
for i = 1:L-1;
b_latitude=pi/L; % L划分层数
latitude(i)= i*b_latitude; %第i层时的纬度
latitude(i)= pi-i*b_latitude; %计算天顶距
b_longitude(i)= acos((cos(b_latitude)-(cos(latitude(i)))^2)/(sin(latitude(i)))^2); %每层上的经度差
ni=floor(2*pi/b_longitude(i)); %每条纬线上的格网点数
for j=1:ni;
longitude(i,j)= (j - 0.5)*2*pi/ni; %经度
end
end