%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第七章例7－3－2程序
% 單自由度阻尼系統的強迫振動,卷積方法
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
m=1;c=4;K=100;dt=0.015;				% 輸入給定的參數
w0=sqrt(K/m);		 				% 求系統固有頻率
zeta=c/sqrt(m*K)/2;					% 求系統固有阻尼系數
a=[1,2*zeta*w0,w0^2];b=1; 			% 求分母分子系數
[r,p]=residue(b,a);					% 求極點、留數
t=0:dt:1.2;
h=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);	% 求出系統的脈沖響應
f=[1:10,10*ones(1,70)];				% 給出外加力的采樣值
x=conv(h,f)*dt;						% 把脈沖響應和外加力作卷積
plot(t(1:80),x(1:80))				% 繪圖
v1=diff(x)/dt;						% 求導得出速度，注意求導后列長少一
[t(1:80)',f(1:80)',x(1:80)',[0,v1(1:79)]']  % 列出結果
set(gcf,'color','w'),grid on