最新课件-机电系统动态仿真matlab电子教案第7章系统频
合集下载
机电系统动态仿真matlabPPT电子教案课件-第七章-SIMULINK仿真
离散滤波器 离散零极点 离散状态方程
零阶保持器 一阶保持器
Simulink的基本模块
5.函数与表模块库
模块
功能
Lookup Table
线性插值查表
Lookup Table（2-D） 二维线性插值
Lookup Table（n-D） N维线性插值
PrelookUp Index S 预查下标
Interpolation（n-D） N维插值
认识Simulink
7.1.2 Simulink的启动与退出 1．Simulink的启动
在MATLAB窗口的工具栏中单击
图标
在命令窗口中输入命令： >>simulink
认识Simulink
7.1.2 Simulink的启动与退出 1．Simulink的启动
认识Simulink
7.1.2 Simulink的启动与退出 1．Simulink的启动
Simulink模块库浏览器
标题 菜单栏 工具栏 模块说明 框
基本模块 库
已安装专 用 模块库
模块查找 框
模块显示 框
认识Simulink
如何打开模型编辑窗口？ 启动方式： (1)模块库浏览器的菜单“File”/“New”/“Model”命令 (2)单击工具栏上的 图标
认识Simulink
方法一：
认识Simulink
7.1.1 Simulink简介 Simulink是MATLAB的重要组成部分，提供建立
系统模型、选择仿真参数和数值算法、启动仿真程 序对该系统进行仿真、设置不同的输出方式来观察 仿真结 启动Simulink并打开模型编辑窗口。 (2) 将所需模块添加到模型中。 (3) 设置模块参数并连接各个模块组成仿真模型。 (4) 设置系统仿真参数。 (5) 开始系统仿真。 (6) 观察仿真结果。
2024版matlab教程(全)资料ppt课件
进行通信系统的建模、仿真和分析。
谢谢聆听
B
C
变量与赋值
在MATLAB中，变量不需要事先声明，可以 直接赋值。变量名以字母开头，可以包含字 母、数字和下划线。
常用函数
MATLAB提供了丰富的内置函数，如sin、 cos、tan等三角函数，以及abs、sqrt等数 学函数。用户可以通过help命令查看函数的
D
使用方法。
02 矩阵运算与数组操作
错误处理
阐述try-catch错误处理机制的语法、 执行流程及应用实例。
04
函数定义与调用
函数概述
阐述函数的概念、作用及分类，包括内置函数和 自定义函数。
函数调用
深入剖析函数的调用方法，包括直接调用、间接 调用及参数传递等技巧。
ABCD
函数定义
详细讲解自定义函数的定义方法，包括函数名、 输入参数、输出参数及函数体等要素。
拟合方法
利用已知数据点构造近似函数，如最小二乘法、多项 式拟合、非线性拟合等。
插值与拟合的比较
插值函数经过所有数据点，而拟合函数则追求整体上 的近似。
数值积分与微分
01
数值积分方法
利用数值技术计算定积分的近似 值，如矩形法、梯形法、辛普森 法等。
02
数值微分方法
通过数值技术求解函数的导数或 微分，如差分法、中心差分法、 五点差分法等。
02
01
矩阵运算
加法与减法
对应元素相加或相减，要求矩阵 大小相同
乘法
使用`*`或`mtimes`函数进行矩阵 乘法，要求内维数相同
点乘与点除
使用`.*`、`./`进行对应元素相乘或 相除，要求矩阵大小相同
特征值与特征向量
matlab教学PPT第7讲MATLAB仿真SIMULINK
>>whos NameSizeBytesClas s ScopeData1x13578structarray Grandtotalis307elementsusing3578bytes
第7讲 SIMULINK 图2-18 通过命令启动Simulink仿真
第7讲 SIMULINK
可见，仿真完成之后，工作空间中出现了“ScopeData” 结构变量，其中包含了示波器显示的全部波形数据。通过 “plot”命令可以作出这些数据对应的波形，
组建用户自定义模块库如果建立了许多自定义的子系统并且已经封装好了而这些已经封装的自定义模块又是会反复使用的就像simulink提供的模块库中那些模块一样在这种情况下就有必要对这些自定义的重用性较好的模块进行建库以方便管理和反复使用同时也可以作为新的专业库提供给其他用户使用
第7讲 SIMULINK
第7讲 MATLAB仿真_SIMULINK
第7讲 SIMULINK
• Simulink全方位地支持动态系统的建模仿真，它支持连 续系统、离散系统、连续离散混合系统、线性系统、非 线性系统、时不变系统、时变系统的建模仿真，也支持 具有多采样速率的多速率系统。可以说，在通用系统仿 真领域，Simulink是无所不包的。
• 结合MATLAB编程和Simulink可视化建模仿真各自的特 点，可以构建更为复杂的系统模型，并进行自动化程度 更高的仿真和仿真结果的数据分析，这是MATLAB的高 级应用方面。
第7讲 SIMULINK 图2-15 仿真结果
第7讲 SIMULINK
• 更换信号源为Sources子模块库中的SignalGenerator，并设置 信号源为0.2Hz的方波，幅度为1，如图2-16左边对话框所示。
• 设置示波器显示窗口的属性（Parameters），使之成为双踪 显示，然后将示波器第二输入节点与信号源输出相连，这 样我们就可以同时观察系统的输入输出波形了。系统建模 如图2-16中间窗口所示。
第7讲 SIMULINK 图2-18 通过命令启动Simulink仿真
第7讲 SIMULINK
可见，仿真完成之后，工作空间中出现了“ScopeData” 结构变量，其中包含了示波器显示的全部波形数据。通过 “plot”命令可以作出这些数据对应的波形，
组建用户自定义模块库如果建立了许多自定义的子系统并且已经封装好了而这些已经封装的自定义模块又是会反复使用的就像simulink提供的模块库中那些模块一样在这种情况下就有必要对这些自定义的重用性较好的模块进行建库以方便管理和反复使用同时也可以作为新的专业库提供给其他用户使用
第7讲 SIMULINK
第7讲 MATLAB仿真_SIMULINK
第7讲 SIMULINK
• Simulink全方位地支持动态系统的建模仿真，它支持连 续系统、离散系统、连续离散混合系统、线性系统、非 线性系统、时不变系统、时变系统的建模仿真，也支持 具有多采样速率的多速率系统。可以说，在通用系统仿 真领域，Simulink是无所不包的。
• 结合MATLAB编程和Simulink可视化建模仿真各自的特 点，可以构建更为复杂的系统模型，并进行自动化程度 更高的仿真和仿真结果的数据分析，这是MATLAB的高 级应用方面。
第7讲 SIMULINK 图2-15 仿真结果
第7讲 SIMULINK
• 更换信号源为Sources子模块库中的SignalGenerator，并设置 信号源为0.2Hz的方波，幅度为1，如图2-16左边对话框所示。
• 设置示波器显示窗口的属性（Parameters），使之成为双踪 显示，然后将示波器第二输入节点与信号源输出相连，这 样我们就可以同时观察系统的输入输出波形了。系统建模 如图2-16中间窗口所示。
《MATLAB Simulink与控制系统仿真(第3版)》的课件 第7章 频域分析法
MATLAB 6.0以上版本还提供了计算系统稳定裕度的函数allmargin，其调用法 则如下。 s=allmargin(sys)：计算幅值裕度、相角裕度以及对应的频率。幅值裕度和相角 裕度是针对开环SISO系统而言，输出s是一个结构体，它包括幅值裕度、相角 裕度以及对应的频率、时滞增益裕度。
闭环系统的等M圆
等M圆轨迹
等M圆和系统开环幅相频率特性
闭环系统的等N圆
7.6 MATLAB在频率法中的应用
MATLAB提供了多种求取并绘制频率响应曲线的函数，如Nyquist曲线绘制函 数nyquist( )、Bode图绘制函数bode( )和Nichols曲线绘制函数nichols( )等。 1．Nyquist曲线绘制函数nyquist( ) 2．Bode图绘制函数bode( ) 3．Nichols曲线绘制函数nichols( ) 4．绘制等M圆和等N圆的函数ngrid ( )
7.7 频率法的稳定性分析
Nyquist稳定判据
7.7.2 稳定裕度
系统的相对稳定性用Nyquist曲线相对点（1, j0）的靠近程度来度量，定量表 示为增益裕度和相角裕度。
7.7.3 MATLAB在稳定性分析中的应用
margin函数可以从频率响应数据中计算出幅值裕度、相角裕度以及对应的频率。 幅值裕度和相角裕度是针对开环SISO系统而言的，它指出了系统在闭环时的 相对稳定性。当不带输出变量引用时，margin可在当前图形窗口中绘出带有裕 量及相应频率显示的Bode图，其中的幅值裕度以分贝为单位。
对数幅相图（Nichols图）
7.3.3 对数幅相图（Nichols图）
对数幅相图也称Nichols图，它是将对数幅频特性和相频特性两张图在角频率 为参变量的情况下合成为一张图。
闭环系统的等M圆
等M圆轨迹
等M圆和系统开环幅相频率特性
闭环系统的等N圆
7.6 MATLAB在频率法中的应用
MATLAB提供了多种求取并绘制频率响应曲线的函数，如Nyquist曲线绘制函 数nyquist( )、Bode图绘制函数bode( )和Nichols曲线绘制函数nichols( )等。 1．Nyquist曲线绘制函数nyquist( ) 2．Bode图绘制函数bode( ) 3．Nichols曲线绘制函数nichols( ) 4．绘制等M圆和等N圆的函数ngrid ( )
7.7 频率法的稳定性分析
Nyquist稳定判据
7.7.2 稳定裕度
系统的相对稳定性用Nyquist曲线相对点（1, j0）的靠近程度来度量，定量表 示为增益裕度和相角裕度。
7.7.3 MATLAB在稳定性分析中的应用
margin函数可以从频率响应数据中计算出幅值裕度、相角裕度以及对应的频率。 幅值裕度和相角裕度是针对开环SISO系统而言的，它指出了系统在闭环时的 相对稳定性。当不带输出变量引用时，margin可在当前图形窗口中绘出带有裕 量及相应频率显示的Bode图，其中的幅值裕度以分贝为单位。
对数幅相图（Nichols图）
7.3.3 对数幅相图（Nichols图）
对数幅相图也称Nichols图，它是将对数幅频特性和相频特性两张图在角频率 为参变量的情况下合成为一张图。
机电系统仿真课案PPT课件
4/48
第4页/共46页
一. 仿真技术概论
2. 仿真的必要性
希望能在较短的时间内观察到系统发展的全 过程，以估计某些参数对系统行为的影响．
难以在实际环境中进行实验和观察时，计算 机仿真是唯一可行的方法， 例如太空飞行的 研究．
需要对系统或过程进行长期运行比较，从大 量方案中寻找最优方案．
5/48
14/48
第14页/共46页
一. 仿真技术概论
5. 系统仿真
（4）系统仿真的一般步骤：
①建立系统的数学模型；
②转换成仿真模型；
③编写仿真程序；
④对仿真模型进行修改校验，看 与实际系统是否一致，确认模型 的正确性。
⑤运行仿真程序，在不同的初始 15/48 条件和参数下，对系统进行第反15页复/共46页
3/48
第3页/共46页
一. 仿真技术概论
2. 仿真的必要性
难以用数学公式表示的系统，或者没有建立 和求解数学模型的有效方法．虽然可以用解 析的方法解决问题，但数学的分析与计算过 于复杂，这时计算机仿真可能提供简单可行 的求解方法
对于随机性系统，可以通过大量的重复试验， 获得其平均意义上的特性指标
具有刚性、阻力及惯性这些性能，在输入力作用下，
19/48
系统能产生一定相应输出，如位移输出。
第19页/共46页
二 . 机电系统建模
2. 机械系统的数学模型
（1）机械移动系统
m
k
c
机械平移系统基本元件
(a)质量；(b)弹簧；(c)阻尼
20/48
第20页/共46页
二 . 机电系统建模
2. 机械系统的数学模型
d 2 x2 dt 2
?
30/48
第4页/共46页
一. 仿真技术概论
2. 仿真的必要性
希望能在较短的时间内观察到系统发展的全 过程，以估计某些参数对系统行为的影响．
难以在实际环境中进行实验和观察时，计算 机仿真是唯一可行的方法， 例如太空飞行的 研究．
需要对系统或过程进行长期运行比较，从大 量方案中寻找最优方案．
5/48
14/48
第14页/共46页
一. 仿真技术概论
5. 系统仿真
（4）系统仿真的一般步骤：
①建立系统的数学模型；
②转换成仿真模型；
③编写仿真程序；
④对仿真模型进行修改校验，看 与实际系统是否一致，确认模型 的正确性。
⑤运行仿真程序，在不同的初始 15/48 条件和参数下，对系统进行第反15页复/共46页
3/48
第3页/共46页
一. 仿真技术概论
2. 仿真的必要性
难以用数学公式表示的系统，或者没有建立 和求解数学模型的有效方法．虽然可以用解 析的方法解决问题，但数学的分析与计算过 于复杂，这时计算机仿真可能提供简单可行 的求解方法
对于随机性系统，可以通过大量的重复试验， 获得其平均意义上的特性指标
具有刚性、阻力及惯性这些性能，在输入力作用下，
19/48
系统能产生一定相应输出，如位移输出。
第19页/共46页
二 . 机电系统建模
2. 机械系统的数学模型
（1）机械移动系统
m
k
c
机械平移系统基本元件
(a)质量；(b)弹簧；(c)阻尼
20/48
第20页/共46页
二 . 机电系统建模
2. 机械系统的数学模型
d 2 x2 dt 2
?
30/48
机电系统动态仿真matlabPPT电子教案课件-第七章-SIMULINK仿真
功能
积分 微分 状态方程 传递函数 零极点 传输延时 可变传输延时
23
Simulink的基本模块
SIMULINK仿真 4.离散系统模块库
模块
Unit Delay Discrete-Time I Discrete Transf F
功能
单位延时采样保持 离散时间积分 离散传递函数
Discrete Filter Discrete Zero-Pole Discrete State-Space Zero-order Hold First-order Hold
45
仿真模型的参数设置
SIMULINK仿真
7.4.3 启动系统仿真与仿真结果分析 设置完仿真参数之后，从Simulation中选择Start菜 单项或单击模型编辑窗口中的Start Simulation命令 按钮，便可启动对当前模型的仿真。
为了观察仿真结果的变化轨迹可以采用3种方法： (1) 把输出结果送给Scope模块或者XY Graph模块。 (2) 把仿真结果送到输出端口并作为返回变量，然后使
34
仿真模型的参数设置
SIMULINK仿真 1．模块的参数设置
35
仿真模型的参数设置
SIMULINK仿真 2．模块的属性设置 ✓ 模块上按鼠标右键并在弹出的快捷菜单中选择Block properties ✓ 在模型编辑窗口的Edit菜单下选择Block properties命令，将打开模块属性对话框。
5
认识Simulink
SIMULINK仿真 7.1.2 Simulink的启动与退出 1．Simulink的启动
在MATLAB窗口的工具栏中单击 在命令窗口中输入命令：
>>simulink
(完整版)matlab第七讲教案
西南科技大学本科生课程备课教案计算机技术在安全工程中的应用——Matlab入门及应用授课教师：徐中慧班级：专业：安全技术及工程第七章逻辑函数与控制结构课型：新授课教具：多媒体教学设备，matlab教学软件一、目标与要求能够正确使用函数find，正解和正确使用if/else系列命令，理解switch/case结构，正确使用for循环和while循环。
二、教学重点与难点本堂课教学的重点在于引导学生在编写matlab程序时能够熟练运用控制结构的相关函数实现相应的功能。
三、教学方法本课程主要通过讲授法、演示法、练习法等相结合的方法来引导学生掌控本堂课的学习内容。
四、教学内容课后习题讲解%% onea=[15 3 22;3 8 5;14 3 82];b=[1 5 6]';c=[12 18 5 2];d=a(:,3)%% twoe=[b,d]%% threef=[b;d]%% fourg=[c(1:3);a]%% five%% Variable DefinitionDensity_mercury=13560;g=9.81;p=0:10000:100000;%% mercuryh_mercury=p./(Density_mercury*g)%% waterh_water=h_mercury.*13560/1000%% outputh=[h_mercury',h_water']（3）已知下列矩阵：1104265878235645913232289x ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦ 1234101272127y ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ []1022513z =①查找各矩阵中所有大于10的元素的序号，输出结果用单序号方法表示。
②查找各矩阵中所有大于10的元素的序号，输出结果用该元素的下标表示。
③查找各矩阵中所有大于10的元素。
④查找各矩阵中所有大于10小于40的元素的序号，输出结果用单序号方法表示。
MATLAB在电类专业课程中的应用第7章课件
电枢电流接近最大值，使得电动机以最快速度
开始上升，最高转速达到230r/min，超调量为
15%，稳态时转速为200r/min；当2秒时给定信
号为-10时，电动机从电动状态变成制动状态，
当转速为零时，电动机开始反转。
7.4 交流调压系统
7.4.1交流调压调速系统原理
异步电机的机械特性方程：
=
3 12 2′ /
换器控制，通过示波器查看电流和转速的波形。
双闭环调速系统的结构在
单闭环调速的基础上增加
了电流环，其余部分与单
闭环调速系统的结构相同，
包括双极式可逆PWM变换
器、直流电机（DC
Machine）、转速反馈
（Gain）等模块。
从仿真结果可以看出，当给定电压为10V时，在
电机起动过程中，电流调节器作用下的电动机
Ug1=Ug4=-Ug2=-Ug3。在一个开关周期内，当0 ≤ t <
时，VT1和VT4导通，UAB=Us，电枢电流经过
VT1MVT4流通；当 ≤ t < T时，驱动电流反
向，但VT2和VT3由于受续流二极管VD2和VD3反向
电压钳住却无法导通，UAB=-Us，电枢电流经过
Vd2MVd3流通，因此在一个周期内UAB具有正
负相间的脉冲波形。
电动机的正反转由驱动电压正负脉冲的宽窄决定，当
正脉冲较宽时， > 2 则UAB平均值为正，电机正转；
反之电机反转，如果正负脉冲宽度相等，则平均电压
UAB为零，电机停止。双极式控制可逆变换器的输出平
均电压为：
=
=
2
− 1 = (2 − 1)
转矩保持不变；当启动1s时，负载由20
开始上升，最高转速达到230r/min，超调量为
15%，稳态时转速为200r/min；当2秒时给定信
号为-10时，电动机从电动状态变成制动状态，
当转速为零时，电动机开始反转。
7.4 交流调压系统
7.4.1交流调压调速系统原理
异步电机的机械特性方程：
=
3 12 2′ /
换器控制，通过示波器查看电流和转速的波形。
双闭环调速系统的结构在
单闭环调速的基础上增加
了电流环，其余部分与单
闭环调速系统的结构相同，
包括双极式可逆PWM变换
器、直流电机（DC
Machine）、转速反馈
（Gain）等模块。
从仿真结果可以看出，当给定电压为10V时，在
电机起动过程中，电流调节器作用下的电动机
Ug1=Ug4=-Ug2=-Ug3。在一个开关周期内，当0 ≤ t <
时，VT1和VT4导通，UAB=Us，电枢电流经过
VT1MVT4流通；当 ≤ t < T时，驱动电流反
向，但VT2和VT3由于受续流二极管VD2和VD3反向
电压钳住却无法导通，UAB=-Us，电枢电流经过
Vd2MVd3流通，因此在一个周期内UAB具有正
负相间的脉冲波形。
电动机的正反转由驱动电压正负脉冲的宽窄决定，当
正脉冲较宽时， > 2 则UAB平均值为正，电机正转；
反之电机反转，如果正负脉冲宽度相等，则平均电压
UAB为零，电机停止。双极式控制可逆变换器的输出平
均电压为：
=
=
2
− 1 = (2 − 1)
转矩保持不变；当启动1s时，负载由20
第7章系统频率响应仿真课件
频率特性的一般概念 连续系统频率特性的MATLAB函数
离散系统的频域仿真 系统分析图形用户界面
33
7.4 系统分析图形用户界面
⑴在指令窗中建立系统模型 ⑵在指令窗中输入：ltiview ⑶点击菜单File，选择Import ⑷在LTI浏览器中，选择sys，调入 ⑸显示调入系统模型的仿真曲线
图7-14 LTI Viewer窗口
u(t)
MATLAB仿真演示pltxys.mdl
u(t) 2sin(5t 30)
红色为系统的输入信号 蓝色曲线为系统的全响应信号， 经过一段时间后与黑色重合
-4
黑色为系统的稳态输出
-6
红 —输 入 ， 蓝 —全 响 应 ， 黑 —稳 态 响 应 2
y(t)
-8
1.5
0
1
2
3
4
5
6
t/s
1
yss(t)
传递函数：G(s)
b0 s m a0 s n
b1sm1 a1sn1
bm an
输入量（正弦信号）： r(t) Ar sin(t)
R(s) Ar
Ar
s2 2 (s j)(s j)
则系统输出为
C(s) G(s)R(s) U (s) A V (s) s2 2
U (s)
(s s1)(s s2 )
u(t) 2sin(20t 30)
0.5
幅值
结论：系统稳定后，输出
0
信号的振幅和相位一般均 -0.5
不同于输入信号，且随着
-1
输入信号频率的变化而变
-1.5
y(t)
化，但输出信号频率=输入 -2
u(t)
信号频率
0
81
离散系统的频域仿真 系统分析图形用户界面
33
7.4 系统分析图形用户界面
⑴在指令窗中建立系统模型 ⑵在指令窗中输入：ltiview ⑶点击菜单File，选择Import ⑷在LTI浏览器中，选择sys，调入 ⑸显示调入系统模型的仿真曲线
图7-14 LTI Viewer窗口
u(t)
MATLAB仿真演示pltxys.mdl
u(t) 2sin(5t 30)
红色为系统的输入信号 蓝色曲线为系统的全响应信号， 经过一段时间后与黑色重合
-4
黑色为系统的稳态输出
-6
红 —输 入 ， 蓝 —全 响 应 ， 黑 —稳 态 响 应 2
y(t)
-8
1.5
0
1
2
3
4
5
6
t/s
1
yss(t)
传递函数：G(s)
b0 s m a0 s n
b1sm1 a1sn1
bm an
输入量（正弦信号）： r(t) Ar sin(t)
R(s) Ar
Ar
s2 2 (s j)(s j)
则系统输出为
C(s) G(s)R(s) U (s) A V (s) s2 2
U (s)
(s s1)(s s2 )
u(t) 2sin(20t 30)
0.5
幅值
结论：系统稳定后，输出
0
信号的振幅和相位一般均 -0.5
不同于输入信号，且随着
-1
输入信号频率的变化而变
-1.5
y(t)
化，但输出信号频率=输入 -2
u(t)
信号频率
0
81
MATLAB与电力系统仿真 ppt课件
(2)利用开始（Start）导航区启动：
单击开始按钮，选择仿真（Simulink）命令，再选择电力 系统仿真命令（SimPowerSystem）,在弹出的对话框中选择 电力系统元件库（Block Library）命令即可
2021/3/26
MATLAB与电力系统仿真 ppt课件
MATLAB应用技术
2
（2）输配电线路元件（Lines） 在电力系统设计和分析中，输配电线路一般用各种类型的等
值电路来进行简化以便于简化分析。输配电线路元件的作用 就是构成各种线路的等值电路，在输配电线路元件中包括3 种元件
2021/3/26
MATLAB与电力系统仿真 ppt课件
MATLAB应用技术
6
（3）断路器元件（Circuit Breakers） 在电力系统中，断路器的作用是通断高压电力线路，可靠地
开始时间（Start time）：0 停止时间（Stop time）：0.4 求解程序类型（Type）选项：可变步长（Variable），
Ode45(Domand-Price) 最大步长（Max step size）选项：自动（Auto） 最小步长（Min step size）选项：自动（Auto） 初始步长（Initial step size）选项：自动（Auto） 相对容差（Relative tolerance）选项：1e-3 绝对容差（Absolute tolerance）选项：1e-6
2．电力系统元件库简介 在电力系统元件库对话框中包含了10类库元件，分别是
电源元件（Electrical Sources）、线路元件（Elements）、 电力电子元件（Power Electronics）、电机元件 （Machines）、连接器元件（Connectors）、电路测量仪器 （Measurements）、附加元件（Extras）、演示教程 （Demos）、电力图形用户接口（Powergui）、电力系统元 件库模型（Powelib_models）。 1）电源元件 电源元件库包括7类元件，分别为：直流电压源元件（DC Voltage Source），交流电压源元件（AC Voltage Source），交流电流源元件（AC Current Source），受控 电压源元件（Controlled Voltage Source），受控电流源 元件（Controlled Current Source），三相电源元件（3Phase Source），三相可编程电压源元件（3-Phase Programmable Voltage Source）。
单击开始按钮，选择仿真（Simulink）命令，再选择电力 系统仿真命令（SimPowerSystem）,在弹出的对话框中选择 电力系统元件库（Block Library）命令即可
2021/3/26
MATLAB与电力系统仿真 ppt课件
MATLAB应用技术
2
（2）输配电线路元件（Lines） 在电力系统设计和分析中，输配电线路一般用各种类型的等
值电路来进行简化以便于简化分析。输配电线路元件的作用 就是构成各种线路的等值电路，在输配电线路元件中包括3 种元件
2021/3/26
MATLAB与电力系统仿真 ppt课件
MATLAB应用技术
6
（3）断路器元件（Circuit Breakers） 在电力系统中，断路器的作用是通断高压电力线路，可靠地
开始时间（Start time）：0 停止时间（Stop time）：0.4 求解程序类型（Type）选项：可变步长（Variable），
Ode45(Domand-Price) 最大步长（Max step size）选项：自动（Auto） 最小步长（Min step size）选项：自动（Auto） 初始步长（Initial step size）选项：自动（Auto） 相对容差（Relative tolerance）选项：1e-3 绝对容差（Absolute tolerance）选项：1e-6
2．电力系统元件库简介 在电力系统元件库对话框中包含了10类库元件，分别是
电源元件（Electrical Sources）、线路元件（Elements）、 电力电子元件（Power Electronics）、电机元件 （Machines）、连接器元件（Connectors）、电路测量仪器 （Measurements）、附加元件（Extras）、演示教程 （Demos）、电力图形用户接口（Powergui）、电力系统元 件库模型（Powelib_models）。 1）电源元件 电源元件库包括7类元件，分别为：直流电压源元件（DC Voltage Source），交流电压源元件（AC Voltage Source），交流电流源元件（AC Current Source），受控 电压源元件（Controlled Voltage Source），受控电流源 元件（Controlled Current Source），三相电源元件（3Phase Source），三相可编程电压源元件（3-Phase Programmable Voltage Source）。
《MATLAB与系统仿真》PPT课件
是系统某种特定性能的一种抽象形式。
模型的表达形式有物理模型和数学模型两类。
数学模型是系统的某种特征本质的数学表达式， 是用数学公式来描述所研究的客观对象或系统中 的某一方面的问题。
数
静态模型
学
连续系统模型
模
动态模型 离散系统模型
型
编辑版ppt
15
三、系统仿真
一个较流行于工程技术界的定义是：仿真是 通过对系统模型的实验去研究一个存在的或设计 中的系统。这种定义适用于概括了所有工程的 (技术的)或非工程的(非技术的)系统。
30
工程实际对仿真技术提出的新需求：
１、减少模型的开发时间； ２、提高模型建立的精度和实验的精度； ３、改进人与人、人与计算机的通信。
编辑版ppt
31
本课程主要讲授内容 １、系统建模的基本方法与模型处理技术 ２、连续系统的数字仿真程序通用算法
－数值积分仿真方法学 ３、连续系统模型的离散化处理技术 4、计算机仿真软件－MATLAB
动而发生变化的系统进行仿真称为离散事件系统 仿真。其数学模型多用流程图或网络图来描述。
（3）混合系统仿真：
当系统的数学模型是由上述两类模型混合 构成时，称为混合系统仿真。其仿真方法是将 上述两类方法综合于一体。
（4）系统动力学仿真： 当对象的数学模型是用系统动力学方程式来
描述时，该系统的仿真称为系统动力学仿真。
法。模拟计算机由一些基本的模拟运算部件组成， 这些运算部件有：积分器、加法器、系数器、函 数发生器、乘法器等。
模拟计算机是并行运算的，运算速度快，但 精度不高，由于它可以实现传递函数为1/s的积 分运算，可以方便地求解微分方程。
编辑版ppt
24
（2）数字计算机仿真：
模型的表达形式有物理模型和数学模型两类。
数学模型是系统的某种特征本质的数学表达式， 是用数学公式来描述所研究的客观对象或系统中 的某一方面的问题。
数
静态模型
学
连续系统模型
模
动态模型 离散系统模型
型
编辑版ppt
15
三、系统仿真
一个较流行于工程技术界的定义是：仿真是 通过对系统模型的实验去研究一个存在的或设计 中的系统。这种定义适用于概括了所有工程的 (技术的)或非工程的(非技术的)系统。
30
工程实际对仿真技术提出的新需求：
１、减少模型的开发时间； ２、提高模型建立的精度和实验的精度； ３、改进人与人、人与计算机的通信。
编辑版ppt
31
本课程主要讲授内容 １、系统建模的基本方法与模型处理技术 ２、连续系统的数字仿真程序通用算法
－数值积分仿真方法学 ３、连续系统模型的离散化处理技术 4、计算机仿真软件－MATLAB
动而发生变化的系统进行仿真称为离散事件系统 仿真。其数学模型多用流程图或网络图来描述。
（3）混合系统仿真：
当系统的数学模型是由上述两类模型混合 构成时，称为混合系统仿真。其仿真方法是将 上述两类方法综合于一体。
（4）系统动力学仿真： 当对象的数学模型是用系统动力学方程式来
描述时，该系统的仿真称为系统动力学仿真。
法。模拟计算机由一些基本的模拟运算部件组成， 这些运算部件有：积分器、加法器、系数器、函 数发生器、乘法器等。
模拟计算机是并行运算的，运算速度快，但 精度不高，由于它可以实现传递函数为1/s的积 分运算，可以方便地求解微分方程。
编辑版ppt
24
（2）数字计算机仿真：
《MATLAB Simulink与控制系统仿真(第4版)》的课件 第7章 频域分析法
MATLAB 6.0以上版本还提供了计算系统稳定裕度的函数allmargin，其调用法 则如下。 s=allmargin(sys)：计算幅值裕度、相角裕度以及对应的频率。幅值裕度和相角 裕度是针对开环SISO系统而言，输出s是一个结构体，它包括幅值裕度、相角 裕度以及对应的频率、时滞增益裕度。
极坐标图（Nyquist图）
7.3.2 对数坐标图（Bode图）
Bode图的优点可概括如下： （1）将幅值相乘化为对数相加运算，大大简化了系统频率特性 的绘制工作。 （2）由于横轴采用了对数分度，缩小了比例尺，从而扩大了频 率视野，可以在较大的频段范围内表示系统频率特性。在一张 Bode图上，既画出了频率特性的中、高频段，又能清楚地画出其 低频段，在分析和设计系统时，低频段特性也是非常重要的。 （3）可以绘制渐近的对数幅频特性；也可以制作标准样板，画 出精确的对数频率特性。
第7章 频域分析法
7.1 引言 7.2 频率特性基本概念 7.3 频率特性的表示方法 7.4 系统开环频率特性作图 7.5 频率响应分析 7.6 MATLAB在频率法中的应用 7.7 频率法的稳定性分析 7.8 综合实例及MATLAB/SIMULINK应用 习题
内容提要
控制系统的频率特性反映的是系统对正弦输入信号的 响应性能。频域分析法是一种图解分析法，它依据系 统频率特性对系统的性能（如稳定性、快速性和准确 性）进行分析。
频域分析法的突出优点是可以通过试验直接求得频率 特性来分析系统的品质，应用频率特性分析系统可以 得出定性和定量的结论，并具有明显的物理含义。
对数幅相图（Nichols图）
7.3.3 对数幅相图（Nichols图）
对数幅相图也称Nichols图，它是将对数幅频特性和相频特性两张图在角频率 为参变量的情况下合成为一张图。
极坐标图（Nyquist图）
7.3.2 对数坐标图（Bode图）
Bode图的优点可概括如下： （1）将幅值相乘化为对数相加运算，大大简化了系统频率特性 的绘制工作。 （2）由于横轴采用了对数分度，缩小了比例尺，从而扩大了频 率视野，可以在较大的频段范围内表示系统频率特性。在一张 Bode图上，既画出了频率特性的中、高频段，又能清楚地画出其 低频段，在分析和设计系统时，低频段特性也是非常重要的。 （3）可以绘制渐近的对数幅频特性；也可以制作标准样板，画 出精确的对数频率特性。
第7章 频域分析法
7.1 引言 7.2 频率特性基本概念 7.3 频率特性的表示方法 7.4 系统开环频率特性作图 7.5 频率响应分析 7.6 MATLAB在频率法中的应用 7.7 频率法的稳定性分析 7.8 综合实例及MATLAB/SIMULINK应用 习题
内容提要
控制系统的频率特性反映的是系统对正弦输入信号的 响应性能。频域分析法是一种图解分析法，它依据系 统频率特性对系统的性能（如稳定性、快速性和准确 性）进行分析。
频域分析法的突出优点是可以通过试验直接求得频率 特性来分析系统的品质，应用频率特性分析系统可以 得出定性和定量的结论，并具有明显的物理含义。
对数幅相图（Nichols图）
7.3.3 对数幅相图（Nichols图）
对数幅相图也称Nichols图，它是将对数幅频特性和相频特性两张图在角频率 为参变量的情况下合成为一张图。
MATLAB 第七章 SIMULINK基础PPT课件
14
SIMULINK仿真的运行
启动仿真
▪ 设置仿真参数和选择解法器之后，就可以启动仿真而运行。 ▪ 选择Simulink菜单下的start选项来启动仿真，如果模型中有些参数没有
定义，则会出现错误信息提示框。如果一切设置无误，则开始仿真运 行，结束时系统会发出一鸣叫声。 ▪ 除了直接在SIMULINK环境下启动仿真外，还可以在MATLAB命令窗 口中通过函数进行，格式如下：
– Real-Imag to Complex：由实部和虚部输入合成复数输出
12
7、Sinks（接收器模块）
– Scope：示波器 – XY Graph：显示二维图形 – To Workspace：将输出写入MATLAB的工作空间 – To File(.mat)：将输出写入数据文件 – Terminator：连接到没有连接到的输出端
▪ SIMULINK是MATLAB软件的扩展，它是实现动态系统建模和仿真的 一个软件包，它与MATLAB语言的主要区别在于，其与用户交互接口 是基于Windows的模型化图形输入，其结果是使得用户可以把更多的 精力投入到系统模型的构建，而非语言的编程上。
▪ 所谓模型化图形输入是指SIMULINK提供了一些按功能分类的基本的 系统模块，用户只需要知道这些模块的输入输出及模块的功能，而不 必考察模块内部是如何实现的，通过对这些基本模块的调用，再将它 们连接起来就可以构成所需要的系统模型（以.mdl文件进行存取）， 进而进行仿真与分析。
① Start time 和 Stop time 栏目分别允许用户填写仿真的起始 时间和结束时间。
② Solver options 的 Type 栏目有两个选项，允许用户选择定 步长和变步长算法。
③仿真精度控制有 Relative Tolerance 选项、Absolute Tolerance 等，其中相对误差限的默认值设置为 1e-3，该值在实 际仿真中显得偏大，建议选择 1e-6 和 1e-7。值得指出的是，由 于采用的变步长仿真算法，所以将误差限设置到这样小的值也 不会增加太大的运算量。
SIMULINK仿真的运行
启动仿真
▪ 设置仿真参数和选择解法器之后，就可以启动仿真而运行。 ▪ 选择Simulink菜单下的start选项来启动仿真，如果模型中有些参数没有
定义，则会出现错误信息提示框。如果一切设置无误，则开始仿真运 行，结束时系统会发出一鸣叫声。 ▪ 除了直接在SIMULINK环境下启动仿真外，还可以在MATLAB命令窗 口中通过函数进行，格式如下：
– Real-Imag to Complex：由实部和虚部输入合成复数输出
12
7、Sinks（接收器模块）
– Scope：示波器 – XY Graph：显示二维图形 – To Workspace：将输出写入MATLAB的工作空间 – To File(.mat)：将输出写入数据文件 – Terminator：连接到没有连接到的输出端
▪ SIMULINK是MATLAB软件的扩展，它是实现动态系统建模和仿真的 一个软件包，它与MATLAB语言的主要区别在于，其与用户交互接口 是基于Windows的模型化图形输入，其结果是使得用户可以把更多的 精力投入到系统模型的构建，而非语言的编程上。
▪ 所谓模型化图形输入是指SIMULINK提供了一些按功能分类的基本的 系统模块，用户只需要知道这些模块的输入输出及模块的功能，而不 必考察模块内部是如何实现的，通过对这些基本模块的调用，再将它 们连接起来就可以构成所需要的系统模型（以.mdl文件进行存取）， 进而进行仿真与分析。
① Start time 和 Stop time 栏目分别允许用户填写仿真的起始 时间和结束时间。
② Solver options 的 Type 栏目有两个选项，允许用户选择定 步长和变步长算法。
③仿真精度控制有 Relative Tolerance 选项、Absolute Tolerance 等，其中相对误差限的默认值设置为 1e-3，该值在实 际仿真中显得偏大，建议选择 1e-6 和 1e-7。值得指出的是，由 于采用的变步长仿真算法，所以将误差限设置到这样小的值也 不会增加太大的运算量。
《机电系统与仿真技术》课件7 MatLab基础
2023/5/13
1.2发展阶段与工程应用的意义
发展阶段
第一阶段 20世纪50年代末到60年，为仿真技术的诞生期（只有大企 业用）； 第二阶段 20世纪70年代末到80年，为仿真技术的成长期（开始出现 研究人员专门研究仿真技术）； 第三阶段 20世纪90年代至今，为仿真技术的成熟期（大量仿真软件 出现并开始应用于科研和工程，如MultiSim, Protel, Tanner, MatLab， SolidWorks等）。
7）Length命令（给出向量长度）
8）Format命令（定义输出格式）
2023/5/13
3.1Matlab数值计算—算术运算
1）加法
>> a=[1 2 3;4 5 6];b=[7 8 9;10 11 12];c=1;
>> a+b
ans =
8
10
12
14
16
18
2）减法
>> a=[1 2 3;4 5 6];b=[7 8 9;10 11 12];
工程应用的意义
1）在经济方面，可以降低成本，而且设备可以重复使用，尤其是对 于大型、复杂系统而言； 2）一些危险的装置如核电站等通常是不允许进行实验的，因此采用 仿真技术可以降低危险程度，对系统研究起到保障作用； 3）提高设计效率，如电路设计、模型设计、控制系统设计等等； 4）具有优化设计和预测性能的特殊功能。
2023/5/13
1.3Matlab的特点
1）编程效率高，因为其编程接近于人们通常进行计算 的思维方式； 2）计算功能强，因为有非常丰富的库函数，矩阵、数 组和向量的计算功能特别强，适用于科学与工程计算； 3）使用方便，MatLab将编绎、链接、执行融为一体， 可以在同一窗口上排除书写、语法错误，加快了用户编 写、修改和调试程序的速度； 4）易于扩充，MatLab可以C、C++、Fortran混合编程。
1.2发展阶段与工程应用的意义
发展阶段
第一阶段 20世纪50年代末到60年，为仿真技术的诞生期（只有大企 业用）； 第二阶段 20世纪70年代末到80年，为仿真技术的成长期（开始出现 研究人员专门研究仿真技术）； 第三阶段 20世纪90年代至今，为仿真技术的成熟期（大量仿真软件 出现并开始应用于科研和工程，如MultiSim, Protel, Tanner, MatLab， SolidWorks等）。
7）Length命令（给出向量长度）
8）Format命令（定义输出格式）
2023/5/13
3.1Matlab数值计算—算术运算
1）加法
>> a=[1 2 3;4 5 6];b=[7 8 9;10 11 12];c=1;
>> a+b
ans =
8
10
12
14
16
18
2）减法
>> a=[1 2 3;4 5 6];b=[7 8 9;10 11 12];
工程应用的意义
1）在经济方面，可以降低成本，而且设备可以重复使用，尤其是对 于大型、复杂系统而言； 2）一些危险的装置如核电站等通常是不允许进行实验的，因此采用 仿真技术可以降低危险程度，对系统研究起到保障作用； 3）提高设计效率，如电路设计、模型设计、控制系统设计等等； 4）具有优化设计和预测性能的特殊功能。
2023/5/13
1.3Matlab的特点
1）编程效率高，因为其编程接近于人们通常进行计算 的思维方式； 2）计算功能强，因为有非常丰富的库函数，矩阵、数 组和向量的计算功能特别强，适用于科学与工程计算； 3）使用方便，MatLab将编绎、链接、执行融为一体， 可以在同一窗口上排除书写、语法错误，加快了用户编 写、修改和调试程序的速度； 4）易于扩充，MatLab可以C、C++、Fortran混合编程。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
离散系统的频域仿真 系统分析图形用户界面
2020/10/24
8
7.1 频率响应与频率特性
频率特性是控制工程最基本的概念之 一
何谓频率响应？ 频率响应有什么特性？
何谓频率特性？
如何求系统的频率响应？
为什么可以用频率响应来研究系统的性能？
演示程序：pltxys.mdl
2020/10/24
9
7.1.1 频率特性的概念与定义
2020/10/24
10
频率特性的定义
线性稳定系统在正弦信号作用下，当频率从零变 化到无穷时，稳态输出与输入的幅值比、相位差随频 率变化的特性，称为频率特性。频率特性由幅频特性 和相频特性两部分组成。
幅频特性
系统稳态正弦输出信号与相应的正弦输人信号的幅值之比随输 入频率的变比而变化的特性称为幅频特性，它描述了系统对输入 信号幅值的放大、衰减特性。
sys2=series(sys1,sys05)
sys3=feedback(sys2,1,-1)
sys4=tf(1,[5.6 0 0])
sys02=tf([0.2 0],1)
sys5=feedback(sys4,sys02,-1)
[n,d]=tfdata(sys1,'v')
2020/10/24
6
m2x2 f fc fk
2020/10/24
s^4 + 0.2679 s^3 + 2.005 2 + 04.25 s
m2 x2 f fc1 fc2 m1x1 f c1 f k fc1 c1 (x2 x1 ) f c2 c2 x2
f k kx1
m2s2 X 2 F Fc1 Fc2 m1s2 X1 Fc1 Fk Fc1 c1s( X 2 X1)
系统建模、模型转换及仿真 ⑴用MATLAB语言表示图5-1所示系统的模型，并给出三 种模型的具体表达；
⑵进行单位阶跃响应和周期为30s的方波响应。
clc,clf sys1=tf([1,0.5],[1,0.1]); sys2=zpk([],[0,-2,-10],20) sys3=series(sys1,sys2); sys=feedback(sys3,1,-1) subplot(121),step(sys) [u,t]=gensig('square',30,60); subplot(122),lsim(sys,'y',u,t)
频率响应的特点：
Xi=Aisinωt
稳定系统
Xo=Ao(ω)sin[ωt+φ(ω)]
幅值比 相位差
A() A0 ()
1
Ar
1 2T 2
它是频率ω的函数， 称为幅频特性；
() x0 (t) xi (t) tg 1T
它是频率ω的函数， 称为相频特性。
结论
幅值比、相位差随频率ω变化 频率响应的这种特性，用频率特性描述。
2020/10/24
1
function y=f(x) sym x if x<5
y=0 elseif x>10
y=2*x else
y=(1-x)/4 end
x=input('x=') if x<5
y=0 elseif x>10
y=2*x else
y=(1-x)/4 end
2020/10/24
2
function[sysm1, sysm2]=modelm4(sysp) m1=sysp(1); m2=sysp(2); k=sysp(3); c=sysp(4); A=[0 0 1 0;0 0 0 1;-k/m1 k/m1 -c/m1 c/m1;k/m2 -k/m2 c/m2 -c/m2 ]; B=[0 0 0 1/m2]'; C=[0 1 0 0 ]; M=[1 0 0 0]; D=0; sys2=ss(A,B,C,D); sys1=ss(A,B,M,D); sysm2=tf(sys2); sysm1=tf(sys1);
c1 0.5N / m s 1 c2 0.2N / m s 1
m1 3.5kg
c2=sysp(3);
m2 5.6kg
m1=sysp(4);
m2=sysp(5);
A=[0 1 0 0;-k/m1 -c1/m1 0 c1/m1; 0 0 0 1;0 c1/m2 0 -(c1+c2)/m2]; B=[0 0 0 1/m2]';
Fc2 c2sX 2 Fk kX1
X2 F
m1m2 s 4
(c1m2
m1s2 c1s k c1m1 c2m1)s3 (c1m2 km2
c1c2 )s2
(c1 c2 )ks
X2 F
19.6s4
3.5s2 0.5s 7 5.25s3 39.3s2
4.9s
2020/10/24
5
F
－Fc1 ╳
－
Fc2
1/(m2s2) x2 ╳ －
c2s
c1s
Fc1 ╳
x1 1/(m1s2)
－
Fk
k
clc sys0=tf(1,[3.5 0 0]) sys1=feedback(sys0,7,-1) sys05=tf([0.5 0],1)
sys7=tf(d,n) sys8=series(sys7,sys3) sys=feedback(sys5,sys8,-1)
m1x1 fc fk
fc c(x2 x1)
fk k(x2 x1)
X2 F
m1s2 cs k cm1s3 (m2 km1)s2
X1 F
cm1s3
cs k (m2 km1)s2
2020/10/24
7
第7章 系统频率响应及其仿真
频率特性的一般概念 频率特性的图解方法
系统的频域分析 频率特性的MATLAB函数
sysp=[7 0.5 0.2 3.5 5.6]; X2_F=modelm3(sysp)
C=[0 0 1 0];
D=0;
Transfer function:
sys1=ss(A,B,C,D);
0.1786 s^2 + 0.02551 s + 0.3571
sysm=tf(sys1);
-------------------------------------
sysp=[12 38 1000 0.1]; [X1_F,X2_F]=modelm4(sysp)
2020/10/24
3
试用MATLAB语言表示如图所示系统。当分别以x2和 f 为系统输出、输入 时的传递函数模型和状态空间模型。
k 7N / m
function[sysm]=modelm3 (sysp) k=sysp(1); c1=sysp(2);
2020/10/24
8
7.1 频率响应与频率特性
频率特性是控制工程最基本的概念之 一
何谓频率响应？ 频率响应有什么特性？
何谓频率特性？
如何求系统的频率响应？
为什么可以用频率响应来研究系统的性能？
演示程序：pltxys.mdl
2020/10/24
9
7.1.1 频率特性的概念与定义
2020/10/24
10
频率特性的定义
线性稳定系统在正弦信号作用下，当频率从零变 化到无穷时，稳态输出与输入的幅值比、相位差随频 率变化的特性，称为频率特性。频率特性由幅频特性 和相频特性两部分组成。
幅频特性
系统稳态正弦输出信号与相应的正弦输人信号的幅值之比随输 入频率的变比而变化的特性称为幅频特性，它描述了系统对输入 信号幅值的放大、衰减特性。
sys2=series(sys1,sys05)
sys3=feedback(sys2,1,-1)
sys4=tf(1,[5.6 0 0])
sys02=tf([0.2 0],1)
sys5=feedback(sys4,sys02,-1)
[n,d]=tfdata(sys1,'v')
2020/10/24
6
m2x2 f fc fk
2020/10/24
s^4 + 0.2679 s^3 + 2.005 2 + 04.25 s
m2 x2 f fc1 fc2 m1x1 f c1 f k fc1 c1 (x2 x1 ) f c2 c2 x2
f k kx1
m2s2 X 2 F Fc1 Fc2 m1s2 X1 Fc1 Fk Fc1 c1s( X 2 X1)
系统建模、模型转换及仿真 ⑴用MATLAB语言表示图5-1所示系统的模型，并给出三 种模型的具体表达；
⑵进行单位阶跃响应和周期为30s的方波响应。
clc,clf sys1=tf([1,0.5],[1,0.1]); sys2=zpk([],[0,-2,-10],20) sys3=series(sys1,sys2); sys=feedback(sys3,1,-1) subplot(121),step(sys) [u,t]=gensig('square',30,60); subplot(122),lsim(sys,'y',u,t)
频率响应的特点：
Xi=Aisinωt
稳定系统
Xo=Ao(ω)sin[ωt+φ(ω)]
幅值比 相位差
A() A0 ()
1
Ar
1 2T 2
它是频率ω的函数， 称为幅频特性；
() x0 (t) xi (t) tg 1T
它是频率ω的函数， 称为相频特性。
结论
幅值比、相位差随频率ω变化 频率响应的这种特性，用频率特性描述。
2020/10/24
1
function y=f(x) sym x if x<5
y=0 elseif x>10
y=2*x else
y=(1-x)/4 end
x=input('x=') if x<5
y=0 elseif x>10
y=2*x else
y=(1-x)/4 end
2020/10/24
2
function[sysm1, sysm2]=modelm4(sysp) m1=sysp(1); m2=sysp(2); k=sysp(3); c=sysp(4); A=[0 0 1 0;0 0 0 1;-k/m1 k/m1 -c/m1 c/m1;k/m2 -k/m2 c/m2 -c/m2 ]; B=[0 0 0 1/m2]'; C=[0 1 0 0 ]; M=[1 0 0 0]; D=0; sys2=ss(A,B,C,D); sys1=ss(A,B,M,D); sysm2=tf(sys2); sysm1=tf(sys1);
c1 0.5N / m s 1 c2 0.2N / m s 1
m1 3.5kg
c2=sysp(3);
m2 5.6kg
m1=sysp(4);
m2=sysp(5);
A=[0 1 0 0;-k/m1 -c1/m1 0 c1/m1; 0 0 0 1;0 c1/m2 0 -(c1+c2)/m2]; B=[0 0 0 1/m2]';
Fc2 c2sX 2 Fk kX1
X2 F
m1m2 s 4
(c1m2
m1s2 c1s k c1m1 c2m1)s3 (c1m2 km2
c1c2 )s2
(c1 c2 )ks
X2 F
19.6s4
3.5s2 0.5s 7 5.25s3 39.3s2
4.9s
2020/10/24
5
F
－Fc1 ╳
－
Fc2
1/(m2s2) x2 ╳ －
c2s
c1s
Fc1 ╳
x1 1/(m1s2)
－
Fk
k
clc sys0=tf(1,[3.5 0 0]) sys1=feedback(sys0,7,-1) sys05=tf([0.5 0],1)
sys7=tf(d,n) sys8=series(sys7,sys3) sys=feedback(sys5,sys8,-1)
m1x1 fc fk
fc c(x2 x1)
fk k(x2 x1)
X2 F
m1s2 cs k cm1s3 (m2 km1)s2
X1 F
cm1s3
cs k (m2 km1)s2
2020/10/24
7
第7章 系统频率响应及其仿真
频率特性的一般概念 频率特性的图解方法
系统的频域分析 频率特性的MATLAB函数
sysp=[7 0.5 0.2 3.5 5.6]; X2_F=modelm3(sysp)
C=[0 0 1 0];
D=0;
Transfer function:
sys1=ss(A,B,C,D);
0.1786 s^2 + 0.02551 s + 0.3571
sysm=tf(sys1);
-------------------------------------
sysp=[12 38 1000 0.1]; [X1_F,X2_F]=modelm4(sysp)
2020/10/24
3
试用MATLAB语言表示如图所示系统。当分别以x2和 f 为系统输出、输入 时的传递函数模型和状态空间模型。
k 7N / m
function[sysm]=modelm3 (sysp) k=sysp(1); c1=sysp(2);