matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
: n+ ^, z* o' L& v' o1 Gn=50;rou=300;lamt=2*pi;
* q0 |* |% p0 q% K) K' _Zz=[300 400 500];
Rou=rou*(0:n)/n;
Lamt=lamt*(0:n)/n;
[Rou1 Lamt1]=meshgrid(Rou,Lamt);
" n& I" D3 M. L: r- M3 DXx=Rou1.*cos(Lamt1);
Yy=Rou1.*sin(Lamt1);
fai1=subs(faii,b,bb);
7 L8 a% q( o: F2 Gfor i=1:3
z=Zz(i);
for j=0:n
3 C! r9 d0 l  t! h# h" dfor k=0:n
1 @# [+ E1 g0 D7 ~x=Xx(j+1,k+1);
y=Yy(j+1,k+1);
sita1=double(subs(sita,{rx,rz},{x,z}));
c21=cos(sita1);
0 j& v2 Y; F0 f7 Q8 u: ]8 x: Os21=sin(sita1);
' x9 H& k: }$ Hfai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
c31=cos(fai2);
! F: G0 f/ }6 g4 Ws31=sin(fai2);
* X8 {) }& d* n3 C. x: S. n8 i' ~gang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
end
! G( C1 x/ n' }" e: q! i- w$ nend
+ K- y  @; [2 ?) M. g  Lend

# Z' H' w+ n- c1 x1 U根据以上程序的计算结果，运行如下程序：
7 d7 @# D+ Q: G6 `2 v& ~# H$ T4 x9 {for i=1:3
/ a" @  S5 O5 ifor j=0:n
for k=0:n
; p3 K/ t; a8 Dgang11=eigs(gang2{i,j+1,k+1});
gang12(j+1,k+1,i)=gang11(6)/gang11(1);
end
end
1 G  C, q9 b+ ^! a8 t- hfigure(i)
gang12=double(gang12);
mesh(Xx,Yy,gang12(:,:,i))
+ p7 r: Z# i! H8 X/ _0 iend

+ }1 Y$ V  Q) n9 E: `该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
* n5 ~$ F4 C4 t6 ?3 Y??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
+ D( P1 G/ q/ B& x' i8 c! h
7 m" a& N2 g  M' o+ PError 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的长度是多少就要有你自己看情况决定了。
# Q$ B& I9 T! l3 Y, |% Z

oewqe

再看看别人是怎么说的

EDA_Mr_PENG

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