实验2离散系统的差分方程冲激

《实验2离散系统的差分方程冲激》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验课程：数字信号处理实验内容：实验2离散系统的差分方程、冲激响应和卷积分析2013年6月3日一、实-验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。二、实验原理：离散系统x[n]Discrctc-timcsystmeyW►其输入、输出关系可用以卜*差分方程描述:NMzdky[n-k}-zpkx[n-k}众=0^=000输入信号分解为冲激信号，记系统单位冲激响应S[n]A[n]，则系统响应为如下的卷积计算式：00)[打]=刺*制=h[雌相当=0，众=1，2，...7V时，h[n]是有限L<:度的(n:[

2、0,M])，称系统为

3、实验代码及结果(lM«]+0J5Xn-l]+0.125y[n-2]=^n]-^n-l]程序代码：%(1)用impz()函数求冲激响应：M=50;num=[l-10];den=[l0.750.125];y=impz(num,den,M);subplot(3,2,1);stem(y);xlabel(’时间’)；ylabel(’振幅’);titleC雄位冲激响应：impzO函数方法’)%(2)用filter()函数求冲激响应：n=0:50;x=[lzeros(1,50)];num=[l-10];den=[l0.750.125

4、];y=filter(num,den,x);%filter函数给出的？^数与输入的x序列？^数一样，所以为了不漏点，输入序列的点数尽量多补0,但用此函数求冲激响应不好subplot(3,2,2);stem(n,y);xlabel(’时间’)；ylabel(’振幅’);title(’单位冲激响应:filter()函数方法’)%(3)用filter()函数求余弦输入响应：n=0:50;x=cos(2*pi*0.47*n);num=[l-10];den=[l0.750.125];y=filter(num,den,x);%注意

5、：“〉〉”后加一个空格再写指令，这样就可以直接COPY到MATLAB的命令窗执行，％是注释符号subplot(3,2,3);stem(n,y)xlabel(’时间’)；ylabel(’振幅’);title(’余弦输入响应：filter()函数方法’)%(4)用conv()函数求余弦响应：n=0:50;xl=cos(2*pi*0.47*n);x2=[lzeros(1,50)];num=[l-10];den=[l0.750.125];y=filter(num,den,x2);z=conv(xl,y);subplot(3,2

6、,4);m=0:100;stem(in,z);axis([050-1010]);%z的点数为xl的加上y(也即X2的)的再减1%convO可以显示201点，而用filter()仅显示101点%为了便于比较两种方法的值，减小本方法的横坐标尺度与filter()方法匹配xlabel(’时间’)；ylabel(’振幅’);title(’余弦输入响应：convO函数方法’)%(5)系统零极点分布图：num=[l-10];den=[l0.750.125];subplot(3,2,5);zplane(num,den);grid图形

7、如下:Figure1FileEditViewInsertToolsDesktopWindowHelp单位冲激响应：impz^函数方法2204060时间余弦输人响应:filt邮函数方法1°•2皇单位冲激响应：filte⑴函数方法20204060时间余弦输入响应：ccnvO函数方法10S2040时间11020304050时间0•2101RealPart2IJECLhJaJCIGBUJI网络Intei(2)y[n]=0.25{x[n-l]+x[n-2]+x[n—3]+x["—4]}程序如下：M=50;num=[00.250.

8、250.250.25];den=[l0000];y二impz(num,den,M);subplot(3,2,1);stem(y);xlabel(’时间’)；ylabel(’振幅’);title(’雄位冲击响应：impz()函数方法’)；n=0:50;x=[lzeros(1,50)];num=[00.250.250.250.25];