ch2.3 连续系统的离散相似法

合集下载

连续系统模型的离散化处理方法

连续系统模型的离散化处理方法
只要T不变,三个系数均不变,可以在仿真前预先计算好,这样就减少了以后的计算工作量。 加入一个理想滤波器,保留输入信号主频段,滤掉附加的频谱分量,不失真
在离散化后,模型精度变差,可能不稳定。
S域到Z域的最基本映射关系是:Z=e (T— TS 数值积分法:将微分方程转换成差分方程,这中间是一步步离散,每一步离散都用到连续系统的原模型,这样的速度就慢了。
TeAT
m T
T eATA Bd
0
xKTTTxKTmTUKT
x(k1) TxkmTUk
B 当输入函数u(KT)在两采样 点间线性变化时(一阶保持)
uuKTukT
p
T
TeATABd
0
xkTTTxkTmTUkTpTUkT
xk1TxkmTUkpTUk
当连续系统状态方程系数A、B已知时,
可求出……
此法相比于数值积分法;只要T不变,三个系 数均不变,可以在仿真前预先计算好,这样 就减少了以后的计算工作量。
2 典型环节的离散状态方程
A 积分环节:G(S)=K/S f1=x2 ; f2=x3 ;
依据各环节的连接关系及外部作用函数 稳定性不及双线性替换法,Ts或信号重构器选择不当,离散模型的稳定性变差
二、Z域离散相似方法
1 基本方法
G z
y z u z
z G h s G s
1
z
s a
z exp( aT )
e TS 1 z
1 z
s
z 1
1
Tz
s* s
( z 1 )( z 1 )
Gz
yz uz
zGh
sGs
Gs k
sa
Gh
s
1

离散相似法

离散相似法

3.1时域离散相似法的基本原理
仿真计算过程:

原连续系统 离散后系统
输入信号 u(t) uh(kT)
2012-12-17
3-0 离散控制系统基本原理

3-0 离散控制系统基本原理

4 离散控制系统:

5 x(t) 和 x(kT) 的数学描述

系统中既含有连续信号[x(t), e(t), u(t)],又含有采样 信号[x*(t)]或者数字信号[x(kT), u(kT)]的系统。 离散控制系统是由连续的控制对象、离散的控制器、 采样器和保持器等几个环节所组成。
控制系统建模与仿真 rjliu@ 19
1.连续系统的离散化:微分方程-》差分方程
Ax Bu, x y Cx Du x (0) x0 , t 0
2012-12-17
2012-12-17
控制系统建模与仿真 rjliu@
20
对系统进行离散化:(T为采样周期)
控制系统建模与仿真 rjliu@ 5 2012-12-17 控制系统建模与仿真 rjliu@ 6

2 采样信号 x*(t)


3 数字信号 x(kT)

2012-12-17
北京交通大学
1
控制系统建模与仿真 rjliu@
3-0 离散控制系统基本原理




数值积分的计算,是根据计算步长h,一步一步的逐步计 算,得到t1,t2, …,tk (采样点)的数值x1, x2,…, xk 。 对原连续的系统进行了离散化处理。求解时,用一个等价 离散模型,代替原连续系统。 在自动控制领域,随着计算机的引入和广泛应用,形成了 一套系统的理论——计算机控制(离散控制,数字控制)。 我们可以从连续系统离散化的角度出发,用控制理论中采 样和信号重构的理论,建立离散模型,并用计算机控制的 理论分析计算。 离散后的系统和原来的连续系统是相似的(输入输出信号 在采样点上是近似相等的),所以叫做离散相似法。

4.3离散相似法

4.3离散相似法

x( n
1)
e aT
x( n
)
k a
(1
e aT
)u( n )
k a2
( eaT
1)
k a
T u( n
)
y(n 1) (b a)x(n1) ku( n 1)
2022/1/12
33
第33页,共40页。
三、状态转移矩阵的计算
(1)解析计算(拉氏逆变换)
利用定义 (T ) eAT L1[(sI A)1] 进行求解
2022/1/12
12
第12页,共40页。
h0(t)
h0(t)
1
1
0
Tt
0
Tt
-1
a
b
零阶保持器的特性:
g0 (t) 1(t) 1(t T )
➢ 对采样值既不放大, 也不衰减
➢只能不增不减地保 持一个采样周期
2022/1/12
保持器的传递函数:
1 eTs 1 eTs G0 (s) s s s
b
k
/
s
系统模拟结构图: u
s
y
b
状态空间模型: X kU Y X bU
A 0, B k
离散状态空间模型:x(n 1) x(n) kTu(n) 0.5kT2u(n)
y(n 1) x(n 1) bku(n 1)
2022/1/12
31
第31页,共40页。
(3)惯性环节
传递函数: G( s ) Y ( s ) k k 1 / s U(s) s a 1 a / s
0 1
0 0
1 1
s 0
1 s 1
(sI
A)1
1 s(s 1)

第6章连续系统的离散化方法及近似解

第6章连续系统的离散化方法及近似解

第6章连续系统的离散化方法及近似解在连续系统中,我们经常需要将其离散化为离散系统以便于分析和求解。

离散化方法能够将连续系统的微分方程转化为差分方程,从而得到近似解。

本章将介绍连续系统的离散化方法及近似解的计算。

连续系统的离散化方法有许多种,常见的有Euler方法、Runge-Kutta方法和有限差分方法等。

其中,Euler方法是最简单和最基础的离散化方法,其基本思想是将连续时间轴划分为若干个小时间间隔,并用差分逼近连续系统的导数。

具体地,对于一阶常微分方程:\[\frac{{dy}}{{dt}} = f(y, t)\]可以使用Euler方法将其离散化为:\[y_{n+1} = y_n + h \cdot f(y_n, t_n)\]其中,\(y_n\)是时间点\(t_n\)的近似解,\(h\)是时间步长。

Runge-Kutta方法是一种更精确的离散化方法,其基本思想是利用多个中间步骤来更准确地逼近连续系统的导数。

常见的是四阶Runge-Kutta 方法,其公式为:\[y_{n+1} = y_n + \frac{h}{6} \cdot (k_1 + 2k_2 + 2k_3 +k_4)\]其中\[k_1=f(y_n,t_n)\]\[k_2 = f(y_n + \frac{h}{2}k_1, t_n + \frac{h}{2})\]\[k_3 = f(y_n + \frac{h}{2}k_2, t_n + \frac{h}{2})\]\[k_4 = f(y_n + hk_3, t_n + h)\]这样可以得到更准确的近似解。

有限差分方法是一种常用的离散化方法,其基本思想是将连续的导数用差分逼近。

以二阶偏微分方程为例,该方程的一般形式为:\[\frac{{\partial^2u}}{{\partial x^2}} +\frac{{\partial^2u}}{{\partial y^2}} = f(x, y)\]可以使用中心差分公式将其离散化为:\[\frac{{u_{i+1,j} - 2u_{i,j} + u_{i-1,j}}}{{\Delta x^2}} + \frac{{u_{i,j+1} - 2u_{i,j} + u_{i,j-1}}}{{\Delta y^2}} =f_{i,j}\]其中,\(u_{i,j}\) 是近似解在网格点 \((i, j)\) 处的值,\(\Delta x\) 和 \(\Delta y\) 分别是网格在 \(x\) 和 \(y\) 方向的步长,\(f_{i,j}\) 是离散化后的右侧函数。

连续控制系统的离散相似-可调参数法仿真

连续控制系统的离散相似-可调参数法仿真

连续控制系统的离散相似-可调参数法仿真
许伟明
【期刊名称】《华东工业大学学报》
【年(卷),期】1995(17)2
【摘要】本文讨论的连续控制系统仿真方法是基于离散相似和可调数值积分法相结合的一种快速数字仿真方法.文中推导了基本仿真模型,并讨论了系统仿真的实现.实例表明,此方法是可行的.
【总页数】7页(P38-44)
【关键词】数字仿真;离散相似;可调参数;连续控制系统
【作者】许伟明
【作者单位】华东工业大学仪器仪表学院
【正文语种】中文
【中图分类】TP15
【相关文献】
1.离散相似法在卫星姿态控制系统设计和数字仿真中的应用 [J], 周文忠;李果
2.连续系统离散相似法仿真的排序算法 [J], 李胜勇;常谦
3.提高连续系统离散相似法仿真精度的实用方法 [J], 王福永
4.用离散相似法仿真控制系统动态响应的研究 [J], 孙光东
5.连续控制系统的离散相似-可调参数法仿真 [J], 许伟明
因版权原因,仅展示原文概要,查看原文内容请购买。

第三章、时域离散相似法

第三章、时域离散相似法
1
0 1 1 1 e T s 1
0 e T
m (T ) (T ) Bd
0
T

T
0 T
1 1 e ( T )
1 d (T ) e 0 0
0
现举例说明如何求连续 系统的离散化模型。 例4.1.1 设有一控制系统的状态 空间表达式为: x1 0 0 x1 1 x 0 u x2 1 1 2 x1 y 0 1 x2 求用零阶保持器采样时 系统的离散化模型。 解:
当初值x(0)和u (t )给定以后,可以利用式 3.1.12 ( ) 递推计算出系统不同采 样点上的状态变量和输 出 变量的数值。 因为零阶保持器的输出 在两个采样时刻之间保 持 常量,又称矩形近似。 由图( )可见,输入函 3.2 数u (kT )的误差为u。为了提高计算精度, 可采用 三角形保持器(一阶保 持器)。根据图( )所 3.2 示,可得两相邻采样点 之间的输出为: u ( ) u (kT ) u[(k 1)T u (kT )] T u (kT ) u (kT )
u(t )
T
u * (t )
保持器
x Ax Bu y Cx Du
y(t )
T
y * (t )
图4.1
连续系统时域离散化
对状态空间表达式进行拉普拉斯变换,得:
sX ( S ) X (0) AX ( s) BU ( s) Y ( s) CX ( s) DXU( s) (3.1.2) (3.1.3)
离散状态方程和输出方 程为: KT 2 x(n 1) x(n) KTu (n) u ( n) 2 y (n 1) x(n 1) (3.2.1)

连续系统的离散化方法课件

连续系统的离散化方法课件

离散化方法的意义
精确性
离散化方法可以提供对连续系统的精 确近似,特别是在计算机仿真和数字 控制系统中。
可计算性
离散化方法可以将不可计算的分析转 化为可计算的形式,便于进行数值计 算和控制器设计。
离散化方法的应用场景
01
02
03
数字控制
在数字控制系统中,连续 系统的离散化是必要的步 骤,以便在数字计算机上 进行数值计算和控制。
小波基选择
常用的小波基包括Haar小波、Daubechies小波、Morlet 小波等。
误差分析
小波变换法的误差主要来自于变换误差和离散化误差。
05
离散化方法的评估与优化
评估离散化方法优劣的标准
01
02
03
04
精度
离散化方法是否能准确代表原 连续系统。
稳定性
离散化方法在一定参数变化范 围内是否能保持稳定。
状态空间模型
用状态变量和输入、输出变量描述连续系统的动态特性。
状态空间模型通常形式为:`x'(t) = Ax(t) + Bu(t)` 和 `y(t) = Cx(t) + Du(t)`,其中 `x(t)` 表 示系统状态,`u(t)` 表示系统输入,`y(t)` 表示系统输出,`A`, `B`, `C`, `D` 是系数矩阵。
化率。
通过求解 ODE,可以得到系统 在任意时刻的状态。
传递函数
表示连续系统在输入和输出之间的传递 特性。
传递函数通常形式为:`G(s) = Y(s) / U(s)`,其中 `Y(s)` 和 `U(s)` 分别是输 出和输入的拉普拉斯变换,`s` 是复变
量。
通过分析传递函数的零点、极点和增益 ,可以得到系统的稳定性和性能特性。

ch2.3 连续系统的离散相似法

ch2.3   连续系统的离散相似法
7
很明显,理想信号保持器是无法实现的。 很明显 理想信号保持器是无法实现的。实际上能实现的 理想信号保持器是无法实现的 是各种近似的保持器。 是各种近似的保持器。
u(kt)
~ (t ) y
保持器将输入的离散采样值转换为连续模拟信号。 保持器将输入的离散采样值转换为连续模拟信号。
At 0
15
t
x(t ) = e x(0) + ∫ e A ( t − τ ) Bu( τ)dτ是连续系统的解析解,
At 0
t
对这个解进行离散化。考察kT和(k + 1)T两个相邻采样时刻 状态向量x(t )的值。将t = (k + 1)T代入上式得, x[(k + 1)T] = e = e [e
10
5、三角形保持器
y[( k + 1)T ] − y ( kT ) y ( t ) = y ( kT ) + ( t − kT ) kT ≤ t ≤ ( k + 1)T T 特点: 特点:要求在计算 kT ~ ( k + 1)T 的 h ( t )时,需已知 y[( k + 1)T ]
三角保持器的 单位脉冲过渡 函数
9
零阶保持器的特点: 零阶保持器的特点:
其特性与低通滤波器相似; ① 其特性与低通滤波器相似; 有相位滞后问题(滞后T/2); ② 有相位滞后问题(滞后 ); 若输入信号为阶跃信号,那么零阶保持器能无失真地恢复信号。 ③ 若输入信号为阶跃信号,那么零阶保持器能无失真地恢复信号。
4、一阶保持器
一阶保持器使用当前时 进行外推 , y ( kT + ∆ t ) = y ( kT ) + y ( kT ) − y[( k − 1)T ] ∆t T

7.离散相似法仿真

7.离散相似法仿真

7.2
典型环节的离散模型
按照前面的讨论,我们将常见的典型环节由传递函数表 达式推导出其离散系数及离散状态方程。
7.2.1 典型线性动态环节 ˆ (T ) Φ Φ m (T ) , Φ(T ) , m 的计算
• 典型线性动态环节有:积分、比例积分、惯性、超前-滞 后、比例五种。均以(Ci+Dis)/(Ai+Bis)形式描述。根据每种 情况的Ai,Bi,Ci,Di取值,确定分别所对应的线性环节 类型,进而计算出该环节的状态转移矩阵的值,包括下面 几个。
离散相似法仿真
本章主要教学内容
离散相似法原理 典型环节的离散模型 线性系统离散相似法仿真 非线性系统离散相似法仿真 采样系统仿真分析
7.1
7.1.1 仿真算法描述
离散相似法原理
所谓离散相似法,就是将一个连续系统进行离散化处理, 从而得到等价的系统离散模型,此种方法按系统的动态结 构图建立仿真模型。 在计算过程中,按各典型环节离散相似模型,根据环节 的输入来计算环节的输出。 1. 环节的离散化模型 将连续系统按图7-1所示对其进行离散化处理,在系统的 输入、输出端加上虚拟采样开关,T为采样周期。为保证输 入信号复现原信号,在输入端加上一个保持器。
Ts = 10 T
如果给定系统开环截止频率时,系统的采样周期也可 以按下式进行选择:
1 Ts = (30 ~ 50)ω c
2.
保持器特性对仿真精度的影响
为使经采样后的信号无失真地复现,要在系统中加入保持器。虽然零
阶保持器比较容易实现,但其精度较低。为了提高控制精度,可以采用 三角保持器,它复现信号的高频部分失真较小,并且无相位滞后,可以 得到比较满意的结果。 此外,为了提高精度还可以采用校正补偿措施,在离散模型中加入 一个确定的校正环节,适当调整参数,可使离散模型尽可能地接近原型。

离散相似法

离散相似法

第五章 连续系统离散相似法数字仿真上一章讲的数值积分法,其结果得到的是递推公式。

其出发点不是建立一个与被仿真的连续系统等价(某种误差意义下)的离散模型,而是建立在数值计算的基础上,即对积分进行数值计算(某种求和计算)。

或从另一角度说,对微分算子s(以一阶微分方程为对象) 找更好的替代方法。

所得的公式,只是这种思路的结果。

也就是说,人们为了计算方便,将其设法整理成了递推公式,递推公式在固定步长的情况下,就变成了通常所用的常系数的差分方程,可用z变换与频域联系起来(也正因为这点,使其结果从形式上与离散相似累积法类似)。

对这种方法精度的分析,是从截断、舍入误差累积角度进行的。

其改进思路,是减小这两种误差。

时域看失真。

可为线性系统,也可为非线性系统。

离散相似法的出发点与数值积分法不同。

它的出发点在于:采样系统可用差分方程 和z变换来描述,而用计算机求解差分方程式很方便的(差分方程可看作多步法递推公式)。

所以,设法找到一个采样系统,使其输入—输出值与被仿真的连续系统在采样时刻的值近似相等。

这个采样系统称为连续系统的离散等价模型(二次模型、仿真模型)。

从此模型,可得出原系统的一系列值。

“等价”程度的高低,就是此法的精度问题。

所以精度分析从保持器引入的幅度与相位失真进行。

频域看失真。

只有线性时,才可列出整个系统的差分方程。

上述二法从结果(即所得递推公式)看有某种相似性。

所以精度稳定性等的解释可互换借鉴。

改进思路的不同,可看作看问题的角度不同(或出发点不同)而已。

因角度等价手段不同,各自能达到的“等价”程度也就有差别。

小结:1.数值积分法:出发点: 对积分进行数值计算,对微分算子s找更好的替代方法。

结果: 递推公式。

定步长时,为差分方程。

精度分析: 从截断、舍入误差累积角度进行。

时域看失真。

2.离散相似法:出发点:出发点与数值积分法不同。

用计算机解差分很方便,找一个采样系统,使其I-O值与被仿真连续系统(在采样时刻)近似相等。

第6章 连续系统的离散化方法及近似解

第6章 连续系统的离散化方法及近似解

K - 2M 0
EI 2 39.4784 Al 4
解得
EI 3 68.9944 Al 4
正则特征向量
a (1) 0.5742 2 0 Al 0.0048
T
a (2)
0 2 (3) 1 a Al 0
0
l
设梁上分别受到分布力f(x,t)和 x xd 处的集中力F(x,t)
当梁上有虚位移
l 0
w( x, t ) i qi 外力虚功为
i 1
n
W f ( x, t ) F (t ) ( x xd ) w( x, t )dx
l f ( x, t )i ( x)dx F (t )i ( xd ) qi 0 i 1 n
得到原来问题的模态向量
3 x ( x) sin 0.0681sin 2l 2l
(1)
x
(2) ( x) 0.1955sin
x
2l
sin
3 x 2l
例:等截面简支梁中部有集中质量,并受有集中力 设集中质量 M a 等于梁的质量
集中力的变化的频率
50 EI / Al 4
Kij EI ( x)i( x) j ( x)dx
0

3 0 2 Al M 0 1 0 2 2 0 3
1 0 0 EI K 0 16 0 2l 3 0 0 81
4
代入本征方程
EI 1 5.6825 Al 4
T
0.5199 2 0 Al 0.7746
T
求梁的响应时,将位移写作假设模态的线性组合

计算机仿真教案04-2-第四章 离散相似法的连续系统仿真

计算机仿真教案04-2-第四章 离散相似法的连续系统仿真
)脉冲传递函数模型
离散系统
连续系统
差分方程 Z变换 脉冲传递函数
微分方程 拉氏变换 传递函数
脉冲传递函数的定义 在连续系统中、应用拉氏变换可将描述系统的
微分方程转化为传递函数。同样,在采样系统中, 利用Z变换可将描述采样系统的差分方程转化为 类似于传递函数的另一种数学模型一脉冲传递函 数,或称Z传递函数。脉冲传递函数的定义如下:
( TeA(Tt)Bt)d u'(tkT ) 0
所以,x(n+1)= Φ (T)x(n)+Φ m(T)u(n)+Φ p(T)u’(n)
(3-2-9)
离散化状态方程-系数计算
(T ) e AT
m (T )
T e A (T t ) Bdt
0
p (T )
T e A(T t ) Btdt
x(n+1)=Φ (T)x(n)+Φ m(T)u(n) 加零阶保持器的离散化状态方程
(3-2-8)
加一阶保持器的离散化状态方程
如果采用一阶保持器,那么,
u(kTt)u(kT)u(kT)u[k(1)T]t T
u(kT)u'(kT)t
代入(3-2-7),
x[k (1)T]eAT x(kT )( TeA(Tt)Bd )u(tkT ) 0
G (Z ) Z { 1 e sT G (s) }(1 z 1 )Z { G (s)}

T(
( AT ) i ) B
0 (i 1)!
p ( T )
T e A ( T t ) Btdt
0
令 =T-t,则有,
P ( T )
T e AB(T )d

连续控制器离散化方法

连续控制器离散化方法

模拟(连续)控制器系统
计算机(离散)控制器系统
离散控制器等效控制系统
采用连续与离散控制器的系统系阶越响应的区别
例:已知某系统被控对象的传递函数为 要求设计控制器,使满足性能指标: ①闭环稳定 ②过渡过程时间Ts≤3s ③阶跃响应超调量δ≤5% 设计满足上述要求的数字控制器D(Z)(取采样周期 T=0.2秒,采用双线性近似法) 解: 模拟控制器设计过程略,得到的模拟控制器为:
D( s ) 16 s 2.1 s 8
( k 1)T kT
u ( ) d
T [u (( k 1)T ) u ( kT )] 2
T z 1 u ( kT ) 2 z 1 1 T z 1 2 z 1 , s s 2 z 1 T z 1 2 z 1 ) C ( Cd ( z ) C ( s ) T z 1
Cd ( z) 的极点为
1 aT 1
稳定
(3)Tustin法
aT aT a z 1 2 Cd ( z ) 2 2 z 1 aT aT aT aT a (1 )z ( 1) 1 1 T z 1 2 2 2 z 2 aT 1 2 ( z 1)
( s b1 )( s b2 ) ( s bm ) C (s) k ( s a1 )( s a2 ) ( s an ) ( z 1)d 1 ( z eb1T )( z eb2T ) ( z ebmT ) Cd ( z ) kd ( z e a1T )( z e a2T ) ( z e anT )
例:分别用前向差分、后向差分、Tustin法对 进行离散化 (1)前向差分
a aT z 1 a z 1 aT T Cd ( z) 的极点为 1 aT Cd ( z )

连续系统离散化

连续系统离散化
连续系统的数字仿真
离散相似法
连续系统的离散化
首先要得到一个与被仿真系统等价的离散 模型。这个模型可以通过对连续系统的离 散化过程来获得。它分成以下五步:
① 首先对输入信号u(t)进行采样,即在输入 端加一个采样开关S1,其采样周期为T。
② 连续变化的信号u(t)经过采样开关后,变 成了一个离散信号u(kT)。为保证模型的等 价性,首先要求信号等价,因此它不能直 接进入原来的连续系统,而必须加上一只 信号重构器,它使信号u(k)重新变成一个 连续信号uh(t), uh(t)u(t)。
环节2是一阶惯性环节,其传递函数为
X (s) k0 k0 T1 U (s) T1s 1 s 1 T1 x(k 1) eT T1 x(k ) k0 (1 eT T1 )u(k )
连续系统按结构图的离散相似法仿真
环节3也是一阶惯性环节,其传递函数为
Y (s) 1 1 T2 X (s) T2s 1 s 1 T2
对上式进行拉氏反变换,求得方程的解为
t
X (t) (t) X (0) 0 (t )BU ( )d
时域解法求取离散系统差分方程
(t) L1 (sI A)1
X (t) eAt X (0) t eA(t )ΒU ( )d 0 t
X (t) (t)X (0) 0 (t )BU ( )d
采用一阶保持器描述系统的差分方程
x(k 1)T (T )x(kT) m (T ) n (T )u(kT) n (T )u(k 1)T
e(k 1) R(k 1) y(k 1)
e(k) R(k) y(k)
x1(k
1)
x1(k
)
3 4
Te(k
)
1 4
Te(k

自动控制系统仿真离散相似法算法的改进

自动控制系统仿真离散相似法算法的改进

自动控制系统仿真离散相似法算法的改进自动控制的发展已经成为现代工程工作的一个必要部分,系统仿真也成为了许多项目的重要环节。

然而,仿真的准确性和效率却成为了人们关注和追求的目标。

离散相似法是目前仿真中常用的一种方法,而其算法的改进也一直是研究的焦点。

本文介绍离散相似法算法改进的相关信息,包括原理、目的、方法和效果等方面。

1. 原理离散相似法是一种将实际工程中的连续系统仿真成离散系统的一种方法,是工程仿真中最常用的方法之一。

离散相似法的基本原理是利用与实际系统的定义特性相同的代理模型以及在系统状态变化点上的数据来重构系统。

其原理如下:(1) 基于实际系统定义特性的代理模型代理模型是实际系统的模拟模型,在仿真中被用来模拟实际系统的动态行为。

它的定义特性应与实际系统的定义特性相同,且其输入输出特性应能仿真实际系统的输入输出行为。

通常,这种输入输出特性可能包括时间延迟、频率响应、非线性特性以及各种噪声特性等内容。

(2)基于系统状态变化点上的数据的重构操作系统状态变化点是指仿真过程中离散系统状态向连续系统状态的转化点。

在这些状态变化点处,数据被采样并分析,被用来确定仿真过程中的离散系统参数,包括时延、频率响应等内容。

这些参数被用来重构仿真过程中的连续系统,并最终用于评价仿真结果。

这种操作可以保证仿真结果与实际系统间的相似性并提高仿真的准确度和效率。

2. 目的离散相似法算法的改进,主要是为了提高其仿真效率和准确性,是一个基于最优控制方法的研究方向。

主要目的包括:(1) 改进仿真精度离散相似法算法需要保证仿真精度,在模型参数的测定以及模型拟合方面采用多种优化控制方法,从而实现最优参数测定和最优模型拟合,减少误差的产生,达到最佳仿真效果。

(2) 优化计算过程随着计算技术的发展,优化计算过程成为了提高仿真速度的重要手段。

离散相似法算法的改进直接关系到仿真计算的速度和效率。

因此,在研究离散相似法的改进时,优化计算过程也成为了关注的焦点。

第四章连续系统的离散化方法

第四章连续系统的离散化方法
Ki hf ( xk ij K j , tk i h), i 1, 2,
j 1 i 1
,6;
j 1, 2,
6
i 1
下一步的状态变量可由下式求出:
xk 1 xk i Ki
i 1
四阶/五阶龙格-库塔法系数表
i* ij i
0 16/135 25/216
一般形式
3、a1 a2 , b1 b2 1 则 2
h x1 x0 ( K1 K 2 ) 2 K1 f (t0 , x0 )
一般形式
1
K 2 f (t0 h, x0 hK1 )
h xk 1 xk ( K1 K 2 ) 2 K1 f (tk , xk ) K 2 f (tk h, xk hK1 )
x1 x0 hf (t0 , x0 )
其一般公式为
xk 1 xk hf (tk , xk )
f1
f
f0
c
0
t0
t1
t
h2 h3 (2) h k ( k 1) Rn f (t0 , x0 ) f (t0 , x0 ) f (t0 , x0 ) 称为截断误差 2! 3! k! 例4-1 用欧拉法求下述微分方程的数值解。
x2 0.9 0.1 0.92 0.819 x3 0.819 0.1 0.8192 0.7519 x10 0.4627810
1 x 上式的精确解是 1 t
t 0 0.1
数值计算与精确解的比较见表
0.2 0.3 1.0
精确解x(t)
1
1
0.9090909
0.9
0.8333333
xi ”的固定顺序输入i 1, 2,

钱雪军-第4章 连续系统仿真——离散相似法2

钱雪军-第4章 连续系统仿真——离散相似法2

4.2.2 保持器
系统辨识与仿真
通过保持器对采样后输入信号的重构,通常可采用三种类型的保持器。
◨ 零阶保持器
◨ 一阶保持器
◨ 三角保持器
▶三种保持器比较:零阶保持器最简单,但重构信号的误差较大; 三角保持器最复杂,恢复信号的失真最小。
▶常用的还是零阶或一阶保持器,特别是零阶。
零阶保持器
u(t)
t0
u(t)
u*(t)
保持器
T
uk(t) x(t) = Ax(t) + Bu(t) x(t)
x*(t)
T
图4.2 状态空间模型离散化原理结构图
4.2.1 采样定理
系统辨识与仿真
采样频率ωs大于或等于两倍的采样器输入连续信号e(t)频谱 中的最高频率ωmax。即ωs>=2ωmax,这就是香农采样定理。
抽样频率小于模拟信号最好频率的2倍会造成频谱混叠 。
比较G(z)与G(s),可以得到置换关系 s ≈ 2 ( z −1) T (z +1)
映射关系
系统辨识与仿真
1+ Ts
考虑到
z

1

2 Ts
将 s = σ + jω 代入上式,可得
2
1+ Ts 1+ T (σ + jω) 1+ T σ + j T ω
(1+ T σ )2 + (T ω)2
z

2 1− Ts
◆为了使输入信号u(t)离散化后仍能保持原来的变化规律, 在输入采样开关后设置信号保持器(亦称为信号重构器), 复现原输入信号u(t),其结构如图4.1。
u(t)
u*(t) 保持器

连续系统的离散化方法及近似解课件

连续系统的离散化方法及近似解课件
差分方程
离散化后的控制系统可以用差分方程来描述,差分方程是连续时间微分方程在离散时间域 上的对应形式。通过求解差分方程,可以得到离散控制系统的输出响应。
Z变换
Z变换是离散时间信号和系统分析的重要工具,它可以将差分方程转换为代数方程,从而 简化离散系统的分析和设计。
电路模拟中的离散化方法及近似解应用
离散系统
离散系统是指系统状态在时间上 是离散的,即系统的状态变量只 在某些特定的时刻有定义,且在 这些时刻间不发生变化。
连续系统与离散系统的区别与联系
区别
连续系统和离散系统最主要的区别在于时间的连续性。连续系统的时间变量是连 续的,而离散系统的时间变量是离散的。
联系
两者之间存在密切的联系。实际上,许多连续系统可以通过离散化方法转化为离 散系统进行处理,这是因为数字计算机在处理问题时,只能处理离散的时间信号 。反之,离散系统的某些理论和方法也可以用来处理连续系统。
连续系统的离散化方法 及近似解课件
目 录
• 连续系统与离散系统概述 • 连续系统的离散化方法 • 离散系统的近似解法 • 连续系统离散化及近似解的应用案例 • 实验与仿真
01
连续系统与离散系统概述
连续系统与离散系统的定义
连续系统
连续系统是指系统状态在时间上 是连续的,即系统的状态变量在 任何时刻都有定义且随时间连续 变化。
感谢观看
前向差分法:前向差分法使用当前时刻及其前一时刻的输入信号来近似 计算下一时刻的输出信号。这种方法简单直观,但离散化误差相对较大 。
后向差分法:后向差分法使用当前时刻及其下一时刻的输入信号来近似 计算当前时刻的输出信号。相比前向差分法,后向差分法具有较小的离
散化误差。
以上内容即为连续系统的离散化方法及近似解课件的部分内容。在实际 应用中,可以根据具体需求和场景,选择合适的离散化方法和参数,以 实现连续系统的高效、准确离散化处理。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

保持器
1、采样周期
的采样开关, 假设有一连续信号 e( t ), 若将它通过一个采样周 期为T的采样开关, 得到的离散信号 e* (t ), 可以证明 1 ∞ e* (t ) = ∑ e(t )e jnωs t , 其中ωs = 2π / T为采样频率 T n = ∞ 1 ∞ * 拉氏变换,并使用平移原理 并使用平移原理: 拉氏变换 并使用平移原理: E ( s ) = ∑ E ( s jn ω s ) T n = ∞
2.3 连续系统的离散相似法
数值积分方法进行系统仿真的不足之处: 数值积分方法进行系统仿真的不足之处: 单步法求解过程中,每计算一个步长h (1) 单步法求解过程中,每计算一个步长h,要多次求取函 数的导数值,步骤烦琐。 数的导数值,步骤烦琐。 多步法求解过程中,又要求存储各状态变量值前r (2) 多步法求解过程中,又要求存储各状态变量值前r次时 刻的数据,系统阶次较高时,存储量相当大, 刻的数据,系统阶次较高时,存储量相当大,而且启动时还需 要其他算法配合。 要其他算法配合。 隐式算法求解必须经若干次迭代, (3) 隐式算法求解必须经若干次迭代,才取得一个时刻的变 量数值,计算速度受影响。 量数值,计算速度受影响。 虽能得到各线性环节的输出响应值, (4) 虽能得到各线性环节的输出响应值,但由于数值积分方 法本身的原因所限,它是统一由状态方程求解变量值, 法本身的原因所限,它是统一由状态方程求解变量值,对单个 环节输出的特殊变化(如非线性变化)难以单独考虑. 环节输出的特殊变化(如非线性变化)难以单独考虑.
-T T
hg(t)
1 t
传递函数为 : e H (s ) = T
Ts
1 e s
Ts
2
滞后一拍的 三角保持器
y ( kT ) y[( k 1)T ] y ( kT + t ) = y[( k 1)T ] + t T
(1 e ) H (s ) =
Ts 2
Ts 2
(试比较一阶保持器)
13
2.3.2 连续系统状态方程的离散化
u(t)
u(kt) 保持器
x(t ) = Ax + Bu y = Cx + Du
y(kt)
对输入为u(t),输出为y(t)的连续系统。使用采样周期为T的 虚拟开关将输入、输出离散化,要求输出y(kt)在采样时刻 的值等同于原输出y(t)在同一时刻的值。
14
对状态方程
Ts 2
刻和前一时刻两个采样
值为基础 0 ≤ t ≤ T
Ts
1 e 传递函数为 H1 (s ) = T (1 + Ts ) Ts
1 + Ts 1 e = T s
2
一阶保持器当满足采样定理时,可无失真地重构斜坡信号。 一阶保持器当满足采样定理时,可无失真地重构斜坡信号。 斜坡信号
称为状态转移矩阵
X ( s ) = L [φ ( t )] X ( 0 ) + L [φ ( t )] BU ( s )
拉氏反变换,并利用卷积公式,可得: 拉氏反变换,并利用卷积公式,可得:
x(t ) = φ (t ) x(0) + ∫ φ (t τ ) Bu (τ )dτ
0
t
即:
x(t ) = e x(0) + ∫ e A(t τ ) Bu (τ )dτ
AT
所以, 所以,有 x[(k + 1)T ] = e =e
AT
x (kT ) + e
AT

( k +1)T
kT
e A( kT τ ) Bu(τ )dτ
x (kT ) + ∫
( k +1)T
kT
e A[( k +1)T τ )] Bu(τ )dτ
16
做不同的处理, 对输入向量 u ( τ )做不同的处理,可以 得到不同的时域等价离 散化模型
用 j ω 代替 s ,可得
1 E ( jω ) = T
*
n = ∞
∑ E [ j (ω n ω

s
)]
的频谱, 若 E ( j ω ) 为 e ( t )的频谱,则离散信号 为一周期性频谱, 为一周期性频谱,其周 期为 ω s
e * ( t )的频谱 E * ( j ω )
(什么是频谱?)
4
为信号频谱最高频率。 ω m 为信号频谱最高频率。
22
r(t)
连续系统模型
y (t )
r(t)
虚拟采样 开关
r * (t )
保持器
r(t) 连续系统
模型
y (t )
y * (t )
保持器的目的是使 输入信号在采样间 隔仍保持连续性
从频域看,保持器的作用就是把因离散化产生的高频分量滤去, 从频域看,保持器的作用就是把因离散化产生的高频分量滤去,保留 主频部分。 主频部分。 但是由于保持器并非是理想滤波器, 但是由于保持器并非是理想滤波器,因此它只能在一定程度上近似 3 还原”连续信号。 “还原”连续信号。
5香农Βιβλιοθήκη 样定理: 香农采样定理:对一个具有有限频谱的连续信号 进行采样, 进行采样,只要选择 ω s
2 ≥ ωm
(ω m < ω < ω m )
理想的低通滤波器 ,通过理想的低通滤波器,就能把 通过理想的低通滤波器,
原来信号毫无失真地提取出来。 原来信号毫无失真地提取出来。 (信号无失真复原的两个条件)
6
2、理想的信号保持器
若 ω s / 2 ≥ ω m , 则离散信号的频谱没有 重叠现象 , 在用一个信号 保持器完全恢复连续信 号时 , 要求该信号保持器的频 率特性是 锐截止的 .
Gh(jw)
理想信号保持器的频率特性
-ws/2 w
ws/2
保持器是一个低通滤波器 从频域观点看,保持器是一个低通滤波器 它用来衰减因离 频域观点看 保持器是一个低通滤波器,它用来衰减因离 散化产生的高频部分,尽可能恢复原来的信号 散化产生的高频部分 尽可能恢复原来的信号 时域观点看 保持器是一种 保持器是一种信号外推器 其任务是解决各 从时域观点看,保持器是一种信号外推器 ,其任务是解决各 采样点之间的插值问题 采样点之间的插值问题
10
5、三角形保持器
y[( k + 1)T ] y ( kT ) y ( t ) = y ( kT ) + ( t kT ) kT ≤ t ≤ ( k + 1)T T 特点: 特点:要求在计算 kT ~ ( k + 1)T 的 h ( t )时,需已知 y[( k + 1)T ]
三角保持器的 单位脉冲过渡 函数
x( 0 ) + ∫
A ( kT τ )
( k +1) T
0
e A ( kT+ T τ ) Bu( τ)dτ
( k +1) T kT
x( 0 ) + ∫ e
0
AkT
Bu( τ)dτ + ∫
kT 0
e A ( kT τ ) Bu( τ)dτ]
t = kT时, x (kT ) = e
x (0) + ∫ e A( kT τ ) Bu(τ )dτ
& x = Ax + Bu
拉氏变换, 拉氏变换, sX ( s ) X ( 0 ) = AX ( s ) + BU ( s )
X ( s ) = ( sI A ) 1 X ( 0 ) + ( sI A ) 1 BU ( s )
令: φ (t ) = L1[( sI A) 1 ] = e At
其特性与低通滤波器相似; ① 其特性与低通滤波器相似; 有相位滞后问题(滞后T/2); ② 有相位滞后问题(滞后 ); 若输入信号为阶跃信号,那么零阶保持器能无失真地恢复信号。 ③ 若输入信号为阶跃信号,那么零阶保持器能无失真地恢复信号。
4、一阶保持器
一阶保持器使用当前时 进行外推 , y ( kT + t ) = y ( kT ) + y ( kT ) y[( k 1)T ] t T
零阶近似(即采用零阶保持器) 1. 零阶近似(即采用零阶保持器)
11
为使采样后的信号能无失真的再现, 为使采样后的信号能无失真的再现,首先要求采样频率比 信号的最高频率高2倍以上,否则会出现信号重叠。 信号的最高频率高 倍以上,否则会出现信号重叠。同时需要 倍以上 加保持器,另外,由于保持器与理想保持器的特性有差别, 加保持器,另外,由于保持器与理想保持器的特性有差别, 要想无失真地再现原信号是不可能的。如果采样频率足够了, 要想无失真地再现原信号是不可能的。如果采样频率足够了, 合适的保持器, 又选择了 合适的保持器,则由保持器引起的失真可以减到足 够小。 够小。 从离散相似算法离散连续系统的角度讲,选择采样周期就 从离散相似算法离散连续系统的角度讲,选择采样周期就 采样周期 相当于数值积分算法的步距,选择保持器的形式相当于选择 相当于数值积分算法的步距,选择保持器的形式相当于选择 数值积分算法的步距 保持器的形式 算法。保持器的阶次越高,相当于数值积分算法的阶次就越 算法。保持器的阶次越高,相当于数值积分算法的阶次就越 越高 数值积分算法的阶次 高。 对于高阶保持器,其重构信号好,对于高频干扰灵敏, 对于高阶保持器,其重构信号好,对于高频干扰灵敏,但 相位滞后严重,硬件实施困难,因此工程上很少使用。 相位滞后严重,硬件实施困难,因此工程上很少使用。
h g ( t )是一个幅值为 1,持续时间为 T 的矩形脉冲,并可表示 的矩形脉冲, 阶跃函数之和: 阶跃函数之和: h g ( t ) = 1t 1t T 取拉氏变换得到零阶保 持器的传递函数: 持器的传递函数: 1 e Ts 1 e Ts H 0 ( s) = = s s s
相关文档
最新文档