|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为
. D5 k0 f2 C$ R" s2 w8 c. ?# f- U8 R6 W; b1 f
y=filter(b,a,f)2 G7 d2 ?( i6 ]9 P$ n
& B4 j# g. w8 u% `' @
$ Q, z: I; [! {3 D6 F2 n- e. e/ \其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。* Q( J( a1 W ~: N; [
# i, H" ]% j* X+ m5 ^4 Q. }! v g
) q7 ^- H1 ]0 I2 ?' i5 R1 P9 @0 m! {
R=51;
: ?! X; K2 U4 [8 @ cd=rand(1,R)-0.5;
" z# q, b# s( x# rk=0:R-1;
( w% l) i9 \8 D: Ls=2*k.*(0.9.^k); X. P/ e7 A2 W' W& r1 \" O3 O
f=s+d;9 o! j" F& L5 _- X) s( Z
figure(1);
/ W. G3 x6 `, x+ ^& K% P* Lstem(k,f,'.');
6 B- T& V- B4 }$ AM=5;$ Y. s' V9 y7 |# T$ k P
b=ones(M,1)/M;
0 x6 G: W4 b7 }7 Q( Qa=1;, O9 p* z$ x5 L4 v3 [0 G* L# P
y=filter(b,a,f);9 p Z8 o1 \/ M) I5 e& m1 ?
figure(2)
. [. Q Y) a( {: Rstem(k,y,'.')" H( c) L; Y4 N+ v, n
* }+ h5 O* `: Q1 j" f5 ~6 I
5 O, l8 w r3 ~8 E. O
1 R% V/ N1 m4 }* K5 ^1 ~3 D1 u; L
5 M& G0 P# _2 t1 k) P |
|