2014matlab设计报告陈志强

合集下载

2014年春 GUI

2014年春    GUI

1. GUIDE的启动
• 以上面任一种方式打开GUIDE后,即下图的启动界面。
GUIDE启动界面
2.GUI设计工作台
主要由4个功能区,六个实用工具组成 菜单条 工具条
对 齐 工 具
菜 单 编 辑
文 件 编 辑
对对 象象 属列 性表
运 行 界 面
设计工作区
工作区 缩放句柄
控件选 择区
3. 用户控件介绍
例4. 基于例3,通过互斥 按钮旋转图片。 左旋90度 右旋90度
竖直反转
左右反转
x=(handles.img); if get(handles.radiobutton1,'value') y=imrotate(x,-90); elseif get(handles.radiobutton2,'value') y=imrotate(x,90); elseif get(handles.radiobutton3,'value') if isrgb(x) for k=1:3 y(:,:,k)=flipud(x(:,:,k)); end else y=flipud(x); end elseif get(handles.radiobutton4,'value') if isrgb(x) for k=1:3 y(:,:,k)=fliplr(x(:,:,k)); end else y=fliplr(x); end else end; axes(handles.axes_result); imshow(y);
<labe l>
<labe l>
<labe l>
<labe l> <labe l> <labe l>

自动化综合设计——基于Matlab仿真实验实现系统校正环节的设计

自动化综合设计——基于Matlab仿真实验实现系统校正环节的设计

·自动化专业综合设计报告一、 设计目的:通过实验设计熟练掌握Matlab 与Simulink 的编程与仿真功能,提高自身的动手能力。

设单位负反馈系统的开环输出为:)12.0)(11.0()(++=S S S Ko s G 用相应的频率校正法对系统进行校正设计,是系统的性能指标达到:1. 相角裕度 r ≥45°2. 在单位斜坡输入下的稳态误差 e ss <0.053. 系统的w c <3rad/s二、 设计要求:要求:1. 手工进行计算,设计校正方法,选择合适的校正装置2. 利用 Matlab 编程实现设计,并提供仿真结果3. 利用Simulink 进行仿真实现单位阶跃响应三、 设计内容:设单位负反馈系统的开环输出为:)12.0)(11.0()(++=S S S Ko s G用相应的频率校正法对系统进行校正设计,是系统的性能指标达到:4. 相角裕度 r ≥45°5. 在单位斜坡输入下的稳态误差 e ss <0.056. 系统的w c <3rad/s要求:4. 手工计算,设计校正方法5. 利用 Matlab 编程实现设计,仿真结果6. 利用Simulink 进行仿真实现单位阶跃响应1、手工设计:解:(1)由系统在单位斜坡输入下的稳态误差ess<0.05可得:速度误差系数Kv=0lim→S it S*G(S)=Ko>05.01=20 故可取Ko=25;计算原系统的相角裕度)(w ϕ和截止频率Wc0:由)(w A =22)2.0(1)1.0(125w w w ++=1可得 Wc0=9;由w w w o2.0arctan 1.0arctan 90)(--=ϕ当Wc0=9时;o o w 9.12)(-=γ 显然此时超前校正已经不能使系统满足要求了,必须使用滞后校正先将系统的相角裕度满足要求。

相角裕度要求o45>γ,加上滞后装置o 6的估值,则要求原系统相角裕度为o o o w 51645)('=+=γ,即o o w w w 512.0arctan 1.0arctan 90)(=--=ϕ可解得Wc1=2.37;取Wc=2.3;由0lg 20)(=+b Wc L ;0lg 203.225lg20=+b ; 092.0253.2==b ; Wc bT 1.01=; T=47; 滞后环节为:s s Ts bTs s Gc 4713.4111)(++=++=相角裕度:oo o Wc Wc Wc Wc Wc 4511.472.0h 4b ?Wc=2.3<3rad/s均满足设计要求,校正之后的传递函数为:)2.01)(1.01)(471()3.41(25)(s s s s s s G ++++=2、Matlab 实验程序:num0=25; %初始化w1=3;r=45;ee=6;pm=r+ee;den0=conv([1 0],conv([0.1 1],[0.2 1])); %输入函数G=tf(num0,den0);[gm1,pm1,wcg1,wcp1]=margin(num0,den0); %原系统的相角裕度pm1,截止频率wcp1for w=wcp1:-0.01:0 %计算原系统中满足γ的截止频率wc 的值gamma=pi/2-atan(0.1*w)-atan(0.2*w);ga=gamma*180/pi;if (ga>pm)wc0=w;break ;endendif (wc0<w1) ; %限制wc 的值使其满足小于3的要求 elseif (wc0>=w1)wc0=w1;endx=wc0;for wc0=x:-0.01:0 %计算校正装置的参数b、T以及验算校正后系统的相角裕度if(wc0<5)b=1/num0*wc0;elseif(wc0<wcp1)b=0.1*wc0^2/num0;endT=10/b/wc0;num2=[b*T 1];den2=[T 1];Gc=tf(num2,den2);[num,den]=series(num0,den0,num2,den2);[gm2,pm2,wcg2,wcp2]=margin(num,den);if(pm2>r);break;endendw=logspace(-3,1);subplot(2,1,1);[mag,phase]=bode(num,den,w); %校正后的Gc*G[mag1,phase1]=bode(num0,den0,w); %初始传函:G[mag2,phase2]=bode(num2,den2,w); %滞后校正装置:Gcsemilogx(w,20*log10(mag),w,20*log10(mag1),'*',w,20*log10(mag2),'-' );ylabel('·幅值(dB)');title('原系统*G,滞后环节Gc,校正后-GGc')xlabel('频率(rad/s)');grid on;subplot(2,1,2);semilogx(w,phase,w,phase1,'*',w,phase2,'-',w,(w-180-w),':');ylabel('相角(度)');xlabel('频率(rad/s)');title('原系统*G,滞后环节Gc,校正后-GGc ')grid on;G %原系统的传递函数Gc %校正系统的传递函数GS=tf(num,den) %校正后系统的传递函数pm2 %校正后系统的相角裕度wcp2 %校正后系统的截止频率3、利用Simulink进行仿真:仿真系统的连接图:在输入端施加阶跃信号,通过Scope检测未加校正环节系统的阶跃响应曲线;通过Scope1来检测加入校正环节后系统的阶跃响应曲线。

MATLAB程序设计实践2014

MATLAB程序设计实践2014

《MATLAB程序设计实践》课程内容一、实践目标1、回顾复习《科学计算与MATLAB》课程内容和基本编程方法;2、编程实现某一科学计算的算法,并举一例应用之和利用MATLAB 自带函数实现之;3、编程解决科学计算和工程实际问题。

二、实践内容1. Matlab基础部分内容:完成一个从文件(data.txt)中读取数据,然后做出相关的图形,每位同学的要求一样。

2. Matlab科学计算部分内容:每位同学根据学号后两位数选取一个内容完成相关程序设计。

三、实践报告的递交1. 实践报告的内容两个实践题目都需提供科学计算的算法说明(流程图)、程序实现代码和最终的结果。

2. 实践报告的递交形式文件包:班级-学号-姓名;包括,M文件、考核题(流程图、程序、运行结果)word文档。

纸质版:考核题、流程图、程序、运行结果;A4纸双面打印、全班按学号顺序订收齐后交指导老师。

一、《MATLAB程序设计实践》Matlab基础班级:学号:姓名:表示多晶体材料织构的三维取向分布函数(f=f(φ1,φ,φ2))是一个非常复杂的函数,难以精确的用解析函数表达,通常采用离散空间函数值来表示取向分布函数,Data.txt是三维取向分布函数的一个实例。

由于数据量非常大,不便于分析,需要借助图形来分析。

请你编写一个matlab程序画出如下的几种图形来分析其取向分布特征:(1)用Slice函数给出其整体分布特征;(2)用pcolor或contour函数分别给出(φ2=0, 5, 10, 15, 20, 25, 30, 35 …90)切面上f分布情况(需要用到subplot函数);(3) 用plot函数给出沿α取向线(φ1=0~90,φ=45,φ2=0)的f 分布情况。

备注:data.txt 数据格式说明数据说明部分,与作图无关此方向表示f 随着φ1从0,5,10,15, 20 …到90的变化而变化此方向表示f 随着φ从0,5,10,15, 20 …到90的变化而变化表示以下数据为φ2=0的数据,即f (φ1,φ,0)二《MATLAB程序设计实践》科学计算(01)班级:学号:姓名:1、编程实现以下科学计算算法,并举一例应用之。

MATLAB程序设计2014春概要

MATLAB程序设计2014春概要

②使用示波器模块观察仿真输出
三个示波器模块
Scope
XYGraph
Display
例:用三种示波器观察正弦波信号的波形曲线与
波形信息。
③使用To Workspace模块将仿真输出信息返回 到MATLAB命令窗口
如果不用示波器直接观察结果,可将控制 系统仿真信息输入“To Workspace”中。该模 块会自动将数据输出到MATLAB命令窗口里, 用变量保存起来,可用绘图命令绘制图形。
主要特点:简单、准确、快捷、形象、灵活等
2、用Simulink建立系统模型
Simulink完全采用标准模块方框图的拷贝方 法来构造动态系统的结构图模型。
创建过程:
模块的查找与选择
模块的操作模块的模拷块贝的、连移接动和删除
模块标题名称的修改
模块内部参数的修改**
例题:绘制一个二阶系统进行阶跃响应仿真的
10节课 2节课
3、 时域分析法的MATLAB实现 4、 根轨迹分析法的MATLAB实现 5、 频域分析法的MATLAB实现 6、 MATLAB在控制系统校正中应用
3节课 2节课
1节课 3节课
7、 MATLAB在非线性系统分析中应用 8、 MATLAB软件在电子电路中的应用
1节课 2节课
5、如何学习?
主讲人:李绍武
课程说明:
1、学习这门课所必须的基础课
基本数学知识、自动控制原理知识、计算机编程 语言(比如C语言)
2、课程安排及考试方式
理论:24 考试方式:
实验:16
指定教材
3、这门课在控制学科中的地位和作用
MATLAB的控制工具箱,已覆盖了控制领域的 每个领域,每一个工具箱都是当今世界上该控制 领域最顶尖、最优秀的计算与仿真软件; MATLAB已经成为控制领域内最流行的、被广泛采 用的控制系统计算、仿真、辅助设计软件。

Matlab实验报告

Matlab实验报告

实验一:Matlab操作环境熟悉一、实验目的1.初步了解Matlab操作环境。

2.学习使用图形函数计算器命令funtool及其环境。

二、实验内容熟悉Matlab操作环境,认识命令窗口、内存工作区窗口、历史命令窗口;学会使用format命令调整命令窗口的数据显示格式;学会使用变量和矩阵的输入,并进行简单的计算;学会使用who和whos命令查看内存变量信息;学会使用图形函数计算器funtool,并进行下列计算:1.单函数运算操作。

求下列函数的符号导数(1)y=sin(x); (2) y=(1+x)^3*(2-x);求下列函数的符号积分(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1-x^2);(4)y=(x1)/(x+1)/(x+2)求反函数(1)y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2));代数式的化简(1)(x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2.函数与参数的运算操作。

从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6)y6=x^2/23.两个函数之间的操作求和(1)sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5乘积(1)exp(-x)*sin(x) (2) sin(x)*x商(1)sin(x)/cos(x); (2) x/(1+x^2); (3) 1/(x-1)/(x-2);求复合函数(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x实验二:MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。

MATLAB与系统仿真实验指导书(2014新模板)_65032(1)

MATLAB与系统仿真实验指导书(2014新模板)_65032(1)

《MATLAB与系统仿真》实验指导书西安电子科技大学机电工程学院2014 年9 月实验一MATLAB基本使用、矩阵操作及编程一、实验目的和要求熟悉MATLAB的界面和基本操作,掌握矩阵的建立方法及各种运算,熟练使用MATLAB编程语言编写程序。

1.熟悉MATLAB软件的界面和帮助系统。

2.掌握MA TLAB软件中关于矩阵建立、矩阵初等变换以及矩阵算术、关系、逻辑运算的各种命令。

3.掌握MA TLAB软件中M函数和M文件的编写,以及程序结构与控制,学会编写一般程序。

二、实验仪器计算机一台,注意正确开、关机及打开软件。

三、实验内容与步骤1.MA TLAB的基本使用(1)启动与退出双击MATLAB图标,进入MA TLAB命令窗口,观察各个窗口以及菜单;单击File菜单中的Exit,或使用MATLAB的Exit命令退出。

(2)变量和矩阵的建立>>a=25>>b=3-9i问题1.1:输入“>>a=25;”,结果有什么区别?>>c=[1,2,3,4;5,6,7,8]>>d=0:1:10>>e=linspace(1,11,6)问题1.2:体会以上输入方法有什么区别和联系?若A为在0~2 之间均匀分布的22个数据,B=(1.3,2.5,7.6,2,-3),C=(23,20,17,14,11,8,5,2),各用何种方法输入较简单?>>A=[1,2,3;4,5,6;7,8,9]问题1.3:输入“A=[1 2 3;4 5 6;7 8 9]”,结果相同吗?MATLAB中对变量名有什么规定?(3)基本命令>>save data a b c>>load data使用↑键和↓键查看或执行以前执行过的命令>>d=1223.58>>format short e %显示格式设定>>d>>format rat>>d>>who a b %工作空间管理>>whos>>workspace>>path>>clear %注意工作空间的变化问题1.4:执行命令“clc”和“clear”结果有何不同?>>help eig %帮助命令>>look for max>>demo2. 矩阵的操作(1)测试矩阵大小>>f=length(A)>>[n,m]=size(A)>>g=ndims(A)问题1.5:以上三个命令有什么区别?(2)矩阵元素的引用>>A (1,:)>>A ([1,3],:)>>A (2:3,1:2)问题1.6:输入“A (2,3)”,结果任何,输入“A (6)”,结果相同吗?>>A ([1,3],:)=A ([3,1],:)>>A (2,:)=4>>A (find(A==4))=0>>A (3,:)=[]>>reshape (A,3,2)问题1.7:如何将A的第1列和第2列互换?输入“reshape (A,2,4)”能否得到正确结果?>>A (4,5)=2>>[A(1:3,2:4),A(2:4,1:3);A,A(:,2)](3)矩阵变换>>diag(A,1)>>tril(A,2)>>triu(A,2)>>flipud(A)>>fliplr(A)>>rot90(A)问题1.8:输入“rot90(A,2)”和“rot90(A,-2)”结果有区别吗?(4)特殊矩阵的产生>>B=eye(5)>>C=ones(2,3)>>D=zeros(3,2)>>E=rand(2,5)>>randn(3,4)>>vander([2 3 9])问题1.9:产生一个在[10,20]内均匀分布的4阶随机矩阵。

MATLAB图像处理技术及用户界面设计毕业论文

MATLAB图像处理技术及用户界面设计毕业论文

2014届毕业生毕业设计说明书题目:Matlab图像处理技术及用户界面设计院系名称:信息科学与工程专业班级:电信1004班学生姓名:瑞学号:************指导教师:教师职称:副教授2014年6月3日摘要随着科学技术、信息技术、人类社会文明水平的不断发展和提高,数字图像处理技术越来越多地被应用在人类生活、生产、发展、延续的各个领域。

然而,随着对图像处理和图像分析要求的提高,现有的图像软件因其采购成本和功能的限制,往往不能满足专业人士的需要,这是就需要专业人士编写自己的图形处理软件。

MATLAB简单的编程、强大的功能使得人们可以少受时间和专业精力的限制来完成对图像处理软件的开发工作。

MATLAB GUI设计师首先利用MATLAB图形用户界面设计了该系统的静态界面。

当静态界面设计完成之后,GUI将自动生成.FIG和.M文件。

其中.FIG文件保存了关于静态窗口界面的所有对象的属性值,.M文件包括GUI设计、控制函数以及为子函数的用户控件回调函数,主要用于控制GUI展开时的各种特征。

这个.M文件可以分为GUI初始化和回调函数两个部分,用户控件的回调函数根据用户与GUI的具体交互方式分别调用。

回调函数就是在调用对象时,该对象所要回应的动作。

因此,如何编写对象的回调函数是该系统的一大难点。

在为编写回调函数时,获得该函数的句柄是实现对象动作功能的关键所在。

句柄实际上就是分配给每个对象的数字标识,每次创建对象时,MATLAB就会自动为它创建一个唯一的句柄,这样只要我们能找到该句柄,就能对该对象进行操作。

本系统支持多种格式的图像文件读操作、写操作、显示操作。

并且包含一些图像处理功能,例如图像增强模块的灰度增强、直方图增强,butterworth滤波、中值滤波等;图像几何变换的图像剪切、旋转、缩放功能;编辑实现了对图像添加噪声的功能;其他的功能则实现了我感兴趣的图像亮度调整、底片处理效果、图像翻转、图像对比度调整效果。

2014《MATLAB及应用》实验报告1

2014《MATLAB及应用》实验报告1

电气工程学院实验报告实验项目名称MATLAB数值计算所属课程名称MATLAB及应用实验类型上机实验实验日期2014-03- 5指导教师盛义发班级XXXXXX学号2012XXXXXXXXX 姓名XXXXXX成绩一、实验名称MATLAB数值计算二、实验目的(1)掌握MATLAB变量的使用(2)掌握MATLAB数组的创建(3)掌握MATLAB数组和矩阵的运算(4)熟悉MATLAB多项式的运用三、实验原理1. 矩阵分析矩阵转置:单引号(’)矩阵的旋转:rot90(A,k),功能是将矩阵A旋转90度的k倍,缺省值是1 矩阵的左右翻转:fliplr(A)矩阵的上下翻转:flipud(A)矩阵的逆:inv(A),与A^(-1)等价矩阵的行列式:det(A)矩阵的秩:rank(A)矩阵的迹:trace(A)将矩阵化为最简式:rref(A)矩阵的特征值与特征向量:(1) E = eig(A);矩阵A的所有特征值构成向量E;(2) [V,D]=eig(A);A的所有特征值构成对角阵D,A的特征向量构成V的列向量;2. 多项式多项式的建立:若多的项的全部根构成的向量为X,则以X为根的多项式为poly(X)多项式的根:roots(p)计算以向量p为系数的多项式的根,包括重根,复根多项式求值:polyval(p,x),p是多项式的系数,x可以是一个数也可以是一个矩阵多项式求拟合次数:polyfit(x,y,n),x可以是一个数也可以是一个矩阵,y 是x对应的数或矩阵多项式的四则运算:(1)P1+P2;(2)P1-P2;(3)conv(P1,P2), (4)deconv(P1,P2)四、实验内容1. 已知矩阵11 12 13 1421 22 23 24A= 31 32 33 344142 43 44(1) A(:,1) (2) A(2,:) (3) A(:,2:3)(4) A(2:3,2:3) (5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:) (9) ones(2,2)(10) eye(2) (11) [A,[ones(2,2);eye(2)]] (12) diag(A)(13) diag(A,1) (14) diag(A,-1) (15) diag(A,2)2. 已知2354127354217623A-⎡⎤⎢⎥--⎢⎥=⎢⎥-⎢⎥--⎣⎦.求(1)A的逆;(2)A的行列式;(3)A的迹;(4)A的所有特征向量和特征值。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《Matlab仿真技术》设计报告题目倒立摆PID控制及其Matlab仿真专业班级电气工程及其自动化1321 学号 201350712103 学生姓名陈志强指导教师薛鹏学院名称电气信息工程学院完成日期: 2014 年 5 月 25日倒立摆PID控制及其Matlab仿真Inverted Pendulum PID Control and ItsMatlab Simulation摘要倒立摆系统是一个典型的快速、多变量、非线性、不稳定系统,对倒立摆的控制研究无论在理论上和方法上都有深远的意义。

本论文以实验室原有的直线一级倒立摆实验装置为平台,重点研究其PID 控制方法,设计出相应的PID控制器,并将控制过程在MATLAB上加以仿真。

本文主要研究内容是:首先概述自动控制的发展和倒立摆系统研究的现状;介绍倒立摆系统硬件组成,对单级倒立摆模型进行建模,并分析其稳定性;研究倒立摆系统的几种控制策略,分别设计了相应的控制器,以MATLAB为基础,做了大量的仿真研究,比较了各种控制方法的效果;借助固高科技MATLAB实时控制软件实验平台;利用设计的控制方法对单级倒立摆系统进行实时控制,通过在线调整参数和突加干扰等,研究其实时性和抗千扰等性能;对本论文进行总结,对下一步研究作一些展望。

关键词:倒立摆;PID控制器;MATLAB仿真设计报告正文1. 简述一级倒立摆系统的工作原理;由轴角编码器测得小车位置和摆杆相对垂直方向的角度,作为系统的两个输出量被反馈至控制计算机。

计算机根据一定的控制算法,计算出控制量,并转化为相应的电压信号提供给驱动电路,以驱动直流力矩电动机的运动,从而通过牵引机构带动小车的位移来控制摆杆和保持平衡。

在忽略了空气阻力,各种摩擦之后,可将直线一级倒立摆系统抽象成小车和匀质杆组成的系统。

2. 依据相关物理定理,列写倒立摆系统的运动方程;小车质量为M ,倒立摆的质量为m ,摆长为2l ,小车的位置为x ,摆的角度为摆杆绕其重心的转动方程:θθθcos sin l F l F J xy -= 摆杆重心的水平运动方程:)sin (22θl x dtd m F X += 摆杆重心的垂直运动方程:)cos (22θl dt d m mg F y =- 小车水平方向运动方程:22dt xd M F F x =- 一级倒立摆系统的动力学模型:()⎪⎪⎩⎪⎪⎨⎧++-+-⋅+⋅=-++-⋅+++=)(cos sin lg )(cos sin cos cos ))((cos sin sin )()(22222222222222m M ml J l m m M l m F ml l m m M ml J g l m ml J lm F ml J xθθθθθθθθθθθθ对系统进行线性化:1cos ,sin ,0≈≈≈θθθθ系统的简化模型:⎪⎪⎩⎪⎪⎨⎧++-+=++-+=θθθθsin )(lg )()()(22222mMl m M J mlF m M mMl m M J g l m F ml J x3. 根据倒立摆的运动方程搭建被控对象在Simulink 环境下的仿真模型;模型中小车质量M=1kg,摆的质量m=1kg ,摆长2l=0.6m ,重力加速度g=10m /s 2,得到小车精确模型:⎪⎪⎩⎪⎪⎨⎧--⋅+⋅=--⋅+=24.0cos 09.0sin 6cos sin 09.0cos 3.0cos 09.024.0cos sin 9.0sin 036.012.022θθθθθθθθθθθθ F F x简化的模型:⎨⎧+=+-=u u x2408.06θθθ4. 结合单位反馈控制系统的控制原理,为被控对象设计PID 控制器。

(1)双闭环PID 控制器设计图4—1 一级倒立摆系统位置伺服控制系统方框图(2)内环控制器的设计内环采用反馈校正进行控制,其方框图如下图4—2图4—2 内环反馈校正方框图一级倒立摆内环外环反馈校正采用PD 控制器212)(K s K s D +=,为了抑制干扰在前向通道上加上一个比例环节K s D =)(2。

(3)控制参数的整定设2()D s 的增益20K =-,则内环控制系统的闭环传递函数为:40646464)()(1)()(2122222-++=+=K s K s s D s G KK s G KK s W s s 令7.0=ξ⎪⎩⎪⎨⎧⨯==-647.0264K1464064K 2x 解得:⎩⎨⎧==1.625K 0.175K 21 内环控制器的传递函数为:625.1175.0)(2+=s s D 内环控制器的闭环传递函数为:642.1164)(22++=s s s W (4)外环控制器的设计外环系统前向通道的传递函数为:)642.11()104.0(64)()(22222+++-=s s s s s G s W 外环系统结构图如下图4—3外环系统结构图系统开环传递函数为一个高阶且带不稳定零点的非最小相为系统,对系统外环模型进行降阶处理,若忽略)(2s W 的高次项,则可近似为一阶传递函数为:64s 2.1164)(2+≈s W()D s对模型1()G s 进行近似处理,则1()G s 的传递函数为:2110)(s s G ≈外环控制器采用PD 形式,其传递函数为: ()1)(31+≈s K s D τ为使系统具有较好的跟随性,采用单位反馈构成外环反馈通道,即1)(1≈s D ,则系统的开环传递函数为:()1)57(57)()()()(321122++==s K s s s D s G s W s W τ 采用基于Bode 图法的希望特性设计方法,得30.12,0.877K τ==,取1τ=,则外环控制器的传递函数为:()112.0)(1+=s s D一阶倒立摆双闭环控制系统方框图如下:图4—4 一阶倒立摆双闭环控制系统方框图(5)对一级倒立摆双闭环控制系统进行仿真。

一级倒立摆双闭环控制系统Simulink 仿真结构图如下图所示:图4—5一级倒立摆双闭环控制系统Simulink 仿真结构图一级倒立摆一级倒立摆双闭环控制系统Simulink仿真结果如下图所示:图4—6一级倒立摆双闭环控制系统Simulink仿真结果图系统仿真程序clear allload PID.mat%ÔËÐÐt14210.mdlt=signals(1,:);q=signals(2,:);x=signals(3,:);figure(1)ht=line(t,q(:));grid on;xlabel('ʱ¼ä(s)');ylabel('°Ú½Ç±ä»¯(rad)');axis([0,10,-0.3,1.2]);axet=axes('Position',get(gca,'Position'),'XAxisLocation','bottom',...'YAxisLocation','right','color','None','XColor','k','YColor','k'); ht=line(t,x,'color','r','parent',axet);ylabel('xλÖñ仯(m)');axis([0,10,-0.3,1.2]);title('\Theta(t) ºÍ x(t)µÄ½×Ô¾ÏìÓ¦ÇúÏß')gtext('\leftarrow x(t)'),gtext('\Theta(t)\uparrow');5. 分析综述比例P、积分I、微分D三个调节参数对系统控制性能的影响。

PID控制器就是根据系统的误差,利用比例、积分、微分计算出控制量进行控制的。

比例(P)控制比例控制是一种最简单的控制方式。

其控制器的输出与输入误差信号成比例关系。

当仅有比例控制时系统输出存在稳态误差(Steady-state error)。

积分(I)控制在积分控制中,控制器的输出与输入误差信号的积分成正比关系。

对一个自动控制系统,如果在进入稳态后存在稳态误差,则称这个控制系统是有稳态误差的或简称有差系统(System with Steady-state Error)。

为了消除稳态误差,在控制器中必须引入“积分项”。

积分项对误差取决于时间的积分,随着时间的增加,积分项会增大。

这样,即便误差很小,积分项也会随着时间的增加而加大,它推动控制器的输出增大使稳态误差进一步减小,直到等于零。

因此,比例+积分(PI)控制器,可以使系统在进入稳态后无稳态误差。

微分(D)控制在微分控制中,控制器的输出与输入误差信号的微分(即误差的变化率)成正比关系。

自动控制系统在克服误差的调节过程中可能会出现振荡甚至失稳。

其原因是由于存在有较大惯性组件(环节)或有滞后(delay)组件,具有抑制误差的作用,其变化总是落后于误差的变化。

解决的办法是使抑制误差的作用的变化“超前”,即在误差接近零时,抑制误差的作用就应该是零。

这就是说,在控制器中仅引入“比例”项往往是不够的,比例项的作用仅是放大误差的幅值,而目前需要增加的是“微分项”,它能预测误差变化的趋势,这样,具有比例+微分的控制器,就能够提前使抑制误差的控制作用等于零,甚至为负值,从而避免了被控量的严重超调。

所以对有较大惯性或滞后的被控对象,比例+微分(PD)控制器能改善系统在调节过程中的动态特性。

相关文档
最新文档