求解离散化状态方程

IBB-EUT

程序如下:
) @3 A0 _! h0 L9 `
3 ?0 h* J. J1 b' I  WA=[0 0 1 0;
4 c: J2 Z- o, M9 ~7 q2 F   0 0 0 1;
4 U+ Z" @) {. K! ]- y5 Y1 q# h. m   1 1 1 2;
- D7 H; E$ v* B+ I4 s4 q. P, b   5 -6 3 -4];
& q# Q$ X; D4 @% I4 h" z9 n% [B=[0 0;
8 v5 N! W7 |: N* A3 q8 |. A   0 0;
   -1 0;
, k/ V" N+ q) r& `: p' l! j: F   0 7];
( f! P( V1 U. V/ |: o- p# |0 VC=[1 1 1 2;
0 ~" m, A! P& E0 h   1 -1 0 0;
" n2 N' @/ J/ \+ H; E0 E+ H1 U   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
+ F  E7 k/ D# {: {- H2 U$ r# f" lD=[-1 0;
    0 0;
   -1 0];
1 _5 h; k4 r/ q1 JN=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);       %%初始化
  @2 Y5 l+ f$ Z8 K# w( m' f%%u1=sinx;
%%u2=cosx;
5 ?- W+ f- |- k( G% l[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
for i=1
    U(1,i)=sin(h*(i-1));
    U(2,i)=cos(h*(i-1));
1 J6 R% B( T- tend      %%%  输入
for j=1
' w: ^$ a7 Q5 s7 F0 B6 J/ E2 s    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
plot(X(3,);
2 H5 D. ~5 f# z  xhold on;
) E% L: ]% i: j' ?; dplot(Y(1,)


按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？
+ K0 ^7 R( k+ ?3 }5 m, C

Uifhjvv

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

小小鲁班

楼主你的代码我跑出来报这个错误

shuddkk

来学习一下