|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
ys=tf(50,[0.125,7, 0]); %tf是传递函数 即被控对象函数G(); & g9 A1 Z# V3 F9 B, n( S& V
dsys=c2d(sys,ts,'z'); %把控制函数离散化取Z变换n阶定常离散系统差分方程2 q3 B# c# A* q, X3 ?! ]9 n3 I4 W
%在零初始条件下取Z变换:
7 P m) d; d" u2 q %dsys即Y(z)/U(z)' M% R6 v2 H1 Q0 r9 A+ n* U$ e6 l
[num,den]=tfdata(dsys,'v');% 离散化后提取分子、分母
2 A* [( u0 @2 ^# u( W9 }5 G b( F) @+ ^( q* ~. `$ z. X
& `( G/ z2 k! C6 W/ b7 Q3 D3 ^
yout(k)=-den(2)*y_1-den(3)*y_2+num(2)*u_1+num(3)*u_2;
4 Q6 [# i' g7 Z
, j: A9 M* k$ c h; m9 \: i: D; E传递函数如下:
! [! S7 G, H' K; t8 J1 Q. V! k7 Z& Hsys =
6 J" o* l& B- Q" }
- J; a/ ~1 L: B2 h: q! ]& Q1 I3 W 50, l+ d& @' \) Q1 Z* z
---------------+ v6 A: V% f/ }7 g" m; b
0.125 s^2 + 7 s
8 C& O5 ~2 S7 ^1 e5 x% c& P8 q
- U, Q2 W4 n! o& S% g+ ~7 } N# A1 _9 W1 E! N: A2 h
z变换以后:- w P' t% o, }: _9 w' t
dsys =
* b$ T* _& Y" R& m: w
5 x# M/ u8 N& B) H! m- r 0.0001963 z + 0.00019275 t& N9 m/ c+ J
-----------------------
4 J2 f. P# b) G1 @* } z^2 - 1.946 z + 0.9455" z' l+ U: ~, k! Z; D# {
. y: ~9 e1 m2 w7 i4 \7 ]7 q
: L8 ~' B9 L8 `
请问这个yout(k)是怎么来的???求助!
1 H: ?6 L. `& y. w( M) w |
|