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

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}

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

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

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

谢大家!!

相关文章:

  • 光谱相机识别瓶子材质的技术原理和应用案例
  • JS超过Number的最大值
  • 如何启用 HTTPS 并配置免费的 SSL 证书
  • hubilder打包ios app, 并上传TestFlight
  • c语言笔记 字符串函数---strstr strlen strtok以及sizeof
  • 1.5.3 掌握Scala内建控制结构 - for循环
  • 告别数据库束缚!用esProc在 csv 文件上执行 SQL
  • 搭建主从服务器
  • K8S之QoS详解
  • oracle 基础知识之 多表查询
  • soul能用ip精准定位吗?ip属地准确吗
  • hot100算法刷题:二叉树的层序遍历
  • Spring Boot 核心知识点深度详解:自动化配置 (Auto-configuration) - 解锁 Spring Boot 的 “魔法”
  • 远程访问家里电脑上部署的Stable diffusion - 免费篇
  • error: Microsoft Visual C++ 14.0 or greater is required解决不需要安装到c盘等
  • Mybatis使用Druid连接池
  • 0x04.若依框架微服务开发(含AI模块运行)
  • Kotlin 中 let 方法的作用和使用场景
  • 人工智能助力家庭机器人:从清洁到陪伴的智能转型
  • ngx_http_conf_ctx_t
  • 刘洪洁已任六安市委副书记、市政府党组书记
  • 国家卫健委对近日肖某引发舆情问题开展调查
  • 停电催生商机,中国品牌 “照亮” 西班牙
  • 经营业绩持续稳中向好,国铁集团2024年度和2025年一季度财务决算公布
  • 庄语乐︱宋代历史是被“塑造”出来的吗?
  • 2025年度中国青年五四奖章暨新时代青年先锋奖评选揭晓