《数字信号处理》课程_实验报告.docVIP

《数字信号处理》课程 实验报告 西安交通大学 电子与信息工程学院 实验一 快速傅立叶变换 一 实验目的 在理论学习的基础上，通过本实验加深对快速傅立叶变换的理解； 熟悉并掌握按时间抽取FFT算法的程序； 了解应用FFT进行信号频谱分析过程中可能出现的问题，例如混淆、泄漏、栅栏效应等，以便在实际中正确应用FFT。 二 实验内容 仔细分析教材第六章‘时间抽取法FFT ’的算法结构，编制出相应的用FFT进行信号分析的C语言（或MATLAB 语言）程序； 用FFT程序计算有限长度正弦信号 分别在以下情况下所得的DFT结果并进行分析和讨论： 信号频率f＝50Hz，采样点数N=32,采样间隔T=0.000625s 信号频率f＝50Hz，采样点数N=32,采样间隔T=0.005s 信号频率f＝50Hz，采样点数N=32,采样间隔T=0.0046875s 信号频率f＝50Hz，采样点数N=32,采样间隔T=0.004s 信号频率f＝50Hz，采样点数N=64,采样间隔T=0.000625s 信号频率f＝250Hz，采样点数N=32,采样间隔T=0.005s 将c） 信号后补32个0，做64点FFT 三 实验步骤 1 根据要求，使用MATLAB软件编写出相应的用FFT进行信号分析的程序。 程序如下： clear; clc; f=input(信号频率f(HZ)\n); N=input(采样点数N\n); T=input(采样间隔T(s)\n); %输入N,T,f的具体数值， n=0:N-1; t=n/f; for j=0:1:N-1 x(j+1)=sin(2*pi*f*j*T); end BL=input(是否需要补0？1(是)/0(否)\n); %判断是否需要补0 if BL==1 gs=input(补0个数\n); %输入补0个数 for j=N:1:N+gs-1 %给采样点后补0 x(j+1)=0; end N=N+gs; %修正采样点数 end M=log2(N); %判断分解级数 for t=1:1:N %码位倒置 s=dec2bin(t-1,M); %把十进制整数t-1转换成2进制形式s, M表示转换成2进制后的数的位数 s=fliplr(s); % 实现矩阵的左右翻转 s=bin2dec(s); %把2进制数s转换成十进制形式 y(s+1)=x(t); end x=y; %按时间抽取的FFT蝶形运算 for L=1:1:M %将FFT分为M级进行 for J=0:1:(2^(L-1)-1) for k=(J+1):2^L:N T=y(k)+y(k+2^(L-1))*exp((-i*2*pi*J*2^(M-L))/N); y(k+2^(L-1))=y(k)-y(k+2^(L-1))*exp((-i*2*pi*J*2^(M-L))/N); y(k)=T; end end end x=abs(y); %将所得频谱取模值 y=max(x); %模值归一化 X=x/y; for j=1:1:N %绘制棒状图 stem(j-1,X(j)); hold on end axis([0 N 0 1]); %设置坐标值范围 2、运行程序，完成实验各项要求，并对所得的DFT结果进行分析和讨论。 （a）信号频率f＝50Hz，采样点数N=32,采样间隔T=0.000625s 分析：对于信号，由于采样点数为N=32，采样间隔为T=0.000625s，而信号的频率为f=50HZ，所以刚好是在信号的一个周期中采样了32个点，在FF