|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
ys=tf(50,[0.125,7, 0]); %tf是传递函数 即被控对象函数G(); 5 p6 t* i* S' g; P3 @$ V5 m4 q0 K
dsys=c2d(sys,ts,'z'); %把控制函数离散化取Z变换n阶定常离散系统差分方程
- Y F. t# S3 K: X9 w %在零初始条件下取Z变换:
7 q7 e) m, R- x4 {9 |% V J* ~3 W* Y0 V %dsys即Y(z)/U(z)5 A3 b& t5 W- A& c
[num,den]=tfdata(dsys,'v');% 离散化后提取分子、分母
3 W8 y) L6 l+ \
0 O. n [/ G+ m, k1 ~) \; s5 w6 H) e) _# ~
yout(k)=-den(2)*y_1-den(3)*y_2+num(2)*u_1+num(3)*u_2;
, t& w& c2 C) ^7 a. ~4 w" {; }$ I1 t( h+ p5 z3 S
传递函数如下:
}/ _0 b8 X, _% E8 Q( Asys =) C. `. }" `7 E4 y3 Q, x
* X5 ~: c1 D4 ^/ }. ] 500 N. R: Z3 X& w0 N ^
---------------" h+ e5 C4 ?. o2 m4 u; V
0.125 s^2 + 7 s
7 B- x6 D* H- P/ _! }+ C/ J4 P5 Q9 p
$ f) p/ H: j# \% G) R- Z
z变换以后:5 \$ D3 t) \6 m% [0 {
dsys =
* ?: }) P, e: \0 X
! ^. ]* @2 h2 Q 0.0001963 z + 0.0001927
, {% [ c2 A' e, L c7 A -----------------------9 ~/ ^5 V- u1 r9 d. ?7 }
z^2 - 1.946 z + 0.9455- Z' s9 r/ I& G7 b, D
/ w7 V- W, x6 M% T8 V
0 Y2 O! }+ N7 c) c' Y请问这个yout(k)是怎么来的???求助! ^: ^6 G! J" d$ K3 U
|
|