ea703.m

matlab 在线性代数中的应用 利用matlab实现可以提高对所学知识的理解 更有助于你考研

% 《工程线性代数（MATLAB版）》第七章例题7.3程序ea703
%
a=input('输入分母系数向量a=[a1,a2,...]=(书上取 [1,2,9,3] ) '); 
n=length(a)-1;
disp('输入初始条件向量(长度应比a小一) Y0=[y0,Dy0,D2y0,...] ')
Y0=input('Y0=(书上分别取[1,0,0],[0,1,0],[0,0,1]) ');
p=roots(a);						% 求系统的极点
V=rot90(vander(p));c= V\Y0';	% 用范德蒙特矩阵求对应于各极点的系数	
dt=input('dt=(书上取0.1) '); 
tf=input('tf=(书上取 5) ')
t=0:dt:tf; y=zeros(1,length(t));
for k=1:n  y= y+c(k)*exp(p(k)*t);end		% 将各分量叠加
plot(t,y),grid on