exn731b.m

MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第七章例7－3－1程序b
% 阻尼系数对单自由度振动的影响
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
wn=10;tf=2;x0=0;v0=1; 
for j=1:10
    zeta(j)=0.1*j;					% 对不同的ζ
    wd(j)=wn*sqrt(1-zeta(j)^2); 		% 求wd
    a=sqrt((wn*x0*zeta(j)+v0)^2+(x0*wd(j))^2)/wd(j); 	% 求振幅A 
    phi=atan2(wd(j)*x0,v0+zeta(j)*wn*x0); 	% 用atan2是为了求四象限相角phi
    t=0:tf/1000:tf;				%	设定自变量数组
    x(j,:)=a*exp(-zeta(j)*wn*t).*sin(wd(j)*t+phi); 	%  公式求过渡过程
end
plot(t,x(1,:), t,x(2,:), t,x(3,:), t,x(4,:), t,x(5,:),...	% 绘图
    t,x(6,:), t,x(7,:), t,x(8,:), t,x(9,:), t,x(10,:))
grid on