ea622.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 14 行

% 《工程线性代数（MATLAB版）》第六章例题6.22程序ea622
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 物料混合问题
% 
clear
A=[-1/10,1/40;1/10,-1/10]
[p,lamda]=eig(A)
for t = 1:100;
    d11=exp(lamda(1,1)*t); d22=exp(lamda(2,2)*t);   % 求两个对角元素对指数
    Y(:,t)=p*[d11,0;0,d22]*inv(p)*[60;0];
end
plot([1:100],Y(1,:),[1:100],Y(2,:),':')
grid on
set(gcf,'color','w')