|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
M文件如下:
function[p,a]=eigfreg(m,k)
p=[];
h=inv(m)*k;
[a,d]=eig(h);
n=size(h);
for i=1:n
p(i)=d(i,i);
a(:,i)=a(:,i)/a(1,i);
end
for j=1:n-1
for i=j+1:n
if p(j)>p(i)
t=p(i);
p(i)=p(j);
p(j)=t;
q=a(:,i);
a(:,i)=a(:,j);
a(:,j)=q;
end
end
end
for i=1:n
subplot(4,1,i)
plot(a(:,i))
grid on
end
执行:m= [25 25 0 0; 25 125 0 0; 0 0 30 0;0 0 0 30];
k= [400000 0 -400000 0; 0 400000 0 -400000;-400000 0 500000 0;0 -400000 0 500000];
eigfreg(m,k)
得出如下图形:
我想要把图片改成和文献中的一样,如何改啊,求大神帮忙!
|
|