当前位置: 首页 > article >正文

MATLAB 控制系统设计与仿真 - 27

状态空间的标准型

传递函数和状态空间可以相互转换,接下来会举例如何有传递函数转成状态空间标准型。

对角标准型

当 G(s)=\frac{B(s)}{A(s)}=\frac{b_2s^2+b_1s+b_0}{s^3+a_2s^2+a_1s+a_0} 

G(s)可以写成:

G(S)=\frac{k_1}{s-p_1}+\frac{k_1}{s-p_2}+\frac{k_3}{s-p_3}, p_1\neq p_2\neq p_3

即:

根据上图可知:

A=\begin{bmatrix} p_1 & 0 & 0\\ 0& p_2 & 0\\ 0 & 0 & p_3 \end{bmatrix} , B=\begin{bmatrix} 1\\ 1\\ 1 \end{bmatrix} , C=\begin{bmatrix} k_1 &k_2 &k_3 \end{bmatrix} , D=0

约当标准型

当 G(s)=\frac{B(s)}{A(s)}=\frac{b_2s^2+b_1s+b_0}{s^3+a_2s^2+a_1s+a_0}

G(s)可以写成:

G(s)=\frac{k_1}{s-p_1}+\frac{k_2}{(s-p_m)^2}+\frac{k_3}{s-p_m}

即:

根据上图可知:

A=\begin{bmatrix} p_1 & 0 & 0\\ 0& p_m & 1\\ 0 & 0 & p_m \end{bmatrix} , B=\begin{bmatrix} 1\\ 0\\ 1 \end{bmatrix} , C=\begin{bmatrix} k_1 &k_2 &k_3 \end{bmatrix} , D=0

能控标准型

当传递函数:Y(s)=\frac{1}{s^3+a_2s^2+a_1s+a_0}U(s)时,

Y(s)(s^3+a_2s^2+a_1s+a_0)=U(s)

\dddot{y}+a_2\ddot{y}+a_1\dot{y}+y=u\\ y=-\dddot{y}-a_2\ddot{y}-a_1\dot{y}+u

如果:x_1=y,x_2=\dot{x_1}=\dot{y},x_3=\dot{x_2}=\ddot{y}

\begin{bmatrix} \dot{x_1}\\ \dot{x_2}\\ \dot{x_3}\end{bmatrix}=\begin{bmatrix} 0 & 1 & 0\\ 0& 0 & 1\\ -a_0& -a_1 & -a_2 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}+\begin{bmatrix} 0\\ 0\\ 1 \end{bmatrix}u

y=\begin{bmatrix} 1 &0 & 0 \end{bmatrix}\begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}

如果传递函数是一般形式如下:

Y(s)=\frac{b_2s^2+b_1s+b_0}{s^3+a_2s^2+a_1s+a_0}U(s),我们可以分解如下:

第一部分我们可以得到

\widetilde{y}=x_1

所以第二部分为:

y=b_2 \ddot{x_1} +b_1 \dot{x_1}+b_0=\begin{bmatrix} b_0 &b_1 &b_2 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}

所以

\begin{bmatrix} \dot{x_1}\\ \dot{x_2}\\ \dot{x_3}\end{bmatrix}=\begin{bmatrix} 0 & 1 & 0\\ 0& 0 & 1\\ -a_0& -a_1 & -a_2 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}+\begin{bmatrix} 0\\ 0\\ 1 \end{bmatrix}u

y=\begin{bmatrix} b_0 &b_1 &b_2 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}

其方框图可表示为:

能观标准型

如果Y(s)=\frac{b_2s^2+b_1s+b_0}{s^3+a_2s^2+a_1s+a_0}U(s),分子分母同时除以s^3可以得到

(1+a_2\frac{1}{s}+a_1\frac{1}{s^2}+a_0\frac{1}{s^3})Y(s)=(b_2\frac{1}{s}+b_1\frac{1}{s^2}+b_0\frac{1}{s^3})U(s)

因此:

Y(s)=-a_2\frac{1}{s}Y(s)-a_1\frac{1}{s^2}Y(s)-a_0\frac{1}{s^3}Y(s)+b_2\frac{1}{s}U(s)+b_1\frac{1}{s^2}U(s)+b_0\frac{1}{s^3}U(s)

其方框图可以表示如下:

其中:

Y(s)=X_1(s);\\ sX_1(s)=-a_2X_1(s)+X_2(s)+b_2U(s) \\ sX_2(s)=-a_1X_1(s)+X_3(s)+b_1U(s) \\ sX_3(s)=-a_0X_1(s)+b_0U(s)

对上式进行Laplace逆变换可得:

y(t)=x_1(t) \\ \dot{x_1}(t)=-a_2x_1(t)+x_2(t)+b_2u(t) \\ \dot{x_2}(t)=-a_1x_1(t)+x_3(t)+b_1u(t) \\ \dot{x_3}(t)=-a_0x_1(t)+b_0u(t) \\

因此:

\begin{bmatrix} \dot{x_1}\\ \dot{x_2}\\ \dot{x_3}\end{bmatrix}= \begin{bmatrix} -a_2 & 1 & 0\\ -a_1& 0 & 1\\ -a_0& 0 & 0 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}+\begin{bmatrix} b_2\\ b_1\\ b_0 \end{bmatrix}u

y=\begin{bmatrix} 1 &0 &0 \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3 \end{bmatrix}

********************************************************************************************************

最后,欢迎大家有问题给我留言。

另外欢迎小伙伴们点赞转发加关注

谢大家!!


http://www.kler.cn/a/590277.html

相关文章:

  • 【机器人-基础知识】标定 - 相机标定全解
  • 【FAQ】HarmonyOS SDK 闭源开放能力 —Map Kit(6)
  • json字符串转对象,对象转JSON
  • c++基础知识-图论进阶
  • ArcGIS10.X影像智能下载!迁移ArcGIS Pro批量智能高清影像下载工具至ArcGIS!
  • Matlab 汽车悬架系统动力学建模与仿真
  • Docker封装镜像、分发、部署实践:nginx
  • 轨道交通CPU+FPGA控制器,支持codesys/vxWorks/翼辉等实时系统
  • 基于eNSP的IPV4和IPV6企业网络规划
  • 商用电脑VS家用电脑,为何企业更需要华硕NUC 14 Pro?
  • 【GIS】重要技术3DGS
  • 计算机网络-网络接入技术
  • 【每日学点HarmonyOS Next知识】状态栏字体、生命周期、自定义对话框屏幕中间、透明度、tab居中
  • 从零开始 | C语言基础刷题DAY3
  • 数据库管理员助理(DP-300)适合什么群体考?
  • DeepSeek:从入门到精通
  • 开发一个go模块并在其他项目中引入
  • 耗时sql分析
  • Webpack 模块联邦 vs npm 私服
  • QT编译器mingw与msvc区别及环境配置