求解离散化状态方程

IBB-EUT

程序如下:

A=[0 0 1 0;
5 n! V: n: o! f; }5 q   0 0 0 1;
, X: O: z. D  L7 x. G   1 1 1 2;
9 S- u8 V' I; D9 F( D   5 -6 3 -4];
3 g* d' k. \9 I; C5 [$ L2 S# L/ rB=[0 0;
   0 0;
   -1 0;
, P5 c+ M3 h* J8 N9 @- `   0 7];
" ]# Z- \. v" N8 G( m! o) }) y- N; {1 ~7 rC=[1 1 1 2;
   1 -1 0 0;
, C, x& @( H3 {  P5 g+ Q! d   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
. O$ U( k  k" C9 f1 M& \6 v: m% \D=[-1 0;
3 D3 A2 P5 i  G    0 0;
0 ?( \+ ~; x: h5 V. s# x   -1 0];
N=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长
: T$ t% ]' O  d& h: {X=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化
6 [% C9 Z* Q7 X5 I( D% A% I7 p%%u1=sinx;
1 Z  k* g' _  N* z# ]! p%%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));
end      %%%  输入
! Z$ j& `2 Z6 O, l) jfor j=1
' F; H8 R; I* R7 t  }) Q    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
plot(X(3,);
hold on;
7 @4 k3 Z5 H# ?( v4 \plot(Y(1,)


按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？

3 {" t* P0 d# Y& ^4 E( d

Uifhjvv

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

小小鲁班

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

shuddkk

来学习一下