求解离散化状态方程

IBB-EUT

程序如下:
" f2 J; ^, _+ d  S5 [2 S
; E0 H+ f$ Y( Y5 C5 NA=[0 0 1 0;
   0 0 0 1;
   1 1 1 2;
   5 -6 3 -4];
B=[0 0;
4 r+ h+ S7 d# Q( L. `( j   0 0;
   -1 0;
   0 7];
$ o( B+ W! g+ V: CC=[1 1 1 2;
   1 -1 0 0;
8 J5 ]% |4 V5 _' K   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
D=[-1 0;
) d1 q7 P& h5 l; q    0 0;
   -1 0];
  R; i6 M8 K) q5 zN=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);       %%初始化
* H# p5 Y( W1 u+ b%%u1=sinx;
%%u2=cosx;
[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
; n. W0 K- H- n- B) gfor i=1
    U(1,i)=sin(h*(i-1));
; @" e0 V9 K( W    U(2,i)=cos(h*(i-1));
( D- c7 {/ S1 ?5 Nend      %%%  输入
for j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
- e1 |, a! X5 r9 z, R% l2 v: n7 aplot(X(3,);
hold on;
plot(Y(1,)
1 u: Y; C# t- O

按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？
) a0 r! t; X$ M& {: K5 Q: E
# C) I4 Y) I3 }& r- p

shuddkk

来学习一下

小小鲁班

楼主你的代码我跑出来报这个错误
4 T. l+ ^1 w9 k3 u* p# {3 w% P

Uifhjvv

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