|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为
1 Z u, t' e* o- Z: e ?9 o K8 U# p! p
y=filter(b,a,f)
' `4 w& R- ]. q, [& ^: C, z, @* u- M- b, t
. R# L. r1 z x% [6 k: i; n. L其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
8 A0 R. @) {- p0 y* U* U7 u# L
$ f1 b0 J, \7 s- }+ T% H
. |2 \: X& K, A9 n( y- d- P. j
R=51;" D2 f% c+ ]. q, I3 g6 K3 m) k
d=rand(1,R)-0.5;
8 ]$ k) i) i# l, nk=0:R-1;* h3 w4 ^3 f: O! C$ H3 ~2 N4 T
s=2*k.*(0.9.^k);: v) q, k2 o8 i
f=s+d;
1 F/ M1 `4 ~5 Cfigure(1);, @3 s" x% D! M. M
stem(k,f,'.');
6 ~: J) n+ ^& O, [$ c# V! O- JM=5;0 M7 N$ J: D7 x7 k0 Z
b=ones(M,1)/M;
8 e, U8 o' |4 ]0 za=1;
! h# X" F7 p0 Z# ?& v8 G; my=filter(b,a,f);( C9 b3 N k) k0 U+ ~
figure(2)( ^ o: V5 z% c. X
stem(k,y,'.')
7 G( @% A$ e' o' N; z9 {5 W' p5 W1 y0 K9 `3 `/ c8 v$ I! b
1 i* X2 b0 h; a e' W# ?! ~. r$ z
$ u4 T; ~$ P! B+ W2 o
|
|