|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为, j" M$ g8 U- C3 k& f3 [
, \& U4 f4 B- b* C F* iy=filter(b,a,f)% m( |6 Z8 f* A& n1 l7 F, ]5 B
' R0 k- N: P. N9 |9 i6 S8 e
9 t" ]0 D! z. i9 I其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
) Z% I1 Y8 J. Q1 R% d$ t' r# r, w
% A- X) g/ ^3 J
1 a$ t/ u! Z+ E6 i, `
. G) M: [0 P( X7 P+ fR=51;; W7 S' O8 P: e( b+ l
d=rand(1,R)-0.5;9 M6 u, v% ]6 F
k=0:R-1;1 |+ J. Y( b+ s
s=2*k.*(0.9.^k);
- e- {% k8 i5 L, J% pf=s+d;
2 d% B: ^' k- S! p9 d1 Ofigure(1);% `) @3 _7 X+ j" E- P0 t0 P
stem(k,f,'.');
8 e" b# B% A. m: f1 m( wM=5;5 e+ u; X( a8 z8 ^. ~
b=ones(M,1)/M;
/ s/ G- J7 J& J2 xa=1;- h! w$ [) r# x& m# c# B0 O
y=filter(b,a,f);1 Q: c3 U# g, ` I+ A8 N$ }
figure(2)! \* ^3 _- G4 C$ p! H7 C
stem(k,y,'.')
S4 P' L, R) v! t H8 F5 b
* r# D9 @, x2 C! t0 @# X
3 ?, a( A |5 Y, x* C
' }) b9 A+ G! J3 {% j
1 f' S$ G8 Q; A! s6 f. v( K
|
|