作者rainfull (小米酒)
看板MATLAB
标题[讨论] 求最小特徵值的程式码
时间Wed May 15 21:26:10 2013
A矩阵是包含变数k的矩阵
想求k=1~10的最小特徵值
然後画出k与最小特徵值的关系图
程式码如下，但是一直无法顺利执行
请问是哪个环节出了问题呢？（错误讯息在最下方）
k=1:10;
A=[-2 0 0 1 -(0.5)*k 0 0 0 0;
0 -2 0 (0.5)*k 1 -(0.5)*k 0 0 0;
0 0 -2 0 (0.5)*k 1 0 0 0;
1 (0.5)*k 0 -2 0 0 1 -(0.5)*k 0;
-(0.5)*k 1 (0.5)*k 0 -2 0 (0.5)*k 1 -(0.5)*k;
0 -(0.5)*k 1 0 0 -2 0 (0.5)*k 1;
0 0 0 1 (0.5)*k 0 -2 0 0;
0 0 0 -(0.5)*k 1 (0.5)*k 0 -2 0;
0 0 0 0 -(0.5)*k 1 0 0 -2];
B=[22 -8 1 -8 2 0 1 0 0;
-8 21 -8 2 -8 2 0 1 0;
1 -8 22 0 2 -8 0 0 1;
-8 2 0 21 -8 -1 -8 2 0;
2 -8 2 -8 20 -8 2 -8 2;
0 2 -8 1 -8 21 0 2 -8;
1 0 0 -8 2 0 22 -8 1;
0 1 0 2 -8 2 -8 21 -8;
0 0 1 0 2 -8 1 -8 22];
C=inv(A);
D=C*B;
y=min(eig(D));
plot(k,y);
xlabel('k=1~10');
ylabel('minimum eigenvalue');
出现以下错误讯息
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
Error in ==> aa at 2
A=[-2 0 0 1 -(0.5)*k 0 0 0 0;
0 -2 0 (0.5)*k 1 -(0.5)*k 0 0 0;
0 0 -2 0 (0.5)*k 1 0 0 0;
1 (0.5)*k 0 -2 0 0 1 -(0.5)*k 0;
-(0.5)*k 1 (0.5)*k 0 -0 (0.5)*k 1 -(0.5)*k;
0 -(0.5)*k 1 0 0 -2 0 (0.5)*k 1;0 0
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.120.109.14