|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
7 \6 q+ G! Q; { b6 w# m" O在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为6 C( O1 f- E- x) A
/ Z) z7 p( {: X4 O7 z- k# A. _- y=filter(b,a,f)
2 p& S6 U) Y& y6 z% c$ w/ R4 ~ ; R# p9 w. l' `! e6 j& ]
" T3 t( J4 U, `# m4 W+ N
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
4 u$ r( v$ h# {$ L( K' i7 m$ ~; y, Z7 x3 E/ e1 ^5 N
) b% @: ?7 a- k. Q/ X
/ F$ n2 V7 d M0 ~& u4 ?
- 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,'.')
, b: C6 j" A' ^9 v) w+ ^+ Q4 ^, Q' E 4 F& B" R8 G: z. u) x8 @% H
3 P r# V% A( {% j# I
|
|