- 1、本文档共6页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
序列的傅里叶变换实验报告.doc
实验一 序列的傅立叶变换
一、实验目的
1.进一步加深理解DFS,DFT算法的原理;
2.研究补零问题;
3.快速傅立叶变换(FFT)的应用。
二、实验步骤
复习DFS和DFT的定义,性质和应用;
熟悉MATLAB语言的命令窗口、编程窗口和图形窗口的使用;
利用提供的程序例子编写实验用程序;
按实验内容上机实验,并进行实验结果分析;
写出完整的实验报告,并将程序附在后面。
三、实验内容
周期方波序列的频谱试画出下面四种情况下的的幅度频谱,
已知
并分析补零后,对信号频谱的影响。
有限长序列x(n)的DFT
取x(n)(n=0:10)时,画出x(n)的频谱X(k) 的幅度;
将(1)中的x(n)以补零的方式,使x(n)加长到(n:0~100)时,画出x(n)的频谱X(k) 的幅度;
取x(n)(n:0~100)时,画出x(n)的频谱X(k) 的幅度。利用FFT进行谱分析
3.已知:模拟信号
以t=0.01n(n=0:N-1)进行采样,求N点DFT的幅值谱。
请分别画出N=45; N=50;N=55;N=60时的幅值曲线。
4. 自己编写基2 DIT-FFT的FFT函数,并用编写MATLAB程序,利用DFT计算所给序列的线性卷积;在程序中利用自己编写的FFT函数。
已知,
求
四、图
1、(1)
(2)
(3)
(4)
2、
3、
4、
附录:程序
%Example1
L=5;N=20;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
Xk=dfs(xn,N);
magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
k=[-N/2:N/2];
figure(1)
subplot(2,1,1);
stem(n,xn);xlabel(n);ylabel(xtide(n));
title(DFS of SQ.wave:L=16,N=64);
subplot(2,1,2);
stem(k,magXk);
axis([-N/2,N/2,0,16]);
xlabel(k);ylabel(Xtide(k));
%Example2
M=100;
N=100;
n=1:M;
xn=cos(0.48*pi*n)+cos(0.52*pi*n);
n1=[0:1:N-1];y1=[xn(1:1:M),zeros(1,N-M)];
figure(1)
subplot(2,1,1);
stem(n1,y1);xlabel(n);
title(signal x(n) ,0=n=100);
axis([0,N,-2.5,2.5]);
Y1=fft(y1);
magY1=abs(Y1(1:1:N/2+1));
k1=0:1:N/2;
w1=2*pi/N*k1;
subplot(2,1,2);
title(Samples of DTFT Magnitude);
stem(w1/pi,magY1);
axis([0,1,0,10]);
xlabel(frequency in pi units);
%example 3
figure(1)
subplot(2,2,1)
N=45;n=0:N-1;t=0.01*n;
q=n*2*pi/N;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
y=fft(x,N);
plot(q,abs(y))
stem(q,abs(y))
title(FFT N=45)
%
subplot(2,2,2)
N=50;n=0:N-1;t=0.01*n;
q=n*2*pi/N;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
y=fft(x,N);
plot(q,abs(y))
title(FFT N=50)
%
subplot(2,2,3)
N=55;n=0:N-1;t=0.01*n;
q=n*2*pi/N;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
y=fft(x,N);
plot(q,abs(y))
title(FFT N=55)
%
subplot(2,2,4)
N=16;n=0:N-1;t=0.01*n;
q=n*2*pi/N;
x=2*sin(4*pi*t)+5*cos(8*pi*t);
y=fft(x,N);
plot(q,abs(y))
title(FFT N=16)
function[Xk]=dfs(xn,N)
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n*k;
WNnk=WN.^nk;
Xk=xn*WNnk;
您可能关注的文档
- 实测与分析LED路灯照度.pdf
- 实测竣工测绘合同范本.doc
- 实验一太阳辐射与日照时数的观测.doc
- 实验三十四反射时的测定与反射弧的分析.doc
- 实验六Linux设备驱动中的阻塞与非阻塞IO.doc
- 实验十一多普勒效应的验证.doc
- 实验四傅立叶光学的空间频谱与空间滤波实验.doc
- 实验室天平的正确称量.ppt
- 实验室常用蒸发浓缩方法.doc
- 实验室生物安全工作及危险防范措施.doc
- 2024年江西省寻乌县九上数学开学复习检测模拟试题【含答案】.doc
- 2024年江西省省宜春市袁州区数学九上开学学业水平测试模拟试题【含答案】.doc
- 《GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语》.pdf
- 中国国家标准 GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语.pdf
- GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- 《GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构》.pdf
- 中国国家标准 GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 中国国家标准 GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 《GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南》.pdf
文档评论(0)