当前位置:文档之家› matlab实验九控制系统的PI校正设计及仿真

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

matlab实验九控制系统的PI校正设计及仿真
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;

滞后:

-----------

s + 1

练习9-2设被控对象开环传递函数为:

系统的相位裕度γ0

40≥

系统斜坡输入稳态误差e ss =0.05;

截止频率为c ω=2 1/ S

要求:

(1)用频率法设计满足上述要求的串联滞后校正控制器;

(2)画出校正前后的Bode 图;

(3)用Simulink 对校正前后的闭环系统进行仿真,求出其阶跃响应;

(4)分析设计效果。

滞后:

5 s + 1

-----------

s + 1

k0=20;n1=10;

d1=conv([1,0],[1,5]);

w=logspace(-1,3,1000);

sope=tf(k0*n1,d1);

figure(1)

margin(sope);grid on

wc=2;

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;

练习9-3 已知单位反馈系统被控对象开环传递函数为:

试用BODE 图设计方法对系统进行滞后串联校正设计,使之满足:

(1)在单位斜坡信号r (t )= t 的作用下,系统的速度误差系数K V 1

30-≥s ;

(1) 系统校正后剪切频率ωC 13.2-≥s ;

(3)系统斜校正后相角裕度γ040≥

要求:

(1)用频率法设计满足上述要求的串联滞后校正控制器;

(2)画出校正前后的Bode 图;

(3)用Simulink 对校正前后的闭环系统进行仿真,求出其阶跃响应;

(4)分析设计效果。

clear

%s=tf('s')

%G=30/s+1)+1)

k0=30;n1=1;

d1=conv(conv([1,0],[,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超前滞后校正概要

课程设计任务书 学生姓名: 张弛 专业班级: 电气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通信原理仿真 学号: 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 仿真实验报告 课题名称:MATLAB 仿真——图像处理 学院:机电与信息工程学院 专业:电子信息科学与技术 年级班级:2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作,熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能,理解多维图像的相关性质 3、熟悉Help 命令的使用,掌握对相关函数的查找,了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用,能够处理多维图像。 二、实验条件

MATLAB调试环境以及相关图像处理的基本MATLAB语句,会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents : ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating(联系起来)A1, A2 ... along the dimension DIM. Calls to CAT can be nested(嵌套). (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts(整数下标). (4)、Manipulating multi-dimensional arrays

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 及仿真实验(控制系统的时域分析) 一、实验目的 学习利用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 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仿真实验解析

实验报告 课程名称:控制理论(乙) 指导老师:林峰 成绩:__________________ 实验名称:MATLAB 仿真实验 实验类型:________________同组学生姓名:__________ 一、实验目的和要求(必填) 二、实验内容和原理(必填) 三、主要仪器设备(必填) 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析(必填) 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的: 1.用计算机辅助分析的办法,掌握系统的时域分析方法。 2.熟悉Simulink 仿真环境。 二、实验原理及方法: 系统仿真实质上就是对系统模型的求解,对控制系统来说,一般模型可转化成某个微分方程或差分方程表示,因此在仿真过程中,一般以某种数值算法从初态出发,逐步计算系统的响应,最后绘制出系统的响应曲线,进而可分析系统的性能。控制系统最常用的时域分析方法是,当输入信号为单位阶跃和单位冲激函数时,求出系统的输出响应,分别称为单位阶跃响应和单位冲激响应。在MATLAB 中,提供了求取连续系统的单位阶跃响应函数step ,单位冲激响应函数impulse ,零输入响应函数initial 等等。 二、实验内容: 二阶系统,其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求: 1.编制MATLAB 程序,画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应; (1)画出系统的单位阶跃响应曲线; A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];

MATLAB仿真实验报告

MATLA仿真实验报告 学院:计算机与信息学院 课程:—随机信号分析 姓名: 学号: 班级: 指导老师: 实验一

题目:编写一个产生均值为1,方差为4的高斯随机分布函数程序, 求最大值,最小值,均值和方差,并于理论值比较。 解:具体的文件如下,相应的绘图结果如下图所示 G仁random( 'Normal' ,0,4,1,1024); y=max(G1) x=mi n(G1) m=mea n(G1) d=var(G1) plot(G1);

实验二 题目:编写一个产生协方差函数为CC)=4e":的平稳高斯过程的程序,产生样本函数。估计所产生样本的时间自相关函数和功率谱密度,并求统计自相关函数和功率谱密度,最后将结果与理论值比较。 解:具体的文件如下,相应的绘图结果如下图所示。 N=10000; Ts=0.001; sigma=2; beta=2; a=exp(-beta*Ts); b=sigma*sqrt(1-a*a); w=normrnd(0,1,[1,N]); x=zeros(1,N); x(1)=sigma*w(1); for i=2:N x(i)=a*x(i-1)+b*w(i); end %polt(x); Rxx=xcorr(x0)/N; m=[-N+1:N-1]; Rxx0=(sigma A2)*exp(-beta*abs(m*Ts)); y=filter(b,a,x) plot(m*Ts,RxxO, 'b.' ,m*Ts,Rxx, 'r');

periodogram(y,[],N,1/Ts); 文件旧硯化)插入(1〕 ZMCD 克闻〔D ]窗口曲) Frequency (Hz) 50 100 150 200 250 300 350 400 450 500 NH---.HP)&UO 二 balj/ 」- □歹

基于MatlabSimulink控制系统校正器的设计与仿真

第3卷第3期2003年9月 南京工业职业技术学院学报 J0u瑚】ofN皿j堍hlsdn】te0f111dus田Tcchnology V0l3.No3 SeD..2003 基于Matlab/Simulink控制系统校正器的设计与仿真 周昱英 (南京工业职业技术学院自动化系,江苏南京210016) 摘要:介绍了№dab软件的主要功能和特.軎'.并利用其提供的siⅡ出nk,讨论了如何进行控制系统校正嚣的设计和仿真。 关键词:M“d日b软件;陆ntdink;模糊控制器;P∞控制器 中图分类号:T挖73文献标识码:A文章编号:1671—46“(2003)03—0017—04 D鹤ignaIldSimIllationofthe C咖pen鼢torinControlSyste】咀 byUsiIlgMauab/SimuliIll【 ZHOUYu—ving 【肫唧曙肺∞i№。,h妇叫z酗∞妇y。讹耐昭210016,蕊iⅥ) Al妞rad:Thispaperi11删ucesthemajor矗mcti蚰s锄dfhturesofMadabsoftwar|e锄ddiscusses110wtodesj即andsir眦latethe∞mpensa时inconnolsystemby毓删b11l【0fMadab. K叮words:mtlab蚓ham;simlllink;血可con嘣1er;PDcon响Uer 引言 作为研究控制系统重要手段,计算机仿真技术随着计算机语言的不断更新,也在不断发展。Madab作为一套集数值计算、符号运算及图形处理等强大功能于一体的科学计算语言,最早于1984年由Mat}1works公司推出,目前已发展成为国际公认的最出色的仿真软件之一。作为强大的科学计算平台,它几乎能够满足所有的计算需求,其强大的扩展功能为各领域的应用提供了基础。面向控制领域,它推出的模糊控制、神经网络、控制系统等工具箱为控制系统的设计和仿真提供了有力的支持、极大地推动了其领域内仿真研究的发展。 №dab的最新版本为MaL】ab6.1,它具有一系列新的优势和特点,其中包括更加友好的工作平台和编程环境,增强的模块集和工具箱,实用的程序接口和发布平台等。特别是它的模块化设计和动态仿真工具simuliI】k4.1,使用的时候,不需要编写任何MaⅡab语句,只要利用鼠标直接在模型窗口上“画”出所需要的控制系统模型,然后直接利用Simulink提供的功能对系统进行仿真和线性分析,从而使一个很复杂的系统仿真变得相当容易且直观,简化了设计过程,减轻了设计负担,提高了仿真的集成化和可视化。鉴于6.1版的这些特点,本文将利用它的SiIllIllink来研究如何进行控制系统校正器设计、仿真和分析。 1基于Simllink的控制系统建模与仿真 假设控制系统的传递函数为:专淄=r墨未b,其中c(s)=溅,求它的阶跃输出 响应。在simdh山环境中,利用鼠标将相应的模块 收稿日期:20∞一晒一26 作者简介:周昱荚(1976一),女,湖北武汉人’南京工业职业技术学院自动化系讲师.工学硕士。

MATLAB程序设计实验报告

MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法; 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法,利用绘图命令绘制出典型信号的波形,了解这些信号的基本特征; 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法,利用绘图命令绘制出典型信号的波形,了解这些信号的基本特征; 4. 通过绘制信号运算结果的波形,了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 : 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘,信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 (1)信号的相加和相乘:已知信号)(1t f 和)(2t f ,信号相加和相乘记为 )()(1t f t f =)(2t f +;)()(1 t f t f =)(2t f *。 (2)信号的微分和积分:对于连续时间信号,其微分运算是用diff 函数来完成的,其语句格式为:diff(function,’variable’,n),其中function 表示需要进行求导运算的信号,或者被赋值的符号表达式;variable 为求导运算的独立变量;n 为求导的阶数,默认值为求一阶导数。连续信号的积分运算用int 函数来完成,语句格式为:diff(function,’variable’,a,b),其中function 表示需要进行被积信号,或者被赋值的符号表达式;variable 为求导运算的独立变量;a,b 为积分上、下限,a 和b 省略时为求不定积分。 (3)信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移,信号的翻转包含信号的倒相与折叠,平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化,可使信号压缩或扩展。)(at f 将原波形压缩a 倍,)/(a t f 将原波形扩大a 倍。 ¥ 对系统频率特性的分析

用MATLAB进行控制系统的超前校正设计汇总

课程设计任务书 学生姓名: ________________ 专业班级:____________ 指导教师:陈启宏_______ 工作单位:自动化学院 题目:用MATLA进行控制系统的超前校正设计。 初始条件:已知一单位反馈系统的开环传递函数是 K G(s) s(1 0.05s)(1 0.5s) 要求系统跟随2r/min的斜坡输入产生的最大稳态误差为2°,45。 要求完成的主要任务:(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求) 1、用MATLAB^出满足初始条件的最小K值的系统伯德图,计算系统的幅值裕量和相 位裕量。 2、在系统前向通路中插入一相位超前校正,确定校正网络的传递函数。 3、用MATLA画出未校正和已校正系统的根轨迹。 4、课程设计说明书中要求写清楚计算分析的过程,列出MATLAE程序和 MATLA输出。说明书的格式按照教务处标准书写。 时间安排: 2012 年12月17 日

系主任(或责任教师)签名:

用MATLA进行控制系统的 超前校正设计 1、超前校正概述 i.i、何谓校正 所谓校正,就是在系统中加入一些其参数可以根据需要而改变的机构或装置,使系统整个特性发生变化,从而满足给定的各项性能指标。校正的目的就是为了当我们在调整放大器增益后仍然不能满足设计所要求的性能指标的情况下,通过加入合适的校正装置,使系统的性能全面满足设计要求。 按照校正装置在控制系统中的连接方式,可以将校正方式分为串联校正和并联校正两种。在用分析法进行串联校正时,校正环节的结构通常采用超前校正、滞后校正、滞后一一超前校正这三种类型,也就是工程上常用的PID调节器。 在实际的分析设计中,具体采用哪种校正方式,取决于系统的校正要求、信号的性质、系统各点的功率、可选元件和经济条件等。 本次课程设计的要求为用MATLA进行控制系统的超前校正设计,已知一单位反馈系统的开环传递函数是: G(s) K s(1 0.05s)(1 0.5s) 要求系统跟随2r/min的斜坡输入产生的最大稳态误差为 2 所以接下来将对超前校正进行相应的介绍。45

《MATLAB与控制系统。。仿真》实验报告

《MATLAB与控制系统仿真》 实验报告 班级: 学号: 姓名: 时间:2013 年 6 月

目录实验一MATLAB环境的熟悉与基本运算(一)实验二MATLAB环境的熟悉与基本运算(二)实验三MATLAB语言的程序设计 实验四MATLAB的图形绘制 实验五基于SIMULINK的系统仿真 实验六控制系统的频域与时域分析 实验七控制系统PID校正器设计法 实验八线性方程组求解及函数求极值

实验一MATLAB环境的熟悉与基本运算(一) 一、实验目的 1.熟悉MATLAB开发环境 2.掌握矩阵、变量、表达式的各种基本运算 二、实验基本原理 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MATLAB常用命令 表1 MATLAB常用命令 变量与运算符 3.1变量命名规则 3.2 MATLAB的各种常用运算符 表3 MATLAB关系运算符 表4 MATLAB逻辑运算符

| Or 逻辑或 ~ Not 逻辑非 Xor逻辑异或 符号功能说明示例符号功能说明示例 :1:1:4;1:2:11 . ;分隔行.. ,分隔列… ()% 注释 [] 构成向量、矩阵!调用操作系统命令 {} 构成单元数组= 用于赋值 的一维、二维数组的寻访 表6 子数组访问与赋值常用的相关指令格式 三、主要仪器设备及耗材 计算机 四.实验程序及结果 1、新建一个文件夹(自己的名字命名,在机器的最后一个盘符) 2、启动MATLAB,将该文件夹添加到MATLAB路径管理器中。 3、学习使用help命令。

自动控制原理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) 1.1设计内容 (1) 1.2 设计要求 (1) 1.3 Matlab软件 (2) 1.3.1基本功能 (2) 1.3.2应用 (3) 第二章控制系统程序设计 (4) 2.1 校正装置计算方法 (4) 2.2 课程设计要求计算 (4) 第三章利用Matlab仿真软件进行辅助分析 (6) 3.1校正系统的传递函数 (6) 3.2用Matlab仿真 (6) 3.3利用Matlab/Simulink求系统单位阶跃响应 (10) 3.2.1原系统单位阶跃响应 (10) 3.2.2校正后系统单位阶跃响应 (11) 3.2.3校正前、后系统单位阶跃响应比较 (12) 3.4硬件设计 (13) 3.4.1在计算机上运行出硬件仿真波形图 (14) 课程设计心得体会 (16) 参考文献 (18)

第一章 课程设计内容与要求分析 1.1设计内容 针对二阶系统 )1()(+= s s K s W , 利用有源串联超前校正网络(如图所示)进行系统校正。当开关S 接通时为超前校正装置,其传递函数 11 )(++-=Ts Ts K s W c c α, 其中 132R 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 ,开环截止频率ωc’≥4.4弧度/秒,相位裕量γ’≥45°; 2)根据性能指标要求,确定串联超前校正装置传递函数; 3)利用对数坐标纸手工绘制校正前、后及校正装置对数频率特性曲线; c R R

线性代数MATLAB仿真实验报告

合肥学院 2018—2019学年第2学期 线性代数及应用 (模块) 实验报告 实验名称:线性代数MATLAB实验 实验类别:综合性 设计性□验证性 专业班级: 17通信工程(2)班 实验时间: 9-12周 组别:第组人数 3人 指导教师:牛欣成绩: 完成时间: 2019年 5 月9日

一. 小组成员 姓名学号具体分工 汪蔚蔚(组长) 1705022025 A报告最后的整合,编写,案例四的计算与应用 以及案例一的计算与证明 陶乐 1 1705022009 C案例二,化学方程式配平问题 程赢妹1505022036 A案例三,应用题灰度值的计算问题 二. 实验目的 1、案例一利用MATLAB进行线性代数计算,求出矩阵B 2、案例二利用MATLAB计算出每一个网格数据的值,然后每一个网格数据的值乘以256以后进行归一化处理,根据每个网格中的灰度值,绘制出灰度图像。 3、案例三利用MATLAB完成对化学方程式进行配平的应用 4、案例四利用MATLAB求极大线性无关组,并表示出其余向量 三. 实验内容 1、案例一: 0,1,0 ,=1,0,0, 0,0,0 A B AB BA A B ?? ?? =?? ?? ?? 已知矩阵和矩阵满足乘法交换律,即且求矩阵。 2、案例二 配平下列化学方程式: 3、案例三: 3*32 0.81.21.70.20.3 0.6021.61.20.6. 1MATLAB 2256MATLAB 给定一个图像的个方向上的灰度叠加值:沿左上方到右 下方的灰度叠加值依次为,,,,;沿右上方到左下 方的灰度叠加值依次为,。,,, )建立可以确定网络数据的线性方程组,并用求解 )将网络数据乘以,再取整,用绘制该灰度图像

matlab实验报告

软件课程设计 实验报告 模块名称:__MATLAB_____________ 专业班级:__ __________ 姓名:__ _____________ 学号:__ ____________ 实验日期:年月日实验报告日期:年月日

一、要求练习的实验部分 1. 题目:在时间区间 [0,10]中,绘出t e y t 2cos 15.0--=曲线。 实现方法: 在[0,10]区间中每隔0.1取一点,构造t 数组,然后利用exp 和cos 函数计算出y ,最后用plot 函数画出图形。 程序: t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果及结果分析: 2. 题目:写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波,要 求它的负半波被置零,且在 2 3 处被削顶。 图 E2-1 题15所要表示的波形(逐段解析函数) 实现方法: 将t 在区间[0,3*pi]用向量表示,分成500点,乘法要用.*,通过y=0和a=sin(pi/3)将函数y=sin(t)截断比较再相加得到f(t); 对于y>0的部分进行截取,代码为z=(y>=0).*y ; 对于y=a).*b+(y

z=(y>=0).*y; z=(y>=a).*a+(ytol) - rank(A) 是多少 ? (2) S(1,1) - norm(A) = 0 是多少 ? (3) sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? (4) S(1,1)/S(3,3) - cond(A) 的结果是什么 ? (5) S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么? (6) V(:,1)'*null(A) 得到什么结果 ? (7) abs(A*null(A)) < tol 得到什么结果 ? (8) U(:,1:2) = = orth(A) 的运行结果是什么 ? 实现方法: 键入 A=[1 2 3;4 5 6;7 8 9] 得 A =1 2 3 4 5 6 7 8 9 键入 [U,S,V]=svd(A) 得 U = -0.2148 0.8872 0.4082 -0.5206 0.2496 -0.8165 -0.8263 -0.3879 0.4082 S =16.8481 0 0 0 1.0684 0 0 0 0 V =-0.4797 -0.7767 -0.4082

相关主题
文本预览
相关文档 最新文档