matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
5 M% V( W* s  c( t- R; m6 v" M* on=50;rou=300;lamt=2*pi;
# I) ^$ j- y! C& }# W5 ^  SZz=[300 400 500];
: j9 A. k4 L9 X9 l4 M( V: v4 }Rou=rou*(0:n)/n;
Lamt=lamt*(0:n)/n;
; ]) S7 s. G) p3 I[Rou1 Lamt1]=meshgrid(Rou,Lamt);
Xx=Rou1.*cos(Lamt1);
! x; q3 W" s. q1 a& ?, W) h9 t0 m4 mYy=Rou1.*sin(Lamt1);
) ~) r5 e" s5 ffai1=subs(faii,b,bb);
for i=1:3
( y( w$ l7 ]5 R% X% R" U- r. g( x$ xz=Zz(i);
for j=0:n
for k=0:n
x=Xx(j+1,k+1);
" k0 u) B' q8 h7 H. T; Z8 k. {1 hy=Yy(j+1,k+1);
sita1=double(subs(sita,{rx,rz},{x,z}));
c21=cos(sita1);
  m/ o3 }# L. f4 y% {s21=sin(sita1);
fai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
c31=cos(fai2);
s31=sin(fai2);
: v) T$ W# e& T% ~2 zgang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
end
end
) b) F7 Y9 x0 s/ x7 j# Send
2 A$ m) J/ M8 F/ X; [
根据以上程序的计算结果，运行如下程序：
' ~& O5 T% V/ Afor i=1:3
for j=0:n
! N( ]5 D# J5 D5 ?( a* V1 q, {) ufor k=0:n
gang11=eigs(gang2{i,j+1,k+1});
# X2 k' }) x, F2 [gang12(j+1,k+1,i)=gang11(6)/gang11(1);
' W8 ?) B. w) ]8 `. m0 n/ A; gend
end
figure(i)
gang12=double(gang12);
mesh(Xx,Yy,gang12(:,:,i))
end

该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
0 h+ W7 U. o. y( n% _; k
Error in ==> gangdutu at 5
gang12(j+1,k+1,i)=gang11(6)/gang11(1);
6 K7 c: F5 ]8 B5 {! u+ {' p

Jame33

你的gang11是一个值，而不是一个向量，所以调用gang11（6）和gang（11）有误，超出了范围。应把上面一句 gang11=eigs(gang2{i,j+1,k+1});改为 gang11（k）=eigs(gang2{i,j+1,k+1});此处gang11的长度是多少就要有你自己看情况决定了。
# }& g  @: I; D: Z. r- G( S9 c3 s' V6 ~

oewqe

再看看别人是怎么说的
& Q& K3 z2 v, y: f" p8 ]8 [

EDA_Mr_PENG

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