06常用数字信号分析算法应用实验.docVIP

常用数字信号分析算法应用实验 实验目的 了解如何用计算机数字计算来实现信号周期、均值、均方值、方差等参数的识别原理。掌握用数学计算和计算机显示代替复杂的测量电路和机械结构的方法。了解信号频谱和相关系数的直接计算方法，加深对FFT等快速算法优点的理解。1.周期信号的周期/频率测量 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差、或过零点的时间差。 图1、信号周期的识别 工程实际中常用信号过零点检测法识别周期信号周期，其计算机编程实现很简单，用Signal VBScript实现的程序代码为：n=0 For K = 3 To 512 If data(k-1)0 And data(k-2)0? And data(k)=0 And data(k+1)0 Then ???? Rem 用线形插值做过零位置检测 ???? If data(k)=0 Then ?????? ti(n)=K*dt ???? ?? n=n+1 ???? Else ?????? x1=k-1 ?????? y1=data(k-1) x2=k ?????? y2=data(k) a=(y2-y1)/(x2-x1) b=a*x1-y1 ti(n)=b/a n=n+1 ??? End If End If Next T=(ti(1)-ti(0))*dt f=1.0/T 2 信号的时域波形参数测量 (1) 均值E[x(t)]计算机编程实现很简单，用Signal VBScript实现的程序代码为： U=0 For K = 0 To N U=U+data(k) Next U=U/N (2) 均方值 信号x(t)的均方值E[x2(t)]，或称为平均功率，其表达式为： 　 (2) 均方值表达了信号的强度，其正平方根值，又称为有效值，也是信号的平均能量的一种表达。在工程信号测量中一般仪器的表头示值就是信号的有效值。有效值计算机编程实现很简单，用Signal VBScript实现的程序代码为 E2=0 For K = 0 To N E2=E2+data(k)*data(k) Next RMS=sqr(E2/N) 3. 信号的相关函数估计 在离散情况下，信号x(n)和y(n)的相关函数定义为： 在不考虑计算效率的情况下编程计算也很容易： For i = 0 To N r(i)=0 For j = 0 To N r(i)=r(i)+x(j)*y(j+i) Next Next For i = 0 To N r(i)=r(i)/N Next 4. 信号频率成分直接估计算法 离散傅里叶变换的计算公式为： (3) 式中N为信号采样长度，Δt为采样间隔，T=N*Δt为采样窗口宽度。如不考虑计算效率，用计算机编程计算信号中的频率成分是很容易实现的。例如，下面是用Signal VBScript编程计算频率点f的幅值与相位的程序代码.。 f=？ //计算的频率点 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=2*sqr(XR*XR+XI*XI)/T Q=atn(XI/XR) 实验内容 自己用Signal VBScript编写数字信号处理小程序，插入DRVI中，识别测量信号周期、均值、均方值。建立实验环境观察。 图1 信号频率过零检测算法实验环境 点击该实验中的“信号频谱直接算法实验”，建立实验环境。与DRVI中的FFT芯片计算效率对比，理解FFT的优点。 图2 信号频谱直接算法参考实验环境 点击该实验附录中的“信号自相关系数直接算法实验”，建立实验环境。与DRVI中的自相关芯片计算效率对比，理解用FFT计算相关系数的优点。 图3 信号自相关系数直接算法实验环境 实验报告要求 简述实验目的和原理，画出上述实验的虚拟仪器装配图。 拷贝实验系统运行界面，插入到Word格式的实验报告中，用Winzip压缩后通过Email上交实验报告。 思考题 对上述实验功能进行扩充，编程计算出信号的强度，也就是信号的有效值。 图4 编程计算信号的强度 深圳市德普施科技有限公司 SHENZHEN DEPUSH CO., LTD. 工程