找回密码
 注册
关于网站域名变更的通知
查看: 458|回复: 3
打印 上一主题 下一主题

求解离散化状态方程

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-12-25 10:40 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
程序如下:" ]0 Q! K8 r1 |! z
2 e  m# }' ?$ w  K
A=[0 0 1 0;
2 Q% P" B0 P' p( Q% r   0 0 0 1;! t2 a3 @( z# _: O3 s  r5 U
   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;( ?; C9 E" T/ ~4 A
   -1 0;
% l; _/ _9 ?  v0 R/ k. ^' |   0 7];
1 [/ Q9 f) v1 ^5 d: ^5 LC=[1 1 1 2;2 c# m4 a4 f( E: z! G$ o+ X$ U
   1 -1 0 0;% X0 `) a1 M+ E6 {  t  s/ ^$ ?
   0 1 0 0];%%  A第三行与C第一行一样;B第三行与D第一行一样6 ~7 p. T% {( R9 Y/ O+ m
D=[-1 0;
/ @+ C8 p. l. o2 @0 X! o& Q7 O    0 0;5 |! y2 F/ U) A5 \* r
   -1 0];  E6 z5 x5 M9 U+ {, k( U, F( F
N=500;T=5;a=0;h=(T-a)/N;L=length(a:h:T);    %%步长' q  Y* e" C. d# f6 ~2 |
X=zeros(4,L);Y=zeros(3,L);U=zeros(2,L);       %%初始化$ d# @( \" J6 {9 O, s/ J# ?, @, h
%%u1=sinx;
$ F2 ~$ n/ P1 G6 O# n" {) i5 q%%u2=cosx;
1 b9 l5 f* ?$ [% l- o6 }" h0 h, A
[Ad,Bd,Cd,Dd]=c2dm(A,B,C,D,h);       %%%状态方程离散化: `+ |3 I* |5 L1 _3 D; j
for i=1( {7 O. J+ V! \" K) r1 V
    U(1,i)=sin(h*(i-1));1 I% F; k# ^; O
    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=15 l2 v; V* l/ e8 y2 Z1 ?- L- {
    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,);' E/ _, q0 K  O
hold on;
- f6 N3 t. u- {* k' l3 eplot(Y(1,)6 @' |+ U+ D) U! O
. Z. D& J# k2 g* P; [) {

5 c# o6 z( o- D/ z) s* [按道理两个曲线应该重合,结果(附件)却不重合,这是什么原因?3 W; G- C8 g9 p& o1 Z5 R: I5 M
/ L/ d; ~& c* H

该用户从未签到

2#
发表于 2020-12-25 11:13 | 只看该作者
从你的C矩阵看,你的Y(1)和X(3)可不相等哦,这两个曲线不重合很正常。

该用户从未签到

3#
发表于 2020-12-25 13:14 | 只看该作者
楼主你的代码我跑出来报这个错误
4 N8 m3 H, t! q

201605k82858yrmtdd9z20.png (27.2 KB, 下载次数: 4)

201605k82858yrmtdd9z20.png

该用户从未签到

4#
发表于 2020-12-28 11:21 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-31 19:41 , Processed in 0.171875 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表