第六章试验数据处理技术.pptVIP

6.5 DFT与FFT 　 1、离散傅立叶变换 离散傅里叶变换（Discrete Fourier Transform）一词是为适应计算机作傅里叶变换运算而引出的一个专用名词。 x(t) 截断、周期延拓 xT(t) 周期信号xT(t)的傅里叶变换： 对周期信号xT(t)采样，得离散序列xT(n),将积分转为集合： 展开，得连续傅立叶变换计算公式： 用计算机编程很容易计算出指定频率点值： f=？ //计算的频率点 Fs=? N=1024 dt=1.0/Fs pi=3.1415926 XR=0 XI=0 For n=0 To N-1 XR=XR+x(n)*cos(2*pi*f*n*dt)*dt XI=XI+x(n)*sin(2*pi*f*n*dt)*dt Next A=sqr(XR*XR+XI*XI) Q=atn(XI/XR) VBScript 样例 6.5 DFT与FFT 　 6.5 DFT与FFT 　 连续傅立叶变换编程计算实验： 　 采样信号频谱是一个连续频谱，不可能计算出所有频率点值，设频率取样间隔为： Δf = fs / N 频率取样点为｛0,Δf,2Δf,3Δf,....｝，有： 6.5 DFT与FFT 　 该公式就是离散傅立叶计算公式(DFT) 6.5 DFT与FFT 　 2、快速傅立叶变换 快速傅立叶变换(FFT)是离散傅立叶变换的一种有效的算法，通过选择和重新排列中间结果，减小运算量。 展开各点的DFT计算公式： XR(1)=x(0).cos(2pi*0*1/N)+x(1).cos(2pi*1*1/N)+x(2).cos(2pi*2*1/N)….. XR(2)=x(0).cos(2pi*0*2/N)+x(1).cos(2pi*1*2/N)+x(2).cos(2pi*2*2 /N)….. 6.5 DFT与FFT 　 有大量重复的cos、sin计算，FFT的作用就是用技巧减少cos、sin项重复计算。 当采样点数为1024点,DFT要求一百万次以上计算量，而FFT则只要求一万次。 6.6 栅栏效应与窗函数 　 1、栅栏效应 为提高效率,通常采用FFT算法计算信号频谱，设数据点数为N，采样频率为Fs。则计算得到的离散频率点为: Xs(Fi) ， Fi = i *Fs / N , i = 0，1，2，.....，N/2 X(f) f 0 Δf Δ 如果信号中的频率分量与频率取样点不重合，则只能按四舍五入的原则，取相邻的频率取样点谱线值代替。 6.5 DFT与FFT 　 栅栏效应误差实验： 　 6.6 栅栏效应与窗函数 　 2 能量泄漏与栅栏效应的关系 频谱的离散取样造成了栅栏效应，谱峰越尖锐，产生误差的可能性就越大。 例如，余弦信号的频谱为线谱。当信号频率与频谱离散取样点不等时，栅栏效应的误差为无穷大。 * 华中科技大学机械学院 * 华中科技大学机械学院 第六章 数字信号处理技术 本章学习要求： 1.了解信号模数转换和数模转换原理　　 2.掌握信号采样定理，能正确选择采样频率 3.了解数字信号处理中信号截断、能量泄露、栅栏效应等现象　 4.掌握常用的数字信号处理方法 6.1 数字信号处理概述 　 1、数字信号处理的主要研究内容 数字信号处理主要研究用数字序列来表示测试信号，并用数学公式和运算来对这些数字序列进行处理。内容包括数字波形分析、幅值分析、频谱分析和数字滤波。 0 A t X(0) X(1) X(2) X(3) X(4) 2、测试信号数字化处理的基本步骤 物理信号 对象 传感器 电信号 放大调制 电信号 A/D 转换 数字信号 计算机 显示 D/A转换 电信号 控制 物理信号 3、数字信号处理的优势 1)用数学计算和计算机显示代替复杂的电路 和机械结构 2)计算机软硬件技术发展的有力推动 a)多种多样的工业用计算机。 ?????????????????????????????????????????????????????????????????????????????????????? b)灵活、方便的计算机虚拟仪器开发系统 案例：铁路机车FSK信号检测与分析 京广线计划提速到200公里/小时 合作任务：机车状态信号识别(频率解调) 虚拟仪器设计方案 6.2 模数(A/D)和数模(D/A) 　 采样――利用采样脉冲序列，从信号中抽取一系列 离散值，使之成为采样信号x(nTs)的过程． 编码――将经过量化的值变为二进制数字的过程。 量化――把采样信号经过舍入变为只有有限个有