求解离散化状态方程

IBB-EUT

程序如下:

A=[0 0 1 0;
" A( M+ j3 i, H& Q- [: |* O   0 0 0 1;
' ?: B) a- Y2 Z7 o$ `0 b   1 1 1 2;
* D/ f2 @. F- a8 ?   5 -6 3 -4];
B=[0 0;
! n( I  a( I0 D/ O- |( v   0 0;
   -1 0;
/ g2 g+ c, Z8 \' a. f7 c   0 7];
6 v" t  o( A+ Q, rC=[1 1 1 2;
   1 -1 0 0;
& G' M9 H# g! v  i+ w. [2 t   0 1 0 0];%%  A第三行与C第一行一样；B第三行与D第一行一样
- w  C' [9 N* F0 E  kD=[-1 0;
    0 0;
   -1 0];
) ^! Y9 j% z- q. S" R1 n1 UN=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;
%%u2=cosx;
[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化
for i=1
$ D. D3 W$ k4 B  D$ c$ Y    U(1,i)=sin(h*(i-1));
    U(2,i)=cos(h*(i-1));
end      %%%  输入
. I/ y, Q. o) |0 F5 kfor j=1
! a1 e9 {& n  V" k    Y(:,j)=Cd*X(:,j)+Dd*U(:,j);
    X(:,j+1)=Ad*X(:,j)+Bd*U(:,j);
end        %%%  求离散化后的状态方程
plot(X(3,);
) t: l  Q  v; }* ~5 Y; N4 m4 ^: Mhold on;
1 J% B8 g- \1 z& E$ E% f. Vplot(Y(1,)
. Z0 f7 h3 D1 N, a& p# `( Y! `
; Y0 ?! w% n5 a
按道理两个曲线应该重合，结果（附件）却不重合，这是什么原因？

/ k3 j+ C) p9 q7 C  i, {/ A

Uifhjvv

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

小小鲁班

楼主你的代码我跑出来报这个错误
/ U: L6 M4 U0 U: c: V

shuddkk

来学习一下