求解离散化状态方程

IBB-EUT

程序如下:

A=[0 0 1 0;
2 Q% P" B0 P' p( Q% r   0 0 0 1;
   1 1 1 2;
( c4 h$ p8 _9 ^; c: `5 Y   5 -6 3 -4];
$ C+ P3 f+ g9 c8 u6 n( S4 f8 NB=[0 0;
3 k; r! Y3 s0 [( n   0 0;
   -1 0;
% l; _/ _9 ?  v0 R/ k. ^' |   0 7];
1 [/ Q9 f) v1 ^5 d: ^5 LC=[1 1 1 2;
   1 -1 0 0;
   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
D=[-1 0;
/ @+ C8 p. l. o2 @0 X! o& Q7 O    0 0;
   -1 0];
N=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;
$ F2 ~$ n/ P1 G6 O# n" {) i5 q%%u2=cosx;
[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));
+ w: t$ ?' L7 N9 {$ Y* x/ Y) c* kend      %%%  输入
- r3 d( j9 h! J0 J, q6 Gfor j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
, T8 }9 D% J& B, r+ k' K    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
. m$ S& [* Q$ Mend        %%%  求离散化后的状态方程
) O7 t* `1 L8 k# b; U5 bplot(X(3,);
hold on;
- f6 N3 t. u- {* k' l3 eplot(Y(1,)


5 c# o6 z( o- D/ z) s* [按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？

Uifhjvv

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

小小鲁班

楼主你的代码我跑出来报这个错误
4 N8 m3 H, t! q

shuddkk

来学习一下