求解离散化状态方程

IBB-EUT

程序如下:

A=[0 0 1 0;
   0 0 0 1;
6 V6 Y' D4 }5 h   1 1 1 2;
   5 -6 3 -4];
2 M$ h, \( Z! C  pB=[0 0;
) r# m0 D% ~1 F, {  A   0 0;
( c6 Z) n) q7 t0 c8 R9 l   -1 0;
1 O) h2 T  b4 a1 l5 ^" s  _   0 7];
' @- V7 c( G- YC=[1 1 1 2;
- x; N/ y7 T8 X   1 -1 0 0;
) p8 Y+ x, O1 A, ^4 z1 ]( @7 A7 J   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
! i+ q/ k: N% K/ h+ G; p' f# F7 g7 AD=[-1 0;
    0 0;
   -1 0];
N=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长
& K$ ^( U- B4 Q  l$ w) Z7 |. uX=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化
. a8 Z( j/ F1 d8 B  e/ o0 c$ a- R. `8 s%%u1=sinx;
4 T  {8 g/ I8 g- L' N- R%%u2=cosx;
8 R1 V% G2 r# k4 @0 A[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
for i=1
7 R/ P# V5 P! y0 v* i    U(1,i)=sin(h*(i-1));
7 T) X) H4 k8 j, s3 D7 I  g4 l    U(2,i)=cos(h*(i-1));
* z; X' a7 X# vend      %%%  输入
for j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
; K& I' ^( m) [8 M, ~' N+ D    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
* b* t: G6 p$ c- Lend        %%%  求离散化后的状态方程
. N9 Z1 L# |" [% k/ D2 M: i, wplot(X(3,);
9 t# o/ h% I. Yhold on;
plot(Y(1,)


  L; b' h/ O# q2 |5 \" P+ _' S按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？

# n# N! ~3 P1 [9 D# c+ k2 R

Uifhjvv

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

小小鲁班

楼主你的代码我跑出来报这个错误
) f/ F) d% s* i) r8 @: ]$ @

shuddkk

来学习一下