原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
.
实验二 用 FFT做谱分析
一、实验目的
1. 进一步加深 DFT 算法原理和基本性质的理解（因为 FFT 只是 DFT 的一种快
速算法，所以 FFT 的运算结果必然满足 DFT 的基本性质）。
2. 熟悉 FFT 算法原理和 FFT 子程序的应用。
3. 学习用 FFT 对连续信号和时域离散信号进行谱分析的方法
二、实验原理
如果给出的是连续信号 xa(t) ，则首先要根据其最高频率确定抽样频率 f s 以及由频
率分辨率选择抽样点数 N ，然后对其进行软件抽样 ( 即计算 x(n)=x a(nT)，0≤n≤
N-1) ，产生对应序列 x(n) 。再利用 MATLAB 所提供的库函数 fft(n,x) 进行 FFT 计
算
三、实验内容
① 实验信号 ：
x1(n) = R4(n)
n 1, 0 n 3
2
x (n) = 8 n, 4 n 7
0 , 其他 n
4 n, 0 n 3
3
x (n) = n 3, 4 n 7
0 , 其他 n
4
x (n) = cos( πn /4)
5
x (n) = sin( πn /8)
6
x (t) = cos8 πt + cos16 πt + cos20 π
FFT 变换区间及 x6 (t)抽样频率 fs
1 2 3 4 5
x (n) , x (n) , x (n) , x (n) , x (n)：N = 8 , 16
6 s
x (t) ：f = 64(Hz) , N = 16 , 32 , 64
② MATLAB 程序代码
N1=8;
N2=16;
x1=ones(1,4);
x2=[1:4,4:-1:1];
x3=[4:-1:1,1:4];
n=0:1:16;
x4=cos(pi*n/4);
..
.
x5=sin(pi*n/8);
X11=fft(x1,N1);
X11=abs(X11);
X21=fft(x2,N1);
X21=abs(X21);
X31=fft(x3,N1);
X31=abs(X31);
X41=fft(x4,N1);
X41=abs(X41);
X51=fft(x5,N1);
X51=abs(X51);
X12=fft(x1,N2);
X12=abs(X12);
X22=fft(x2,N2);
X22=abs(X22);
X32=fft(x3,N2);
X32=abs(X32);
X42=fft(x4,N2);
X42=abs(X42);
X52=ff
文档评论（0）