《离散富立叶变换(DFT)和频谱分析.doc》由会员分享，可在线阅读，更多相关《离散富立叶变换(DFT)和频谱分析.doc（7页珍藏版）》请在金锄头文库上搜索。
1、 实 验 报 告姓 名曾霞学 号2200551070140班 级2006050702成 绩实验名称实验三 离散富立叶变换（DFT）和频谱分析实验日期2008/12/1一实验内容信号函数h(t)= e-a*t*tsin(2ft)编程计算：1. 画出信号函数h(t)的图像。2. 编程计算的离散富立叶变换（DFT）。3. 画出的IDFT图像。4. 计算振幅谱，相位谱以及检测信号主频二算画出信号函数法分析和描述1. 周期序列的离散富立叶变换周期序列的离散富立叶变换公式如下： X(k)=x(n)Wnk, 0kN-1dft运算的方法和步骤如下:初始化X(n=0) 第一重循环：对于k的一个取值ko，计算X(
2、k)=x(n)Wnko,其中旋转因子W= e-j*2*pi/N，循环控制变量n从0变化到N+1第二重循环：对于k=0,1,N-1, 依次计算出X(k),循环变量k从0变化到N-12. 离散富立叶反变换计算x(n)=X(k)W-nk, 0nN-1编程实现计算IDFT计算时，只需要把输入序列换成，输出序列换成，并把旋转因子指数取负号，即，最后将输出序列乘以，计算公式如下：x=1/N*X*H-1=1/N*X*W-G3. 振幅谱，相位谱振幅谱：X(k)=Xr2(k)+ Xi2(k)相位谱: (k)=arctan(Xi(k)/Xr(k)先在振幅谱中搜索振幅最大值，检测到的振幅最大值所对应的频率即为信号主
3、频。for i=1:length(h)-n0 h1(i+n0)=h(i)endsubplot(3,2,5);stem(t2,h1,b.); xlabel (h1(n)的图形)grid on;%h(n)时移n0n3=fix(length(h1)/2)-1;t5=-n3*dt:dt:n3*dt;rhh1=corrl(h,h1);subplot(3,2,6);stem(t4,rhh1,b.); xlabel (h(n)与h1(n)相关函数的图形)grid on%互相关函数rhh(m)计算卷积的 Conv函数：function y=conv(x,h)nx=length(x);nh=length(h);
4、for n1=1:nx+nh-1; sum=0; for m1=1:nh t=n1-m1 if(t0&t0&t=nx) sum=sum+x(m2)*h(m2-n2); end end r(n2+w)=sum;end四实验结果分析（含执行结果验证、输出显示信息和图形等）五结论通过这次试验，更加深了我对卷积和自相关的理解，更加熟悉了matlab的使用，并理解到函数调用的使用方法。通过实验，理解到相关与卷积的运算是类似的，基本运算过程都包含平移，相乘，求和。二者的差异是相关运算不对y(n)反褶，卷积运算则需对y(n)反褶。因此若要利用卷积运算取相关函数，应在卷积前先对y(n)进行一次反褶，其卷积结果既为相关函数，同样也可利用相关计算实现卷积。当x(n) 或y(n)是实偶函数时。卷积与相关运算就完全相同，因为实偶函数反褶与否完全相同。报告提交日期20081118信号分析实验报告 姓名：曾霞 学号：2200551070140 班级：06级软件工程2班
