|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为 }* m& m/ x0 V1 ^+ _* N6 W+ a
m4 c% t* N( P. c
y=filter(b,a,f)% g; J0 V# l4 {' |
7 U: ?$ x0 O: N. f& e7 m. `6 m
4 c; }) Q- `. _5 M: U其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
, |8 ~+ ~0 V! a* m* F3 s! [2 N/ o" |& Z( R/ q- N
E; W% Q/ \5 e) m+ P2 ~8 E; r9 @3 B. d/ z6 X3 x
R=51;
x- w; b0 B6 Z2 o$ {8 ^- Qd=rand(1,R)-0.5;" c# J6 f# j0 o/ |. h9 m7 l
k=0:R-1;
N$ b' U& G) x. q8 j9 is=2*k.*(0.9.^k);
! V8 ]2 l* i8 A$ f! cf=s+d;
- N0 y% R- t5 O/ j4 Q0 Kfigure(1);
+ N9 t) L/ i# r; q# B7 ], fstem(k,f,'.');
; R! g8 w- ~" ~* w4 z' u. AM=5;4 X. b, a! N# l
b=ones(M,1)/M;
% V% l7 F. |: D6 t2 E) Ya=1;
* I- v/ b$ j! s1 h/ xy=filter(b,a,f);8 g' ~1 k5 R1 V* O- F# H) j
figure(2)
. D7 Z# ?1 a8 a4 estem(k,y,'.'): M- R- ^3 ^+ f
5 G1 I. s6 y+ I8 w5 p
: U" z+ D3 F. [6 f8 T5 k& q2 D0 L2 a0 H4 r3 p, F- n" t% ^( ]
$ V$ E' b$ ]) V |
|