EDA365电子论坛网

标题: 求解离散化状态方程 [打印本页]

---

作者: IBB-EUT    时间: 2020-12-25 10:40
标题: 求解离散化状态方程
程序如下:
9 B" E( y- W0 v4 U
- S7 S) T: [2 ~3 ~! k% Z5 ]  _A=[0 0 1 0;
3 H; }% V( i& Y& U7 k2 w4 T$ z$ X   0 0 0 1;
; E6 f# s0 `% w  d- D. a6 t% h   1 1 1 2;
   5 -6 3 -4];
3 v2 w7 F: m# `B=[0 0;
/ E+ U/ v. r* O+ w( n! O; P, Y) f0 s   0 0;
   -1 0;
   0 7];
& `8 f. \& k5 \. [: hC=[1 1 1 2;
9 a. i9 F1 j- F' Q. \5 [   1 -1 0 0;
   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
D=[-1 0;
    0 0;
2 R1 h  T$ I. M   -1 0];
  M( `( ?+ s$ Q/ z5 vN=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长
6 d; @8 u. u4 _. GX=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化
%%u1=sinx;
' y) ^% i  M( h* o. W%%u2=cosx;
[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
* Y3 E1 A+ ~& A0 Lfor i=1
    U(1,i)=sin(h*(i-1));
    U(2,i)=cos(h*(i-1));
& ]: o/ _+ V5 Z. O3 f, Hend      %%%  输入
for j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
  A0 g: }/ ^: |) ?' a8 r3 G" Kplot(X(3,);
hold on;
( ~( Z* J  c% ^8 Zplot(Y(1,)


# W. @0 j/ v6 R按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？
3 q- B8 r1 m: A" I8 p

---

作者: Uifhjvv    时间: 2020-12-25 11:13
从你的C矩阵看，你的Y(1)和X(3)可不相等哦，这两个曲线不重合很正常。

---

作者: 小小鲁班    时间: 2020-12-25 13:14
楼主你的代码我跑出来报这个错误

201605k82858yrmtdd9z20.png (27.2 KB, 下载次数: 6)

下载附件  保存到相册

2020-12-25 13:14 上传

---

作者: shuddkk    时间: 2020-12-28 11:21
来学习一下

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2