DSP数字信号处理实验报告材料优秀.docx
- 文档编号:6556459
- 上传时间:2023-01-07
- 格式:DOCX
- 页数:24
- 大小:609.82KB
DSP数字信号处理实验报告材料优秀.docx
《DSP数字信号处理实验报告材料优秀.docx》由会员分享,可在线阅读,更多相关《DSP数字信号处理实验报告材料优秀.docx(24页珍藏版)》请在冰豆网上搜索。
DSP数字信号处理实验报告材料优秀
FFT频谱分析
一、实验目的
a)进一步加深DFT算法原理和基本性质的理解
b)熟悉FFT算法原理和FFT程序的应用
c)学习用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确的应用FFT
二、实验原理
a)离散傅里叶变换(DFT):
离散傅里叶变换在作为有限长序列的傅里叶变换表示法在理论上相当重要;由于存在着计算离散傅里叶变换的快速算法(FFT),从而离散傅里叶变换在各种数字信号处理的算法中起到了核心的作用。
其对应的离散傅里叶变换对为:
需要注意:
有限长序列的离散傅里叶变换及周期序列的离散傅里叶级数之间的关系是:
它们仅仅是n、k的取值不同,DFT只取主值区间。
X(n)、X(k)是一个有限长序列的离散傅里叶变换对,已知其中一个序列,就可以唯一确定另一个序列,这是因为x(n)、X(k)都是长为N的序列,都有N个独立值,所以信息量相同。
b)DFT计算量:
长度为N的DFT的计算量是N个复数乘法和N-1个复数加法(4N个实数乘法和4N-2个复数加法)
c)FFT的计算量:
长度为N的FFT的计算量是
个复数乘法和
个复数加法。
三、实验步骤
a)复习DFT的定义、性质和用DFT做谱分析的有关内容
b)复习FFT算法原理与编程思想,熟悉DIT-FFT运算流图
c)编制信号产生程序,产生典型信号尽心谱分析。
d)进行以下几个信号的谱分析
i.
ii.
iii.
iv.
v.
vi.
vii.令
N=8,16
viii.令
N=8,16
针对上述信号进行逐一的谱分析,下面给出针对各个信号的FFT点数N及对连续信号
的采样频率
,供实验时参考:
N=8,16
四、实验内容
a)对
进行谱分析
1.编辑代码
x1=[1111];
y11=fft(x1,8);
y12=fft(x1,16);
subplot(2,2,1);
stem(0:
3,x1);
title('函数X1的图像');
subplot(2,2,2);
stem(0:
7,abs(y11));
title('N=8的DFT');
subplot(2,2,4);
stem(0:
15,abs(y12));
title('N=16的DFT');
2.谱分析图片
b)对
进行谱分析
i.编辑代码
x2=[12344321];
y11=fft(x2,8);
y12=fft(x2,16);
subplot(2,2,1);
stem(0:
7,x2);
title('函数X2的图像');
subplot(2,2,2);
stem(0:
7,abs(y11));
title('N=8的DFT');
subplot(2,2,4);
stem(0:
15,abs(y12));
title('N=16的DFT');
ii.谱分析图像
c)对
进行谱分析
i.谱分析程序
x3=[43211234];
y11=fft(x2,8);
y12=fft(x2,16);
subplot(2,2,1);
stem(0:
7,x2);
title('函数X2的图像');
subplot(2,2,2);
stem(0:
7,abs(y11));
title('N=8的DFT');
subplot(2,2,4);
stem(0:
15,abs(y12));
title('N=16的DFT');
ii.谱分析图片
d)对
进行谱分析
i.谱分析程序
n=0:
1:
19;
x2=cos(0.25*pi*n);
y11=fft(x2,32);
y12=fft(x2,64);
subplot(2,2,1);
stem(0:
19,x2);
title('函数X4的图像');
subplot(2,2,2);
stem(0:
31,abs(y11));
title('N=32的DFT');
subplot(2,2,4);
stem(0:
63,abs(y12));
title('N=64的DFT');
ii.谱分析图片
e)对
进行谱分析
i.谱分析程序
n=0:
1:
19;
x2=sin(0.125*pi*n);
y11=fft(x2,32);
y12=fft(x2,64);
subplot(2,2,1);
stem(0:
19,x2);
title('函数X5的图像');
subplot(2,2,2);
stem(0:
31,abs(y11));
title('N=32的DFT');
subplot(2,2,4);
stem(0:
63,abs(y12));
title('N=64的DFT');
ii.谱分析图片
f)对
进行谱分析
i.谱分析程序
n=0:
1:
15;
x1=cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);
n=0:
1:
31;
x2=cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);
n=0:
1:
63;
x3=cos(8*pi*n/64)+cos(16*pi*n/64)+cos(20*pi*n/64);
y1=fft(x1,16);
y2=fft(x2,32);
y3=fft(x3,64);
subplot(3,2,1);
stem(0:
15,x1);
title('函数X6N=16的图像');
subplot(3,2,2);
stem(0:
15,abs(y1));
title('N=16的DFT');
subplot(3,2,3);
stem(0:
31,x2);
title('函数X6N=32的图像');
subplot(3,2,4);
stem(0:
31,abs(y2));
title('N=32的DFT');
subplot(3,2,5);
stem(0:
63,x3);
title('函数X6N=64的图像');
subplot(3,2,6);
stem(0:
63,abs(y3));
title('N=64的DFT');
ii.谱分析图片
g)对
进行谱分析
i.谱分析程序
n=0:
1:
19;
x2=j*sin(0.125*pi*n)+cos(0.25*pi*n);
y11=fft(x2,8);
y12=fft(x2,16);
subplot(2,2,1);
stem(0:
19,x2);
title('函数X=X4+X5的图像');
subplot(2,2,2);
stem(0:
7,abs(y11));
title('N=8的DFT');
subplot(2,2,4);
stem(0:
15,abs(y12));
title('N=16的DFT');
ii.谱分析图片
h)对
进行谱分析
i.谱分析程序
n=0:
1:
19;
x2=j*sin(0.125*pi*n)+cos(0.25*pi*n);
y11=fft(x2,8);
y12=fft(x2,16);
subplot(2,2,1);
stem(0:
19,x2);
title('函数X=X4+X5的图像');
subplot(2,2,2);
stem(0:
7,abs(y11));
title('N=8的DFT');
subplot(2,2,4);
stem(0:
15,abs(y12));
title('N=16的DFT');
ii.谱分析图片
五、实验分析
六、实验结论
通过这次利用FFT对信号进行频谱分析的实验,更加深刻的理解了DFT算法的理解和性质的理解;同时也更为熟悉了FFT算法的原理和应用;学会使用FFT对离散信号和连续信号进行频谱分析,了解了可能出现的分析误差和原因。
七、思考题
a)在N=8时,
和
的幅频特性会相同吗?
为什么?
N=16呢?
答:
在N=8的时候,幅频特性相同。
1.知X(k)=
2.由表达式和图像可以知道
3.
4.由离散傅里叶变换的性质:
序列的循环位移性质
可以得到:
5.可以知道:
a.当N=8:
不改变幅度特性,只改变相位特性,则
b.当N=16:
由于进行了补零处理,
和
序列就无法满足
,即不满足序列的循环位移,导致
b)对一个周期序列分别取一个周期或多个周期所做的频谱分析结果是否有区别?
答:
有区别。
1.知道DFT的计算为X(k)=
2.当N越大,所计算得到的在X(k)在0~2π范围内有更多的N的点数,这意味着频域的频率分辨率得到提高,使得频谱信息更为精细丰富。
窗函数设计FIR滤波器
一、实验目的
a)熟悉熟悉矩形窗、汉宁窗、海明窗等常用窗函数
b)掌握用上述窗函数法设计FIR数字滤波器的原理和方法
c)熟悉线性相位FIR数字滤波器特性
d)了解各种窗函数对滤波特性的影响
二、实验原理
a)如果所希望的滤波器的理想频率响应函数为
,则其对应的单位冲激响应为
用窗函数(n)将
截断,得到:
h(n)就作为实际设计的FIR数字滤波器的单位冲激响应序列,其频率响应函数为
如果要求线性相位特性,则h(n)还必须满足:
根据上式中的正、负号和长度N取为奇数或偶数又将线性相位FIR滤波器分成四类。
要根据所设计的滤波特性正确选择其中一类。
三、实验步骤
a)复习用窗函数法设计FIR数字滤波器一节内容,阅读本实验原理,掌握设计步骤
b)编写程序,画出滤波器的幅频特性和相频特性。
其中幅频特性要求用分贝(dB)表示
设
画图时,用
画出幅度特性
注意:
补零的作用
N=15,N=33,
用三种窗函数设计线性相位低通滤波器。
要求在两种窗口长度下,绘制相应的幅频特性和相频特性曲线,比较三种窗函数对滤波器特性的影响。
四、实验内容
a)使用矩形窗设计
i.设计代码
clc;
clear;
wc=pi/4;
N=33;
N1=33;
a=(N-1)/2;
n=0:
(N-1);
hdn=sin(wc*(n-a+eps))./(pi*(n-a+eps));
wn=rectwin(N);
hn=hdn.*wn';
fn=1024;
H=fft(hn,fn);
Hdb=20*log10((abs(H)+eps
(1))./(max(abs(H)+eps
(1))));
w=0:
2*pi/fn:
2*pi*(fn-1)/fn;
subplot(2,3,1);
stem(n,hdn);
xlabel('n');
ylabel('hd(n)');
axis([0,N,-0.1,0.3]);
title('hd(n)');
subplot(2,3,2);
stem(n,wn);
xlabel('n');
ylabel('w(n)');
axis([0,N,-0.2,1.2]);
title('w(n)');
subplot(2,3,3);
stem(n,hn);
xlabel('n');
ylabel('h(n)');
axis([0,N,-0.1,0.3]);
title('h(n)');
subplot(2,3,4);
plot(w/pi,Hdb);
xlabel('w/pi');
ylabel('20log');
%axis([0,1,-20,0]);
title('幅度');
subplot(2,3,5);
plot(w,unwrap(angle(H)));
xlabel('w/rad');
ylabel('相位/rad');
axis([0,pi,-20,0]);
title('相位');
ii.图像分析
矩形窗N=16
矩形窗N=33
b)使用汉宁窗设计
i.代码设计
clc;
clear;
wc=pi/4;
N=33;
N1=33;
a=(N-1)/2;
n=0:
(N-1);
hdn=sin(wc*(n-a+eps))./(pi*(n-a+eps));
wn=hanning(N);
hn=hdn.*wn';
fn=1024;
H=fft(hn,fn);
Hdb=20*log10((abs(H)+eps
(1))./(max(abs(H)+eps
(1))));
w=0:
2*pi/fn:
2*pi*(fn-1)/fn;
subplot(2,3,1);
stem(n,hdn);
xlabel('n');
ylabel('hd(n)');
axis([0,N,-0.1,0.3]);
title('hd(n)');
subplot(2,3,2);
stem(n,wn);
xlabel('n');
ylabel('w(n)');
axis([0,N,-0.2,1.2]);
title('w(n)');
subplot(2,3,3);
stem(n,hn);
xlabel('n');
ylabel('h(n)');
axis([0,N,-0.1,0.3]);
title('h(n)');
subplot(2,3,4);
plot(w/pi,Hdb);
xlabel('w/pi');
ylabel('20log');
%axis([0,1,-20,0]);
title('幅度');
subplot(2,3,5);
plot(w,unwrap(angle(H)));
xlabel('w/rad');
ylabel('相位/rad');
axis([0,pi,-20,0]);
title('相位');
i.图像分析
汉宁窗N=15
汉宁窗N=33
c)使用海明窗设计
i.代码设计
clc;
clear;
wc=pi/4;
N=33;
N1=33;
a=(N-1)/2;
n=0:
(N-1);
hdn=sin(wc*(n-a+eps))./(pi*(n-a+eps));
wn=hamming(N);
hn=hdn.*wn';
fn=1024;
H=fft(hn,fn);
Hdb=20*log10((abs(H)+eps
(1))./(max(abs(H)+eps
(1))));
w=0:
2*pi/fn:
2*pi*(fn-1)/fn;
subplot(2,3,1);
stem(n,hdn);
xlabel('n');
ylabel('hd(n)');
axis([0,N,-0.1,0.3]);
title('hd(n)');
subplot(2,3,2);
stem(n,wn);
xlabel('n');
ylabel('w(n)');
axis([0,N,-0.2,1.2]);
title('w(n)');
subplot(2,3,3);
stem(n,hn);
xlabel('n');
ylabel('h(n)');
axis([0,N,-0.1,0.3]);
title('h(n)');
subplot(2,3,4);
plot(w/pi,Hdb);
xlabel('w/pi');
ylabel('20log');
%axis([0,1,-20,0]);
title('幅度');
subplot(2,3,5);
plot(w,unwrap(angle(H)));
xlabel('w/rad');
ylabel('相位/rad');
axis([0,pi,-20,0]);
title('相位');
ii.图像分析
海明窗N=15
海明窗N=33
五、实验分析
通过对于三种滤波器的设计,通过其对应的频谱分析图像,可以比较容易地看出,
六、实验结论
通过窗函数设计FIR数字滤波器的过程,熟悉了矩形窗、汉宁窗、海明窗常用的窗函数;熟悉了窗函数设计法的基本原理和方法;熟悉了线性相位数字滤波器的特性;通过FFT对窗函数的频谱分析,了解到了各种窗函数对滤波特性的影响。
总的来说,加深了对书本中关于FIR数字滤波器设计部分的知识的感性理解,提高了相关的实验技能,提高了我们的专业知识的实践能力。
七、思考题
a)如果给定通带截至频率和阻带截止频率以及阻带最小衰减,如何用窗函数设计线性相位低通滤波器?
写出设计步骤?
答:
1.通过阻带最小衰减指标选择合适的窗函数
2.由给定的通带截止频率和阻带截止频率得出过渡带宽
3.通过过渡带宽选择合适的窗长度N(N取奇数)
4.得到理想低通滤波器的单位冲激响应
5.通过查表得到所需要的窗函数的表达式
,然后与理想低通滤波器相乘,得到所想设计的实际滤波器:
b)如果要求用窗函数法设计带通滤波器,且给定上、下带截止频率为
与
试求理想带通的单位冲激响应。
答:
1.知
2.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- DSP 数字信号 处理 实验 报告 材料 优秀
![提示](https://static.bdocx.com/images/bang_tan.gif)