例6-8.m

MATLAB 7基础与提高书的配套源程序

%设定公共参数
wn=10;                         				 
x0=1;             
v0=1;		
tf=2;     					                    	%计算终点时间

for j=1:10
     zeta(j)=0.1*j;              				 	%设定不同的zeta
     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是为了求四象限相角
     t=0:tf/1000:tf                                	%设定自变量数组
     x(j,:)=a*exp(-zeta(j)*wn*t).*sin(wd(j)*t+phi);      	%求过渡过程
end

figure

subplot(1,2,1)
for j=1:2:10
      plot(t,x(j,:));
      text(0.3,x(j,151),['zeta=',num2str(zeta(j))]);  
      %用于显示zeta（阻尼系数）值，针对第2种情况，此句可以修改为
      %text(0.5,x(j,151),...),以适应图形的调整
      hold on
end

xlabel('t'),ylabel('x')
grid
subplot(1,2,2), mesh(x)                   			%画出三维图形