大作业-----双闭环直流调速系统的MATLAB仿真

大作业 双闭环直流调速系统的MATLAB 仿真

（设计性实验）

一、实验目的

1、了解Simulink 下数学模型的仿真方法。

2、掌握数学模型的仿真建模方法和仿真参数设置要求

3、进一步掌握双闭环反馈控制系统的基本特性。 二、实验原理

双闭环控制系统的结构原理框图如下所示：

要求：在实验前需下根据双闭环控制原理计算出各环的PI 参数，写出正确的系统控制数学模型。

三、实验内容

1、直流电机双闭环控制系统SIMULINK 仿真模块建立。

2、仿真参数的设置。

3、各控制点的波形分析。

4、改变给定, 观察输出响应的变化。

5、用传递函数方程仿真。 四、实验步骤

1、按上述结构和参数建立仿真模型。

2、设定输入为单位阶跃信号，用scope 观察系统输出响应。

3、将改变给定，其余参数不变，观察系统输出响应。

4、突加负载，观察系统的抗扰性能。改变开环放大系数，观察系统的稳定性能。 5．用建立的传递函数方程仿真。 五、实验报告

按规定的实验报告要求写出实验报告，报告的内容有实验目的、建模名称，参数设置，实验电路，仿真模型结构图、仿真结果波形，结果分析。

名称

路径

参数设置 AC V oltage Source Simpowersystem/electrical/ Ground(output) Simpowersystem/ connectors Thyristor Simpowersystem/Power electronics

Series RLC branch Simpowersystem/Elenemts C=inf L=0纯电阻 Vlotage Measurement Simpowersystem/ Measurements Current Measurement Simpowersystem/ Measurements Multimeter Simpowersystem/ Measurements T Connector Simpowersystem/ connectors Scope

Simulink/sinks

α

1/C e

U *n

n U d0

U n

+

- ASR

1/R

T l s+1

R T m s

K s T s s+1

ACR

β

U

*

i

U i

- -

E

I d

±?I dL

负载扰

Mux Simulink/Signal Routing

Demux Simulink/Signal Routing

Pulse Generator Simulink/ Sources

step Simulink/sources ‘0’‘0’‘10’‘0’Integrator Simulink/ continuous/ ‘0’‘10’‘-10’Transfer Fcn Simulink/ continuous/

gain Simulink/ math operations

sum Simulink/ math operations +-

讨论：计算电流环设计成典型一型系统选用PI调节器

KT Kp τi ki=Kp/τi 0.25 0.5067 0.03 16.89

0.5 1.013 0.03 34.77

1 2.027 0.03 67.567

MATLAB大作业

选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c

(2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵：

（1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2

时域有限差分法的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仿真期末大作业

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

实验一 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大作业

MATLAB大作业 作业要求： （1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。 （2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电子稿要包含MATLAB 程序。 （3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。 （4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。 （5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。 第一类：绘制图形。（B级） 问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据 斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。 问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢 尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。

问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮 亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如https://www.docsj.com/doc/7213839873.html,/content/16/0103/14/5315_525141100.shtml），请分析构图原理并用MATLAB实现。 问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。 点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。 提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。 第二类：插值与拟合。（B级） 问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未 （1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。 （2）绘制插值图形并标注样本点。 问题二：估算矩形平板各个位置的温度。已知平板长为5m，宽为3m，平板上3×5栅格 点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。 （1）分别使用最近点插值、线性插值和三次样条插值进行计算。 （2）用杆图标注样本点。 （3）绘制平板温度分布图。 问题三：自行车道的设计。对9条道路上的自行车道宽度以及自行车与过往机动车之间 （1）对数据进行线性拟合。 （2）绘制拟合曲线和样本点。 （3）如果自行车与过往机动车之间安全距离的最小距离是1.8m，试计算相应的自行车道宽度的最小值。 问题四：在水资源工程学中，水库的大小与为了蓄水而拦截的河道中的水流速度密切相 关。对于某些河流来说，这种长时间的历史水流记录很难获得。然而通常容易得到过去若干年间关于降水量的气象资料。鉴于此，推导出流速与降水量之间的关系式往往特别有用。只

各种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仿真大作业报告

. 基于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 时，磁路饱和(为比较不同饱和电流对应的效果，饱和电流可以自己改变）。

增量式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大作业

Matlab 大作业 （组内成员：彭超杰、南彦东、江明伟） 一、研究模型 （电车）通过控制油门（保持一定角度）来调节电动机能输出稳定的转速，从而控制车速稳定。 数学依据说明如下： 由图可知存在以下关系：a d a a u w k R i dt di L =++ (w k e d d =) L M M dt dw J -= a m i k M = L a m M i k dt dw J -=

k为反电势常数，m k为电动机电磁力矩常数，这里忽略阻尼力矩。d

二、数学模型 再看整个研究对象，示意图以课本为依据，不同点是这里将数控的进给运动，转换为汽车行驶所需要的扭矩。（这里不说明扭矩的具体产生过程，仅仅说明输出车轮旋转的角速度w ） 对照课本不同，() s θ变为()s N ，1 221z z w w =，1w 为电动机的转速，2w 为轮胎的转速，1z 为电动机的光轴齿轮的齿数，2z 为与轮胎相连光轴的 齿轮齿数。 )(*10110w x w k x ==，1 21z z k = ()c a m m d b a m x K K K k s k k JRs JLs K K K k s G i 1231+++= () c a m m d M K K K k s k k JRs JLs R Ls K s G L 1231)(++++-= 同理，忽略电枢绕组的电感L ，简化系统传递函数方框图如下

()JR K K K k JR s k k s JR K K K k s G c a m m d b a m x i 121++= ()JR K K K k JR s k k s K K K K k s k k Rs R K s G c a m m d c a m m d M L 121121++-=++-=

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控制系统与仿真设计

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与数学实验大作业

《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日

一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。

（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：

内点法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);

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仿真（题目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仿真-电机特性曲线)

电机大作业 专业班级：电气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');

复杂过程控制系统设计与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，这为过程控制系统设计与参数整定的计算与仿真提供了一个强有力的工具，使过程控制系统的设计与整定发生了革命性的变化。