|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为# k! s4 C! n1 S) f5 `3 e
/ L6 J0 D7 } u$ l% Q+ T9 V
y=filter(b,a,f)! ]6 {' p2 T9 M
$ w6 v1 Z2 {! ]" D1 N
1 u% o0 v5 {, ~4 |6 Y- g$ O其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。2 B9 H* R5 V5 W, p+ {
! B2 b7 y. G% ^
; q# q1 Z/ L5 N: r" M9 E, J- l- j9 v4 \6 m: P: ]
R=51;
/ P/ b6 w3 {" @d=rand(1,R)-0.5;2 l7 P* ^& c2 C
k=0:R-1;; Q7 \& U% d$ }7 {, q
s=2*k.*(0.9.^k);- q6 q m1 h! C4 s# p o) C4 C" s5 b
f=s+d;
# {/ j, \- ]/ \5 gfigure(1);
) Z; o3 u# q; o' o* H( x/ cstem(k,f,'.');
- ?+ W, N, K7 H8 R) v8 _1 s7 D" GM=5;
# q, Z0 | t6 M. {' d4 j' Jb=ones(M,1)/M;
0 T6 m2 R+ N1 {2 Q3 u# J- N& Xa=1;
$ c. r8 E% ^" N: G* q7 xy=filter(b,a,f);* ^( Q7 Q$ s5 E4 L
figure(2)9 d& {% O, s2 \4 D) c
stem(k,y,'.')& H0 j( T9 t, K& I+ {0 i
1 c' m9 s j1 a
6 K4 y; p% b3 _
* j9 M& C! L8 Y9 ~3 Y* {# \% Z1 h& ?6 h7 |6 c3 E1 H0 z: C
|
|