|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为* L, F" V' s6 r9 y" _+ n5 z+ G
. k) ^$ p- K3 b- m) f% Iy=filter(b,a,f)
7 |1 T+ D/ f+ g0 B# k$ I, K( Z* H( X" M4 P; H; l% ~( l' i
* D& l+ K& ?& M, r$ d# O
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。# t! ?& x* ?& g8 G2 B9 ]
* I4 h* B7 X1 E) D
$ T6 k% J! I7 p: g& G/ S s
. m* o4 b# j% s. X; {R=51;
* O. ^& _. C- @& ? h% s+ Z* Wd=rand(1,R)-0.5;
# V) H: ?# U* Z) f; a2 ?; hk=0:R-1;
0 X6 S( ?& [0 b% F: Xs=2*k.*(0.9.^k);$ P5 K6 i" @1 R' y
f=s+d;
! S) K+ C' d3 cfigure(1);) K8 V. u9 A# |2 q. p: m
stem(k,f,'.');1 X& r8 Z! f7 X% f7 _2 l* y
M=5;
, B* E% G3 c( r# m& T& ib=ones(M,1)/M;4 z2 j8 i I$ Z# Z* M9 o
a=1;
7 i8 m& }! J, V9 Ty=filter(b,a,f);
% d b& Y, \% Pfigure(2)6 T9 ?+ {$ h) [! F
stem(k,y,'.')
2 p' h# J- u% {3 o7 o2 f |) R2 [) T |
6 Z9 ]! z' B- H
' R8 s' y& C8 f% {
6 \4 N6 l# @2 C: y$ Y5 c9 z& H |
|