|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
' X6 J# c* O. j! `
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为
- w5 x# r/ d- `- c) r$ n* j0 B8 m E* {; r; u, u
- y=filter(b,a,f)* Z( D7 A* x' z8 `
0 Y7 J9 Y+ U8 S) H. b
3 ~$ ?) F) o* ~3 w+ d1 B7 t8 T
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
7 m2 @+ x% |: L% @9 o: [1 c" N B2 X2 J8 R9 M4 C0 V8 \
# y* M# B) c, H6 f, P' k0 |& d! r+ U! p: L# i5 }
- 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,'.')
: `1 z. ~" @1 h% Z' h, c2 z4 W. Q
6 U6 E" g7 e0 ^3 b* m, p
) E! l8 y: {/ a9 q- c2 F, t |
|