数字信号处理-实验一离散系统的时域分析

《数字信号处理-实验一离散系统的时域分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验一离散系统的时域分析实验目的:（1）熟悉并掌握离散系统的差分方程表示法；（2）加深对冲激响应和卷积分析方法的理解。二、实验原理与方法离散系统x[n>Discrete-timesystme其输入、输出关系可用以下差分方程描述：/VMy[n-k]=^pkx[n-k]々=()k=G输入信号分解力冲激信号，OOx[n]=x[m]8[n-m]o州=一OO记系统单位冲激响应S[n]h[n],则系统响应为如下的卷积计算式：OOy[n]=x[n]^h[n]=^x[m]h[n-m]当<=0，々=l，2，..JV时，h[nl是有限长度的（n:[0,Mp,称系统为FIR系统;反之，称

2、系统为IIR系统。三、实验内容及要求1、要求分别用filter、conv、impz三种函数完成。（1）用conV（x，y）函数计算序列｛-201-13｝和序列｛120-1｝的离散卷积程序如下：xn=[-201-13];hn=[l20-1];yn=conv(xn,hn);结果如下:1-3CDACommandWindowFileEditDebugDesktopWindowHelpyn-2-4A»(2)用yn=filter(B,A,xn)函数计算差分方程沖)+0.7Xi-Q-⑽沖-2)-0.6X*-X=aSz(R)-0.44冲-1)*0.36场-2)+102jc(n-3)当

3、输入序列为《(

4、<

5、=逐(#时的输出结果，C3x340。程序如下：xn=[l，zeros(l，40)];B=[l0.7-0.45-0.6];A=[0.8-0.440.360.22];yn=filter(A,B,xn)n=O:length(yn)-l;subplot(1,1,1);stem(n,yn,7)title('(a),);xlable(,n,);ylable('y(n)')结果如下：FileEditDebugDesktopWindowhdelpyn=Columns1through50.8000-1.00001.4200-0.94400.6998Columns6t

6、hrough10-0.0627-0.20760.5370-0.50690.4719Colujnns11through15-0.23630.07360.1253-0.19640.2380Columns16through20-0.17980.1151-0.0187-0.04300.0908Columns21through25-0.09410.0809-0.04450.01110.0207Columns26through30-0.03620.0414-0.03280.0198-0.0038Columns31through35-0.00810.0158-0.01700.014

7、2-0.0081Coliunns36through400.00180.0036-0.00650.0073-0.0059Column410.0035A»ACommandWindow口回图如下:(3)用impz(B,A,N)绘制出差分方程^»)4-0.7X>-Q-QL45Xk-2)-0.6X*-X=Q.Sz(ji)-0.44j(M-D4-O.Kz(«-2)4-102je(«-3)所对应的单位脉冲响应/?(/!),0

8、th(hn)-1;subplot(1,1,1);stem(n,hn,’.’)title(’⑹’);xlabel(’n’)；ylabel(’y(n)’)结果如下：人CommandWindowFileEditDebugDesktopWindowHelp0.8000-1.00001.4200-0.94400.6998-0.0627-0.20760.5370-0.50690.4719-0.23630.07360.1253-0.19640.2380-0.17980.1151-0.0187-0.04300.0908图如下:四、实验报告要求(1)简述实验原理，画出程序框图，列出实验

9、程序淸单，并附上必要的程序说明。(2)记录调试运行情况及所遇问题的解决方法。(3)画出实验结果图，实验后，对结果进行分析。