求解离散化状态方程

IBB-EUT

程序如下:
) P" S( V' ?" V( Y. E: k; w# i( i
A=[0 0 1 0;
   0 0 0 1;
   1 1 1 2;
, b0 }% ?1 X& G% O, I: K1 U   5 -6 3 -4];
B=[0 0;
   0 0;
- W: s$ S3 z; y0 ^/ s8 k) S   -1 0;
   0 7];
C=[1 1 1 2;
4 k3 \* L& o: \/ g7 t) ^5 A- W6 W   1 -1 0 0;
   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
D=[-1 0;
    0 0;
) X; v8 i" Z; V: u. m" F  l   -1 0];
) y  X. I; V' h2 Z. RN=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长
X=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化
%%u1=sinx;
7 S+ c" N( ?! B. s%%u2=cosx;
& o3 b: X# R. e- V9 a; f% x[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
# N! |6 [) C' r. {, tfor i=1
$ Z: ^$ V$ ?6 \0 ^; T+ M    U(1,i)=sin(h*(i-1));
    U(2,i)=cos(h*(i-1));
end      %%%  输入
for j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
3 J! G, a0 k, F4 p& h6 F3 D" h    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
plot(X(3,);
hold on;
plot(Y(1,)
; P2 S/ S1 y: _# e, K1 T' f1 Y

9 G$ k1 p, s1 T$ W按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？
' h* {+ |+ p  j2 ]

Uifhjvv

从你的C矩阵看，你的Y(1)和X(3)可不相等哦，这两个曲线不重合很正常。

小小鲁班

楼主你的代码我跑出来报这个错误
& g3 C6 K, A3 Q

shuddkk

来学习一下