%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第七章例7－3－1程序c
% 阻尼系數對單自由度振動的影響的直接求解
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
wn=10;tf=2;x0=0;v0=1;zeta=0.3,t=0:tf/1000:tf;	%	設定自變量數組
p=roots([1,2*zeta*wn,wn^2])
C=inv([1,1;p(1),p(2)])*[x0;v0]
x=real(C(1)*exp(p(1)*t)+C(2)*exp(p(2)*t)); 	%  公式求過渡過程
plot(t,x),grid on