MATLAB上机实验报告.docVIP

平顶山学院计算机语言类课程实验报告（4） 课程名称 MATLAB语言及应用 实验机房 3305 院系 电气信息工程学院 专业 电气工程及其自动化 班级 电气二班 学号 姓名 机器号 A4 实验日期 2013/5/3 任课教师 王凯 实验学时 2 实验成绩 一．实验名称：MATLAB求解线性方程组和矩阵的初等计算 二．实验目的和要求 1、掌握利用MATLAB程序编辑器编写应用程序的方法； 2、掌握MATLAB求解线性方程组的方法； 3、掌握MATLAB进行矩阵的初等计算的方法 三．实验内容 1、给出一个信号t=0:0.001:3，其正弦信号频率响应特性u＝sin(300t)+2cos(200t)，求其幅频特性曲线。 2、已知某系统的闭环传递函数为 试求该系统的特征根，并判别其稳定性。 3、已知某系统的开环传递函数为 ，画出该系统的开环Nyquist图和Bode图，并判断系统的稳定性。 4、已知某系统的开环传递函数为 ，计算该系统的相位裕度和幅值裕度。 5、已知某系统的闭环传递函数为，设ζ＝0.707,ωn＝1:1:5，画出该系统的单位阶跃响应曲线。 四．实验设计方案（实验步骤或开发过程） 1、clear all close all clc t=0:0.001:3; u=sin(300*t)+2*cos(200*t); y=fft(u); yy=abs(y); plot(t,yy) 2、d=[1 10 20 30 41 54]; r=roots(d) 3、clear all close all num=1; den=[12 7 1 0]; bode(num,den) figure(2) nyquist(num,den) 4、num=5; den=conv([1,0],conv([1,4],[1,5])); sys=tf(num,den); w=logspace(0,4,50); bode(sys,w); grid; [Gm,Pm,Wg,Wc]=margin(sys) 5、 i=0; t=0:0.1:10; for wn=1:5 i=i+1; zeta=[0.707]; num=[wn^2]; den=[1,2*zeta*wn,wn^2]; sys=tf(num,den); y(:,i)=step(sys,t); end plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4),t,y(:,5)); grid xlabel(time(s)) ylabel(step response) 五．实验中存在问题及解决办法 这一次上机实验遇到很多问题，首先关于用MATLAB求自动控制相关问题以前接触到一点， 但是现在完全不会，相关函数都不清楚，在上网查阅资料，仿照网友给出的程序最终得除了正确结果，但是对于怎么用程序函数判断系统的稳定性还是没有掌握。最大的问题还是对于MATLAB软件不熟悉。 五．实验中存在问题及解决办法 六．实验结果 五．实验中存在问题及解决办法 六．实验结果 五．实验中存在问题及解决办法 六．实验结果 五．实验中存在问题及解决办法 六．实验结果 六．实验结果 1、 5、 2、r = -7.8739 -1.4201 + 0.9942i -1.4201 - 0.9942i 0.3570 + 1.4679i 0.3570 - 1.4679i 由此看出系统不稳定 3、 4、Gm = 36（幅值裕度） Pm =83.5835（相角裕度） 七．附录（源程序清单） 1、clear all close all clc t=0:0.001:3; u=sin(300*t)+2*cos(200*t); y=fft(u); yy=abs(y); plot(t,yy) 2、d=[1 10 20 30 41 54]; r=roots(d) 3、clear all close all num=1; den=[12 7 1 0]; bode(num,den) figure(2) nyquist(num,den) 4、num=5; den=conv([1,0],conv([1,4],[1,5])); sys=tf(num,den); w=logspace(0,4,50); bode(sys,w); grid; [Gm,Pm,Wg,Wc]=margin(sys) 5、i=0; t=0:0.1:10; for wn=1:5 i=i+1; zeta=[