Matlab课程设计5.docVIP

设计题目：信号的频域分析及MATLAB实现 班级： 姓名： 同组 2012年 6 月 4 日 设计目的 通过该设计，理解傅里叶变换的定义及含义，掌握对信号进行频域分析的方法。 设计内容 1.信号的离散傅里叶变换。 有限长序列的离散傅里叶变换公式为： MATLAB函数：fft功能是实现快速傅里叶变换，fft函数的格式为：y=fft(x),返回向量x的不连续fourier变换。 若x(n)=cos()是一个N=12的有限序列，利用MATLAB计算x(n)它的离散傅里叶变换X(k)并画出图形，然后再对X(k)进行离散傅里叶反变换，求出x(n)并画出其波形。 2.频率分辨率与DFT参数的选择 在DFT问题中，频率分辨率是指在频率轴上所能得到的最小频率间隔,即最小频率间隔反比于数据的长度N。若在x(n)中有两个频率分别和的信号，对x(n)用矩形窗截断时，要分辨出这两个频率，N必须满足 通过下面实验，验证上面的结论： 设一序列中含有两种频率成分，，采样频率取为，表示为 根据上面的结论，要区分出这两种频率成分，必须满足N400. 取x（n）（）时，计算x（n）的DFT X(k)； 取x（n）（）时，计算x（n）的DFT X(k)。 三总体设计方案 按题目要求进行设计，写出设计报告，给出源程序。 画出x(n)=cos()和它的离散傅里叶变换X（k）的图形，然后再对X（k）进行离散傅里叶反变换求出x（n）并画出X（k）和x（n）的波形，观察两次实验结果是否吻合。 ，取x（n）（）时，画出x(n)和x(n)的DFT X(k)的图形；并取x(n)（）时，画出x(n)和x(n)的DFT X(k)的图形，比较两次实验是否能区分出两种频率成分。 主要参数 f1=2Hz f2=2.05Hz 采样频取为 fs=10Hz 五、源程序 1，画出想x(n)=cos(nπ/6)和它的离散傅里叶变换X（k）的图形，然后再对X（k）进行离散傅里叶反变换求出x(n)并画出X（k）和x(n)的波形，观察两次实验结果是否吻合。 （1） N=12,n=1:N-1; N = 12 Xn=cos(pi*n/6); K=0:N-1; Wn=exp(-j*2*pi*n*K/N); Xk=Xn*Wn; figure,subplot(1,2,1),stem(n,Xn) subplot(1,2,2),stem(K,abs(Xk)) y=fft(Xk); figure,stem(K,abs(y)) （2） Xn=ifft(Xk); figure,subplot(1,2,1),stem(n,Xn) subplot(1,2,2),stem(K,abs(Xk)) 2.，取x（n）（）时，画出x(n)和x(n)的DFT X(k)的图形；并取x(n)（）时，画出x(n)和x(n)的DFT X(k)的图形，比较两次实验是否能区分出两种频率成分。 f1=2;f2=2.05;fs=10; N=500; n=0:128;k=n;nk=n*k; Xn=sin(2*pi*f1*n/fs)+sin(2*pi*f2*n/fs); WN=exp(-j*2*pi/N); Wnk=WN.^nk; Xk=Xn*Wnk; subplot(1,2,1),plot(n,Xn) subplot(1,2,2),plot(k,Xk) Xn=sin(2*pi*f1*n/fs)+sin(2*pi*f2*n/fs); N=length(Xn); n=0:N-1; k=n;nk=n*k; Wnk=WN.^nk; Xk=Xn*Wnk; subplot(1,2,1),plot(n,Xn) subplot(1,2,2),plot(k,Xk) 六、仿真结果 讨论、心得 短短的课程设计已经结束了，虽然时间很短暂，但是学到的东西却是只有课堂教学所无法比的。这样的课程设计，让我们每个人都能亲自动手去学习知识，去运用知识，通过实践来加深对课堂知识的理解。而且可以弥补课堂教