1、>> wn=0.4wn = 0.4000>> t=0.01:80,num=wn2t = Columns 1 through 10 0.0100 1.0100 2.0100 3.0100 4.0100 5.0100 6.0100 7.0100 8.0100 9.0100 Columns 11 through 20 10.0100 11.0100 12.0100 13.0100 14.0100 15.0100 16.0100 17.0100 18.0100 19.0100 Columns 21 through 30 20.0100 21.0100 22.0100 23.01

2、00 24.0100 25.0100 26.0100 27.0100 28.0100 29.0100 Columns 31 through 40 30.0100 31.0100 32.0100 33.0100 34.0100 35.0100 36.0100 37.0100 38.0100 39.0100 Columns 41 through 50 40.0100 41.0100 42.0100 43.0100 44.0100 45.0100 46.0100 47.0100 48.0100 49.0100 Columns 51 through 60 50.0100 51.0100 52.0100

3、 53.0100 54.0100 55.0100 56.0100 57.0100 58.0100 59.0100 Columns 61 through 70 60.0100 61.0100 62.0100 63.0100 64.0100 65.0100 66.0100 67.0100 68.0100 69.0100 Columns 71 through 80 70.0100 71.0100 72.0100 73.0100 74.0100 75.0100 76.0100 77.0100 78.0100 79.0100num = 0.1600>> zeta1=0;den1=1 2*ze

4、ta1*wn wn2den1 = 1.0000 0 0.1600>> zeta3=0.3;den3=1 2*zeta1*wn wn2den3 = 1.0000 0 0.1600>> zeta5=0.5;den5=1 2*zeta*wn wn2den5 = 1.0000 0 0.1600>> y1,x,t=step(num,den1,t)y1 = 0 0.0789 0.3033 0.6376 1.0292 1.4161 1.7374 1.9422 1.9983 1.8968 1.6536 1.3073 0.9125 0.5315 0.2244 0.0398 0

5、.0068 0.1306 0.3916 0.7487 1.1455 1.5193 1.8111 1.9748 1.9847 1.8391 1.5610 1.1943 0.7970 0.4317 0.1561 0.0138 0.0272 0.1941 0.4883 0.8633 1.2598 1.6154 1.8737 1.9942 1.9577 1.7699 1.4607 1.0787 0.6843 0.3397 0.0994 0.0012 0.0608 0.2686 0.5919 0.9796 1.3706 1.7030 1.9245 2.0000 1.9176 1.6903 1.3541

6、0.9620 0.5758 0.2567 0.0548 0.0023 0.1072 0.3531 0.7011 1.0963 1.4763 1.7811 1.9626 1.9921 1.8650 1.6013 1.2427 0.8457 0.4732 0.1838 0.0232 0.0169x = t = Columns 1 through 16 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Columns 17 through 32 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 Columns 33 throug

7、h 48 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 Columns 49 through 64 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 Columns 65 through 80 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79>> y3,x,t=step(num,den3,t);>> y5,x,t=step(num,den5,t);>> plot(t,y1,t,y3,t,y5)>> grid on>> wn=0.4;>> z=0:0.1:2;>> t=0:0.1:80;>> figure(1)hold onfor z1=znum=wn2;den=1,2*z1*wn,wn2;y,x,t=step(num,den);plot(t,y)endtitle('step response')hold off>> wn=0.4;z=0:0.1:2;t=0:0.1:80;figure(1