%《MATLAB在電子信息課程中的應用》第六章例6.2程序q602
% 求任意線性系統的零輸入分量：
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

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