ea621.m

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

% 《工程线性代数（MATLAB版）》第六章例题6.21程序ea621
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 人口迁移模型
% 
A=[0.94,0.02;0.06,0.98]
x0=[0.3;0.7], x1=A*x0,
x10=A^10*x0, 
x30=A^30*x0, 
x50=A^50*x0
pause,[p,lamda]=eig(A)
%pause,q=inv(p),q*x0
syms k
xk=p*diag([lamda(1,1)^k,lamda(2,2)^k])*inv(p)*x0