q605.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第六章例6.5程序q605
% 求任意线性系统的冲击函数和卷积
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clf,clear
a=input('多项式分母系数向量a=(书上取 [1,2,8] ）  ');
b=input('多项式分子系数向量b=(书上取 1 ）  ');
t=input('输入时间序列t=[0:dt:tf](书上取 [0:0.1:5]）  ');
u=input('输入序列u=(书上取 3*t+cos(0.1*t)）  ');
tf=t(end);dt=tf/(length(t)-1);
[r,p,k] = residue( b, a);				% 用极点留数法求冲击函数
h = r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);	% 冲击函数h(t)
subplot(2,1,1),plot(t,h);grid
set(gcf,'color','w')                    % 设置图形背景色为白色
x=conv(u,h)*dt;						    % 求u和h的卷积
subplot(2,1,2),plot(t,x(1:length(t)));grid