自动控制原理Matlab程序作业(精)
自动控制原理(第2版)(余成波_张莲_胡晓倩)习题全解及MATLAB实验第5章习题解答
第5章频率特性法频域分析法是一种图解分析法，可以根据系统的开环频率特性去判断闭环系统的性能，并能较方便地分析系统参量对系统性能的影响，从而指出改善系统性能的途径，已经发展成为一种实用的工程方法，其主要内容是：1）频率特性是线性定常系统在正弦函数作用下，稳态输出与输入的复数之比对频率的函数关系。
频率特性是传递函数的一种特殊形式，也是频域中的数学模型。
频率特性既可以根据系统的工作原理，应用机理分析法建立起来，也可以由系统的其它数学模型（传递函数、微分方程等）转换得到，或用实验法来确定。
2）在工程分析和设计中，通常把频率特性画成一些曲线。
频率特性图形因其采用的坐标不同而分为幅相特性(Nyquist图)、对数频率特性(Bode图)和对数幅相特性(Nichols图)等形式。
各种形式之间是互通的，每种形式有其特定的适用场合。
开环幅相特性在分析闭环系统的稳定性时比较直观，理论分析时经常采用；波德图可用渐近线近似地绘制，计算简单，绘图容易，在分析典型环节参数变化对系统性能的影响时最方便；由开环频率特性获取闭环频率指标时，则用对数幅相特性最直接。
3）开环对数频率特性曲线(波德图)是控制系统分析和设计的主要工具。
开环对数幅频特性L(ω)低频段的斜率表征了系统的型别(v)，其高度则表征了开环传递系数的大小，因而低频段表征系统稳态性能；L(ω)中频段的斜率、宽度以及幅值穿越频率，表征着系统的动态性能；高频段则表征了系统抗高频干扰的能力。
对于最小相位系统，幅频特性和相频特性之间存在着唯一的对应关系，根据对数幅频特性，可以唯一地确定相应的相频特性和传递函数。
4）奈奎斯特稳定性判据是利用系统的开环幅相频率特性G(jω)H(jω)曲线，又称奈氏曲线，是否包围GH平面中的(－l，j0)点来判断闭环系统的稳定性。
利用奈奎斯特稳定判据，可根据系统的开环频率特性来判断闭环系统的稳定性，并可定量地反映系统的相对稳定性，即稳定裕度。
稳定裕度通常用相角裕量和幅值裕量来表示。
自动控制原理实验用Matlab软件编制劳斯判据程序并解题(【学习辅导】例4.3.5)
上海电力学院实验报告自动控制原理实验课程题目：用Matlab软件编制劳斯判据程序并解题(《学习辅导》例4.3.5) 班级：姓名：学号：时间：2012年11月4日自动化工程学院自动化（电站自动化）专业实验报告目录一、问题描述 (3)二、理论方法分析 (3)三、实验设计与实现 (3)四、实验结果与分析 (5)五、结论与讨论 (6)六、实验心得体会 (6)七、参考文献 (7)八附录 (7)一、问题描述用MATLAB编制劳斯判据列出其劳斯矩阵并判断相对应系统的稳定性二、理论方法分析采用M文件实现Matlab编程。
1) M文件的建立与调用➢从Matlab操作桌面的“File”菜单中选择“New”菜单项，再选择“M-file”命令，屏幕将出现Matlab文本编辑器的窗口。
➢在Matlab命令窗口的“File”菜单中选择“Open”命令，则屏幕出现“Open”对话框，在文件名对话框中选中所需打开的M文件名。
2) M文件的调试在文件编辑器窗口菜单栏和工具栏的下面有三个区域，右侧的大区域是程序窗口，用于编写程序；最左面区域显示的是行号，每行都有数字，包括空行，行号是自动出现的，随着命令行的增加而增加；在行号和程序窗口之间的区域上有一些小横线，这些横线只有在可执行行上才有，而空行、注释行、函数定义行等非执行行的前面都没有。
在进行程序调试时，可以直接在这些程序上点击鼠标以设置或去掉断点。
三、实验设计与实现（1）程序%RouthMatrix**劳斯矩阵(带参数的特征多项式)并判断对应系统稳定性** clear;syms k z q %定义变量k z qp=input('请输入特征多项式的参数='); %提示输入参数n=length(p); %得到p的长度for i=0:ceil(n/2)-1 %将多项式进行劳斯矩阵排序a(1,i+1)=p(2*i+1);if 2*(i+1)>na(2,i+1)=0;breakenda(2,i+1)=p(2*(i+1));endfor k=3:n %计算从第三行开始劳斯矩阵内容for j=1:ceil((n-k+1)/2)if a(k-1,1)==0 %判断是否有共轭虚根disp('系统有共轭虚根')breaksign=1;breakenda(k,j)=(a(k-1,1)*a(k-2,j+1)-a(k-1,j+1)*a(k-2,1))/a(k-1,1);endenddisp('劳斯矩阵') %输出对应的劳斯矩阵disp(double(a))for i=3:k %用劳斯判据判断系统的稳定性if a(i-1,1)<=0 %判断第一列元素是否不大于0 q=1;breakendendif q==1disp('系统不稳定')elsedisp('系统稳定') %输出系统稳定性判定结果end（2）设计劳斯表的内容根据书上的原理进行设计，将多项式进行劳斯矩阵排序，可得其前两行的内容，余下几行则要根据书上原理写出。
自动控制原理中MATLAB常用编程语言学习及例题示例
系统数学模型建立num=[1 2];den=[4 3 2 1];model=tf(num,den)model=s + 2-----------------------4 s^3 + 3 s^2 + 2 s + 1num=7*[2 3];den=conv([1 0],[3 2 1]);model=tf(num,den)model =14 s + 21-----------------3 s^3 + 2 s^2 + s系统连接串联：函数调用:series(sys1,sys2)并联：函数调用:parallel(sys1,sys2)反馈：函数调用:feedback(sys1,sys2,sign)sys1前向通道sys2反馈通道sign缺省为负sign=1为正num1=[1];num2=[2];den1=[2 1];den2=[1 2];sys1=tf(num1,den1);sys2=tf(num2,den2);syss=series(sys1,sys2)syss =2---------------2 s^2 + 5 s + 2sysp=parallel(sys1,sys2)sysp =5 s + 4---------------2 s^2 + 5 s + 2sysf=feedback(sys1,sys2)sysf =s + 2---------------2 s^2 + 5 s + 4模型转换传递函数模型—>零极点增益模型函数调用:[z,p,k]=tf2zp(num,den) num=[1 1];den=[1 5 6];[z p k]=tf2zp(num,den);sys=zpk(z,p,k)sys =(s+1)-----------(s+3) (s+2)零极点增益模型—>传递函数模型函数调用:[num,den]=zp2tf(z,p,k) z=[-1];p=[-2 -3];k=1;[num,den]=zp2tf(z,p,k);sys=tf(num,den)sys =s + 1-------------s^2 + 5 s + 6线性系统的阶跃响应函数调用:[Y,T]=step(num,den) [Y ,T]=step(num,den,t) [Y ,T]=step(G) wn=1;zetas=[0:0.2:1,2,3,5];t=0:0.1:12;hold onfor i=1:length(zetas)G=tf(wn^2,[1,2*zetas(i)*wn,wn^2]);step(G,t),grid onendhold offtitle('当\omega_n=1,\zeta=0,0.2,...1,2,3,5时系统的单位阶跃响应曲线')024681012当ωn =1,ζ=0,0.2,...1,2,3,5时系统的单位阶跃响应曲线Time (seconds)A m p l i t u d e控制系统性能指标计算上升时间tr 峰值时间tp 超调量pos 调整时间ts线性系统的单位脉冲响应函数调用:[Y,T]=impulse(num,den) [Y ,T]=impulse(num,den,t) [Y,T]=impulse(G) G=tf(1,[1 0.4 1]);impulse(G);grid ontitle('G(s)=1/(s^2+0.4s+1)的单位脉冲响应曲线')05101520253035G(s)=1/(s 2+0.4s+1)的单位脉冲响应曲线Time (seconds)A m p l i t u d e线性系统的任意输入响应函数调用:[Y,T]=lsim(num,den,u,t) [Y ,T,X]=lsim(num,den,u,t,x0) u 指定时间序列上的输入值；x0指定初始状态；G=tf(1,[1 0.4 1]);t=0:0.1:10;r=t;y=lsim(G,r,t);plot(t,r,t,y,'o');grid on;title('G(s)=1/(s^2+0.4s+1)的单位斜坡响应曲线')G(s)=1/(s2+0.4s+1)的单位斜坡响应曲线线性系统的稳定性分析函数调用:roots(den)计算特征方程den全部极点den=conv([1 3],[1 2 2]);roots(den)ans =-3.0000 + 0.0000i-1.0000 + 1.0000i-1.0000 - 1.0000i绘制系统开环零极点函数调用:[p,z]=pzmap(num,den) [p,z]=pzmap(p,z) num=[1 2];den=conv([1 3],[1 2 2]);pzmap(num,den)title('G(s)=(s+2)/((s+3)(s^2+2s+2))开环零极点');-1.5-1-0.500.511.5G(s)=(s+2)/((s+3)(s 2+2s+2))开环零极点Real Axis (seconds -1)I m a g i n a r y A x i s (s e c o n d s -1)绘制系统根轨迹函数调用:[R,K]=rlocus(num,den) [R,K]=rlocus(num,den,K) num=[1 2];den=conv([1 3],[1 2 2]);rlocus(num,den);title('G(s)=(s+2)/((s+3)(s^2+2s+2))180度根轨迹曲线');G(s)=(s+2)/((s+3)(s 2+2s+2))180度根轨迹曲线Real Axis (seconds -1)I m a g i n a r y A x i s (s e c o n d s -1)计算给定一组特征根的根轨迹增益ki函数调用:[K,poles]=rlocfind(num,den) [K,poles]=rlocfind(num,den,P) num=[1 2];den=conv([1 3],[1 2 2]);rlocfind(num,den)title('G(s)=(s+2)/((s+3)(s^2+2s+2))给定一组特征根的根轨迹增益ki');根轨迹设计工具函数调用:rltool() rltool(num,den) rltool(G)绘制系统的奈式图函数调用:nyquist(sys) nyquist(sys,w)[re,im,w]=nyquist(sys) [re,im,w]=nyquist(sys,w) num=[1 2];den=conv([1 3],[1 2 2]);nyquist(num,den);title('G(s)=(s+2)/((s+3)(s^2+2s+2))奈式图曲线');-1-0.8-0.6-0.4-0.200.20.4G(s)=(s+2)/((s+3)(s 2+2s+2))奈式图曲线Real Axis I m a g i n a r y A x i s绘制系统的伯德图函数调用:bode(sys) bode(sys,w)[mag,phase,w]=bode(sys) [mag,phase,w]=bode(sys,w) num=10*[1 2];den=conv([1 3],[1 2 2]);bode(num,den);title('G(s)=(s+2)/((s+3)(s^2+2s+2))伯德曲线');M a g n i t u d e (d B )1010101010P h a s e (d e g )G(s)=10(s+2)/((s+3)(s 2+2s+2))伯德曲线Frequency (rad/s)系统的相对稳定性函数调用:margin(sys) [Gm,Pm,Wcg,Wcp]=margin(sys) num=10*[1 2];den=conv([1 3],[1 2 2]);[Gm,Pm,Wcg,Wcp]=margin(num,den)Gm =Inf 无穷大Pm =54.1970 相角裕度Wcg =InfWcp =2.8597基于MATLAB 的控制系统设计串联超前校正 leadc();I 型系统的渐进稳定性分析i ）绘制I 型系统根轨迹Real Axis (seconds -1)I m a g i n a r y A x i s (s e c o n d s -1)num1=500;den1=conv(conv([1,0],[0.007,1]),[0.9,1]);rlocus(num1,den1);grid onhold offtitle('根轨迹曲线');ii ）求解I 型系统特征根Gp=tf([500],conv([1,0],conv([0.007,1],[0.9,1]))); H=tf([1]);G=feedback(Gp,H);p=eig(G)p =1.0e+02 *-1.4658 + 0.0000i0.0131 + 0.2323i0.0131 - 0.2323i结论：根据根轨迹与求解其特征根的形式判断该I 型系统不稳定。
自动控制原理习题全解及MATLAB实验 第6章习题解答
系统开环传递函数为 G0 s
s0.1s
K
10.2s
1
，要求：
（1）系统响应斜坡信号 r(t)=t 时，稳态误差 ess 0.01 ；
(2) 系统相位裕量 ' 40 。
试用分析法设计一个串联滞后-超前校正装置。
解：（1）系统为Ⅰ型系统，在单位斜坡信号下
分稳态误差为
essr
1 k
令 essr
稳态性能与动态性能？ 答：PID 兼有 PI、PD 控制的特点，它相当于提供了一个积分环节与两个一阶微分环节。
积分环节改善稳态性能，两个一阶微分环节改善动态性能。 试分别叙述利用比例负反馈和微分负反馈包围振荡环节所起到的作用。
答：二阶振荡环节的频率特性为
1
T 2S 2 2 S 1
用比例负反馈 H(s)=h
0.2s 1 0.0143s 1
（5） Gc (s)
s 1 14s 1
0.2s 1 0.0143s
(6)
G
k
(s)
s(14s
100(s 1) 1)(0.1s 1)(0.0143s
1)
' 180 [90 arctan 7 arctan(14 7) arctan(0.1 7) arctan(0.0143 7)] 41.9 40
10lg( 12)
6dB
,
最后得出 c' m 4.47rad/s>4.4rad/s
(4) 确定校正装置的转折频率
1 m
2.2rad/s ，2 m
8.8rad/s ，T 1 0.45s ， 1
G(s)=
s
2.2 s
1 1
0.45s 0.11s
1 1
(最新版)自动控制原理MATLAB仿真实验报告
实验一 MATLAB及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性；二、预习要点1、系统的典型响应有哪些？2、如何判断系统稳定性？3、系统的动态性能指标有哪些？三、实验方法（一）四种典型响应1、阶跃响应：阶跃响应常用格式：1、；其中可以为连续系统，也可为离散系统。
2、；表示时间范围0---Tn。
3、；表示时间范围向量T指定。
4、；可详细了解某段时间的输入、输出情况。
2、脉冲响应：脉冲函数在数学上的精确定义：其拉氏变换为：所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式：①；②③（二）分析系统稳定性有以下三种方法：1、利用pzmap绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den) 运行结果： p =-1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
自动控制原理MATLAB仿真实验报告
实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点1、 系统的典型响应有哪些2、 如何判断系统稳定性3、 系统的动态性能指标有哪些 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den)运行结果： p =+ - + -P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
matlab课程设计自动控制原理
matlab课程设计自动控制原理一、教学目标本课程的目标是使学生掌握自动控制原理的基本概念和MATLAB在自动控制领域的应用。
通过本课程的学习，学生应能理解自动控制系统的组成、工作原理和设计方法，熟练运用MATLAB进行自动控制系统的分析和仿真。
知识目标：学生通过本课程的学习，应掌握自动控制基本理论、MATLAB基本操作和自动控制系统仿真方法。
技能目标：学生应能熟练使用MATLAB进行自动控制系统的建模、仿真和分析，具备一定的实际问题解决能力。
情感态度价值观目标：培养学生对自动控制技术的兴趣和热情，提高学生运用现代技术手段进行科学研究的能力，培养学生的创新精神和团队合作意识。
二、教学内容本课程的教学内容主要包括自动控制原理的基本概念、MATLAB的基本操作和自动控制系统的仿真方法。
1.自动控制原理：包括自动控制系统的组成、数学模型、稳定性分析、控制器设计和校正方法等。
2.MATLAB基本操作：包括MATLAB的安装和启动、变量和数据类型、矩阵运算、编程和函数的使用等。
3.自动控制系统仿真：包括MATLAB仿真环境的设置、Simulink的介绍和应用、控制系统仿真的方法和步骤等。
三、教学方法本课程采用讲授法、案例分析法和实验法相结合的教学方法。
1.讲授法：通过教师的讲解，使学生掌握自动控制原理的基本概念和MATLAB的基本操作。
2.案例分析法：通过分析实际案例，使学生理解和掌握自动控制系统的建模和仿真方法。
3.实验法：通过上机实验，使学生熟练掌握MATLAB自动控制系统仿真工具的使用，提高学生的实际操作能力。
四、教学资源本课程的教学资源包括教材、多媒体资料和实验室设备。
1.教材：选用《自动控制原理》和《MATLAB基础教程》作为主要教材，为学生提供系统的理论知识和实践指导。
2.多媒体资料：制作课件、教学视频等，以图文并茂的形式展示自动控制原理和MATLAB的操作方法。
3.实验室设备：提供计算机和MATLAB软件，供学生进行自动控制系统的仿真实验。
MATLAB自动控制原理实验
实验七 控制系统的MATLAB 分析一、 实验目的1)、掌握如何使用Matlab 进行系统的时域分析 2)、掌握如何使用Matlab 进行系统的频域分析 3)、掌握如何使用Matlab 进行系统的根轨迹分析 4)、掌握如何使用Matlab 进行系统的稳定性分析 5)、掌握使用Bode 图法进行控制系统设计的方法 二、 实验内容 1、时域分析法根据下面传递函数模型：绘制其单位阶跃响应曲线并从图上读取最大超调量，绘制系统的单位脉冲响应、零输入响应曲线。
1)、某单位负反馈系统传递函数为：8106)65(5)(232+++++=s s s s s s Gt (seconds)c (t )t (seconds)c (t )结论：2)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ζ=0.7,ωn 取2、4、6、8、10、12的单位阶跃响应。
Step ResponseTime (seconds)00.51 1.52 2.53 3.54结论：3)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。
Time (seconds)结论：2、频率分析法根据下面传递函数模型，绘制出系统的频率响应曲线，包括Bode 图和Nyquist 图，并从图上读取相角交接频率、截止频率，并求出幅值裕度和相角裕度。
1)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ζ=0.7,ωn 取2)、4)、6)、8、1)0、1)2)的伯德图和奈奎斯特图。
Wn=2M a g n i t u d e (d B )10101010P h a s e (d e g )Bode DiagramGm = Inf dB (at Inf rad/s) , Pm = 164 deg (at 0.4 rad/s)Frequency (rad/s)Real AxisI m a g i n a r y A x i sWn=4M a g n i t u d e (d B )10101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=6M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=8M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=10M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=12M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i s2)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的伯德图和奈奎斯特图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
e显示综合后系统的阶跃响应
二、运行结果
输入传递函数分子分母
生成根轨迹图
选择点并得到该点各项参数在下方输出面板输出
获得阶跃响应图
用rltool(辅助,选择合适的插入零点
输入零点,并得到根轨迹图
% varargout cell array for returning output args (see VARARGOUT;
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @Liushuai20122510_OpeningFcn, ...
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before Liushuai20122510_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
function edit1_Callback(hObject, eventdata, handles
% hObject handle to edit1 (see GCBO
% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA
选择根轨迹图上的任一点,得到数据,在下方输出面板输出
得到阶跃响应图
三、运行说明
第一步,在请输入分子后的输入框输入传递函数分子的矩阵,在下一输入框输入传递函数分母并按“生成根轨迹图”按钮获得根轨迹
第二步,按选择点并显示各参数获得根轨迹图上任一点的各项数据,数据全部输出在下方输出面板
第三步,按“生成阶跃响应图”按钮可以获得该函数的阶跃响应
四、附录............................................ 9代码. ............................................. 9参考文献......................................... 17
一、题目
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles
num = str2num(get(handles.edit1,'String' ;
den = str2num(get(handles.edit2,'String' ;
G = tf(num,den;
rlocus(G
hold on ;
sgrid(0.707,[]
hold off ;
% --- Executes on button press in pushbutton2.
% varargin command line arguments to Liushuai20122510 (see VARARGIN
% Choose default command line output for Liushuai20122510
handles.output = hObject;
% Update handles structure
第四步,在“请输入插入零点”后的输入框中输入参数,并按“生成综合后根轨迹图”按钮产生根轨迹(可以通过点击“根轨迹校正”按钮,调用工具箱拖动零点进行快速查看根轨迹图,选择合适的根轨迹再在输入框中输入零点的值
第五步,按“选择点并显示各参数(综合后系统”选取各点,查阅参数,数据输出在下方输出面板上
第六步,按“生成阶跃响应图(综合后系统”可以得到综合后系统的阶跃响应
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA
最后,点击“退出”结束程序
四、附录
代码:
function varargout = Liushuai20122510(varargin
% LIUSHUAI20122510 MATLAB code for Liushuai20122510.fig
% LIUSHUAI20122510, by itself, creates a new LIUSHUAI20122510 or raises the existing
自控控制原理MATLAB程序设计作业
指导老师:汪晓宁
一、题目............................................ 2
二、运行结果........................................ 3
三、程序说明........................................ 8
function pushbutton3_Callback(hObject, eventdata, handles
% hObject handle to pushbutton3 (see GCBO
% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA close(gcbf;
% singleton*.
%
% H = LIUSHUAI20122510 returns the handle to a new LIUSHUAI20122510 or the handle to
% the existing singleton*.
%
% LIUSHUAI20122510('CALLBACK',hObject,eventData,handles,...
'gui_OutputFcn', @Liushuai20122510_OutputFcn, ... 'gui_LayoutFcn', [] , ...
'gui_Callback', [];
if nargin && ischar(varargin{1}
gui_State.gui_Callback = str2func(varargin{1};
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:};
else
gui_mainfcn(gui_State, varargin{:};
end
% End initialization code - DO NOT EDIT
function pushbutton2_Callback(hObject, eventdata, handles
% hObject handle to pushbutton2 (see GCBO
% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA
% --- Executes just before Liushuai20122510 is made visible.
function Liushuai20122510_OpeningFcn(hObject, eventdata, handles, varargin
% This function has no output args, see OutputFcn.
num = str2num(get(handles.edit1,'String' ;
den = str2num(get(handles.edit2,'String' ;
sys=tf(num,den;
[y,t]=step(sys;
plot(t,y;
grid on ;
% --- Executes on button press in pushbutton3.