clear
close all
el(1)=64.035821/180.*pi;%仰角
az(1)=180.563562/180.*pi;%方位角
el(2)=42.215059/180.*pi;
az(2)=55.293444/180.*pi;
el(3)=16.531392/180.*pi;
az(3)=305.909139/180.*pi;
el(4)=38.166797/180.*pi;
az(4)=303.774466/180.*pi;
el(5)=60.448430/180.*pi;
az(5)=39.094687/180.*pi;
el(6)=33.807673/180.*pi;
az(6)=192.801643/180.*pi;
el(7)=13.587918/180.*pi;
az(7)=66.028819/180.*pi;
el(8)=53.779274/180.*pi;
az(8)=246.708353/180.*pi;
el(9)=31.749929/180.*pi;
az(9)=310.181418/180.*pi;
el(10)=7.725467/180.*pi;
az(10)=306.914246/180.*pi;
el(11)=21.615636/180.*pi;
az(11)=117.716822/180.*pi;
r=6400; %单位：千米
lamda=116.5326272014/180.*pi; %用户位置经度
faim=40.0139164358/180.*pi; %用户位置纬度
hi=350;%电离层高度
for (i=1:1:11 )
pefai(i)=0.5*pi-el(i)-asin([r/(r+350)]*cos(el(i))); %观测站与电离层穿透点的地心夹角
faim(i)=asin(sin(faim)*cos(pefai(i))+cos(faim)*sin(pefai(i))*cos(az(i)));%穿透点纬度
lamda(i)=lamda+asin(sin(faim)*sin(az(i))/cos(faim(i)));%穿透点经度
%F=sqrt(1-[r*cos(el)/(r+hi)]).^(-1);
end
报错 ???  In an assignment  A(I) = B, the number of elements in B and I must be the same.


Error in ==> ipp at 31
faim(i)=asin(sin(faim)*cos(pefai(i))+cos(faim)*sin(pefai(i))*cos(az(i)));%穿透点纬度
这个循环应该怎么改才能正确？