求解离散化状态方程

IBB-EUT

程序如下:

. c7 I9 i# S+ X4 j4 ]3 a; xA=[0 0 1 0;
- o' l9 {* C: L1 p$ Y   0 0 0 1;
6 |8 A0 v, w7 d5 t- `  _3 m- J8 D- L  ]   1 1 1 2;
6 \) d6 F6 v1 L  c7 K   5 -6 3 -4];
B=[0 0;
+ |# m" {5 M% I0 B; p; s9 o   0 0;
9 m3 m0 K" }- c" Q) i# a   -1 0;
   0 7];
C=[1 1 1 2;
   1 -1 0 0;
   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
" w, D0 ]8 }& X! y* Q8 mD=[-1 0;
$ \1 B2 x$ c& I& F( ?, P    0 0;
& H% S/ V3 \3 {5 L% ^   -1 0];
N=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长
$ N% w) k! Z9 H" P' ?# i, hX=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化
4 ]8 R" Q8 N* N%%u1=sinx;
6 b. A' @( I; d0 X. ?3 F%%u2=cosx;
[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
! F; Z. `: v; E+ rfor i=1
% X" O  o5 u7 f9 Y, v% D3 [5 Y    U(1,i)=sin(h*(i-1));
* s) Y7 b- J. `    U(2,i)=cos(h*(i-1));
end      %%%  输入
for j=1
    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
0 i; S7 ^6 ~# D( x: z9 t" W3 G9 u    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
. ^( y& z3 ]/ x8 c6 Cend        %%%  求离散化后的状态方程
1 J" W9 Y& `: p* ^/ qplot(X(3,);
hold on;
" z7 a" ]; R( I  k) {) S4 ]* \plot(Y(1,)
3 j& r) b- z* v) R1 k- @3 C
1 X' C4 t4 V! Z- P
+ K3 b8 p1 u& e- i" h按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？
) M8 G! h0 s" T

Uifhjvv

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

小小鲁班

楼主你的代码我跑出来报这个错误
! s% V( j3 |- K/ Z+ O. T; f

shuddkk

来学习一下