文档介绍：数字信号处理
实验报告
学生姓名
孟令江
班级
电子信息工程1301班
学号
1312140103
指导教师
郑伟华
实验一快速傅里叶变换及其应用
一、实验目的
在理论学习的基础上,通过本实验,加强对FFT的理解,熟悉MATLAB中有关函数。
应用FFT对典型信号进行频谱分析。
了解应用FFT进行信号频谱分析过程中可能出现的问题,以便在实际中正确应用FFT.
应用相关FFT实现序列的线性卷积和相关。
二、实验内容
实验用到的信号序列:
高斯序列
衰减正弦序列
三角波序列
反三角波序列
上机实验内容:
观察高斯序列的时域和幅频特性,固定信号中参数p=8,改变q的值,使q分别等于2、4、8,观察他们的时域和幅频特性,了解当q去不同值时,对信号序列的时域和频域特性的影响;固定q=8,改变p,使p分别为等于8、13、14,观察参数p变化对信号序列的时域以及频域特性的影响,注意p等于多少时,会发生明显的泄露现象,混叠是否也随之出现?记录实验中观察到的现象,绘制相应的时域序列和幅频特性曲线。
代码如下:
clear all;
%p=8,q=2%
[Xa1,Fa1]=gauss(8,2);
k=0:15;
subplot(5,2,1);
plot(k,Xa1);
Xlabel('n');
Ylabel('时域特性');
text(10,,'p=8,q=2');
subplot(5,2,2);
plot(k,Fa1);
Xlabel('n');
Ylabel('幅频特性');
text(8,3,'p=8,q=2');
%p=8,q=4%
[Xa1,Fa1]=gauss(8,4);
k=0:15;
subplot(5,2,3);
plot(k,Xa1);
Xlabel('n');
Ylabel('时域特性');
text(10,,'p=8,q=4');
subplot(5,2,4);
plot(k,Fa1);
Xlabel('n');
Ylabel('幅频特性');
text(8,3,'p=8,q=4');
%p=8,q=8%
[Xa1,Fa1]=gauss(8,8);
k=0:15;
subplot(5,2,5);
plot(k,Xa1);
Xlabel('n');
Ylabel('时域特性');
text(10,,'p=8,q=8');
subplot(5,2,6);
plot(k,Fa1);
Xlabel('n');
Ylabel('幅频特性');
text(8,3,'p=8,q=8');
%p=8,q=8%
[Xa1,Fa1]=gauss(13,8);
k=0:15;
subplot(5,2,7);
plot(k,Xa1);
Xlabel('n');
Ylabel('时域特性');
text(10,,'p=13,q=8');
subplot(5,2,8);
plot(k,Fa1);
Xlabel('n');
Ylabel('幅频特性');
text(8,3,'p=13,q=8');
%p=8,q=8%
[Xa1,Fa1]=gauss(14,8);
k=0:15;
subplot(5,2,9);
plot(k,Xa1);
Xlabel('n');
Ylabel('时域特性');
text(10,,'p=14,q=8');
subplot(5,2,10);
plot(k,Fa1);
Xlabel('n');
Ylabel('幅频特性');
text(8,3,'p=14,q=8');
观察衰减正弦序列的时域和幅频特性,a=,f=,检查峰值出现的位置是否正确,注意频谱的形状,绘出频谱特性曲线,改变发,,观察这两种情况下,频谱的形状和频谱出现的位置,有无混叠、泄露现象?说明产生现象的原因。
代码如下:
clear all;
a=;
f=;
i=0:160
x=exp(-a*(i)).*sin(2*pi*f*(i));
subplot(3,2,1);
plot(i,x);
xlabel('n');
title('a=,f= 时域特性')
subplot(3,2,2);
G=fft(x);
plot(i,abs(G));
xlabel('k');
title('a=,f= 幅频特性')
f=;
x=exp(-a*(i)).*sin(2*pi*f*(i));
subplot(