6武建明数字信号处理.docx

数字信号处理实验报告学生姓名 武建明 专业班级 信息工程11-01 学号 541107070136院 （系） 计算机与通信工程学院实验一信号、系统及系统响应实验目的：1、熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对时域采样定理的理解。2、熟悉离散信号和系统的时域特性。3、熟悉线性卷积的计算编程方法；利用卷积方法观察分析系统的时域特性。4、掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析。1、理想采样信号序列(1)首先产生信号X(n)，0=n=50n=0:50;A=444.128;a=50*sqrt(2.0)*pi;T=0.001;w0=50*sqrt(2.0)*pi;x=A*exp(-a*n*T).*sin(w0*n*T);subplot(3,1,1);stem(x);title(理想采样信号序列);(2)绘制信号x(n)的幅度谱和相位谱k=-25:25;W=(pi/12.5)*k;f=(1/25)*k*1000;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,1,2);stem(magX);title(理想采样信号序列的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(理想采样信号序列的相位谱);(3)改参数为：A=1，a=0.4，Ω0=2.0734，T=1n=0:50;A=1;a=0.4;w0=2.0734;T=1;x=A*exp(-a*n*T).*sin(w0*n*T);subplot(3,1,1);stem(n,x);title(理想采样信号序列);k=-25:25;W=(pi/12.5)*k;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,1,2);stem(magX); title(理想采样信号序列的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(理想采样信号序列的相位谱);2、单位脉冲序列n=0:50;x=[1 zeros(1,50)];subplot(3,1,1);stem(n,x);title(单位脉冲信号序列);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,1,2);stem(magX);title(单位脉冲信号序列的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(单位脉冲信号序列的相位谱);3、矩形序列n=0:50;x=[ones(1,10) zeros(1,41)];subplot(3,1,1);stem(n,x);title(矩形序列);axis([0 50 0 1.2]);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n*k);magX=abs(X);subplot(3,1,2);stem(magX);title(矩形序列的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(矩形序列的相位谱);实验总结：在分析理想采样序列特性的实验中，采样频率不同，则其周期T不同，相应理想采样序列的傅立叶变换频谱的数字频率也不相同，而因为是同一信号，故它们所对应的模拟频率是相同的。在卷积定理验证的实验中，如果选用不同的频域采样点数M值，所得结果之间有差异。因为所得图形由其采样点数唯一确定，由频域采样定理可知，若M小于采样序列的长度N，则恢复原序列时会发生时域混叠现象。实验二应用FFT对信号进行频谱分析实验目的：1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉FFT算法及其程序的编写。2、熟悉应用FFT对典型信号进行频谱分析的方法。3、了解应用FFT进行信号分析的过程中可能出现的问题，以便在实际中正确应用FFT。1、观察高斯序列的时域和频域特性n=0:15;p=8;q=2;x=exp(-1*(n-p).^2/q);close all;subplot(3,2,1);stem(x);subplot(3,2,2);stem(abs(fft(x)))p=8;q=4;x=exp(-1*(n-p).^2/q);subplot(3,2,3);stem(x);subplot(3,2,4);stem(abs(fft(x)))p=8;q=8;x=exp(-1*(n-p).^2/q);subplot(3,2,5);stem(x);subplot(3,2,6);stem(abs(fft(x)))2、衰减正弦序列n=0:15;a=0.1;f=0.0625