matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
( z! c4 X2 c5 P( ^5 on=50;rou=300;lamt=2*pi;
Zz=[300 400 500];
Rou=rou*(0:n)/n;
Lamt=lamt*(0:n)/n;
; \  R6 W& y) J[Rou1 Lamt1]=meshgrid(Rou,Lamt);
' a% R- t) ]- N4 |! ]: SXx=Rou1.*cos(Lamt1);
3 o$ m8 ]2 l' w3 x9 wYy=Rou1.*sin(Lamt1);
; S0 Q! S& ]2 }" X" hfai1=subs(faii,b,bb);
$ H% T% @0 r- y- Kfor i=1:3
$ |5 t' l; x8 V% J$ M8 bz=Zz(i);
for j=0:n
# X; j" y* B/ A/ w8 A. |7 m, }for k=0:n
, p/ J! W* y) v: h3 N0 K, p# v! wx=Xx(j+1,k+1);
y=Yy(j+1,k+1);
% C' d& `- q4 X1 t3 hsita1=double(subs(sita,{rx,rz},{x,z}));
c21=cos(sita1);
s21=sin(sita1);
fai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
3 X/ w: k/ F& o) J. jc31=cos(fai2);
s31=sin(fai2);
& m6 M2 A% @! U" H9 e2 igang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
end
! P2 i% A6 S5 y9 k4 C2 \% M7 Mend
: a; `6 Q$ b( pend
8 \4 U& F, i3 O, @" O7 w0 o# v0 s
7 q8 H; r/ x1 G* R根据以上程序的计算结果，运行如下程序：
( t% q$ q7 a9 {2 Vfor i=1:3
. k7 K) e* S7 ?" v* {3 q9 u' z& L9 ifor j=0:n
for k=0:n
gang11=eigs(gang2{i,j+1,k+1});
: z; Y- R9 p. v$ \6 Igang12(j+1,k+1,i)=gang11(6)/gang11(1);
3 C1 D. Q; A# l" z, send
end
1 p$ K9 W( u' I0 Cfigure(i)
* x- \  I9 d) w' b) Fgang12=double(gang12);
mesh(Xx,Yy,gang12(:,:,i))
end
8 a6 _9 u* ]- l7 a, F0 z8 I( L& b
该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
/ I, [# Y( S9 Q: ~: o2 L; c$ @??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
% ?; y. L* e7 O1 V& g3 _
Error in ==> gangdutu at 5
gang12(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的长度是多少就要有你自己看情况决定了。
( V7 H" Z- x4 }/ N5 O5 F; L# A9 n

oewqe

再看看别人是怎么说的
" v9 x: V9 A& v' I/ W

EDA_Mr_PENG

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