作业帮 > 数学 > 作业

关于matlab矩阵特征向量的问题

来源:学生作业帮 编辑:灵鹊做题网作业帮 分类:数学作业 时间:2024/05/02 17:13:07
关于matlab矩阵特征向量的问题
要求一个含参数的矩阵的特征值和特征向量,如国结果是以虚数表示时,我不想把所有的结果都显示出来,只想把存在实部大于0的那部分显示出来,该如何实现?
例如 矩阵for Xc=0.85:0.4575:10
A =[-4179/14317 38009/28634 -6558011107873165*2^(-61) 0 6474480396118621*2^(-59) 4744231324742093*2^(-58) -9950/42951 0;-1 -42/191 50/191 0 8538633403822592*2^(-51) 6256757235559658*2^(-50) 0 -100/573;-7875/28634 71625/57268 -6262900608018873*2^(-59) 0 6100327006393174*2^(-59) 8940134405921030*2^(-59) -3125/14317 0;0 0 0 0 1 0 0 0;0 -7684061700217884*2^(-56) -5083102811422187*2^(-66) -5083/15660 -6761890687512570*2^(-58) 0 0 0;0 -4881639197785480*2^(-57) 6458530630983485*2^(-68) 0 -6420807358500570*2^(-62) 0 0 0;Xc 0 0 0 21721/1000 12733/400 0 1;0 Xc 0 0 -26883/53125 -927/1250 -1 0]
s=eig(A)
if real(s)>0
disp(s)
end
end
这是我编写的,可不能达到目的,
关于matlab矩阵特征向量的问题
k=0;B=[];C=[];m=0;
for Xc=0.85:0.4575:10
k=k+1;
A =[-4179/14317 38009/28634 -6558011107873165*2^(-61) 0 6474480396118621*2^(-59) 4744231324742093*2^(-58) -9950/42951 0;-1 -42/191 50/191 0 8538633403822592*2^(-51) 6256757235559658*2^(-50) 0 -100/573;-7875/28634 71625/57268 -6262900608018873*2^(-59) 0 6100327006393174*2^(-59) 8940134405921030*2^(-59) -3125/14317 0;0 0 0 0 1 0 0 0;0 -7684061700217884*2^(-56) -5083102811422187*2^(-66) -5083/15660 -6761890687512570*2^(-58) 0 0 0;0 -4881639197785480*2^(-57) 6458530630983485*2^(-68) 0 -6420807358500570*2^(-62) 0 0 0;Xc 0 0 0 21721/1000 12733/400 0 1;0 Xc 0 0 -26883/53125 -927/1250 -1 0];
B(:,k)=eig(A);
end
[a,b]=size(A);
for i=1:k
for j=1:a
if real(B(i))>0
m=m+1;
C(m)=real(B(i));
end
end
end