实验二 应用FFT对信号进行频谱分析Word文档格式.docx
- 文档编号:15918710
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:17
- 大小:283.10KB
实验二 应用FFT对信号进行频谱分析Word文档格式.docx
《实验二 应用FFT对信号进行频谱分析Word文档格式.docx》由会员分享,可在线阅读,更多相关《实验二 应用FFT对信号进行频谱分析Word文档格式.docx(17页珍藏版)》请在冰豆网上搜索。
二、实验原理与方法
(1)复习DFT的定义、性质和用DFT作谱分析的有关内容。
(2)复习FFT算法原理与编程思想,并对照DIT-FFT运算流图和程序框图,读懂本实验提供的FFT子程序。
(3)编制信号产生子程序,产生以下典型信号供谱分析用:
应当注意,如果给出的是连续信号xa(t),则首先要根据其最高频率确定采样速率fs以及由频率选择采样点数N,然后对其进行软件采样(即计算择要以能分辨开其中的三个频率对应的谱线为准则。
对周期序列,最好截取周期的整数倍进行分析,否则有可能产生较大的分析误差。
请实验者根据DFT的隐含周期性思考这个问题。
(4)编写主程序。
三、上机实验内容及实验结果图
(1)对2中所给出的信号逐个进行谱分析。
下面给出针对各信号的FFT变换区间N以及对连续信号x6(t)的采样频率fs。
x1(n),x2(n),x3(n),x4(n),x5(n),:
N=8,16
x6(t):
fs=64(hz),N=16,32,64
x1(n)=R4(n)
对应x1(n)的matlab源程序代码:
n=[0:
7];
x=[11110000]
f1=fft(x,8)
f2=fft(x,16)
subplot(2,2,1)
stem(n,x);
axis([0802])
xlabel('
n'
)
ylabel('
xl(n)'
title('
xl的波形'
subplot(2,2,4)
k=[0:
15]
stem(k,abs(f2));
axis([01605])
k'
|xl(k)|'
xl(n)的16点fft'
subplot(2,2,3)
7]
stem(k,abs(f1));
axis([01005])
xl(n)的8点fft'
得到如下波形图:
2.
对应x2(n)的matlab源程序代码:
x=[12344321]
axis([0804])
x2(n)'
x2的波形'
axis([016020])
|x2(k)|'
x2(n)的16点fft'
axis([010020])
x2(n)的8点fft'
3.
对应x3(n)的matlab源程序代码:
x=[43211234]
x3(n)'
x3的波形'
|x3(k)|'
x3的16点fft'
axis([08020])
x3(n)的8点fft'
4.x4(n)=cos(pi/4*n)
对应x4(n)的matlab源程序代码:
x=cos(0.25*pi*n)
axis([08-44])
x4(n)'
x4的波形'
subplot(2,2,2)
axis([016-44])
|x4(k)|'
x4(n)的16点fft'
|x4(k)'
x4(n)的8点fft'
X5(n)=sin(pi/8*n)
对应x5(n)的matlab源程序代码:
clear
x=sin((pi*n)/8)
x5(n)'
x5的波形图'
x=sin(0.125*pi*n)
|x5(k)|'
x5的16点fft'
x5的8点fft'
X6(t)=cos(pi*8*n)+sin(pi*16*t)+cos(20*pi*t)
对应x6(n)的matlab源程序代码:
Ts=1/16;
n=0:
15;
Xa=cos(8*n*Ts*pi)+cos(16*n*Ts*pi)+cos(20*n*Ts*pi);
f1=fft(Xa,16);
subplot(3,2,1);
stem(n,Xa);
axis([015-23]);
);
x6(n)'
x6(n)N=16'
)%Ï
Ô
Ê
¾
x6£
¨
n£
©
N=16;
k=0:
15
subplot(3,2,2);
axis([016015])
|x6(k)|'
x6(n)N=16的16点fft'
31;
Xb=cos(8*n*Ts*pi)+cos(16*n*Ts*pi)+cos(20*n*Ts*pi);
f2=fft(Xb,32);
subplot(3,2,3);
stem(n,Xb);
axis([032-23]);
x6(n)N=32'
N=32;
subplot(3,2,4);
axis([032020])
x6(n)N=32的32点fft'
63;
Xc=cos(8*n*Ts*pi)+cos(16*n*Ts*pi)+cos(20*n*Ts*pi);
f3=fft(Xc,64);
subplot(3,2,5);
stem(n,Xc);
axis([064-23]);
N=64;
subplot(3,2,6);
stem(k,abs(f3));
axis([064040])
x6(n)N=64的64点fft'
(2)令x7(n)=x4(n)+x5(n),用FFT计算8点和16点离散傅里叶变换,
X(k)=DFT[x(n)]
x=cos(0.25*pi*n)+sin(0.125*pi*n)
x7(n)'
x7的波形图'
|x7(k)|'
x7(n)的16点fft'
x7(n)的8点fft'
(3)令x8(n)=x4(n)+jx5(n),重复
(2)。
对应x8(n)的matlab源程序代码:
15];
x=cos(0.25*pi*n)+j*sin(0.125*pi*n)
f1=fft(x,16)
stem(n,abs(f1));
axis([016010])
|x8(k)|'
x8(n)的16点fft
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验二 应用FFT对信号进行频谱分析 实验 应用 FFT 信号 进行 频谱 分析