新能源汽车性能仿真第3章动态系统模型及表示课件.ppt

《新能源汽车性能仿真第3章动态系统模型及表示课件.ppt》由会员分享，可在线阅读，更多相关《新能源汽车性能仿真第3章动态系统模型及表示课件.ppt（85页珍藏版）》请在三一办公上搜索。

1、新能源汽车性能仿真第3章动态系统模型及表示,课程主要内容及章节安排,第1章SIMULINK与系统仿真第2章 MATLAB数值计算分析第3章 动态系统模型及表示第4章 创建SIMULINK模型第5章 动态系统SIMULINK仿真及原理第6章 SIMULINK子系统与S函数第7章 电动汽车电池特性及仿真第8章 电动汽车用电机性能及仿真第9章 电动汽车控制策略及仿真第10章 电动汽车整车性能仿真,2023/3/20,第 3 页,第3章 动态系统模型及表示,3.1 简单系统模型及表示3.2 离散系统模型及表示3.3 连续系统模型及表示3.4 混合系统模型及表示,2023/3/20,第 4 页,3.1

2、简单系统模型及表示,3.1.1 简单系统的基本概念不同系统具有不同数量的输入与输出；一般来说，输入输出数目越多，系统越复杂。最简单的系统一般只有一个输入与输出，而且任一时刻的输出只与当前时刻的输入有关。本节首先介绍简单系统的基本概念以及简单系统的simulink表示。,2023/3/20,第 5 页,【定义】简单系统：(1)系统某一时刻的输出直接且唯一依赖于该时刻的输入量。(2)系统对同样的输入，其输出影响不随时间的变化而变化。(3)系统中不存在输入的状态量，所谓的状态量是指系统输入的微分项（及输入的导数项）。如果一个系统满足上述的条件，则称之为简单系统。,2023/3/20,第 6 页,设简

3、单系统的输入为x，系统的输出为y，x可以具有不同的物理含义。对于任何系统，都可以将它视为对输入量x的某种变换，因此可以用T 表示任意一个系统，即Y=Tx对于简单系统，x一般为时间变量或其它的物理变量，并具有一定的输入范围。系统输出变量y仅与x的当前值相关，从数学的角度来看，y是x的一个函数，给出一个x值，便有一个y值与之对应。,2023/3/20,第 7 页,实例3-1,其中u为系统的输入变量，t为时间变量，y为系统的输出变量。显然，此系统满足简单系统的条件，为一简单系统。系统输出仅由系统当前时刻的输入决定。,对下式所描述的一个系统,2023/3/20,第 8 页,3.1.2 简单系统的描述方

4、式一般来讲，简单系统都可以采用代数方程与逻辑结构相结合的方式进行描述。1、代数方程采用数学方程对简单系统进行描述，可以很容易由系统输入求出系统输出，并且由此可方便地对系统进行定量分析。,2023/3/20,第 9 页,2、逻辑结构一般来说，系统输入都有一定的范围。对于不同范围的输入，系统输出与输入之间遵从不同的关系。由系统的逻辑结构可以很容易了解系统的基本概况。,2023/3/20,第 10 页,3.1.3 简单系统的simulink描述本章主要介绍动态系统的基本知识，为使用simulink进行系统仿真打下基础。因此这里并不准备建立系统的simulink模型，而是采用编写M脚本文件的方式对系统

5、进行描述并进行简单的仿真。下面以【实例3-1】中的简单系统为例，说明在simulink中如何对简单系统进行描述。,2023/3/20,第 11 页,程序实例3-1,【例3.1】中的简单系统，编写如下的systemdemo1.m脚本文件进行描述与分析。u=0:0.1:10;%设定系统输出范围与仿真步长leng=length(u);%计算系统输入序列长度for i=1:leng%计算系统输出序列 if u(i)=1%逻辑判断 y(i)=u(i).2;else y(i)=sqrt(u(i);endendplot(u,y);grid on;,2023/3/20,第 12 页,实例3-1结果,简单系统的

6、输入输出关系图,2023/3/20,第 13 页,3.2 离散系统模型及表示,3.2.1 前面所涉及到的系统中，无论是系统的输入还是系统的输出均是连续的变量，在这里连续指的是系统的输入与输出均在时间变量上连续取值（与数学上函数连续概念并不相同）。本节将简单介绍离散系统的基本概念，系统的描述与简单仿真。所谓离散系统，是指系统的输入与输出仅在离散的时间上取值，而且离散的时间具有相同的时间间隔。下面给出离散系统更全面的定义。,2023/3/20,第 14 页,【定义】离散系统：(1)系统每隔固定的时间间隔才“更新”一次，即系统的输入与输出每个固定的时间间隔便改变一次。固定的时间间隔称为系统的“采样”

7、时间。(2)系统的输出依赖于系统当前的输入、以往的输入与输出，即系统的输出是它们的某种函数。(3)离散系统具有离散的状态。其中状态指的是系统迁移时刻的输出量。满足上述条件的系统称为离散系统。,2023/3/20,第 15 页,3.2.2 离散系统的数学描述前面给出了离散系统的定义，这里给出离散系统的数学描述。设系统输入变量为，其中为系统的采样时间，为采样时刻。显然，系统的输入变量每隔固定的时间间隔改变一次。由于为一固定的值，因而系统输入常被简记为。设系统输出为，同样也可以简记为。由离散系统的定义可知，其数学描述应为,2023/3/20,第 16 页,实例3-2,对如下的离散系统模型,其中系统的

8、初始状态为y(0)=3，系统输入为，则系统在时刻0,1,2 的输出分别为,2023/3/20,第 17 页,离散系统除了采用一般的数学描述方式之外，还可以采用差分方程进行描述。使用差分方程描述方程形式如下：设系统的状态变量为，离散系统差分方程由以下两个方程构成：,状态更新方程：系统输出方程：,2023/3/20,第 18 页,3.2.3 离散系统的simulink描述,这里以【实例3-2】中的离散系统为例，说明如何利用simulink对离散系统进行描述，并在此基础上对系统进行简单的分析。与前面类似，此处并不建立系统的simulink模型进行仿真，二是编写M脚本文件从原理上对离散系统进行说明，并

9、说明离散系统与连续系统的区别之处。,2023/3/20,第 19 页,程序实例3-2,编写脚本文件systemdemo2.m对【例3.2】中的离散系统进行描述分析。y(1)=3;u(1)=0;for i=2:11 u(i)=2*i;y(i)=u(i).2+2*u(i-1)+3*y(i-1);endplot(u,y);grid on;,2023/3/20,第 20 页,系统从时刻0到时刻10的输入与输出的关系如图3.2所示。其中横坐标表示离散系统的输入向量，而纵坐标表示离散系统的输出向量。说明：这里并没有指定离散系统的采样时间，而仅仅举例说明离散系统的求解分析。在实际的系统中，必须指定系统的采样

10、时间，只有这样才能获得离散系统真正的动态性能。,2023/3/20,第 21 页,实例3-2结果图,2023/3/20,第 22 页,3.2.4 线性离散系统,对于任何系统而言，系统的描述都可以采用抽象的数学形式来描述。这是因为任何系统可以被看做是输入到输出的某种变换。例如，离散系统可以由下述的变换进行描述,2023/3/20,第 23 页,在离散系统之中，线性离散系统具有重要的地位。下面对线性离散系统进行简单的介绍。首先给出如下的两个概念：(1)齐次性：若对于离散系统，如果对任意的输入与给定的任意常数，恒有,则称该系统满足齐次性。,2023/3/20,第 24 页,(2)叠加性：如果系统对于

11、输入和，输出分别为和，恒有,则称系统满足叠加性。,2023/3/20,第 25 页,【定义】线性离散系统：当离散系统同时满足如下式所示的齐次性与叠加性时，称该离散系统为线性离散系统。,2023/3/20,第 26 页,例如对如下的离散系统,2023/3/20,第 27 页,3.2.5 线性离散系统的数学描述,对于线性离散系统来说，可以使用一般的方式对其进行描述，如采用如下所示的数学方程进行描述：,或采用差分方程进行描述状态更新方程：系统输出方程：,2023/3/20,第 28 页,除了使用一般的方式描述离散系统外，针对线性离散系统本身的特点，经常使用Z变换来描述线性离散系统。Z变换是对离散信号

12、进行分析的一个强有力的工具，尤其是对线性离散系统。Z变换有丰富的内容，此处只简单介绍线性离散系统的Z变换域描述以及MATLAB中一些比较常用的对线性离散系统进行分析的函数。,2023/3/20,第 29 页,Z变换具有多种不同的性质，此处介绍几个：(1)线性性。即对于离散信号和，设它们的Z变换分别为与，所谓Z变换的线性性指的是Z变换满足下面的关系：,2023/3/20,第 30 页,(2)位移性1、双边Z变换的位移性质若序列x(n)的双边Z变换为,则其右移位后的Z变换为,左移位后的Z变换为,2023/3/20,第 31 页,2、单边z变换的位移性左移位性质若则其中m为正整数。,对于m=1,2的

13、情况,2023/3/20,第 32 页,右移位性质若则其中m为正整数。,对于m=1,2的情况,2023/3/20,第 33 页,初值定理：若x(n)为因果序列，已知,则有,2023/3/20,第 34 页,终值定理：若x(n)为因果序列，已知,则有,此时需要n，x(n)收敛。,2023/3/20,第 35 页,实例3-3,对于如下的线性离散系统,同时对等式两边进行Z变换，则有。一般在系统分析中，往往对系统输出与系统输入的比值比较关心，将此式化成分式的形式，有,2023/3/20,第 36 页,在对系统进行描述分析时，此种形式的描述称之为滤波器描述。对上式进行等价变换，可以得到系统的传递函数描述

14、线性系统最常见的一种描述方式：,2023/3/20,第 37 页,还可以得到系统的零极点描述：,2023/3/20,第 38 页,3.2.6 线性离散系统的simulink描述,线性离散系统的描述方式有如下四种方式：(1)线性离散系统的滤波器模型：在simulink中，滤波器表示为num=n0 n1 n2；den=d0 d1；其中num表示Z变换分式的分子系数向量，den为分母系数向量。(2)线性离散系统的传递函数模型：在simulink中，系统的传递函数表示为num=n0 n1 n2；den=d0 d1；,2023/3/20,第 39 页,(3)线性离散系统的零极点模型：在simulink中

15、，系统零极点表示为gain=K；zeros=z1 z2；poles=0 p1(4)线性离散系统的状态空间模型：在simulink中，设系统差分方程为如下形式：x(n+1)=Fx(n)+Gu(n)；y(n)=Cx(n)+Du(n)。其中x(n)，u(n)，y(n)分别为线性离散系统的状态变量、输入向量、输出向量。F，G，C，D分别为变换矩阵。在simulink中，其表示很简单，只需要输入相应的变换矩阵F，G，C，D即可。,2023/3/20,第 40 页,实例3-4,对如下所示的线性离散系统：,2023/3/20,第 41 页,在MATLAB中输入下面的语句，可以绘制出此系统的Bode图：num

16、=2-1-5;den=1 3 6 2;dbode(num,den,1)grid on此离散系统的Bode图如下所示。,2023/3/20,第 42 页,离散系统Bode图,2023/3/20,第 43 页,也可以用下面的语句起初系统的幅值与相位而不绘制图形：mag,phase=dbode(num,den,1);mag=0.3334 1.0000phase=-181.1458-360.0000,2023/3/20,第 44 页,此外，在MATLAB中，离散系统的不同描述模型之间可以相互转化。这里给出几个比较常用的函数：zeros,poles,k=tf2zp(num,den)将系统传递函数模型转化

17、为零极点模型。num,den=zp2tf(zeros,poles,k)将系统零极点模型转化传递函数模型。其中num，den分别为系统的传递函数表示；zeros，poles，k为系统的零极点模型。,2023/3/20,第 45 页,至于线性离散系统的状态空间模型描述，只给出它与传递函数模型、零极点模型相互转化的函数命令：zeros,poles,k=ss2zp(F,G,C,D)将系统状态空间模型转化为零极点模型。F,G,C,D=zp2ss(zeros,poles,k)将系统零极点模型转化为状态空间模型。num,den=ss2tf(F,G,C,D)将系统状态空间模型转化为传递函数模型F,G,C,D=

18、tf2ss(num,den)将系统传递函数模型转化为状态空间模型。,2023/3/20,第 46 页,实例3-5,以下面的线性离散系统为例说明系统模型的转化。,解：将传递函数模型转化为零极点模型：num=2 1 5;den=1 3 6 2;zeros,poles,k=tf2zp(num,den)结果为,2023/3/20,第 47 页,zeros=1.8508-1.3508poles=-1.2980+1.8073i-1.2980-1.8073i-0.4039 k=2,2023/3/20,第 48 页,将传递函数模型转化为状态空间模型：num=2-1-5;den=1 3 6 2;F,G,C,D=

19、tf2ss(num,den)结果为,2023/3/20,第 49 页,F=-3-6-2 1 0 0 0 1 0G=1 0 0C=2-1-5D=0,2023/3/20,第 50 页,3.3 连续系统模型及表示,3.3.1 连续系统的基本概念与离散系统不同，连续系统是指系统输出在时间上连续变化，而非仅在离散的采样时刻取值。连续系统的应用非常广泛，下面给出连续系统的基本概念。,2023/3/20,第 51 页,【定义】连续系统：(1)系统输出连续变化。变化的间隔为无穷小量。(2)对系统的数学描述来说，存在系统输入或输出的微分项（导数项）。(3)系统具有连续的状态。在离散系统中，系统的状态为时间的离散

20、函数，而连续系统的状态为连续量。满足上述条件的系统为连续系统。,2023/3/20,第 52 页,3.3.2 连续系统的数学描述,设连续系统的输入变量为u(t)，其中t为连续取值的时间变量，设系统的输出为y；由连续系统的基本概念可以写出连续系统的最一般数学描述，即,系统的实质为输入变量到输出变量的变换，注意这里系统的输入变量与输出变量既可以是标量（单输入单输出系统），也可以是向量（多输入多输出系统）；而且在系统的数学描述中含有系统输入或输出的导数。,2023/3/20,第 53 页,除了采用最一般的数学方程描述连续系统外，还可以使用连续系统的微分方程形式对连续系统进行描述，即,这里x(t)、x

21、(t)分别为连续系统的状态变量、状态变量的微分。,2023/3/20,第 54 页,对于线性连续系统来说，由连续系统的微分方程描述可以容易的推导出连续系统的状态空间模型。这与使用差分方程对离散系统进行描述类似。下面举例说明连续系统的数学描述。,2023/3/20,第 55 页,实例3-6,对于如下所示的连续系统：,显然此系统为单输入单输出连续系统，且含有输入变量的微分项。由此方程可以容易得出系统的输出变量为,2023/3/20,第 56 页,3.3.3 连续系统的simulink描述,前面给出了连续系统的基本概念与系统的基本描述方法：数学方程描述与微分方程描述。本部分使用例3-6给出的连续系统

22、,说明如何利用simulink对连续系统进行描述，并在此基础上对连续系统进行简单分析。与前面类似，在此并不建立系统的simulink模型进行仿真，而是采用编写M脚本文件从原理上对连续系统进行说明，并进行简单的仿真。,2023/3/20,第 57 页,实例3-7,编写脚本文件systemdemo3.m，对【例3.6】中的连续系统进行分析。%systemdemo3.m脚本文件t=0:0.1:5;%系统仿真范围，时间间隔为0.1 sut=t+sin(t);%系统输入变量utdot=1+cos(t);%系统输入变量的导数yt=ut+utdot;%系统输出plot(yt);grid;%绘制系统输出曲线下

23、图为此连续系统在时间0,5内的输出曲线。由此可见，使用简单的MATLAB语句可对系统性能进行简单的分析。,2023/3/20,第 58 页,连续系统输入输出关系图,2023/3/20,第 59 页,3.3.4 线性连续系统,在介绍线性离散系统时，已经给出线性系统的基本概念，这里做一个简单的回归并介绍线性连续系统的概念。连续系统可以用如下的方式来表达：,2023/3/20,第 60 页,【定义3.5】线性连续系统：如果一个连续系统能够同时满足如下所示的齐次性和叠加性，则称该系统为连续线性系统。(1)齐次性。对于任意的参数，系统满足,(2)叠加性。对于任意输入变量与参数，系统满足,2023/3/2

24、0,第 61 页,3.3.5 线性连续系统的数学描述,线性连续系统最一般的描述为连续系统的输入输出方程形式，也可以使用连续系统的微分方程模型进行描述：,2023/3/20,第 62 页,除了使用这两种连续系统通用的形式描述线性连续系统之外，还可以直接用传递函数、零极点模型与状态空间模型对其进行描述。与线性离散系统相类似，线性连续系统的传递函数模型与零极点模型采用连续信号的拉氏变换来实现。,2023/3/20,第 63 页,拉氏变换具有如下两个性质：(1)线性性。即对于连续信号和，设它们的拉氏变换分别为Lu1与Lu2，则拉氏变换的线性性是指拉氏变换满足下面的关系：,(2)设连续信号u(t)的拉氏

25、变换为U(s)，则u(t)的拉氏变换为sU(s)，u(t)的拉氏变换为s2U(s)。,2023/3/20,第 64 页,对于如下的信号,同时对等式的两边进行拉氏变换，并将其化为分式的形式，可得到如下所示的系统的传递函数模型。,一般来说，线性连续系统的拉氏变换总可以写成如下传递函数的形式。,2023/3/20,第 65 页,将其进行一定的等价变换，可以得出线性连续系统的零极点模型：,其中z1为线性连续系统的零点，p1、p2为系统的极点，k为系统的增益。,2023/3/20,第 66 页,线性连续系统的另外一种模型为状态空间模型。前面已经提到，对于线性连续系统，使用其微分方程很容易推导出系统的状态

26、空间模型。这里给出线性连续系统用状态空间模型进行描述的一般方式：,其中，x(t)为线性连续系统的状态变量，u(t)、y(t)分别为系统的输入与输出变量，可以为标量，也可以为向量。,2023/3/20,第 67 页,3.3.6 线性连续系统的simulink描述,现在介绍如何在Simulink中实现对线性连续系统的描述。一般来说，在simulink中对线性连续系统的描述方式有以下三种：(1)线性连续系统的传递函数模型描述：在simulink中，传递函数表示为num=n0,n1；den=d0,d1,d2；其中num表示传递函数的分子系数向量，den为分母系数向量。,2023/3/20,第 68 页

27、,(2)线性连续系统的零极点模型描述：在simulink中，零极点模型表示为gain=k；zeros=z1；poles=p1,p2；其中gain表示系统增益，zeros表示系统零点，poles表示系统极点。(3)线性连续系统的状态空间模型描述：如果系统的状态空间表示为则在simulink中直接输入变换矩阵ABCD即可。,2023/3/20,第 69 页,一般来说，线性连续系统的不同模型之间可以相互转化，MATLAB中有内置的函数可以完成线性连续系统间的转化。我们在线性离散系统模型间转化中已经做了介绍，这里仅列出这些函数原型：zeros,poles,k=tf2zp(num,den)；num,de

28、n=zp2tf(zeros,poles,k)；zeros,poles,k=ss2zp(A,B,C,D)；A,B,C,D=zp2ss(zeros,poles,k)；num,den=ss2tf(A,B,C,D)；A,B,C,D=tf2ss(num,den)。,2023/3/20,第 70 页,实例3-8,对于如下采用传递函数模型进行描述的线性连续系统,要求绘制此系统的Bode图、Nyquist图，并求取系统的零极点模型与状态空间模型描述。,2023/3/20,第 71 页,实例3-8程序,解：在MATLAB中输入下面的语句：num=1-3;den=2-3-5;w=logspace(-1,1);su

29、bplot(2,1,1)bode(num,den,w)grid on subplot(2,1,2)nyquist(num,den,w)grid on,2023/3/20,第 72 页,线性连续系统的Bode图与Nyquist图,2023/3/20,第 73 页,再输入以下语句求解零点、极点以及增益。zeros,poles,k=tf2zp(num,den)zeros=3poles=2.5000-1.0000k=0.5000,2023/3/20,第 74 页,输入下面语句求解状态空间向量。A,B,C,D=tf2ss(num,den)A=1.5000 2.5000 1.0000 0B=1 0C=0.

30、5000-1.5000D=0,2023/3/20,第 75 页,3.4 混合系统模型及表示,3.4.1 混合系统的数学描述混合系统是由不同类型的系统共同构成的，因此混合系统的数学描述可以由不同类型系统描述共同构成。但是由于混合系统的复杂性，一般难以用单独的数学模型进行描述或表达，因此混合系统一般都是由系统各部分输入与输出间的数学方程所共同描述的，下面举例说明。,2023/3/20,第 76 页,实例3-9,对于如下的一个混合系统：设系统的输入为一离散变量，系统由离散系统与连续系统串联构成，其中离散系统输出经过一个零阶保持器后作为连续系统的输入。其中离散系统的输入输出方程为且，系统采样时间为Ts

31、=1s。,2023/3/20,第 77 页,连续系统的输入输出方程为,由于此混合系统中离散系统的输出经过一零阶保持器后作为连续系统的输入，因此“与”的数学关系为,2023/3/20,第 78 页,其中Ts=1s为离散系统的采样时间。故此混合系统的输入与输出之间的关系可以由下面的方程来描述：,2023/3/20,第 79 页,3.4.2 混合系统的simulink描述与简单分析,在对单独离散系统或连续系统进行描述时，由于系统一般比较简单，因而可以采用诸如差分方程、传递函数、状态空间等模型表示。但对于混合系统，由于系统本身的复杂性，即使是很简单的混合系统，如【例3-9】给出的例子，都难以用一个简单

32、的模型进行描述。因此，这里采用简单的数学方式对系统进行描述与分析。,2023/3/20,第 80 页,例3-10,编写M脚本文件systemdemo4.m，对【例3-9】中的混合系统进行分析。t=1:0.1:99.9;n=1:100;un=0.5*n;yn=un+1;for i=1:length(n)-1 for j=1:length(t),2023/3/20,第 81 页,if t(j)=n(i),2023/3/20,第 82 页,混合系统输入与输出关系,2023/3/20,第 83 页,从系统的输出曲线图中可以看出：由于系统中离散部分的输出经过零阶保持器后作为连续部分的输入，而零阶保持器具有阶跃的特性，在系统仿真结果中出现阶跃现象。另外，系统呈现类似正弦发散的特征表明系统为一发散不稳定系统。,2023/3/20,第 84 页,本章对动态系统做了简单的介绍，其中涉及到简单系统、离散系统、线性离散系统、连续系统、线性连续系统、混合系统等系统的概念及数学描述和simulink描述，并且使用简单的MATLAB语句对不同的系统进行了简单的仿真与分析。至于使用simulink进行动态系统建模、仿真与分析将在下面章节部分进行详细介绍。,2023/3/20,第 85 页,谢谢！,