matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
/ g8 X* d( l$ c, vn=50;rou=300;lamt=2*pi;
2 }1 t, |0 _7 u$ Q7 \; cZz=[300 400 500];
Rou=rou*(0:n)/n;
Lamt=lamt*(0:n)/n;
[Rou1 Lamt1]=meshgrid(Rou,Lamt);
$ @9 n  e; l9 K. a1 i* v/ H9 T; n  uXx=Rou1.*cos(Lamt1);
" }# _; _# g* k/ t% v& {Yy=Rou1.*sin(Lamt1);
fai1=subs(faii,b,bb);
for i=1:3
4 `2 K6 b/ b+ a; {4 C$ lz=Zz(i);
for j=0:n
9 P- k/ \+ J% w& D7 n1 afor k=0:n
: Q$ g; s# l9 p% [( ox=Xx(j+1,k+1);
y=Yy(j+1,k+1);
sita1=double(subs(sita,{rx,rz},{x,z}));
( O4 j% |! x; G9 O+ J8 Ic21=cos(sita1);
8 @/ I: K  U; C7 Ds21=sin(sita1);
3 z) \0 e6 P7 V6 {& J! O. Ofai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
c31=cos(fai2);
s31=sin(fai2);
gang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
end
! w$ p) [# n; w7 L* xend
end

1 ^5 a2 U) |4 v; b- _1 U/ X; f根据以上程序的计算结果，运行如下程序：
3 {3 e! b0 d: |* J" g0 u$ G5 ufor i=1:3
1 f5 e& K+ a, I1 R1 hfor j=0:n
for k=0:n
gang11=eigs(gang2{i,j+1,k+1});
$ ^1 J/ O$ ~! t/ A2 d9 L, j% F9 Xgang12(j+1,k+1,i)=gang11(6)/gang11(1);
end
$ x7 w$ s8 `" `. l/ ?" Vend
figure(i)
# {. K& r  g! E9 k! ?* tgang12=double(gang12);
mesh(Xx,Yy,gang12(:,:,i))
end

该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
9 U# }6 @% F$ c. l/ [' q& G??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
. `( q9 {3 R% Y' l' Y* Y/ M
Error in ==> gangdutu at 5
  f8 o% a1 d) R- @4 r- k: cgang12(j+1,k+1,i)=gang11(6)/gang11(1);
4 Q! I) Y, v! T. v

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

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