|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
9 L& u% q+ m% O0 {% d/ o在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为, h8 g' P% q0 J: A3 f( O0 A
x5 L" L7 f3 p2 J5 Q& I" v- y=filter(b,a,f)
! Y- U& S4 G( g( q) E& W
8 r+ y' U+ D d4 O+ @* p; `7 l0 J
5 C' p D% Y$ n& h0 G, i% n其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。4 K6 b% z* j* M& D$ p+ L
5 @/ a0 n- R2 }: F
- m& f5 C' `- m& i4 y2 f5 F! j
$ r* y. z; N- o9 ^0 ^! t4 j' n- R=51;
- d=rand(1,R)-0.5;
- k=0:R-1;
- s=2*k.*(0.9.^k);
- f=s+d;
- figure(1);
- stem(k,f,'.');
- M=5;
- b=ones(M,1)/M;
- a=1;
- y=filter(b,a,f);
- figure(2)
- stem(k,y,'.')
6 W% K6 o, S$ e* o6 W( o! Y 6 R" |% }- p( G$ ?- E
, i( S) ~* y9 n" v9 Q& y/ o |
|