ea622.m

来自「matlab 在线性代数中的应用 利用matlab实现可以提高对所学知识的理解 」· M 代码 · 共 13 行

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