原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
北邮dsp实验二实验报告
实验二：数字信号的 FFT 分析题目1假设信号 x(n) 由下述信号组成：这个信号有两根主谱线 0.3pi 和 0.302pi 靠的非常近，而另一根谱线 0.45pi 的幅度很小，请选择合适的长度 N 和窗函数，用 DFT 分析其频谱，得到清楚的三根谱线。步骤：1.编写离散傅里叶变换DFT函数：function [Xk] = dft(xn,N)% Computes Discrete Fourier Transform Coefficients% [Xk] = dft(xn,N)% Xk = DFT coeff. array over 0 = k = N-1% xn = input signal% N = length of DFTn = [0:1:N-1]; % row vector for nk = [0:1:N-1]; % row vecor for kWN = exp(-j*2*pi/N); % Wn factornk = n*k; % creates a N by N matrix of nk valuesWNnk = WN .^ nk; % DFT matrixXk = xn * WNnk; % row vector for DFT coefficients2.代码实现：n=0:1:999;x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-0.25*pi);stem(n,x);title(signal x(n), 0=n=999);xlabel(n);X=dft(x,1000);% 计算1000点DFTmagX=abs(X(1:1:501));% 镜像对称，只画出一半k=0:1:500;w=2*pi*k/1000;stem(w/pi,magX);title(DTFT Magnitude);xlabel(frequency in pi units);axis([0.29,0.31,0,500]);xlabel(frequency between 0.29pi and 0.31pi);axis([0.44,0.46,0,0.5]);xlabel(frequency between 0.44pi and 0.46pi);3.图片：4.分析：x(n)由3个正弦函数叠加而成，周期分别是40, 20, 1000。这里序列长度选择n=1000，为周期的最小公倍数。从频谱分析中可以看出，此时的谱线中有3条清晰的谱线，所在位置分别是0.3pi, 0.302pi和0.45pi。这说明当取样点数为函数周期的整倍数时，其频谱中科出现3条清晰的谱线。题目2已知信号这里，N=25，Q= 0.9+j0.3。可以推导出，首先根据这个式子计算X(k) 的理论值，然后计算输入序列x(n) 的32个值，再利用基2时间抽选的FFT算法，计算x(n) 的DFT X(k)，与X(k) 的理论值比较（要求计算结果最少6位有效数字）。步骤：1.计算X(k) 的理论值写入代码：Q= 0.9+j*0.3;N=25;n=0:1:N-1;k=0:1:N-1;WN=exp(-j*2*pi/N);X=(1-Q.^N)./(1-Q*WN.^k);X=vpa(X,6)得到结果：X =[ 1.83992 + 2.88851*i, 10.0793 + 7.6341*i, 0.751948 - 5.86055*i, 0.231478 - 2.56282*i, 0.28018 - 1.62973*i, 0.333453 - 1.18569*i, 0.374736 - 0.920169*i, 0.407012 - 0.73897*i, 0.433275 - 0.603825*i, 0.455542 - 0.496182*i, 0.475155 - 0.405862*i, 0.493047 - 0.326701*i, 0.50991 - 0.254627*i, 0.526307 - 0.18669*i, 0.542748 - 0.120528*i, 0.559746 - 0.0540151*i, 0.577892 + 0.0150219*i, 0.597934 + 0.0890957*i, 0.620917 + 0.171456*i, 0.64843 + 0.266736*i, 0.683099 + 0.382148*i, 0.729676 + 0.529983*i, 0.797838 + 0.733482*i, 0.91
文档评论（0）