jacky
问题补充说明：clear all,close allfid1=fopen('A619','wt');fid2=fopen('B619','wt');M(1,1)='m';M(2,2)='m';M(3,3)='m';K(1,1)='k',K(2,2)='k';K(3,3)='k';D=inv(K+1000*M)*M;A=ones(3,1);for i=1:3 pp0=0; i B=D*A; pp=1.0/B(1); A=B/B(1);while abs((pp-pp0)/pp)>1.e-12 pp0=pp; B=D*A; pp=1.0/B(1); A=B*pp;endif(pp-1000)>0 f=sqrt(pp-1000)/2/pielse f=0endfprintf(fid1,'%20.5f',A);fprintf(fid2,'%20.5f',A); D=D-A*A'*M/(A'*M*A*pp);endfid1=fopen('A619','rt');A=fscanf(fid1,'%f',[3,3]);fid2=fopen('B619','rt');A=fscanf(fid2,'%f',[3,1]);t=1:3;h1=figure('numbertitle','off','name','1','pos',[50 200 420 420]);bar(t,f(t,1)),xlabel('频率阶级次'),ylabel('Hz'),title('固有频率'),hold on,gridh1=figure('numbertitle','off','name','2','pos',[50 200 420 420]);bar(t,A(t,1)),xlabel(-'自由度（质量块）'),ylabel('振型向量'),title('0阶主振型'),hold on,gridpause(0.1)h1=figure('numbertitle','off','name','3','pos',[50 200 420 420]);bar(t,A(t,2)),xlabel('自由度（质量块）'),ylabel('振型向量'),title('1阶主振型'),hold on,gridpause(0.1)h1=figure('numbertitle','off','name','4','pos',[50 200 420 420]);bar(t,A(t,3)),xlabel('自由度（质量块）'),ylabel('振型向量'),title('2阶主振型'),hold on,grid索引超出矩阵维度。出错 Untitled (line 35)bar(t,f(t,1)),xlabel('频率阶级次'),ylabel('Hz'),求解决
你查查a1矩阵的大小。
出现这种提示一般就是矩阵中缺少你要的值。也就是说900或者386比矩阵的宽或高还360问答大。
