实验技术与信号处理作业.docVIP

2011-2012学年第二学期研究生 课程考核 （读书报告、研究报告） 考核 科目：实验技术与信号处理 学生所在院（系）：机电工程学院 学生所在学科：机械工程 姓 名： 学 号： 题 目:数字信号处理的步骤与注意事项，并编写1024个采样点的FFT C语言程序 习题:试说明数字信号处理的步骤与注意事项，并用C语言编写1024 个采样点的FFT程序。FFT是什么？大概有哪些方法，各有什么优 点，自己选一种方法编程(方法，背景发展) 数字信号处理步骤： 数字信号处理的步骤主要有预处理、采样与A/D转换和计算机处理三人步骤。预处理 包括解调、滤波、电平转换。采样与A/D转换包括采样、保持、量化。计算机处理包括截 断、加权、数字滤波和计算机信号分析四个部分。 注意事项： 1、 采样间隔的选择，采样间隔太小(采样频率高)，则对定长的时间记录來说其数字 序列就很长，计算工作景迅速增大；如果数字序列长度较短，则只能处理较短的时间历程, 可能会产生较大的谋差。若采样间隔过大(采样频率低)，则可能丢失有川的信息。 2、 要尽最减小混叠现彖的产生。采样频率人需满足：fs 2fh ,几为信号的最高频 率。在实际工作中,常选人=(3?4)办。 3、 由于计算机只能进行有限长序列的运算，所以必须从采样厉信号的时间序列截取有 限长的一?段来计算，其余部分视为零而不了考虑。 4、 截断后会产生泄漏，为了减小或抑制泄漏，需要用窗函数对时域信号进行加权处理, 以改善时域截断处的不连续状况。所选的窗函数应尽最使其频谱的主瓣宽度窄些、旁瓣幅度 小些。窄的主瓣可以提高频率分辨能力；小的旁瓣可以减小泄漏。 快速傅立叶变换(FFT)介绍: 一、基本原理 在信号处理中，DFT的计算具有举足轻重的地位，有限长序列可以通过离散傅立叶变 换(DFT)将其频域离散化，信号的相关、滤波、谱估计等等都可通过DFT來实现。对于N 点序列x (n),其离散傅立叶变换对为： .V-1 X⑹=工X(/?W (0 W k W N小 I .V-l x()=万工X伙)w /V 厶_“ 式中W二LW 由上式，求一个N点的DFT要N?次复数乘法和N (N-1)次负数加法。 当N很大时，其计算量是相当大。如果将(1)、(2)两式写成矩阵形式，如门 x(0) ■ w° … ? x(0)■ 41) ■ ■ ■ — ■ w,xl ? ■ ■ w2xl ■ ■ ■ …W(NT)xl ■ ■ x(l) ■ ? ? ⑶ 兀(N —1) yylMN-l) W2x(N_l) ??? W(NT)x(N-l) X(N —1) x(0) _ w。 IV0 … w° ■ x(0)■ XI) ■ ■ ■ 1 ? ■ ■ IV,xI ? ■ ■ ? ■ ? ... ivGV_,)x, ? ■ ■ x(l) ? ■ ■ ⑷ x(N-V)_ _w° 闪|心-|) ”2心-1) ...g) X(N — 1) 经过观察，在DFT运算屮，包含人量重复运算。观察W矩阵，虽然其屮有N?个元素，但 由于W的周期性，其屮只有N个独立的值，即W°,W〔…,W^,且在这N个值中有一些 W的収值很简单，简言之，W的収值有如下特点： 1、 W°=l , wy/2= -1 2、 = (周期性) 3、 W 2 二一W* (对称性) 利用W因子的周期性和对称性，便可对W矩阵进行简化，根据这种原理，库利(J.W.Cooley) 和图基(J.W.Tukey)在1965年提出了高效的快速算法。 二、发展背景 数字信号处理技术的发展能有效地促进各个工程技术领域的技术改造和学科发展，应用 领域也更加广泛，深入。在数字信号处理屮，离散傅立叶变换(Discrete Fourier Transform, DFT)是常用的变换方法，我们知道频域分析常常比时域分析更优越，不仅简单，且易于分 析复杂信号。但用较精确的数字方法，即DFT进行频谱分析，是不切实际的，因为其计算 量过于庞犬。直到1965年出现FFT。 库利(J.W.Cooley)和图基(J.W.Tukey)于1965年在《计算数学》杂志上发表了 “机 器计算傅立叶级数的一种算法”的文章，这是一篇关于计算DFT的一种快速有效的计算方 法的文章。其他快速算法也如雨后春笋般迅速被发现，这些算法统称为快速傅立叶变换 (FFT)o 1967年至1968年间FFT的数字换件就制成了。至此DFT的运算大为简化。因而各 个科学领域广泛地采用了 FFT技术，它人大推动了近30年来信号处理技术的发展，成为数 字信号处理应用领域强有力的工具，广泛应用于雷达、声纳、通信、地质勘探、图像处理、 生物医学等领域中，在以后的儿十年屮，FFT算法有了进-?步的发展。 三、各种FFT算法比较 IzICooly与Tukey捉出基2 FFT算法以后，新的算