matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
n=50;rou=300;lamt=2*pi;
$ k( c, D2 g0 s) M8 v/ `1 UZz=[300 400 500];
Rou=rou*(0:n)/n;
Lamt=lamt*(0:n)/n;
7 `* T4 W6 n( X) u7 R$ q[Rou1 Lamt1]=meshgrid(Rou,Lamt);
Xx=Rou1.*cos(Lamt1);
+ \; b! C- i, N3 @' D9 `Yy=Rou1.*sin(Lamt1);
fai1=subs(faii,b,bb);
for i=1:3
z=Zz(i);
5 R- ?: [( I: Qfor j=0:n
for k=0:n
x=Xx(j+1,k+1);
2 T: T! t- x$ j6 v3 _y=Yy(j+1,k+1);
9 W: S; `4 t( i: `# Y4 w/ [7 Zsita1=double(subs(sita,{rx,rz},{x,z}));
c21=cos(sita1);
4 P) L8 y1 m; n) H8 ks21=sin(sita1);
1 }' G& [1 b5 ffai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
0 _3 e. A! q* f& Y4 ic31=cos(fai2);
  }$ o% Z# s% g, R4 Ds31=sin(fai2);
6 L: ]  q2 q- T& R+ C( Fgang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
& l' N# m$ G: }/ Y; R0 Tend
end
4 P* t1 ?3 M  o& J, ]9 z/ \) Xend

根据以上程序的计算结果，运行如下程序：
( ?1 p$ s7 c5 K( o7 d; zfor i=1:3
4 T, p" F; N! T. x" I4 ofor j=0:n
4 S5 u9 ^. u" C  o5 zfor k=0:n
& Z  j; ^/ }% Zgang11=eigs(gang2{i,j+1,k+1});
$ Z- h+ H7 e& {) Sgang12(j+1,k+1,i)=gang11(6)/gang11(1);
end
end
1 H9 r" Y7 x; J; S3 Q  _figure(i)
gang12=double(gang12);
mesh(Xx,Yy,gang12(:,:,i))
end

该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
( V5 V0 I" s3 L??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
* E+ h9 i+ J/ b. |0 T  S
4 f' a- V, x2 j* w  uError in ==> gangdutu at 5
+ \3 u6 S7 @* f7 a, Fgang12(j+1,k+1,i)=gang11(6)/gang11(1);

Jame33

你的gang11是一个值，而不是一个向量，所以调用gang11（6）和gang（11）有误，超出了范围。应把上面一句 gang11=eigs(gang2{i,j+1,k+1});改为 gang11（k）=eigs(gang2{i,j+1,k+1});此处gang11的长度是多少就要有你自己看情况决定了。

oewqe

再看看别人是怎么说的

EDA_Mr_PENG

% b  ^2 Q0 \4 b4 ?你的gang11是一个值，而不是一个向量，所以调用gang11（6）和gang（11）有误，超出了范围。应把上面一句 gang11=eigs(gang2{i,j+1,k+1});改为 gang11（k）=eigs(gang2{i,j+1,k+1});此处gang11的长度是多少就要有你自己看情况决定了