njupt
[ -1/2*(-1/2*5^(1/2)-1/2*i*3^(1/2))^3-1/4*5^(1/2)-1/4*i*3^(1/2)] y =
[ 1/2*5^(1/2)+1/2*i*3^(1/2)] [ 1/2*5^(1/2)-1/2*i*3^(1/2)] [ -1/2*5^(1/2)+1/2*i*3^(1/2)] [ -1/2*5^(1/2)-1/2*i*3^(1/2)]
6. 在某激励条件下，二阶系统归一化响应可表示为y(t) 1
1
其中 e tsin( t )，
为阻尼系数，
2
， arctg(
2
)。请用不同的颜色或线型，在同一张
图上，绘制 0.2:0.2:2.0取值下系统在t [0,18]区间内的响应曲线，并要求用
0.2和 2.0对它们相应的两条曲线进行醒目的文字标志。
程序：
b=0.2:0.2:2.0; t=0:0.1:18;
color=['r','g.','bo','cx','m+','y*','kv','r--','g:','b-.']; for i=1:10
p(i)=sqrt(1-b(i).*b(i));
q(i)=atan(sqrt(1-b(i).*b(i))./b(i));
y=1-1./p(i).*exp(-b(i).*t).*sin(p(i).*t+q(i)); if i==1
plot(t,y,color(i))
text(t(30),y(30),' \leftarrow ξ=0.2','FontSize',18) hold on
else if i==10
plot(t,y,color(i))
text(t(50),y(50),' \leftarrow ξ=2.0','FontSize',18) hold on else
plot(t,y,color(i)) hold on end end end 结果：