信号实验报告三.docxVIP

下载本文档

19
0
约2.68千字
约 6页
2019-09-14 发布于广东
举报
版权申诉

信号实验报告三.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

1、比较DFT和FFT的运算时间。自行选择进行计算的序列(或产生随机序列)。利用计吋函数tic、toe,计算点数N=64、128、256、 tic Nl=64 xl=[ones(l,64)J Xkl=dft(xl,Nl) toe %Elapscd time is 0.016000 seconds. tic N2=512 x2=fones( 1,512)] Xk2=dft(x2,N2) toe %Elapsed time is 0.687000 seconds. tic N3=1024 x3=[ones(l,1024)] Xk3=dft(x3,N3) toe %Ebpsed time is 3.016000 seconds. tic Nll=64 xll=[ones(l,64)] Xkll=fTt(xll,Nll) toe %Elapsed time is 0.000000 seconds. tic N22=512 x22=Lones(l,512)] Xk22=fft(x22,N22) toe % Elapsed time is 0.015000 seconds. tic N33=1024 x33=Lones(l,1024)J Xk33=fft(x33,N33) toe % Elapsed time is 0.047000 seconds. 2、利用FFT计算信号功率谱。自行设计信号，利用周期图法实现信号功率谱并慚图。 figure subplot(2J J) N=20;n=0:N-l ;t=0.01 *n; q=n*2*pi/N; x=3*cos(8*pi*t)+6*cos(20*pi*t); y=fft(x,N); plot(q/pi,abs(y));title(FFT N=45) FFT N=50 FFT N=50 3、利用FFT求信号频谱及分析采样频率对频谱的影响。设模拟信号 x（/） = 3*cos（8加）+ 6cos（20加）,以 f = 0.0In 0 /? - 1 进行取样，求N点DFT的幅值谱（N分别取50, 256）,利用FFT实现。分析信号频谱所对应频率轴的数字频率和频率Z间的关系。在截取长度不变的条件下改变采样频率，观察信号频谱的变化，分析产生变化的原因。 N=50 n=0:49 t=0.01*n %采样间隔 0.01s q 二 n*2*pi/N 1; x=3*cos(8*pi*t)+6*cos(20*pi*t) y=fft(x,N) subplot(2,2,l) plot(q/pi,abs(y)); titleCFFT N=50r) %N=50采样 N=256 n=0:255 t=0.01*n %采样间隔 0.01s q=n*2*pi/N; x=3*cos(8*pi*t)+6*cos(20*pi*t) y=fft(x,N)subplot(2,2,l)plot(q/pi,abs(y));titleCFFT N=256‘) %N=256采样N=50n=0:49t=0.1*n %采样间隔0.1q=n*2*pi/Nl;x=3*cos(8*pi*t)+6*cos(20*pi*t)y=fft(x,N)subplot(2,2,l)plot(q/pi,abs(y));title(!FFT N=50f) %N=50采样N=256n=0:255t=0.1*n %采样间隔0.1s y=fft(x,N) subplot(2,2,l) plot(q/pi,abs(y)); titleCFFT N=256‘) %N=256采样 N=50 n=0:49 t=0.1*n %采样间隔0.1 q=n*2*pi/Nl; x=3*cos(8*pi*t)+6*cos(20*pi*t) y=fft(x,N) subplot(2,2,l) plot(q/pi,abs(y)); title(!FFT N=50f) %N=50采样 N=256 n=0:255 t=0.1*n %采样间隔0.1s q=n*2*pi/N; x=3*cos(8*pi*t)+6*cos(20*pi*t) y=fft(x,N) subplot(2,2,l) plot(q/pi,abs(y)); title(fFFT N=256f) %N=256采样 4、利用FFT计算快速卷积。利用FFT计算2个序列的快速卷积, 并与直接计算线性卷积的结杲及时间进行比较。 %线性卷积 tic n=0:10 xl=cos(2*n) x2=2.An y 1 =conv_m(x 1 ,n,x2,n) toe %yl= 1.0000 1.5839 2.5141 5.9883 11.8311 22.8231 46.4901 93.1169 185.2761 371.2126 74