wn=sqrt(10); kosi=1;
G=tf([wn*wn],[1 2*kosi*wn wn*wn]); step(G);
title('Step Response of kosi=1');
Step Response of kosi=110.90.80.70.6Amplitude0.50.40.30.20.1000.511.5Time (sec)22.53
%kosi=2的阶跃响应曲线 wn=sqrt(10);kosi=2;
G=tf([wn*wn],[1 2*kosi*wn wn*wn]);step(G); title('Step Response of kosi=2');
6 Step Response of kosi=210.90.80.70.6Amplitude0.50.40.30.20.100123456789Time (sec)
当wn不变时，由和的响应曲线可归纳：
①平稳性，由曲线看出，阻尼系数ζ ↑，超调量↓，响应的振荡↓，平稳性好；反之， ζ ↓，振荡↑，平稳性差。
②快速性，ζ↑，ts↑，快速性差；反之， ζ ↓， ts ↓；但ζ过小，系统响应的起始速度较快，但振荡强烈，影响系统稳定。
第5）题：
%wn1=0.5w0的阶跃响应曲线
w0=sqrt(10);kosi=1sqrt(10);wn1=0.5*w0;
G=tf([wn1*wn1],[1 2*kosi*wn1 wn1*wn1]);step(G); title('Step Response of wn1=0.5w0');
7 Step Response of wn1=0.5w01.41.21Amplitude0.80.60.40.200246Time (sec)81012
图1-6 wn1=0.5w0的阶跃响应曲线
%wn2=2w0的阶跃响应曲线
w0=sqrt(10);kosi=1sqrt(10);wn2=2*w0; G=tf([wn2*wn2],[1 2*kosi*wn2 wn2*wn2]); step(G);
title('Step Response of wn2=2w0');
8 Step Response of wn2=2w01.41.21Amplitude0.80.60.40.2000.511.5Time (sec)22.53
图1-7 wn2=2w0的阶跃响应曲线
由图1-6和图1-7得：
当ζ一定时，ωn↑，ts↓，所以当ζ一定时，ωn越大，快速性越好。
2. 作出以下系统的阶跃响应，并与原系统响应曲线进行比较，作出相应的实验分析结果
（1），有系统零点的情况 （2），分子、分母多项式阶数相等 （3），分子多项式零次项为零 （4），原响应的微分，微分系数为110 %各系统阶跃响应曲线比较
G0=tf([10],[1 2 10]);G1=tf([2 10],[1 2 10]);G2=tf([1 0.5 10],[1 2 10]); G3=tf([1 0.5 0],[1 2 10]);G4=tf([1 0 ],[1 2 10]); step(G0,G1,G2,G3,G4); grid on;
title('实验1.2 Step Response 曲线比较');
9 实验1.2 Step Response 曲线比较1.61.41.210.8Amplitude0.60.40.20-0.2-0.40123Time (sec)4567
图1-8 各系统的阶跃响应曲线比较
3. 单位阶跃响应：
求该系统单位阶跃响应曲线,并在所得图形上加网格线和标题 %单位阶跃响应 G=tf([25],[1 4 25]); step(G); grid on;
title('实验1.3 Step Response of G(s)=25(s^2+4s+25)');
10