转速反馈单闭环直流调速系统Matlab仿真

计算机仿真技术作业一——转速反馈单闭环直流调速系统仿真

姓名：朱龙胜

班级：电气1102

学号： 11291065

日期： 2014年5月4日

指导老师：郭希铮

北京交通大学

计算机仿真技术作业一

题目：转速反馈单闭环直流调速系统仿真

直流电机模型框图如下图所示，仿真参数为R=0.6Ω，T l=0.00833，T m=0.045，Ce=0.2。本次仿真采用算法为ode45，仿真时间5s。

图1 直流电机模型

1、开环仿真：

用simulink实现上述直流电机模型，直流电压U d0取220V，

time n 转速

time

s

静差率

仿真算法为ode23

0.5

1

1.5

2

2.53

3.5

4

4.5

5

time n 转速

time

s

静差率

选择仿真算法ode23时，在2.5s 之后转速发生小幅度震荡，后来趋于稳定,运算时间和精度相差不大 2、闭环仿真：

在上述仿真基础上，添加转速闭环控制器，转速指令为1130rpm ，

0~2.5s ，电机空载，即I d =0；比例环节：Kp=4 2.5s~5s ，电机满载，即I d =55A 。比例环节：Kp=4

(1)控制器为比例环节：试取不同k p 值，画出转速波形，求稳态时n 和s 并进行比较。 在Simulink 中建立系统模型如下图所示

time time

s

静差率Kp=1

Kp=2时转速n 稳态值为5784rpm ，静差率s 稳态值为-4.0609

time n

转速Kp=2

time

s

静差率Kp=2

Kp=4时转速n 稳态值为5827rpm ，静差率s 稳态值为-4.0982

time n

转速Kp=4

time

s

静差率Kp=4

Kp=5时转速n 稳态值为5835rpm ，静差率s 稳态值为-4.1054

time n

转速Kp=5

time

s

静差率Kp=5

记录数据如下表所示

结论:由图像可得，Kp 越大，转速n 超调量越大

(2)控制器为比例积分环节，设计恰当的k p 和k I 值，并与其它不同的k p 和k I 值比较，画出不同控制参数下的转速波形，比较静差率、超调量、响应时间和抗扰性。K I =14.83

Kp overshoot risetime 抗干扰性 Kp=0.62，46.3235 0.0050 0.0064

time

n

转速Kp=0.62 Ki=14.83

time

s

静差率Kp=0.62 Ki=14.83

Kp overshoot risetime 抗干扰性 Kp=0.7, 42.1429 0.0486 0.0112

time

n

转速Kp=0.7 Ki=14.83

time

s

静差率Kp=0.7 Ki=14.83

图2 转速闭环直流电机调速控制框图

控制K I =14.83不变，调节Kp

3、分析

结合《自动控制系统》相关知识，对上述结果进行分析。

PID控制器参数整定的方法很多，概括起来有两大类：

一是理论计算整定法。它主要是依据系统的数学模型，经过理论计算确定控制器参数。这种方法所得到的计算数据未必可以直接用，还必须通过工程实际进行调整和修改。

二是工程整定方法，它主要依赖工程经验，直接在控制系统的试验中进行，且方法简单、易于掌握，在工程实际中被广泛采用。

PID控制器参数的工程整定方法，主要有临界比例法、反应曲线法和衰减法。现在一般采用的是临界比例法。

利用该方法进行PID控制器参数的整定步骤如下：(1)首先预选择一个足够短的采样周期让系统工作；(2)仅加入比例控制环节，直到系统对输入的阶跃响应出现临界振荡，记下这时的比例放大系数和临界振荡周期；(3)在一定的控制度下通过公式计算得到PID控制器的参数。

PID参数的设定：是靠经验及工艺的熟悉，参考测量值跟踪与设定值曲线，从而调整P\I\D的大小。

常用口诀：

参数整定找最佳，从小到大顺序查先是比例后积分，最后再把微分加曲线振荡很频繁，比例度盘要放大曲线漂浮绕大湾，比例度盘往小扳曲线偏离回复慢，积分时间往下降曲线波动周期长，积分时间再加长曲线振荡频率快，先把微分降下来动差大来波动慢。微分时间应加长理想曲线两个波，前高后低4比1 一看二调多分析，调节质量不会低

直流调速系统的MATLAB仿真(参考程序)汇总.

直流调速系统的MATLAB 仿真 一、开环直流速系统的仿真 开环直流调速系统的电气原理如图1所示。直流电动机的电枢由三相晶闸管整流电路经平波电抗器L 供电，通过改变触发器移相控制信号c U 调节晶闸管的控制角α，从而改变整流器的输出电压，实现直流电动机的调速。该系统的仿真模型如图2所示。 图1 开环直流调速系统电气原理图 图2 直流开环调速系统的仿真模型 为了减小整流器谐波对同步信号的影响，宜设三相交流电源电感s 0L =，直流电动机励磁由直流电源直接供电。触发器（6-Pulse ）的控制角（alpha_deg ）由移相控制信号c U 决定，移相特性的数学表达式为 min c cmax 9090U U αα?-=?-

在本模型中取min 30α=?，cmax 10V U =，所以c 906U α=-。在直流电动机的负载转矩输入端L T 用Step 模块设定加载时刻和加载转矩。 仿真算例1 已知一台四极直流电动机额定参数为N 220V U =，N 136A I =， N 1460r /min n =，a 0.2R =Ω，2222.5N m GD =?。励磁电压f 220V U =，励磁电流f 1.5A I =。采用三相桥式整流电路，设整流器内阻rec 0.3R =Ω。平波电抗器 d 20mH L =。仿真该晶闸管-直流电动机开环调速系统，观察电动机在全压起动和起动后加额定负载时的电机转速n 、电磁转矩 e T 、电枢电流d i 及电枢电压d u 的变化情况。N 220V U = 仿真步骤： 1）绘制系统的仿真模型（图2）。 2）设置模块参数（表1） ① 供电电源电压 N rec N 2min 2200.3136 130(V)2.34cos 2.34cos30U R I U α++?= =≈?? ② 电动机参数 励磁电阻： f f f 220146.7()1.5 U R I = ==Ω 励磁电感在恒定磁场控制时可取“0”。 电枢电阻： a 0.2R =Ω 电枢电感由下式估算： N a N N 0.422019.1 19.10.0021(H)2221460136 CU L pn I ?==?≈??? 电枢绕组和励磁绕组间的互感af L ： N a N e N 2200.2136 0.132(V min/r)1460 U R I K n --?= =≈?

信号与系统的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仿真 关键词: Matlab 矩形波导时域有限差分法 摘要：介绍了时域有限差分法的基本原理，并利用Matlab仿真，对矩形波导谐振腔中的电磁场作了模拟和分析。 关键词：时域有限差分法；Matlab；矩形波导；谐振腔 目前，电磁场的时域计算方法越来越引人注目。时域有限差分（Finite Difference Time Domain，FDTD）法［1］作为一种主要的电磁场时域计算方法，最早是在1966年由K. S. Yee提出的。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解，通过建立时间离散的递进序列，在相互交织的网格空间中交替计算电场和磁场。经过三十多年的发展，这种方法已经广泛应用到各种电磁问题的分析之中。 Matlab作为一种工程仿真工具得到了广泛应用［2］。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。 下面将采用FDTD法，利用Matlab仿真来分析矩形波导谐振腔的电磁场，说明了将二者结合起来的优越性。 1FDTD法基本原理 时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。FDTD 空间网格单元上电场和磁场各分量的分布如图1所示。 电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元每个面的中心，每个磁场（电场）分量都有4个电场（磁场）分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足，而且

还允许旋度方程在空间上进行中心差分运算，同时也满足了法拉第电磁感应定律和安培环路积分定律，也可以很恰当地模拟电磁波的实际传播过程。 1.1Maxwell方程的差分形式 旋度方程为： 将其标量化，并将问题空间沿3个轴向分成若干网格单元，用Δx，Δy和Δz 分别表示每个网格单元沿3个轴向的长度，用Δt表示时间步长。网格单元顶点的坐标（x，y，z）可记为： 其中：i，j，k和n为整数。 同时利用二阶精度的中心有限差分式来表示函数对空间和时间的偏导数，即可得到如下FDTD基本差分式： 由于方程式里出现了半个网格和半个时间步，为了便于编程，将上面的差分式改写成如下形式：

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仿真实验报告

实验一 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.

各种BP学习算法MATLAB仿真

3.3.2 各种BP学习算法MATLAB仿真 根据上面一节对BP神经网络的MATLAB设计，可以得出下面的通用的MATLAB程序段，由于各种BP学习算法采用了不同的学习函数，所以只需要更改学习函数即可。 MATLAB程序段如下： x=-4:0.01:4; y1=sin((1/2)*pi*x)+sin(pi*x); %trainlm函数可以选择替换 net=newff(minmax(x),[1,15,1],{'tansig','tansig','purelin'},'trainlm'); net.trainparam.epochs=2000; net.trainparam.goal=0.00001; net=train(net,x,y1); y2=sim(net,x); err=y2-y1; res=norm(err); %暂停，按任意键继续 Pause %绘图，原图（蓝色光滑线）和仿真效果图（红色+号点线） plot(x,y1); hold on plot(x,y2,'r+'); 注意：由于各种不确定因素，可能对网络训练有不同程度的影响，产生不同的效果。如图3-8。 标准BP算法（traingd）

图3-8 标准BP算法的训练过程以及结果（原图蓝色线，仿真图+号线）增加动量法(traingdm) 如图3-9。 图3-9 增加动量法的训练过程以及结果（原图蓝色线，仿真图+号线）弹性BP算法（trainrp）如图3-10 图3-10 弹性BP算法的训练过程以及结果（原图蓝色线，仿真图+号线）

动量及自适应学习速率法（traingdx）如图3-11。 图3-11 动量及自适应学习速率法的训练过程以及结果（原图蓝色线，仿真图+号线）共轭梯度法(traincgf)如图3-12。

直流调速系统的MATLAB仿真参考程序汇总

直流调速系统的MA TLAB仿真 一、开环直流速系统的仿真 开环直流调速系统的电气原理如图1所示。直流电动机的电枢由三相晶闸管整流电路经平波电抗器供电，通过改变触发器移相控制信号调节晶闸管的U L c ，从 而改变整流器的输出电压，实现直流电动机的调速。该系统的仿真控制角模型如图2所示。 L + + GT UCR E d- - 开环直流调速系统电气原理图图1 图2 直流开环调速系统的仿真模型

L?0，直为了减小整流器谐波对同步信号的影响，宜设三相交流电源电感s流电动机励磁由直流电源直接供电。触发器（6-Pulse）的控制角（alpha_deg）由U 决定，移相特性的数学表达式为移相控制信号c???90?min U?90??c U cmax 1 ??。在直流电动机的负载，所以，在本模型中取U?10V6??30?90U?ccmaxmin转矩输 入端用Step模块设定加载时刻和加载转矩。T L仿真算例1已知一台四极直流电动机额定参数为，，136AIU?220V?NN22。励磁电压，励磁电，， 220VminUR?0.2??1460rn?/m?22.5NGD?fNa流。采用三相桥式整流电路，设整流器内阻。平波电抗器??1.5A0.3RI?recf。仿真该晶闸管-直流电动机开环调速系统，观察电动机在全压起动20mHL?d n、电磁转矩、电枢电流和起动后加额定负载时的电机转速及电枢电压的uTi ded变化情况。220V?U N仿真步骤： 1）绘制系统的仿真模型（图2）。 2）设置模块参数（表1） ①供电电源电压 U?RI220?0.3?136NNrec130(V)U??? 2?2.34?cos302.34cos?min②电动机参数 励磁电阻： U220f)146.7(?R???f I1.5f励磁电感在恒定磁场控制时可取“0”。 电枢电阻： ?0.2R?a电枢电感由下式估算： CU0.4?220N?19.1?L?19.1?0.0021(H) a2pnI2?2?1460?136NN L：电枢绕组和励磁绕组间的互感 af U?RI220?0.2?136NNa?K?0.132(V?min/r)?e n1460N 2 6060K??0.132?K?1.26eTπ2π2K1.26T0.84(H)??L?af1.5I f电机转动惯量 222.5GD2 )??0.57(kg?mJ?9.814?4g 额定负载转矩③ 模块参数名参数

增量式PID控制算法的MATLAB仿真

增量式PID 控制算法的MATLAB 仿真 PID 控制的原理 在工程实际中，应用最为广泛的调节器控制规律为比例、积分、微分控制，简称PID 控制，又称PID 调节。PID 控制器问世至今已有近70年历史，它以其结构简单、稳定性好、工作可靠、调整方便而成为工业控制的主要技术之一。当被控对象的结构和参数不能完全掌握，或得不到精确的数学模型时，控制理论的其它技术难以采用时，系统控制器的结构和参数必须依靠经验和现场调试来确定，这时应用PID 控制技术最为方便。即当我们不完全了解一个系统和被控对象,或不能通过有效的测量手段来获得系统参数时，最适合用PID 控制技术。PID 控制，实际中也有PI 和PD 控制。PID 控制器就是根据系统的误差，利用比例、积分、微分计算出控制量进行控制的。 一、 题目：用增量式PID 控制传递函数为G(s)的被控对象 G （s ）=5/(s^2+2s+10), 用增量式PID 控制算法编写仿真程序（输入分别为单位阶跃、正弦信号，采样时间为1ms ，控制器输出限幅：[-5,5],仿真曲线包括系统输出及误差曲线，并加上注释、图例）。程序如下 二、 增量式PID 原理 { U(k)= ?u(k)+ U(k-1) 或 { U(k)= ?u(k)+ U(k-1) 注：U(k)才是PID 控制器的输出 三、 分析过程 1、对G(s)进行离散化即进行Z 变换得到Z 传递函数G(Z)； 2、分子分母除以z 的最高次数即除以z 的最高次得到； )]}2()1(2)([)()]1()({[)(-+--++ --=?n n n T T n T T n n K n U D I P O εεεεεε)] 2()1(2)([)(i )]1()([)(-+--++--=?n n n Kd n K n n K n U P O εεεεεε

单闭环直流调速系统的MATLAB计算与仿真

1 绪论 直流电动机具有良好的起、制动性能，宜于在广泛范围内平滑调速，在轧钢机、矿井卷扬机、挖掘机、海洋钻机、金属切割机床、造纸机、高层电梯等需要高性能可控电力拖动的领域中得到了广泛的应用。近年来直流调速系统发展很快，然而直流拖动控制系统毕竟在理论上和实践上都比较成熟，而且从反馈闭环控制的角度来看，它又是交流拖动控制系统的基础，所以首先应该很好的掌握直流系统。我们可以首先从单闭环转速负反馈直流调速系统来研究。由于系统需要观察较多的性能，计算参数较多，而MATLAB中的Simulink实用工具可直接构建其动态模型，省去大量的计算，通过修改动态模型可完善系统性能。 直流调速系统概述 从生产机械要求控制的物理量来看，电力传动自动控制系统有调速系统、位置伺服系统、张力控制系统等其他多种类型，各种系统往往是通过控制转速来实现的，因此调速系统是最基本的驱动控制系统。调速系统目前分为交流和直流调速控制系统，由于直流调速系统的调速范围广，静差率小、稳定性好并且具有良好的动态性能。因此在相当长的时期内，高性能的调速系统几乎都采用了直流调速系统。相比于交流调速系统，直流调速系统在理论上和实践上更加成熟。 直流调速是现代电力拖动自动控制系统中发展较早的自动控制系统。在20世纪60年代发展起来的电力电子技术，使电能可以转换和控制，产生了现代各种高效、节能的新型电源和交直流调速装置，为工业生产，交通运输，建筑、办公、家庭自动化控制设备提供了现代化的高新技术，提高了生产效率和人们的生活质量，因此，人类社会的生产、生活发生了巨大变化。随着新型电力电子器件的研究和开发，先进控制技术的发展，电力电子和电力传动控制装置的性能也不断优化和提高，这一变化的影响将越来越大。 单闭环直流电机调速系统在现代日常生活中的应用越来越广泛，其良好的调速性能、低廉的价格越来越被大众接受。 单闭环直流调速系统由整流变压器、平波电抗器、晶闸管整流调速装置、电动机-发电机、闭环控制系统组成。我们可以通过调节晶闸管的控制角来调节转速，非常方便，高效。

MATLAB控制系统与仿真设计

MATLAB控制系统与仿真 课 程 设 计 报 告 院（系）：电气与控制工程学院 专业班级：测控技术与仪器1301班 姓名：吴凯 学号：1306070127

指导教师：杨洁昝宏洋 基于MATLAB的PID恒温控制器 本论文以温度控制系统为研究对象设计一个PID控制器。PID控制是迄今为止最通用的控制方法，大多数反馈回路用该方法或其较小的变形来控制。PID控制器（亦称调节器）及其改进型因此成为工业过程控制中最常见的控制器(至今在全世界过程控制中用的84%仍是纯PID调节器，若改进型包含在内则超过90%)。在PID控制器的设计中，参数整定是最为重要的,随着计算机技术的迅速发展，对PID参数的整定大多借助于一些先进的软件，例如目前得到广泛应用的MATLAB仿真系统。本设计就是借助此软件主要运用Relay-feedback法，线上综合法和系统辨识法来研究PID控制器的设计方法，设计一个温控系统的PID控制器，并通过MATLAB中的虚拟示波器观察系统完善后在阶跃信号下的输出波形。 关键词：PID参数整定；PID控制器；MATLAB仿真。 Design of PID Controller based on MATLAB Abstract This paper regards temperature control system as the research object to design a pid controller. Pid control is the most common control method up until now; the great majority feedback loop is controlled by this method or its small deformation. Pid controller (claim regulator also) and its second generation so become the most common controllers in the industry process control (so far, about 84% of the controller being used is the pure pid controller, it’ll exceed 90% if the second generation included). Pid parameter setting is most important in pid controller designing, and with the rapid development of the computer technology, it mostly recurs to some advanced software, for example, mat lab simulation software widely used now. this design is to apply that soft mainly use Relay feedback law and synthetic method on the line to study pid

内点法matlab仿真doc资料

编程方式实现： 1.惩罚函数 function f=fun(x,r) f=x(1,1)^2+x(2,1)^2-r*log(x(1,1)-1); 2.步长的函数 function f=fh(x0,h,s,r) %h为步长 %s为方向 %r为惩罚因子 x1=x0+h*s; f=fun(x1,r); 3. 步长寻优函数 function h=fsearchh(x0,r,s) %利用进退法确定高低高区间，利用黄金分割法进行求解h1=0;%步长的初始点 st=0.001; %步长的步长 h2=h1+st; f1=fh(x0,h1,s,r); f2=fh(x0,h2,s,r); if f1>f2 h3=h2+st; f3=fh(x0,h3,s,r); while f2>f3 h1=h2; h2=h3; h3=h3+st; f2=f3; f3=fh(x0,h3,s,r); end else st=-st; v=h1; h1=h2; h2=v; v=f1; f1=f2; f2=v; h3=h2+st; f3=fh(x0,h3,s,r); while f2>f3 h1=h2; h2=h3; h3=h3+st; f2=f3;

f3=fh(x0,h3,s,r); end end %得到高低高的区间 a=min(h1,h3); b=max(h1,h3); %利用黄金分割点法进行求解 h1=1+0.382*(b-a); h2=1+0.618*(b-a); f1=fh(x0,h1,s,r); f2=fh(x0,h2,s,r); while abs(a-b)>0.0001 if f1>f2 a=h1; h1=h2; f1=f2; h2=a+0.618*(b-a); f2=fh(x0,h2,s,r); else b=h2; h2=h1; f2=f1; h1=a+0.382*(b-a); f1=fh(x0,h1,s,r); end end h=0.5*(a+b); 4. 迭代点的寻优函数 function f=fsearchx(x0,r,epson) x00=x0; m=length(x0); s=zeros(m,1); for i=1:m s(i)=1; h=fsearchh(x0,r,s); x1=x0+h*s; s(i)=0; x0=x1; end while norm(x1-x00)>epson x00=x1; for i=1:m s(i)=1; h=fsearchh(x0,r,s);

09级系统仿真与MATLAB语言实验

系统仿真与MATLAB语言 实验指导书

对参加实验学生的总要求 1、认真复习有关理论知识，明确每次实验目的，了解实验相关软件操作，熟悉实验内容和方法。 2、实验过程中注意仔细观察，认真记录有关数据和图像，并经由指导教师查验后方可结束实验。 3、应严格遵守实验室规章制度，服从实验室教师的安排和管理。 4、对实验仪器的操作使用严格按照实验室要求进行。

实验总要求 1、封面：注明实验名称、实验人员班级、学号（全号）和姓名等。 2、内容方面：注明实验所用设备、仪器及实验步骤方法；记录清楚实验所得的原始数据和图像，并按实验要求绘制相关图表、曲线或计算相关数据；认真分析所得实验结果，得出明确实验结论。并注明该结论所依据的原理和理论；对实验进行反馈回顾，总结出实验方法要领和注意事项，对实验失败的原因进行分析剖解，总结出实验的经验和教训。 3、文字方面，撰写规范，杜绝错别字。 4、杜绝抄袭，杜绝提供不真实的实验内容。

实验一 MATLAB 语言工作环境和基本操作 1 实验目的 1）．熟悉MATLAB 的开发环境； 2）．掌握MATLAB 的一些常用命令； 3）．掌握矩阵、变量、表达式的输入方法及各种基本运算。 2 实验器材 计算机WinXP 、Matlab7.0软件 3 实验内容 (1). 输入 A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]， 在命令窗口中执行下列表达式，掌握其含义： A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A (2).输入 C=1:2:20，则 C （i ）表示什么？其中 i=1,2,3,…,10； （3）掌握MA TLAB 常用命令 >> who %列出工作空间中变量 >> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量 >>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称 如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。 （4） 在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ (5). 试用 help 命令理解下面程序各指令的含义： clear t =0:0.001:2*pi; subplot(2,2,1); polar(t, 1+cos(t)) subplot(2,2,2); plot(cos(t).^3,sin(t).^3) subplot(2,2,3); polar(t,abs(sin(t).*cos(t))) subplot(2,2,4); polar(t,(cos(2*t)).^0.5) （6）（选做）设计M 文件计算: x=0:0.1:10 当sum>1000时停止运算，并显示求和结果及计算次数。 i i i x x sum 2100 2 -= ∑ =

某温度控制系统的MATLAB仿真

课程设计报告 题目某温度控制系统的MATLAB仿真（题目C）

过程控制课程设计任务书 题目C ：某温度控制系统的MATLAB 仿真 一、 系统概况： 设某温度控制系统方块图如图： 图中G c (s)、G v (s)、G o (s)、G m (s)、分别为调节器、执行器、过程对象及温度变送器的传递函数；，且电动温度变送器测量范围（量程）为50～100O C 、输出信号为4～20mA 。G f (s)为干扰通道的传递函数。 二、系统参数 二、 要求： 1、分别建立仿真结构图，进行以下仿真，并求出主要性能指标： (1)控制器为比例控制，其比例度分别为δ＝10％、20％、50％、100％、200％时，系统广义对象输出z(t)的过渡过程； (2)控制器为比例积分控制，其比例度δ＝20％，积分时间分别为T I ＝1min 、3min 、5min 、10min 时，z(t)的过渡过程； 0m v o 0f o o =5min =2.5min =1.5(kg/min)/mA =5.4C/(kg/min) =0.8 C C T T K K K x(t)=80f(t)=10； ；；； ；给定值； 阶跃扰动

(3)控制器为比例积分微分控制，其比例度δ＝10％，积分时间T I＝5min，微分时间T D = 0.2min时，z(t)的过渡过程。 2、对以上仿真结果进行分析比对，得出结论。 3、撰写设计报告。 注：调节器比例带δ的说明 比例控制规律的输出p(t)与输入偏差信号e(t)之间的关系为 式中，K c叫作控制器的比例系数。 在过程控制仪表中，一般用比例度δ来表示比例控制作用的强弱。比例度δ定义为 式中，(z max-z min)为控制器输入信号的变化范围，即量程；(p max-p min)为控制器输出信号的变化范围。 = c p(t)K e(t) max min ( ) =100% ) max min e z z p(p-p δ - ?

PID控制算法的matlab仿真

PID 控制算法的matlab 仿真 PID 控制算法就是实际工业控制中应用最为广泛的控制算法,它具有控制器设计简单,控制效果好等优点。PID 控制器参数的设置就是否合适对其控制效果具有很大的影响,在本课程设计中一具有较大惯性时间常数与纯滞后的一阶惯性环节作为被控对象的模型对PID 控制算法进行研究。被控对象的传递函数如下: ()1d s f Ke G s T s τ-= + 其中各参数分别为30,630,60f d K T τ===。MATLAB 仿真框图如图1所示。 图1 2 具体内容及实现功能 2、1 PID 参数整定 PID 控制器的控制参数对其控制效果起着决定性的作用,合理设置控制参数就是取得较好的控制效果的先决条件。常用的PID 参数整定方法有理论整定法与实验整定法两类,其中常用的实验整定法由扩充临界比例度法、试凑法等。在此处选用扩充临界比例度法对PID 进行整定,其过程如下: 1) 选择采样周期 由于被控对象中含有纯滞后,且其滞后时间常数为 60d τ=,故可选择采样周期1s T =。 2) 令积分时间常数i T =∞,微分时间常数0d T =,从小到大调节比例系数K , 使得系统发生等幅震荡,记下此时的比例系数k K 与振荡周期k T 。 3) 选择控制度为 1.05Q =,按下面公式计算各参数:

0.630.490.140.014p k i k d k s k K K T T T T T T ==== 通过仿真可得在1s T =时,0.567,233k k K T ==,故可得: 0.357,114.17,32.62, 3.262p i d s K T T T ==== 0.0053.57 p s i i p d d s K T K T K T K T === = 按此组控制参数得到的系统阶跃响应曲线如图2所示。 01002003004005006007008009001000 0.20.40.60.811.21.41.6 1.8 图2 由响应曲线可知,此时系统虽然稳定,但就是暂态性能较差,超调量过大,且响应曲线不平滑。根据以下原则对控制器参数进行调整以改善系统的暂态过程: 1) 通过减小采样周期,使响应曲线平滑。 2) 减小采样周期后,通过增大积分时间常数来保证系统稳定。 3) 减小比例系数与微分时间常数,以减小系统的超调。 改变控制器参数后得到系统的阶跃响应曲线如图3所示,系统的暂态性能得到明显改善、

自控-二阶系统Matlab仿真

自动控制原理 二阶系统性能分析Matlab 仿真大作业附题目+ 完整报告内容

设二阶控制系统如图 1所示，其中开环传递函数 ) 1(10 )2()(2+=+=s s s s s G n n ξωω 图1 图2 图3 要求： 1、分别用如图2和图3所示的测速反馈控制和比例微分控制两种方式改善系统的性能，如果要求改善后系统的阻尼比ξ =0.707，则t K 和 d T 分别取多少？ 解： 由)1(10 )2()(2 += +=s s s s s G n n ξωω得10 21,10,102===ξωωn 22n n () s s ωξω+R (s )C (s ) -

对于测速反馈控制，其开环传递函数为：) 2()s (2 2n t n n K s s G ωξωω++=； 闭环传递函数为：2 2 2)2 1(2)(n n n t n s K s s ωωωξωφ+++= ； 所以当n t K ωξ2 1+=0.707时，347.02)707.0(t =÷?-=n K ωξ； 对于比例微分控制，其开环传递函数为：)2()1()(2 n n d s s s T s G ξωω++=； 闭环传递函数为：) )2 1(2)1()(2 22 n n n d n d s T s s T s ωωωξωφ++++=； 所以当n d T ωξ2 1 +=0.707时，347.02)707.0(=÷?-=n d T ωξ； 2、请用MATLAB 分别画出第1小题中的3个系统对单位阶跃输入的响应图； 解： ①图一的闭环传递函数为： 2 22 2)(n n n s s s ωξωωφ++=，10 21 ,10n ==ξω Matlab 代码如下： clc clear wn=sqrt(10); zeta=1/(2*sqrt(10)); t=0:0.1:12; Gs=tf(wn^2,[1,2*zeta*wn,wn^2]); step(Gs,t)

复杂过程控制系统设计与Simulink仿真

银河航空航天大学 课程设计 （论文） 题目复杂过程控制系统设计与Simulink仿 真 班级 学号 学生姓名 指导教师

目录 0. 前言 (1) 1. 总体方案设计 (2) 2. 三种系统结构和原理 (3) 2.1 串级控制系统 (3) 2.2 前馈控制系统 (3) 2.3 解耦控制系统 (4) 3. 建立Simulink模型 (5) 3.1 串级 (5) 3.2 前馈 (5) 3.3 解耦 (7) 4. 课设小结及进一步思想 (15) 参考文献 (15) 附录设备清单 (16)

复杂过程控制系统设计与Simulink仿真 姬晓龙银河航空航天大学自动化分校 摘要：本文主要针对串级、前馈、解耦三种复杂过程控制系统进行设计，以此来深化对复杂过程控制系统的理解，体会复杂过程控制系统在工业生产中对提高产品产量、质量和生产效率的重要作用。建立Simulink模型，学习在工业过程中进行系统分析和参数整定的方法，为毕业设计对模型进行仿真分析及过程参数整定做准备。 关键字：串级；前馈；解耦；建模；Simulink。 0.前言 单回路控制系统解决了工业过程自动化中的大量的参数定制控制问题，在大多数情况下这种简单系统能满足生产工艺的要求。但随着现代工业生产过程的发展，对产品的产量、质量，对提高生产效率、降耗节能以及环境保护提出了更高的要求，这便使工业生产过程对操作条件要求更加严格、对工艺参数要求更加苛刻，从而对控制系统的精度和功能要求更高。为此，需要在单回路的基础上，采取其它措施，组成比单回路系统“复杂”一些的控制系统，如串级控制（双闭环控制）、前馈控制大滞后系统控制（补偿控制）、比值控制（特殊的多变量控制）、分程与选择控制（非线性切换控制）、多变量解耦控制（多输入多输出解耦控制）等等。从结构上看，这些控制系统由两个以上的回路构成，相比单回路系统要多一个以上的测量变送器或调节器，以便完成复杂的或特殊的控制任务。这类控制系统就称为“复杂过程控制系统”，以区别于单回路系统这样简单的过程控制系统。 计算机仿真是在计算机上建立仿真模型，模拟实际系统随时间变化的过程。通过对过程仿真的分析，得到被仿真系统的动态特性。过程控制系统计算机仿真，为流程工业控制系统的分析、设计、控制、优化和决策提供了依据。同时作为对先进控制策略的一种检验，仿真研究也是必不可少的步骤。控制系统的计算机仿真是一门涉及到控制理论、计算机数学与计算机技术的综合性学科。控制系统仿真是以控制系统的模型为基础，主要用数学模型代替实际控制系统，以计算机为工具，对控制系统进行实验和研究的一种方法。在进行计算机仿真时，十分耗费时间与精力的是编制与修改仿真程序。随着系统规模的越来越大，先进过程控制的出现，就需要行的功能强大的仿真平台Math Works公司为MATLAB提供了控制系统模型图形输入与仿真工具Simulink，这为过程控制系统设计与参数整定的计算与仿真提供了一个强有力的工具，使过程控制系统的设计与整定发生了革命性的变化。

实验一 典型环节的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仿真)

任务书 1.设计题目 转速、电流双闭环直流调速系统的设计 2.设计任务 某晶闸管供电的双闭环直流调速系统，整流装置采用三相桥式电路，基本数据为： 直流电动机：U n =440V，I n =365A，n N =950r/min，R a =, 电枢电路总电阻R=, 电枢电路总电感L=, 电流允许过载倍数=， 折算到电动机飞轮惯量GD2=20Nm2。 晶闸管整流装置放大倍数K s =40，滞后时间常数T s = 电流反馈系数=A (10V/ 转速反馈系数= min/r (10V/nN) 滤波时间常数取T oi =，T on = ===15V；调节器输入电阻R a =40k 3.设计要求 (1)稳态指标:无静差

(2)动态指标:电流超调量5%；采用转速微分负反馈使转速超调量等于0。 目录 任务书.......................................................... I 目录........................................................... II 前言 (1) 第一章双闭环直流调速系统的工作原理 (2) 双闭环直流调速系统的介绍 (2) 双闭环直流调速系统的组成 (3) 双闭环直流调速系统的稳态结构图和静特性 (4) 双闭环直流调速系统的数学模型 (5) 双闭环直流调速系统的动态数学模型 (5) 起动过程分析 (6) 第二章调节器的工程设计 (9) 调节器的设计原则 (9) Ⅰ型系统与Ⅱ型系统的性能比较 (10) 电流调节器的设计 (11) 结构框图的化简和结构的选择 (11) 时间常数的计算 (12) 选择电流调节器的结构 (13) 计算电流调节器的参数 (13) 校验近似条件 (14) 计算调节器的电阻和电容 (15) 转速调节器的设计 (15) 转速环结构框图的化简 (15)

LMMSE算法信道均衡MATLAB仿真

一．信道均衡的概念 实际的基带传输系统不可能完全满足无码间串扰传输条件，因而码间串扰是不可避免的。当串扰严重时，必须对系统的传输函数 进行校正，使其达到或接近无码间串扰要求的特性。理论和实践表明，在基带系统中插入一种可调滤波器就可以补偿整个系统的幅频，和相频特性从而减小码间串扰的影响这个对系统校正的过程称为均衡，实现均衡的滤波器称为均衡器。 均衡分为频域均衡和时域均衡。频域均衡是从频率响应考虑，使包括均衡器在内的整个系统的总传输函数满足无失真传输条件。而时域均衡，则是直接从时间响应考虑，使包括均衡器在内的整个系统的冲激响应满足无码间串扰条件。 频域均衡在信道特性不变，且传输低速率数据时是适用的，而时域均衡可以根据信道特性的变化进行调整，能够有效地减小码间串扰，故在高速数据传输中得以广泛应用。 时域均衡的实现方法有多种，但从实现的原理上看，大致可分为预置式自动均衡和自适应式自动均衡。预置式均衡是在实际传数之前先传输预先规定的测试脉冲（如重复频率很低的周期性的单脉冲波形），然后按“迫零调整原理”自动或手动调整抽头增益；自适应式均衡是在传数过程中连续测出距最佳调整值的误差电压，并据此电压去调整各抽头增益。一般地，自适应均衡不仅可以使调整精度提高，而且当信道特性随时间变化时又能有一定的自适应性，因此很受重视。这种均衡器过去实现起来比较复杂，但随着大规模、超大规模集成电路和微处理机的应用，其发展十分迅速。 二．信道均衡的应用 1.考虑如图所示的基带等效数据传输系统，发送信号k x 经过ISI 失真信道传输，叠加高斯加性噪声。 图1基带等效数据传输模型 设发送信号采用QPSK 调制，即(1)k x j =±±ISI 信道的冲击响应以向量的形式表示为h 2211[,,,]T L L L h h h --+=???。典型的ISI 信道响应向量有三种： h [0.04,0.05,0.07,0.21,0.5,0.72,0.36,0,0.21,0.03,0.07]T A =--- h [0.407,0.815,0.407]T B = h [0.227,0.46,0.6888,0.46,0.227]T C = k ω为实部与虚部独立的复高斯白噪声，其均值为零，方差为2 ωσ。 2.实现目的