q603.m

来自「MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全」· M 代码 · 共 18 行

%《MATLAB及在电子信息课程中的应用》第六章例6.3程序q603
% 求任意线性系统的冲击函数
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
a=input('多项式分母系数向量a=(书上取 poly([0,-1+2i,-1-2i,-2,-5]) ）  ');
b=input('多项式分子系数向量b=(书上取 [8,3,1] ）  ');
[r,p] = residue(b, a), 	            % 求留数
disp('解析式h(t)=Σr(i)*exp(p(i)*t)')
disp('给出时间数组t=[0:dt:tf]')
dt=input('dt= ');					% 给定时间数组
tf=input('tf= ');	
t=0:dt:tf;h=zeros(1,length(t)); 	%  h的初始化
for i=1:length(a)-1 				% 根数为a的长度减1
   h = h+ r(i)*exp(p(i)*t); 	    % 叠加各根分量
end
plot(t,h),grid
set(gcf,'color','w')                % 设置图形背景色为白色