随机模拟系统仿真与MATLAB实现

自控课设MATLAB超前滞后校正概要

课程设计任务书 学生姓名： 张弛 专业班级： 电气1002班 指导教师： 刘志立 工作单位： 自动化学院 题 目: 用MATLAB 进行控制系统的滞后－超前校正设计 初始条件：已知一单位反馈系统的开环传递函数是 ) 2)(1()(++= s s s K s G 要求系统的静态速度误差系数110-≥S K v ， 45≥γ。 要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要 求） 1、 M ATLAB 作出满足初始条件的最小K 值的系统伯德图，计算系统的幅值裕量和相位裕量。 2、前向通路中插入一相位滞后－超前校正，确定校正网络的传递函数。 3、用MATLAB 画出未校正和已校正系统的根轨迹。 4、用Matlab 对校正前后的系统进行仿真分析，画出阶跃响应曲线，计算其时域性能指标。 5、课程设计说明书中要求写清楚计算分析的过程，列出MATLAB 程序和MATLAB 输出。说明书的格式按照教务处标准书写。 时间安排: 指导教师签名: 年 月 日 系主任（或责任教师）签名: 年 月 日

摘要 (3) 1基于频率响应法校正设计概述 (4) 2串联滞后-超前校正原理及步骤 (5) 2.1滞后超前校正原理 (5) 2.2滞后-超前校正的适用范围 (6) 2.3串联滞后-超前校正的设计步骤 (6) 3串联滞后-超前校正的设计 (7) 3.1待校正系统相关参数计算及稳定性判别 (7) 3.1.1判断待校正系统稳定性 (7) 3.1.2绘制待校正系统的伯德图 (8) 3.1.3绘制待校正系统的根轨迹图 (9) 3.1.4绘制待校正系统的单位阶跃响应曲线 (10) 3.1.5利用SIMULINK进行控制系统建模仿真 (11) 3.2滞后超前-网络相关参数的计算 (12) 3.3对已校正系统的验证及稳定性分析 (15) 3.3.1绘制已校正系统的伯德图 (15) 3.3.2判断已校正系统的稳定性 (16) 3.3.3绘制已校正系统的根轨迹图 (17) 3.3.4绘制已校正系统的单位阶跃响应曲线 (18) 3.3.5利用SIMULINK进行控制系统建模仿真 (19) 3.3.6串联滞后-超前校正设计小结 (20) 4心得体会 (21) 参考文献 (21) 附录 (22)

信号与系统的MATLAB仿真

信号与系统的MATLAB 仿真 一、信号生成与运算的实现 1.1 实现)3(sin )()(π±== =t t t t S t f a )(sin )sin()sin(sin )()(t c t t t t t t t S t f a '=' '== ==πππ π ππ m11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果： 1.2 实现)10() sin()(sin )(±== =t t t t c t f ππ m12.m t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数 plot(t,f); % 绘制sinc(t)的波形 运行结果： 1.3 信号相加：t t t f ππ20cos 18cos )(+= m13.m syms t; % 定义符号变量t f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：

1.4 信号的调制：t t t f ππ50cos )4sin 22()(+= m14.m syms t; % 定义符号变量t f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.5 信号相乘：)20cos()(sin )(t t c t f π?= m15.m t=-5:0.01:5; % 定义时间范围向量 f=sinc(t).*cos(20*pi*t); % 计算函数f(t)=sinc(t)*cos(20*pi*t) plot(t,f); % 绘制f(t)的波形 title('sinc(t)*cos(20*pi*t)'); % 加注波形标题 运行结果：

Matlab通信系统仿真实验报告

Matlab通信原理仿真 学号： 2142402 姓名：圣斌

实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);

基于Matlab的自动控制系统设计与校正

自动控制原理课程设计设计题目：基于Matlab的自动控制系统设计与校正

目录 第一章课程设计内容与要求分析.................................................... 错误！未定义书签。 1.1设计内容 (1) 1.2 设计要求 (1) 1.3 Matlab软件 (2) 1.3.1基本功能 (2) 1.3.2应用 (2) 第二章控制系统程序设计 (4) 2.1 校正装置计算方法 (4) 2.2 课程设计要求计算 (4) 第三章利用Matlab仿真软件进行辅助分析 (6) 3.1校正系统的传递函数 (6) 3.2用Matlab仿真 (6) 3.3利用Matlab/Simulink求系统单位阶跃响应 (8) 3.2.1原系统单位阶跃响应 (8) 3.2.2校正后系统单位阶跃响应 (8) 3.2.3校正前、后系统单位阶跃响应比较 (8) 3.4硬件设计 (8) 3.4.1在计算机上运行出硬件仿真波形图 (9) 课程设计心得体会 (10) 参考文献 (12)

1 第一章 课程设计内容与要求分析 1.1设计内容 针对二阶系统 )1()(+= s s K s W ， 利用有源串联超前校正网络（如图所示）进行系统校正。当开关S 接通时为超前校正装置，其传递函数 11 )(++-=Ts Ts K s W c c α， 其中 1 3 2R R R K c += ， 1 ) (13243 2>++ =αR R R R R ，C R T 4=， “-”号表示反向输入端。若Kc=1，且开关S 断开，该装置相当于一个放大系数为1的放大器（对原系统没有校正作用）。 1.2 设计要求 1 1.0)(≤∞e ，开环截止频率ω’≥45°； 2 3） 4）设校正装置网络元件参数R4、5R=100K ，C=1μF 、10μF 若干个）； 6）利用Matlab 仿真软件辅助分析，绘制校正前、后及校正装置对数频率特性曲线,并验算设计结果； 7）在Matlab-Simulink 下建立系统仿真模型，求校正前、后系 统单位阶跃响应特性，并进行系统性能比较； 8）利用自动控制原理实验箱完成硬件设计过程，包括：搭建校正前后 c R R

春MATLAB仿真期末大作业

MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：

2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。

仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。

matlab控制系统仿真.

课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称龙蟠学院 专业自动化 班级M10自动化 学生姓名 学号 课程设计地点 C208 课程设计学时一周 指导教师应明峰 金陵科技学院教务处制成绩

一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。

（d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统PID控制器取参数分别为：50 2 5 有积分作用单回路控制系统PID控制器取参数分别为：50 0 5

大比例作用单回路控制系统PID控制器取参数分别为：50 0 0 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响应，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏；

matlab实验九控制系统的PI校正设计及仿真

实验九控制系统的PI校正设计及仿真 一、实验目的 1. 应用频率综合法对系统进行PI校正综合设计； 2. 学习用MATLAB对系统性能进行仿真设计、分析； 二、实验设计原理与步骤 1?设计原理 滞后校正(亦称PI校正)的传递函数为： 其对数频率特性如图9-1所示，参数表征滞后校正的强度。 2 ?设计步骤 基于频率法的综合滞后校正的步骤是： (1)根据静态指标要求，确定开环比例系数K，并按已确定的K画出系统固有部分的Bode 图； (2)根据动态指标要求试选c，从Bode图上求出试选的c 点的相角，判断是否满足相 位裕度的要求(注意计入滞后校正带来的50 ~ 120的滞后量)，如果满足，转下一步。否则， 如果允许降低 c ,就适当重选较低的 C ； (3)从图上求出系统固有部分在c点的开环增益Lg ( c) o如果Lg ( c) >0令Lg ( c) =20lg ，求出，就是滞后校正的强度，如果Lg (c) < 0，则无须校正，且可将开环比例系数提高。 (4)选择21 1 1 (~ ) C，进而确定 1 1 o T 5 10 T (5)画出校正后系统的Bode图，校核相位裕量。 滞后校正的主要作用是降低中频段和高频段的开环增益,但同时使低频段的开环增益 不受影响，从而达到兼顾静态性能与稳定性。它的副作用是会在c点产生一定的相角滞后。 三、实验内容 练习9-1设系统原有开环传递函数为： 系统的相位裕度丫40° 系统的开环比例系数K= 5 S-1

截止频率为c=0 . 5 S1 要求： (1)用频率法设计满足上述要求的串联滞后校正控制器； ( 2)画出校正前后的Bode 图 ( 3)用Simulink 对校正前后的闭环系统进行仿真，求出其阶跃响应； (4)分析设计效果。 k0=5;n1=1; d1=conv(conv([1,0],[1,1]),[0.5 1]); w=logspace(-1,3,1000); sope=tf(k0*n1,d1); figure(1) margin(sope);grid on wc=0.5; num=sope.num{1}; den=sope.den{1}; na=polyval(num,j*wc); da=polyval(den,j*wc); g=na/da; g1=abs(g); h=20*log10(g1); beta=10A(h/20); t=10/wc; bt=beta*t; gc=tf([t,1],[bt,1]) sys1=sope*gc [mag1,phase1,w]=bode(sys1,w); [gm1,pm1,wcg1,wcp1]=margin(mag1,phase1,w); figure(2) margin(sys1); grid on; 滞后： 20 s + 1 173.5 s + 1 练习9-2 设被控对象开环传递函数为： 系统的相位裕度丫400

电机学matlab仿真大作业报告

. 基于MATLAB的电机学计算机辅助分析与仿真 实验报告

一、实验内容及目的 1.1 单相变压器的效率和外特性曲线 1.1.1 实验内容 一台单相变压器，N S =2000kVA, kV kV U U N N 11/127/21=，50Hz ，变压器的参数 和损耗为008.0* ) 75(=C k o R ，0725.0*=k X ，kW P 470=，kW P C KN o 160)75(=。 (1)求此变压器带上额定负载、)(8.0cos 2滞后=?时的额定电压调整率和额定效率。 (2)分别求出当0.1,8.0,6.0,4.0,2.0cos 2=?时变压器的效率曲线，并确定最大效率和达到负载效率时的负载电流。 (3)分析不同性质的负载（),(8.0cos 0.1cos ),(8.0cos 222超前，滞后===???）对变压器输出特性的影响。 1.1.2 实验目的 （1）计算此变压器在已知负载下的额定电压调整率和额定效率 （2）了解变压器效率曲线的变化规律 （3）了解负载功率因数对效率曲线的影响 （4）了解变压器电压变化率的变化规律 （5）了解负载性质对电压变化率特性的影响 1.1.3 实验用到的基本知识和理论 （1）标幺值、效率区间、空载损耗、短路损耗等概念 （2）效率和效率特性的知识 （3）电压调整率的相关知识 1.2串励直流电动机的运行特性 1.2.1实验内容 一台16kw 、220V 的串励直流电动机，串励绕组电阻为0.12Ω，电枢总电阻为0.2Ω。电动势常数为.电机的磁化曲线近似的为直线。其中为比例常数。假设电枢电流85A 时，磁路饱和(为比较不同饱和电流对应的效果，饱和电流可以自己改变）。

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.

基于Matlab仿真PID校正

基于matlab 仿真的pid 校正总结 PID 控制器是目前在过程控制中应用最为普遍的控制器，它通常可以采用以下几种形式：比例控制器，0;D I K K ==比例微分控制器，0;I K =比例积分控制器，0; D K =标准控制器。 下面通过一个例子来介绍PID 控制器的设计过程。 假设某弹簧（阻尼系统）如图1所示，1,10/,20/M kg f N s m k N m ==?=。让 我们来设计不同的P 、PD 、PI 、PID 校正装置，构成反馈系统。来比较其优略。 系统需要满足： （1） 较快的上升时间和过渡过程时间； （2） 较小的超调； （3） 无静差。 图1 弹簧阻尼系统 系统的模型可描述如下：

控制系统建模与仿真论文( 2011) () 2 ()1()X s G s F s M s fs k = = ++ （1）、绘制未加入校正装置的系统开环阶跃响应曲线。 根据系统的开环传递函数，程序如下： clear; t=0:0.01:2; num=1; den=[1 10 20]; c=step(num,den,t); plot(t,c); xlabel('Time-Sec'); ylabel('y'); title('Step Response'); grid; 系统的阶跃响应曲线如图2 图2 未加入校正时系统的开环阶跃响应曲线

（2）、加入P 校正装置 我们知道，增加p K 可以降低静态误差，减少上升时间和过渡时间，因此首先选择P 校 正，也就是加入一个比例放大器。此时，系统的闭环传递函数为： 2 ()10(20) p c p K G s s s K = +++ 此时系统的静态误差为 120p p K K - +。所以为了减少静差，可以选择系统的比例增益 为300p K =。这样就可以把静差缩小到0.0625。虽然系统的比例系数越大，静差越小，但是比例系数也不能没有限制地增大，它会受到实际物理条件和放大器实际条件的限制。一般取几十到几百即可。增大比例增益还可以提高系统的快速性。 加入P 校正后，程序如下： clear; t=0:0.01:2; Kp=300; num=[Kp]; den=[1 10 (20+Kp)]; c=step(num,den,t); plot(t,c); xlabel('Time-Sec'); ylabel('y'); title('Step Response'); gird; 加入P 校正后系统的闭环阶跃响应曲线如图3

matlab 大作业

上海电力学院 通信原理Matlab仿真 实验报告 实验名称： 8QAM误码率仿真 试验日期： 2014年 6月3日 专业：通信工程 姓名：罗侃鸣 班级： 2011112班 学号： 20112272

一、实验要求 写MATLAB程序，对图示的信号星座图完成M=8的QAM通信系统Monte Carlo仿真，在不同SNRindB=0：15时，对N=10000（3比特）个符号进行仿真。画出该QAM系统的符号误码率。 二、实验原理 1 QAM调制原理 QAM(Quadrature Amplitude Modulation)正交幅度调制技术，是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽内的正交性，实现两路并行的数字信息的传输。该调制方式通常有8QAM，16QAM，64QAM。 QAM调制实际上就是幅度调制和相位调制的组合，相位+ 幅度状态定义了一个数字或数字的组合。QAM的优点是具有更大的符号率，从而可获得更高的系统效率。通常由符号率确定占用带宽。因此每个符号的比特（基本信息单位）越多，频带效率就越高。 调制时，将输入信息分成两部分：一部分进行幅度调制；另一部分进行相位调制。对于星型8QAM信号，每个码元由3个比特组成，可将它分成第一个比特和后两个个比特两部分。前者用于改变信号矢量的振幅，后者用于差分相位调制，通过格雷编码来改变当前码元信号矢量相位与前一码元信号矢量相位之间的相位差。 QAM是一种高效的线性调制方式，常用的是8QAM，16QAM，64QAM等。当随着M 的增大，相应的误码率增高，抗干扰性能下降。 2 QAM星座图 QAM调制技术对应的空间信号矢量端点分布图称为星座图。QAM的星座图呈现星状分

MATLAB Simulink系统建模与仿真 实验报告

MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************

实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：

1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：

matlab实验九控制系统的PI校正设计及仿真

实验九 控制系统的PI 校正设计及仿真 一、 实验目的 1. 应用频率综合法对系统进行PI 校正综合设计； 2．学习用MA TLAB 对系统性能进行仿真设计、分析； 二、实验设计原理与步骤 1．设计原理 滞后校正（亦称PI 校正）的传递函数为： 其对数频率特性如图9-1所示，参数β表征滞后校正的强度。 2．设计步骤 基于频率法的综合滞后校正的步骤是： （1）根据静态指标要求，确定开环比例系数K ，并按已确定的K 画出系统固有部分的Bode 图； （2）根据动态指标要求试选c ω，从Bode 图上求出试选的c ω点的相角，判断是否满足相位裕度的要求（注意计入滞后校正带来的0012~5的滞后量），如果满足，转下一步。否则，如果允许降低c ω，就适当重选较低的c ω； （3）从图上求出系统固有部分在c ω点的开环增益Lg （c ω）。如果Lg （c ω）>0令Lg （c ω）=20lg β，求出β，就是滞后校正的强度，如果Lg （c ω）〈0，则无须校正，且可将开环比例系数提高。 （4）选择C T ωω)10 1~51(12==，进而确定T βω11=。 （5）画出校正后系统的Bode 图，校核相位裕量。 滞后校正的主要作用是降低中频段和高频段的开环增益，但同时使低频段的开环增益不受影响，从而达到兼顾静态性能与稳定性。它的副作用是会在c ω点产生一定的相角滞后。 三、实验内容 练习9-1 设系统原有开环传递函数为： 系统的相位裕度γ0 40≥ 系统的开环比例系数K= 5 S -1 截止频率为c ω=0．5 S -1

（1）用频率法设计满足上述要求的串联滞后校正控制器； （2）画出校正前后的Bode图 （3）用Simulink对校正前后的闭环系统进行仿真，求出其阶跃响应；（4）分析设计效果。 k0=5;n1=1; d1=conv(conv([1,0],[1,1]),[ 1]); w=logspace(-1,3,1000); sope=tf(k0*n1,d1); figure(1) margin(sope);grid on wc=; num={1}; den={1}; na=polyval(num,j*wc); da=polyval(den,j*wc); g=na/da; g1=abs(g); h=20*log10(g1); beta=10^(h/20); t=10/wc; bt=beta*t; gc=tf([t,1],[bt,1]) sys1=sope*gc [mag1,phase1,w]=bode(sys1,w); [gm1,pm1,wcg1,wcp1]=margin(mag1,phase1,w); figure(2) margin(sys1); grid on; 滞后：

matlab机电系统仿真大作业

一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环

3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置

M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置

实验一 典型环节的MATLAB仿真汇总

实验一 典型环节的MATLAB 仿真 一、实验目的 1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。 2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。 3．定性了解各参数变化对典型环节动态特性的影响。 二、SIMULINK 的使用 MATLAB 中SIMULINK 是一个用来对动态系统进行建模、仿真和分析的软件包。利用SIMULINK 功能模块可以快速的建立控制系统的模型，进行仿真和调试。 1．运行MATLAB 软件，在命令窗口栏“>>”提示符下键入simulink 命令，按Enter 键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK 仿真 环境下。 2．选择File 菜单下New 下的Model 命令，新建一个simulink 仿真环境常规模板。 3．在simulink 仿真环境下，创建所需要的系统 三、实验内容 按下列各典型环节的传递函数，建立相应的SIMULINK 仿真模型，观察并记录其单位阶跃响应波形。 ① 比例环节1)(1=s G 和2)(1=s G 实验处理：1)(1=s G SIMULINK 仿真模型

波形图为： 实验处理：2)(1=s G SIMULINK 仿真模型 波形图为： 实验结果分析：增加比例函数环节以后，系统的输出型号将输入信号成倍数放大. ② 惯性环节11)(1+= s s G 和15.01)(2+=s s G 实验处理：1 1 )(1+=s s G SIMULINK 仿真模型

波形图为： 实验处理：1 5.01 )(2+= s s G SIMULINK 仿真模型 波形图为： 实验结果分析：当1 1 )(1+= s s G 时，系统达到稳定需要时间接近5s,当

电机大作业(MATLAB仿真-电机特性曲线)

电机大作业 专业班级：电气XXXX 姓名：XXX 学号：XXX 指导老师：张威

一、研究课题（来源：教材习题 4-18 ） 1. 74 、R 2 0.416 、X 2 3.03 、R m 6. 2 X m 75 。电动机的机械损耗p 139W,额定负载时杂散损耗p 320W, 试求额定负载时的转差率、定子电流、定子功率因数、电磁转矩、输出转矩和效 率。 二、编程仿真 根据T 形等效电路: 3D - R Q 运用MATLAB 进行绘图。MATLAB 文本中，P N PN ，U N UN ，尺 R 1， X 1 X1 , R 2 R 2，X 2 X 2，R m Rm, X m Xm ，p pjixiesunh ao ， p pzasansunhao 。定子电流I11，定子功率因数 Cosangle1，电磁转矩Te ， 效率 Xiaolv 。 1.工作特性曲线绘制 MATLA 文本： R1=0.715;X 仁1.74;Rm=6.2;Xm=75;R2=0.416;X2=3.03;pjixiesu nhao=139; pzasa nsu nhao=320;p=2;m 仁 3; ns=1500;PN=17000;UN=380;fN=50; Z1=R1+j*X1; Zm=Rm+j*Xm; for i=1:2500 s=i/2500; nO=n s*(1-s); Z2=R2/s+j*X2; Z=Z1+Zm*Z2/(Zm+Z2); 有一台三相四极的笼形感应电动机， 参数为P N 17kW 、U N 380V (△联 Rm 结）、尺 0. 715 、X j lcr S

U1=UN; I1=U1/Z; l110=abs(l1); An gle 仁an gle(ll); Cosa ngle10=cos(A ngle1); P仁3*U1*l110*Cosa ngle10; l2=l1*Zm/(Zm+Z2); Pjixie=m1*(abs(I2))A2*(1-s)/s*R2; V=(1-s)*pi*fN; Te0=Pjixie/V; P20=Pjixie-pjixies un hao-pzasa nsun hao; Xiaolv0=P20/P1; P2(i)=P20; n (i)=n0; l11(i)=l110; Cosa ngle1(i)=Cosa ngle10; Te(i)=Te0; Xiaolv(i)=Xiaolv0; hold on; end figure(1) plot(P2, n); xlabel('P2[W]');ylabel(' n[rpm]'); figure(2) plot(P2,l11); xlabel('P2[W]');ylabel('l1[A]'); figure(3) plot(P2,Cosa nglel); xlabel('P2[W]');ylabel('go nglvyi nshu'); figure(4) plot(P2,Te); xlabel('P2[W]');ylabel('Te[Nm]'); figure(5) plot(P2,Xiaolv); xlabel('P2[W]');ylabel('xiaolv');

MATLAB与系统仿真实验指导书

实验指导书（标准格式） 《MATLAB与系统仿真》 实验指导书 课程编号：ME2121025课程名称：MATLAB与系统仿真适应专业：电气工程及其自动化课程类别：限选 实验教学种类：上机课程总学时：30＋6 实验学时：6 执笔人：勾燕洁 西安电子科技大学机电工程学院 2006 年9 月

目录 实验一MA TLAB的基本使用及矩阵操作------------------- 3 实验二MA TLAB编程与图形处理---------------------------- 9 实验三MA TLAB中的数值运算与系统仿真---------------- 15

实验一MA TLAB的基本使用及矩阵操作 一、实验目的和要求 熟悉MA TLAB的界面和基本操作，掌握矩阵的建立方法及各种运算。 1．熟悉MA TLAB软件的界面和帮助系统。 2．掌握MA TLAB软件中关于矩阵建立、矩阵初等变换以及矩阵算术、关系、逻辑运算的各种命令。 3．掌握MA TLAB软件中M函数和M文件的编写，以及程序结构与控制，学会编写一般程序。 二、实验内容 1．启动与退出 2．数、数组、矩阵的输入 3．MA TLAB的基本命令 4．矩阵大小的测试 5．矩阵元素的操作 6．特殊矩阵的产生 7．矩阵的算术运算 8．矩阵的关系运算 9．矩阵的逻辑运算 10．常用函数 三、实验仪器、设备（软、硬件）及仪器使用说明 PC机一台，注意正确开、关机及打开软件。 四、实验原理 无 五、实验方法与步骤 1．启动与退出 双击MA TLAB图标，进入MA TLAB命令窗口，即可输入命令，开始运算； 观察各个窗口以及菜单； 单击File菜单中的Exit，或使用MATLAB的Exit命令退出。

运动控制MATLAB仿真

大作业: 直流双闭环调速MATLAB仿真 运动控制技术课程名称： 名：姓电气学院院：学 自动化业：专 号：学 孟濬指导教师： 2012年6月2日

------------------------------------- -------------学浙大江 李超 一、Matlab仿真截图及模块功能描述 Matlab仿真截图如下，使用Matlab自带的直流电机模型： 模块功能描述： ⑴电机模块（Discrete DC_Machine）：模拟直流电机 ⑵负载转矩给定（Load Torque）：为直流电机添加负载转矩 ⑶Demux：将向量信号分离出输出信号 ⑷转速给定（Speed Reference）：给定转速 ⑸转速PI调节（Speed Controller）：转速PI调节器，对输入给定信号与实际信号

的差值进行比例和积分运算，得到的输出值作为电流给定信号。改变比例和积分运算系数可以得到不同的PI控制效果。 ⑹电流采样环节（1/z）：对电流进行采样，并保持一个采样周期 ⑺电流滞环调节（Current Controller）：规定一个滞环宽度，将电流采样值与给定值进行对比，若：采样值＞给定值+0.5*滞环宽度，则输出0； 若：采样值＜给定值—0.5*滞环宽度，则输出1； 若：给定值—0.5*滞环宽度＜采样值＜给定值+0.5*滞环宽度，则输出不变 输出值作为移相电压输入晶闸管斩波器控制晶闸管触发角 ：根据输入电压改变晶闸管触发角，从而改变电机端电压。GTO⑻晶闸管斩波．⑼续流二极管D1：在晶闸管关断时为电机续流。 ⑽电压传感器Vd：测量电机端电压 ⑾示波器scope：观察电压、电流、转速波形 系统功能概括如下：直流电源通过带GTO的斩波器对直流电机进行供电，输出量电枢电流ia和转速wm通过电流环和转速环对GTO的通断进行控制，从而达到对整个电机较为精确的控制。 下面对各个部分的功能加以详细说明： （1）直流电机 双击电动机模块，察看其参数：

一阶、二阶系统时域和频域仿真

西安交通大学 基于MATLAB/Simulink 的一阶、二阶系统的时域和频 域仿真 ——以单位阶跃信号为输入信号 日期：2013年4月 一阶系统时域和频域仿真 1、建立一阶系统典型数学模型 ()1 1 G s Ts =+ 2、建立simulink 仿真方框图

1T.s+1 Transfer Fcn Step Scope ① 时间常数T=1时，一阶系统时域响应为 12345678 910 00.5 1 一阶系统时域相应（T=1） Matlab 程序： %一阶系统仿真编程 num=[1]; den=[1 1]; bode(num,den); grid on ; gtext('低频段频率-20dB/dec'); 运行程序，有时间常数T=1时，一阶系统的频域响应为

10 -210 -1 10 10 1 10 2 -90-45 一阶系统频域响应 P h a s e (d e g ) Bode Di a gram Frequency (rad/s) -40-30-20-100 低频段斜率-20dB/dec System: sys Frequency (rad/s): 1.01Magni t ude (dB): -3.07 M a g n i t u d e (d B ) ② 时间常数T=3时，一阶系统单位阶跃时域响应 12345678910 00.5 1 一阶系统单位阶跃响应（T=3） Matlab 程序： %一阶系统仿真编程 num=[1]; den=[3 1]; bode(num,den);

grid on ; gtext('低频段频率-20dB/dec'); 运行程序，有时间常数T=3时，一阶系统的频域响应为 10 -210 -1 10 10 1 -90-45 P h a s e (d e g ) Bode Di a gram Frequency (rad/s) -30-20-100 低频段频率-20dB/dec System: sys Frequency (rad/s): 0.334Magni t ude (dB): -3.03 M a g n i t u d e (d B ) 3、分析以上一阶系统在不同时间常数下的单位阶跃响应，可以看出时间常数越小，系统响应越快；而且一阶系统的转角频率为1/T ，在转角频率以上时，幅频特性曲线以-20dB/dec 下降，而相频特性以0°和90°为渐近线。

基于MATLAB的系统仿真实验

基于MATLAB的系统仿真实验 实验指导书 新乡学院 二○○八年三月

说明 MATLAB是MATrix LABoratory的缩写，是一种基于矩阵的数学与工程计算系统，可以用作动态系统的建模与仿真。研究系统的结构和参数的变化对系统性能的影响可利用MATLAB强大的计算和作图功能，因此本实验采用MATLAB仿真研究连续控制系统和离散控制系统的性能分析过程。通过该实验提高学生对控制系统的分析与设计能力，加深对《自动控制原理》课程内容的理解。 由于学时有限，该实验由同学们参考有关MATLAB应用的书籍利用课后时间完成。

第一部分 基于MATLAB 连续系统的仿真 一、 实验目的 在研究系统的结构和参数的变化对系统性能的影响时，采用解析和作图的方法比较麻烦，而且误差也大，用MATLAB 仿真实现则简单方便，精度高。本实验采用MATLAB 实现控制系统的数学描述、控制系统的时域分析及根轨迹和频率特性分析。通过该实验，加深学生对系统阶次，型号，参数与系统性能的关系的理解。 二、实验环境 在计算机Windows 环境下安装好MATLAB6.3以上版本后，双击MATLAB 图标或成“开始”菜单打开MATLAB ，即可进入MATLAB 集成环境。 三、MATLAB 应用实例 1．拉氏变换和反变换 例 求22)(2++=t t t f 的拉氏变换 解 键入 syms s t; ft=t^2+2*t+2; st=laplace(ft,t,s) 运行结果为 st= 2/s^3+2/s^2+2/s 例 求) 2)(34(6 )(2 ++++= s s s s s F 的拉氏反变换 解 键入 syms s t; Fs=(s+6)/(s^2+4*s+3)/(s+2); ft=ilaplace(Fs,s,t) 运行结果为 ft= 3/2*exp(-3*t)+5/2*exp(-t)-4*exp(-2*t) 2。求根运算 例 求多项式 43)(2 3 ++=s s s p 的根，再由根建多项式。 解 键入 p=[1 3 0 4]; r=root(p) 运行结果为 r= -3.3553 0.1777+1.0773i 1.7777-1.0773i