zy63.m

电路与系统专业的matlab程序。有很多经典的实例。

%《MATLAB在电子信息课程中的应用》第六章例6.2程序q602
% 求任意线性系统的零输入分量：
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著  2001年10月

a=input('输入分母系数向量a=[a1,a2,...] '); 
n=length(a)-1;
disp('输入初始条件向量(长度应比a小1) Y0=[y0,Dy0,D2y0,...] ')
Y0=input('Y0= ');
p=roots(a);						% 求系统的极点
V=rot90(vander(p));c= V\Y0';	% 用范德蒙特矩阵求对应于各极点的系数	
dt=input('dt= '); 
tf=input('tf= ')
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