matlab仿真中一个运行错误

abc66

在仿真中，首先进行求出反解，然后运行如下程序：
0 a; }  S$ K' H/ Un=50;rou=300;lamt=2*pi;
6 F, |; U! G  o9 b1 ?( ~Zz=[300 400 500];
Rou=rou*(0:n)/n;
) `) b' C) x0 I. C5 e3 o% S. i; B# aLamt=lamt*(0:n)/n;
[Rou1 Lamt1]=meshgrid(Rou,Lamt);
8 B- O. n8 L) r$ p5 A; U6 [0 zXx=Rou1.*cos(Lamt1);
Yy=Rou1.*sin(Lamt1);
" z3 T/ e& K* L, Pfai1=subs(faii,b,bb);
for i=1:3
) I# [% P, I4 `4 e% ^" s$ Hz=Zz(i);
for j=0:n
for k=0:n
9 z: k, y0 W" z0 X' Cx=Xx(j+1,k+1);
y=Yy(j+1,k+1);
sita1=double(subs(sita,{rx,rz},{x,z}));
c21=cos(sita1);
$ r* Y$ H  s, p. C: c/ K6 q, Fs21=sin(sita1);
4 y7 R- G( C+ G- h. |, b1 Pfai2=double(subs(fai1,{rx,ry,rz,s2,c2},{x,y,z,s21,c21}));
c31=cos(fai2);
# j' K6 `4 R* z1 Ds31=sin(fai2);
gang2{i,j+1,k+1}=subs(gang1,{rx,ry,rz,s2,c2,s3,c3},{x,y,z,s21,c21,s31,c31});
* H4 D. \! E' W. W, O3 z$ u# rend
end
. U3 f- N) |& A. {0 Kend
! Z# y" [) q6 x) S+ d9 c
根据以上程序的计算结果，运行如下程序：
for i=1:3
' G: E& q7 W6 D1 |for j=0:n
4 U. Y8 {6 C' {% K/ Kfor k=0:n
5 M( A6 I* ^. W1 ^. e/ }7 Fgang11=eigs(gang2{i,j+1,k+1});
gang12(j+1,k+1,i)=gang11(6)/gang11(1);
end
6 J; `5 }0 y5 j+ D! ?$ A+ M. C8 nend
figure(i)
' M& K% k) d: ]0 kgang12=double(gang12);
4 N3 ~, d! z% V4 qmesh(Xx,Yy,gang12(:,:,i))
end
5 H# X, ~/ p" B: g9 i: n* Q
该程序应该是绘制出三个图形，但是绘制一个图形以后就出现错误，错误如下：
/ ~- |9 X: S/ @' A0 |' y+ p' r/ |4 @??? Attempted to access gang11(6); index out of bounds because numel(gang11)=0.
8 A% i, |0 Q) G$ M! B
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的长度是多少就要有你自己看情况决定了。

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的长度是多少就要有你自己看情况决定了