《数字信号处理的实习和实验报告》.docVIP

数字信号处理实习报告指导老师：姓名：班级：学号：实验一 离散卷积的计算 一、实验内容 设线性时不变（LTI）系统的冲激响应为h(n)，输入序列为x(n) 1、h(n)=(0.8)n，0≤n≤4； x(n)=u(n)-u(n-4) 2、h(n)=(0.8)nu(n), x(n)=u(n)-u(n-4) 3、h(n)=(0.8)nu(n), x(n)=u(n) 求以上三种情况下系统的输出y(n)。 二、实验程序及实验结果 1．1 matlab程序 x = [ones(1,4)]; N1 = length(x); n1 = 0:N1-1; N2 = 5; n2 = 0:N2-1; h = 0.8.^n2; y = conv(x,h); N = N1+N2-1; n = 0:N-1; subplot(2,2,1);stem(n1,x);title(序列x); xlabel(n);ylabel(x(n)); subplot(2,2,2);stem(n2,h);title(序列h); xlabel(n);ylabel(h(n)); subplot(2,2,3);stem(n,y);title(两序列的卷积); xlabel(n);ylabel(y(n)); 1.2实验结果 2.1 matlab程序 x = [ones(1,4)]; N1 = length(x); n1 = 0:N1-1; N2 = 100; n2 = 0:N2-1; h = 0.8.^n2; y = conv(x,h); N = N1+N2-1; n = 0:N-1; subplot(1,2,1);stem(n1,x);title(序列x); xlabel(n);ylabel(x(n)); subplot(1,2,2);stem(n2,h);title(序列h); xlabel(n);ylabel(h(n)); figure,stem(n,y);title(两序列的卷积); xlabel(n);ylabel(y(n)); 2.2实验结果3.1 matlab程序 x = [ones(1,100)]; N1 = length(x); n1 = 0:N1-1; N2 = 100; n2 = 0:N2-1; h = 0.8.^n2; y = conv(x,h); N = N1+N2-1; n = 0:N-1; subplot(2,2,1);stem(n1,x);title(序列x); xlabel(n);ylabel(x(n)); subplot(2,2,2);stem(n2,h);title(序列h); xlabel(n);ylabel(h(n)); figure,stem(n,y);title(两序列的卷积); xlabel(n);ylabel(y(n)); 3.2实验结果实验结果分析 实验程序通过直接调用卷积函数实现卷积运算，题目中h(n)=(0.8)nu(n)理论上为一个无限长序列，但在matlab的编程中，只有有限长的序列才可以参与运算，因此只是选取了有限点（100点）的h（n）。 实验二 离散傅里叶变换及其应用 一．实验内容 设有离散序列 x(n)=cos(0.48πn)+cos(0.52πn) 分析下列三种情况下的幅频特性。 采集数据长度N=16，分析16点的频谱，并画出幅频特性。采集数据长度N=16，并补零到64点，分析其频谱，并画出幅频特性。 采集数据长度N=64，分析46点的频谱，并画出幅频特性。 观察三幅不同的幅频特性图，分析和比较它们的特点及形成原因。 二．实验程序及实验结果 1.1matlab程序 clc; N1 = 16; n1 = 0:15; x = cos(0.48*pi*n1)+cos(0.52*pi*n1); stem(n1,x);title(16采样序列x); xlabel(n);ylabel(x(n)); y = fft(x); Y = abs(y); figure,subplot(121);stem(n1,Y);title(16点频谱); xlabel(\omega);ylabel(X); y1 = fft(x,64); Y1 = abs(y1); subplot(122);stem(0:63,Y1);title(64点频谱); xlabel(\omega);ylabel(X); N2 = 64; n2 = 0:N2-1; x2 = cos(0.48*pi*n2)+cos(0.52*pi*n2); figure,subplot(121);stem(n2,x2);title(64点采样序列x); xlabel(n);ylabel(x(n)); y2 = fft(x,46); Y2 = abs(y2); subplot(122