状态方程离散化处理

连续状态空间方程:

\dot{x}(t)=Ax(t)+Bu(t)

y(t)=Cx(t)+Du(t)

1欧拉法

1.1 状态方程

\dot{x}(t)=\frac{x(k+1)-x(k)}{T}代入状态方程得

\frac{x(k+1)-x(k)}{T}=Ax(k)+Bu(k)

化简为x(k+1)=(AT+I)x(k)+BTu(k)

1.2 输出方程

输出方程为y(k+1)=Cx(k+1)+Du(k+1)

x(k+1)表达式代入上式得

y(k+1)=C(AT+I))x(k)+CBTu(k)+Du(k+1)

1.3 离散化后的状态空间方程

x(k+1)=(AT+I)x(k)+BTu(k)

y(k+1)=C(AT+I))x(k)+CBTu(k)+Du(k+1)

其实欧拉法属于不精准的离散化。

2 公式法

2.1 状态方程

\dot{x}(t)=Ax(t)+Bu(t)的离散公式为

x(k+1)=A_{d}x(k)+B_{d}u(k),其中A_{d}=e^{AT}B_{d}=\int_{0}^{T}e^{As}Bds,即

x(k+1)=e^{AT}x(k)+\int_{0}^{T}e^{As}Bdsu(k)

2.2 输出方程

y(k)=Cx(k)+Du(k)

  • 3
    点赞
  • 42
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在MATLAB中,偏微分方程的离散化主要是通过使用pdepe求解器的算法来实现的。该算法将原始的椭圆型和抛物线型偏微分方程转化为一组常微分方程。这个转换的过程是基于用户指定的网格点,并使用二阶空间离散化技术进行的。然后,该算法使用ode15s指令来求解。ode15s是一种ode求解器,主要用于处理离散化后的代数方程组和微分方程组。因此,通过使用ode15s求解器,我们可以顺利地求解离散化后的偏微分方程。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *3* [偏微分方程的数值解(二): 一维状态空间的偏微分方程的 MATLAB 解法](https://blog.csdn.net/qq_29831163/article/details/89706692)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *2* [MATLAB求解偏微分方程【PDE和差分法】](https://blog.csdn.net/zhounei/article/details/132310247)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值