当前位置：文档之家 › 状态方程求传递函数

状态方程求传递函数

  • 格式：docx
  • 大小：36.99 KB
  • 文档页数：3

状态方程求传递函数

一、前言

在控制系统中，传递函数是一个非常重要的概念。它描述了输入信号和输出信号之间的关系，可以用来分析和设计控制系统。状态方程是另一种描述控制系统的方法，它将系统的状态表示为一组变量，并用微分方程描述它们之间的关系。在本文中，我们将介绍如何使用状态方程求传递函数。

二、什么是传递函数

传递函数是一个数学表达式，用于描述输入信号和输出信号之间的关系。它通常表示为H(s)，其中s是复变量（Laplace变换中使用的变量）。传递函数可以从系统的输入-输出特性中导出，因此它对于分析和设计控制系统非常有用。

三、什么是状态方程

状态方程是一组微分方程，用于描述控制系统中各个部分之间的关系。它将系统的状态表示为一组变量，并根据这些变量之间的关系来定义微分方程。如果我们知道了初始条件和输入信号，就可以使用状态方程来预测系统在任意时间点上的行为。

四、如何从状态方程求解传递函数

要从状态方程求解传递函数，需要进行以下步骤：

1. 将状态方程转换为拉普拉斯域方程

首先，我们需要将状态方程转换为拉普拉斯域方程。这可以通过将每个状态变量的导数替换为s乘以该变量的拉普拉斯变换来实现。例如，如果我们有一个状态方程x' = Ax + Bu，其中x是状态向量，u是输入向量，A和B是常数矩阵，则可以得到以下拉普拉斯域方程：

sX(s) - x(0) = AX(s) + BU(s)

其中x(0)是初始条件。

2. 解出传递函数H(s)

接下来，我们需要解出传递函数H(s)。这可以通过将输出变量表示为输入变量和状态变量的线性组合来实现。例如，如果我们有一个输出变量y = Cx + Du，则可以得到以下传递函数：

H(s) = Y(s)/U(s) = C(sI - A)^(-1)B + D

其中I是单位矩阵。

3. 化简传递函数

最后，我们需要化简传递函数。这可以通过因式分解和合并项来实现。例如，如果我们有一个传递函数H(s) = (s+1)/(s^2+2s+1)，则可以将其化简为H(s) = 1/(s+1)，因为分母可以因式分解为(s+1)^2。

五、总结

在本文中，我们介绍了如何使用状态方程求解传递函数。首先，我们将状态方程转换为拉普拉斯域方程，然后解出传递函数，并最终化简传递函数。这个过程对于分析和设计控制系统非常有用，因为它可以帮助我们理解输入信号和输出信号之间的关系。

合集下载

传递函数的求取 (2)

传递函数的求取 (2)

传递函数的求取

一、 实验内容及目的

本次实验要求如下：

○1用足够多的方法求得以下电路系统的传递函数。

○2当在Ui上加入一个1V的输入电压时仿真出系统的输出曲线

其中Ui是输入，Uo是输出。

本次实验共用了4种方法求得传递函数，分别是利用微分方程求解、利用阻抗法求解、利用方框图化简求解、利用流图与梅森公式求解。之后用了两种方法求得输出曲线，分别是matlab程序仿真和simulink图形仿真。

实验目的是通过实践分析不同求传递函数方法的需求条件，加深对各种工具的熟练程度。

一、 实验方案及内容

1、利用微分方程直接求传递函数

根据电路理论可列得下列等式：

-----------------------------------------○1

-----------------------------------------○2

-----------------------------------------○3

-----------------------------------------○4

-----------------------------------------○5

利用拉布拉斯变换将其转化为频域下的方程：

------------------------------------------○6

------------------------------------------○7

------------------------------------------○8

------------------------------------------○9

传递函数到状态空间的实现

传递函数到状态空间的实现

学生姓名： 刘吕 学号： 20121562

实验题目： 传递函数到状态空间的实现

课程名称： 计算机仿真

一、实验目的：

 理解并掌握传递函数转换为状态空间方程的方法

 理解状态初值的计算方法

二、实验内容：

 应用MATLAB编写一个可以实现传递函数到状态空间方程的可控可观规范型的m文件。并用相应例题验证程序的正确性。

 完善该程序使其可以用来计算状态初值。并用相应的例题验证程序的正确性。

 程序中需要考虑分子分母同阶以及分母首系数不为1的两种情况。

三、报告内容：

（1） 给出m文件的程序框图，及验证结果，并记录出现的错误，并给出解决的方案。若没有得到解决，请说清楚你的问题

（2） 状态初值的求解，请给出相应的验证结果，并计算与精确解之间的误差。

四、实验原理：

1、传递函数为10111011()nnnnnnnnbsbsbsbGsasasasa……

其状态空间模型能控标准型为：

12101000010nnnAaaaa……………………

001B…，0101101nnnnCbbabbabba…，nDb

能观标准型为：

000,,TTTAABCCB，ToDD

2、计算状态变量初值： 11211221'22321''333(2)11(1)(0)1(0)-(0)10(0)0(0)(0)(0)1(0)(0)100(0)nnnnnnnnnnnnnxaaayccccxaayccxaycxayxy………………………………………………………'''(3)1(2)(0)(0)(0)(0)00(00nnuuucuu……………………………………………五、程序流图

传递函数到状态空间的实现 2

传递函数到状态空间的实现 2

计算机仿真实验报告

1

实验题目： 传递函数到状态空间的实现

课程名称： 计算机仿真

一、实验目的

1、理解并掌握传递函数转换为状态空间方程的方法

2、理解状态初值的计算方法

二、实验内容

1、应用MATLAB编写一个可以实现传递函数到状态空间方程的可控可观规范型的m文件。并用相应例题验证程序的正确性。

2、完善该程序使其可以用来计算状态初值。并用相应的例题验证程序的正确性。

3、程序中需要考虑分子分母同阶以及分母首系数不为1的两种情况。

三、报告内容

1、给出m文件的程序框图，及验证结果，并记录出现的错误，并给出解决的方案。若没有得到解决，请说清楚你的问题

2、如果做了程序的状态初值得求解，请给出相应的验证结果，及程序编写过程中出现的问题，若已经解决，给出具体方法。

四、 实验理论

1、 传递函数为G(s)=𝑏0𝑠𝑛+b1𝑠𝑛−1+⋯+𝑏𝑛−1𝑠+𝑐𝑛𝑎0𝑠𝑛+𝑎1𝑠𝑛−1+⋯+𝑎𝑛−1𝑠+𝑎𝑛

那么其状态空间模型能控标准型为：

A=[010⋯0001⋯0⋮⋮⋮⋱⋮−𝑎𝑛−𝑎𝑛−1−𝑎𝑛−1⋯−𝑎1] B=[00⋮1]

C=[(𝑏𝑛−𝑏0𝑎𝑛) (𝑏𝑛−1−𝑏0𝑎𝑛−1)⋯(𝑏1−𝑏0𝑎1)] D=𝑏𝑛

能观标准型为：

𝐴能观=𝐴能控𝑇 𝐵能观=𝐶能控 𝑇 𝐶能观=𝐵能控𝑇 𝐷能观= 𝐷𝑇能控

2、计算状态变量初值：

（1）不含u的导数项时，则有： 计算机仿真实验报告

2

（2）系统微分方程不仅包含u的输入项，而且包含u的导数项，则：

五、程序检验

（1）输入一个分母首系数为1且分子分母不同阶传递函数：

G=2S3+4S2+3S+5S4+2S3+5S2+4S+2

程序运行结果：

能控标准型：

利用bode图求传递函数例题

利用bode图求传递函数例题

例题：已知最小相位系统开环对数频率特性曲线如图所示。试写出开环传递函数)(sGk 。

解：

1) ω<ω1的低频段斜率为[-20]，故低频段为K/ｓ。

ω增至ω1，斜率由[-20]转为[-40]，增加[-20]，所以ω1应为惯性环节的转折频率，该环节为1111s 。

ω增至ω2，斜率由[–40]转为[–20]，增加[+20]，所以ω2应为一阶微分环节的转折频率，该环节为112s 。

ω增到ω3，斜率由[-20]转为[-40]，该环节为1113s,ω>ω3，斜率保持不变。

故系统开环传递函数应由上述各典型环节串联组成，即

)11)(11()11()(312ssssKsGk

2) 确定开环增益K

当ω=ωc时，A(ωc)=1 。

所以 1111)1()1(1)1()(12232122ccccccccKKA

故

12cK 所以，)11)(11()11()(31212sssssGck

练习：

最小相位系统的对数幅频特性如下图所示，试分别确定各系统的传递函数。

(a) (b)

（c）

a：)1(10)(sssG

b：)1)(110(100)(sssG

c)12.0)(15.0(100)(sssG