|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为6 Y) I4 m& N6 f6 J1 K
8 j# e) U0 h! H2 |0 ^
y=filter(b,a,f)
7 [; ~4 q$ x, M" k {" ?1 c5 _: F3 w$ o
* T1 Q) X p% |
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
/ R; F. O) v- S @2 @& m8 v0 y1 F; a( [& t& n* m1 ^" b% M
' V& ~. L# v% O8 H2 S! ]; e
& s# \4 Q2 X) ]" p1 e* tR=51;
7 @( C4 k8 J2 q1 Q2 Ad=rand(1,R)-0.5;) t$ F2 A9 B+ M, [& |, Q3 B
k=0:R-1;
9 w& Z0 P5 A: y0 |& Ms=2*k.*(0.9.^k);
# r2 J2 I% B8 B6 H8 C$ j8 O& X$ mf=s+d;
+ ~( y" _2 l4 r7 _/ V# [0 I. Pfigure(1);
* q7 C* H" {4 k" A- Jstem(k,f,'.');; Z$ O* U* }2 J9 l+ l6 E
M=5;) S- A/ W, E( ~6 V& _
b=ones(M,1)/M;6 B: F% D" k1 E1 B0 G) E% p7 G
a=1;4 U' _8 s$ m M# m* W
y=filter(b,a,f);6 Z# w* Q# I0 G" Y z& G
figure(2)
% x( u% J1 T: F. P: F5 hstem(k,y,'.')
- x- D2 j! ^+ B/ H4 u4 _* {% r& z
! h% y: }0 y# k; P* P* M
~/ E1 U/ N/ `. W: D' A
+ J! ^4 l, n- Y$ c* U0 G$ k' } |
|