如何解决索引超出数组范围

Lx00591

求助，

相关代码：

q=(r.*cos(e)./(r+h)).^2;

' f9 ?. \  c. A9 D- v9 r

am=1-q;

; a1 c/ J4 _! L4 y9 F. m% c

obliquity_factor=sqrt(am);%MF=stec/vtec 此时MF=1/obliquity_factor

: x# ^- l+ m( R# g* q0 B/ \

MF=1./obliquity_factor;

$ B) T: }$ F* @% H

* s& g7 H# m7 G: a/ k" m  k4 u

Orbit_15degree=[];

for i=1:length(Orbit_15min)

    for j=1:length(Sat_15degree)

, z! J0 r9 i7 F% l5 }0 r5 A

       if Orbit_15min(i,:)==Sat_15degree(j,:)

9 |: u& c6 t& }; M8 g. {

           Orbit_15degree_temp=Sat_15degree(j,:);

4 {4 p* j$ X% d: u4 t+ c

           Orbit_15degree=[Orbit_15degree;Orbit_15degree_temp];

6 G. g8 e6 b5 c% e# m+ h# K! {

       end

    end

5 U" \7 h* ^( R/ \. P7 d- C

end

%15degree时的STEC

: {4 w; `: }- v3 ?! D; t$ v9 q4 S

KSsat_15deg=[];

for j=1:length(Sat_15degree)

& t: E8 a3 ~2 `2 |

    for i = 1:size(KESATzong,1)

       if Orbit_15degree(j,:)==Orbit(KESATzong(i,3)).Coord(KESATzong(i,2),:)

7 ]3 L" O, q) x" {2 ]

         KSsat_15deg_temp=KESATzong(i,:);

         KSsat_15deg= [KSsat_15deg;KSsat_15deg_temp];

; _' i4 A  F# m, F

       end

( e2 @, y0 s) F3 x9 d% ~; _7 k* @) Y

    end

end

VTEC=KSsat_15deg(:,4).*obliquity_factor;%每个可视卫星对应时刻IPP点处的VTEC

kssatzong=[KSsat_15deg, Orbit_15degree,e,B,CarA,LLipp,VTEC,MF];

. t% d% d! O) _: j9 L

报错：

位置 2 的索引超出数组范围。

( p+ u' W7 r4 E/ ^

7 q, D# r7 A6 T# x/ |

出错 IPPzuobiao (line 136)

VTEC=KSsat_15deg(:,4)*obliquity_factor;%每个可视卫星对应时刻IPP点处的VTEC

% q  s3 H: f8 G: h8 S9 J# y  p0 v+ K

ExxNEN

如何解决索引超出数组范围