MATLAB上机实验报告_2.docx
平顶山学院计算机语言类课程实验报告(4)课程名称MAT1AB语言及应用实验机房3305院系电气信息工程学院专业电气工程及其自动化班级电气二班学号姓名机器号A4实验日期2013/5/3任课教师王凯实验学时2实验成绩一.实验名称：MAT1AB求解线性方程组和矩阵的初等计算二.实验目的和要求1、掌握利用MAT1AB程序编辑器编写应用程序的方法；2、掌握MAT1AB求解线性方程组的方法；3、掌握MAT1AB进行矩阵的初等计算的方法三.实验内容1、给出一个信号t=0:0.001:3,其正弦信号频率响应特性U=Sin(300t)+2cos(200t),求其幅频特性曲线。2、已知某系统的闭环传递函数为OG)=试求该系统的特征55+154+20/+352+405+50根，并判别其稳定性。3、已知某系统的开环传递函数为G(S)=!，画出该系统的开环NyqUiSt图和BOde5(35+1)(45+1)图，并判断系统的稳定性。4、已知某系统的开环传递函数为G(S)=，计算该系统的相位裕度和幅值裕度。s(s+4)(s+5)2一5、已知某系统的闭环传递函数为¢(S)=-j,设=0.707,3“=1:1:5,回出该s2+2wns+;系统的单位阶跃响应曲线。四.实验设计方案(实验步骤或开发过程)1、 c1eara11c1osea11c1ct=0:0.OO1:3;u=sin(300*t)+2*cos(200*t);y=fft(u);yy=abs(y);p1ot(t,yy)2、 d=11020304154;r=roots(d)3、 c1eara11c1osea11num=1;den=12710;bode(num,den)figure(2)nyquist(num,den)4、 num=5;den=conv(1,0,conv(1,4z1,5);sys=tf(num,den);W=Iogspace(0,4,50);bode(syszw);grid;GmzPmzWgzWc=margin(sys)5、i=0;t=0:0.1:10;forwn=1:5i=i+1;zeta=0.707;num=wn2;den=1z2*zeta*wn,wn2;sys=tf(num,den);y(：,i)=ste(syszt);endp1ot(tzy(:,1),t,y(:z2)zt,y(:,3)zt,y(:,4)zt,y(:z5);gridx1abe1(,time(s),)y1abe1(,stepresponse,)五.实验中存在问题及解决办法这一次上机实验遇到很多问题,首先关于用MAT1AB求自动控制相关问题以前接触到一点，但是现在完全不会，相关函数都不清楚，在上网查阅资料，仿照网友给出的程序最终得除了正确结果，但是对于怎么用程序函数判断系统的稳定性还是没有掌握。最大的问题还是对于MAT1AB软件不熟悉。六.实验结果1、2、r二-7.8739-1.4201+0.9942i-1.4201-0.9942i0.3570+1.4679i由此看出系统不稳定0.3570-1.4679i3、NyquetDiagramA»1Axis508o8DeQramo,o°Requency"adS)5OPnUJesS244、Gm=36（幅值裕度）Pm=83.5835（相角裕度）七.附录(源程序清单)1、 c1eara11c1osea11c1ct=0:0.001:3;u=sin(300*t)+2*cos(200*t);y=fft(u);yy=abs(y);p1ot(tzyy)2、 d=11020304154);r=roots(d)3、 >c1eara11c1osea11num=1;den=12710;bode(num,den)figure(2)nyquist(num,den)4、 >num=5;den=conv(1z0,conv(1,4,1,5);sys=tf(num,den);W=Iogspace(0z4z50);bode(syszw);grid;Gm,PmzWgzWc=margin(sys)5、 i=0;t=0:0.1:10;forwn=1:5i=i+1;zeta=0.707;num=wn2;den=1,2*zeta*wn,w2;sys=tf(num,den);y(：/i)=step(syszt);endp1ot(t,y(:,1),t,y(:,2)rt,y(:,3)zt,y(:,4),t,y(:,5);gridx1abe1(,time(s)')y1abe1(,stepresponse')
