LabVIEW编程及虚拟仪器设计(第八讲)

1、LabVIEW编程及虚拟编程及虚拟 仪器设计仪器设计第八讲第八讲：测量信号的分析与处理测量信号的分析与处理测量信号的分析与处理（测量信号的分析与处理（1 1） 因测量环境中常存在各种可能的干扰量，故实际测量中采因测量环境中常存在各种可能的干扰量，故实际测量中采集到的信号往往不仅是所关心的信号本身，而还伴有不希望集到的信号往往不仅是所关心的信号本身，而还伴有不希望的其它成分。对采集信号的分析与处理，即指如何有效地将的其它成分。对采集信号的分析与处理，即指如何有效地将采得的数字化信号中的各种干扰消除掉，或最大程度地将其采得的数字化信号中的各种干扰消除掉，或最大程度地将其削弱，从而得到（提取出）所关

2、心的有用信号。削弱，从而得到（提取出）所关心的有用信号。 测量信号的分析与处理（测量信号的分析与处理（2 2） 经经“函数函数”选板选板 -Express - - “信号分析信号分析”途径，可找途径，可找到到13个用于信号分析与处理的用个用于信号分析与处理的用Express VI。 由此途径可找到的用于信号分析与处理的快速由此途径可找到的用于信号分析与处理的快速VIVI即即Express VIExpress VI，大多由基本函数构成。，大多由基本函数构成。 频谱测量：频谱测量：进行基于进行基于FFT的频谱测量，如信号的平均幅的频谱测量，如信号的平均幅度频谱、功率谱和相位谱。度频谱、功率谱和相位

3、谱。 失真测量：失真测量：对信号进行失真测量，如音频分析、总谐波对信号进行失真测量，如音频分析、总谐波失真（失真（THD)、信号与噪声失真比（、信号与噪声失真比（SINAD)。 单频测量：单频测量：输出混频信号中具有最高幅值的频率成分输出混频信号中具有最高幅值的频率成分的信息（幅值、频率、相位），或在指定频率范围内查找混的信息（幅值、频率、相位），或在指定频率范围内查找混频信号中具有最高幅值的频率成分的信息。频信号中具有最高幅值的频率成分的信息。 双通道谱测量：双通道谱测量：根据根据A通道（作为激励）和通道（作为激励）和B通道通道(作作为响应）的输入信号，测量其对应系统的频率响应。为响应）的输

4、入信号，测量其对应系统的频率响应。 幅值和电平测量幅值和电平测量：测量信号的电压：测量信号的电压 信号的时间与瞬态特性测量：信号的时间与瞬态特性测量：测量信号（通常是脉冲）测量信号（通常是脉冲）的瞬态特性，如频率、周期、占空比。的瞬态特性，如频率、周期、占空比。 曲线拟合：曲线拟合：根据所选的模型类型，计算最能代表输入数据的数根据所选的模型类型，计算最能代表输入数据的数学模型系数。学模型系数。 滤波器：滤波器：通过滤波器和窗对信号进行处理。通过滤波器和窗对信号进行处理。 统计：统计：返回波形中第一个通道信号的平均值、均方根、方差等返回波形中第一个通道信号的平均值、均方根、方差等选中参数。选中参

5、数。 卷积和相关：卷积和相关：在输入信号上进行卷积、反卷积、相关等操作。在输入信号上进行卷积、反卷积、相关等操作。 仿真信号：仿真信号：产生正弦、方波、三角波、锯齿波和噪声等仿真信产生正弦、方波、三角波、锯齿波和噪声等仿真信号。号。 信号掩区和边界测试：信号掩区和边界测试：在信号上进行边界测试、该在信号上进行边界测试、该Express VI根据用户设定的上下边界比较信号。根据用户设定的上下边界比较信号。 创建直方图：创建直方图：创建信号的直方图。创建信号的直方图。 随着随着LabVIEW版本的不断翻新，为更便于用户版本的不断翻新，为更便于用户的使用，的使用，LabVIEW将函数类型丰富为了三种

6、，即：将函数类型丰富为了三种，即：基本函数、波形函数以及基本函数、波形函数以及Express VI。 经经 “函数函数”选板选板 - “信号处理信号处理”途径途径，还可找到更，还可找到更多用于信号分析与处理的功能函数。多用于信号分析与处理的功能函数。 基本平均直流基本平均直流-均方根：按照选定的加窗和平均类均方根：按照选定的加窗和平均类型计算波形的直流分量（型计算波形的直流分量（DC)和均方极值（和均方极值（RMS)平均直流平均直流-均方根：按照设定平均类型和时间长度计算均方根：按照设定平均类型和时间长度计算波形的波形的DC和和RMS，可输出，可输出DC和和RMS数值的波形。数值的波形。周期平

7、均值和均方根：计算信号的周期平均值和均方根周期平均值和均方根：计算信号的周期平均值和均方根瞬态特性测量：瞬态测量包括持续时间（上升或下降时瞬态特性测量：瞬态测量包括持续时间（上升或下降时间），摆动率、前冲及过冲等。间），摆动率、前冲及过冲等。脉冲测量：脉冲测量包括周期、脉宽和占空比等脉冲测量：脉冲测量包括周期、脉宽和占空比等幅值和电平：返回波形的幅值、最高电平和最低电平幅值和电平：返回波形的幅值、最高电平和最低电平提取单频信息：在指定频段提取幅值最大的成分频率信号、提取单频信息：在指定频段提取幅值最大的成分频率信号、可输出抽取出的信号的波形、频率、幅值及相位。可输出抽取出的信号的波形、频率、幅

8、值及相位。提取混合单频信息：返回幅值超过提取混合单频信息：返回幅值超过“阈值阈值”的信号的频的信号的频率、幅值和相位。率、幅值和相位。谐波失真分析：谐波分析谐波失真分析：谐波分析SINAD分析：含噪声信号的谐波分析分析：含噪声信号的谐波分析FFT功率谱：计算信号的功率谱功率谱：计算信号的功率谱FFT功率谱密度：计算信号的功率谱密度功率谱密度：计算信号的功率谱密度FFT频谱（幅度频谱（幅度-相位）：相位）：FFT谱分析（幅度谱分析（幅度-相位形式）相位形式）例题：信号谱分析例题：信号谱分析.VIFFT频谱（实部频谱（实部-虚部）：虚部）：FFT谱分析（实部谱分析（实部-虚部形式）虚部形式）频率响

9、应函数（幅度频率响应函数（幅度-相位）：计算系统的频率响应，相位）：计算系统的频率响应，输入为系统的激励和响应的时域信号，输出是系统输入为系统的激励和响应的时域信号，输出是系统的幅频特性和相频特性。的幅频特性和相频特性。频率响应函数（实部频率响应函数（实部-虚部）：计算系统的频率响应，虚部）：计算系统的频率响应，输入为系统的激励和响应的时域信号，输出是以实输入为系统的激励和响应的时域信号，输出是以实部和虚部的形式反映系统的频率特性。部和虚部的形式反映系统的频率特性。交叉谱（幅度交叉谱（幅度-相位）：互功率谱分析（幅度相位）：互功率谱分析（幅度-相位形式）相位形式）交叉谱（实部交叉谱（实部-虚部

10、）：互功率谱分析（实部虚部）：互功率谱分析（实部-虚部形式）虚部形式）1 1、信号的产生、信号的产生编程编程信号处理信号处理波形产生波形产生例题例题 函数发生器函数发生器.VI例题：带幅值的基本混合单频信号发生器例题：带幅值的基本混合单频信号发生器.VIdelta频率频率是介于相邻单是介于相邻单频频率间的幅值。频频率间的幅值。delta频率必须是频率必须是Fs/#s的倍数。的倍数。默认值为默认值为100。 如起始如起始频率为频率为100 Hz，delta频频率为率为10，单频幅值数组，单频幅值数组将包含三个元素，即将包含三个元素，即100Hz、110Hz和和120Hz。 起始频率起始频率是生成

11、的最低单是生成的最低单频频率。值必须是频频率。值必须是Fs/#s的倍数。的倍数。默认值为默认值为10。单频幅值单频幅值该数组的元素为单频的该数组的元素为单频的幅值。数组的大小确定生成的单幅值。数组的大小确定生成的单频数量。频数量。 幅值幅值是所有单频的缩放标准，即是所有单频的缩放标准，即波形的最大绝对值。波形的最大绝对值。默认值为默认值为-1。 将波形输出至模拟输出通道时，将波形输出至模拟输出通道时，可使用幅值。如硬件可输出的最可使用幅值。如硬件可输出的最大值为大值为5伏，可将幅值设置为伏，可将幅值设置为5。如幅值如幅值 为为0，将不进行缩放。，将不进行缩放。 测量信号的分析与处理（测量信号的

12、分析与处理（4 4） “函数函数”选板选板 - “信号处理信号处理” - - “波形生成波形生成”-“仿真信号仿真信号” 测量信号的分析与处理（测量信号的分析与处理（5 5） 对被测信号进行对被测信号进行FFTFFT变换举例（变换举例（A A） 例如，产生一个周期信号，然后对它做例如，产生一个周期信号，然后对它做FFT，可以得到其，可以得到其幅度和相位频谱。幅度和相位频谱。 实施测量任务时，经数据采集卡采到的是时域波形的抽样实施测量任务时，经数据采集卡采到的是时域波形的抽样信号信号采样信号。采样信号。在时域可分析处理信号的时域特征在时域可分析处理信号的时域特征，如信，如信号随时间的号随时间的变

13、化趋势、幅值大小，极值变化趋势、幅值大小，极值等等。等等。 被测信号也有频率特征，如其中含有不同频率成分等。为被测信号也有频率特征，如其中含有不同频率成分等。为分析信号的频率特征，需要将获得的时域波形的采样信号转换分析信号的频率特征，需要将获得的时域波形的采样信号转换到频域去分析处理。到频域去分析处理。在频域分析信号的基本方法是傅立叶变换在频域分析信号的基本方法是傅立叶变换法，其快速实现方法又称法，其快速实现方法又称FFT，且由它逐渐派生出很多种用于且由它逐渐派生出很多种用于在频域分析处理信号的功能函数。在频域分析处理信号的功能函数。对信号进行频域分析（包括对信号进行频域分析（包括频谱分析和扫

14、频分析），能获得在时域测量中得不到的频率特频谱分析和扫频分析），能获得在时域测量中得不到的频率特征信息，如谐波分量、寄生、交调、噪声边带等。征信息，如谐波分量、寄生、交调、噪声边带等。 “函数函数”选板选板 - “信号处理信号处理” - “波形测量波形测量”-“频谱测量频谱测量” 现象：幅度谱存在有现象：幅度谱存在有“泄漏泄漏”现象。这是因为现象。这是因为采集到的被测信号的样本采集到的被测信号的样本数有限。增加样本数，情数有限。增加样本数，情况会有所改善，但不能完况会有所改善，但不能完全消除。全消除。 对被测信号进行对被测信号进行FFTFFT变换举例（变换举例（B B）NI帮助帮助- 分析与处

15、理信号分析与处理信号-电平测量电平测量测量信号的分析与处理（测量信号的分析与处理（6 6） 给被测信号给被测信号“加窗加窗”（a a） 基于计算机构建的虚拟仪器只能分析基于计算机构建的虚拟仪器只能分析/处理有限长的数字信处理有限长的数字信号，被测信号号，被测信号x(t)以以T（采样时间或采样长度）被截取一段送入（采样时间或采样长度）被截取一段送入计算机，就称为被截断，相当于加了计算机，就称为被截断，相当于加了“矩形窗矩形窗”将信号突将信号突然截断，致使在很宽频率范围内向被测信号中添加了另外的成然截断，致使在很宽频率范围内向被测信号中添加了另外的成分。附加频率成分不属于分。附加频率成分不属于x(

16、t)，被称为，被称为“频谱泄漏频谱泄漏”。频谱泄漏。频谱泄漏会带来如下问题：会带来如下问题： 为减少频谱泄漏，可采用变化相对缓慢的非矩形窗函数对被为减少频谱泄漏，可采用变化相对缓慢的非矩形窗函数对被测信号进行截断；这被称为给被测信号测信号进行截断；这被称为给被测信号“加窗加窗”。 使信号的频域曲线上产生许多使信号的频域曲线上产生许多“皱纹皱纹”、频率分辨率降低；、频率分辨率降低； 如果信号为幅值一大一小、频率很接近的两正弦波信号如果信号为幅值一大一小、频率很接近的两正弦波信号 的合成，因泄漏，在频域，幅值小信号可能被的合成，因泄漏，在频域，幅值小信号可能被“淹没淹没”； 所关注信号的频率所关注

17、信号的频率f0附近的频率特性曲线部分可能过于平附近的频率特性曲线部分可能过于平 缓，致使无法准确确定缓，致使无法准确确定f0的值。的值。测量信号的分析与处理（测量信号的分析与处理（7 7） 给被测信号给被测信号“加窗加窗”（b b） 在实际应用中，如何选择窗函数？在实际应用中，如何选择窗函数？ 要仔细分析被测信号的特征，再结合希望达到的目的，并要仔细分析被测信号的特征，再结合希望达到的目的，并可能要经过反复试验。窗函数有多种，如汉宁窗、海宁窗、余可能要经过反复试验。窗函数有多种，如汉宁窗、海宁窗、余弦窗（多种）、布莱克曼窗，等等，它们各有特点；弦窗（多种）、布莱克曼窗，等等，它们各有特点；若使

18、用不若使用不当，甚至可能带来负面效应当，甚至可能带来负面效应。 使用窗函数的原由很多，例如使用窗函数的原由很多，例如：（：（1）限）限定测量的持续时间；定测量的持续时间；（2）减少频谱泄漏；（）减少频谱泄漏；（3）将）将频率接近、幅值不同的信号成分频率接近、幅值不同的信号成分分离出来分离出来。 举例举例：把频率接近、幅值不同的两个正弦信号相互分离提把频率接近、幅值不同的两个正弦信号相互分离提取出来。经过途径取出来。经过途径 “帮助帮助”“查找范例查找范例”“信号分析和处信号分析和处理理”“FFT“FFT和频率分析和频率分析”Window Comparison.viWindow Comparis

19、on.vi 可找到这个例子。可找到这个例子。测量信号的分析与处理（测量信号的分析与处理（8 8） 给被测信号给被测信号“加窗加窗”（c c） “帮助帮助”“查找范例查找范例”“信号分析和处理信号分析和处理”“FFT和频率分析和频率分析”Window Comparison.vi测量信号的分析与处理（测量信号的分析与处理（9 9） ChirpChirp信号及其频谱（一）信号及其频谱（一） Chirp信号又称线性调频信号，其数学表达式为：信号又称线性调频信号，其数学表达式为： F(t)=sin(a*t2+b*t) 其中，设其中，设t为时间变量；为时间变量； 调节参数调节参数a和和b，可得到幅值恒定、

20、频率连续变化（在一定，可得到幅值恒定、频率连续变化（在一定频率范围内）的一个频率范围内）的一个“正弦正弦”信号。信号。LabVIEW提供有此函数。提供有此函数。 从途径从途径 “函数函数”选板选板 “ “信号处理信号处理” “ “信号生成信号生成” Chirp Chirp 信号信号.vi .vi 可得到这个函数。可得到这个函数。 注意注意：利用已给定的缺省值（默认值）为宜，观察效果好：利用已给定的缺省值（默认值）为宜，观察效果好。举例举例：调用调用Chirp函数，并观察它的时域和频域特性。函数，并观察它的时域和频域特性。 Chirp信信号及其频号及其频谱（二）谱（二）F(t)=sin(a*t2

21、+b*t)2、频域分析、频域分析 例例11-2 求两个波形的频率和相位差：信求两个波形的频率和相位差：信号谱分析号谱分析.VI 例例11-3 从频率接近的信号中分离出幅从频率接近的信号中分离出幅值不同的信号值不同的信号 窗函数比较窗函数比较.VI测量信号的分析与处理（测量信号的分析与处理（1010）3、 数字滤波（数字滤波（1 1） 数字滤波是信号分析与处理的重要步骤，在某些应用数字滤波是信号分析与处理的重要步骤，在某些应用领域，尤其是在一些需要灵活性和编程能力的领域它已取代领域，尤其是在一些需要灵活性和编程能力的领域它已取代了模拟滤波器。了模拟滤波器。 与模拟滤波器相比，数字滤波器有下列优点

22、：与模拟滤波器相比，数字滤波器有下列优点： 1）可用软件编程，设计、修改、优化均十分方便、快捷；）可用软件编程，设计、修改、优化均十分方便、快捷； 2）稳定性好，滤波性能可预测；）稳定性好，滤波性能可预测； 3）不因温度、湿度等影响产生误差，不需高精度元器件）不因温度、湿度等影响产生误差，不需高精度元器件； 4）具有较高的性能价格比。）具有较高的性能价格比。 举举 例：例：用用“仿真信号仿真信号”函数（函数（Express VI）发生一个幅值为）发生一个幅值为1、频率为、频率为10Hz的方波信号，并且叠加幅值为的方波信号，并且叠加幅值为0.1的白噪声。把的白噪声。把此信号送给数字滤波器处理，希

23、望仅得到方波信号中的正弦基此信号送给数字滤波器处理，希望仅得到方波信号中的正弦基波信号；并显示出滤波前后的信号（波信号；并显示出滤波前后的信号（查看滤波效果查看滤波效果）。）。测量信号的分析与处理（测量信号的分析与处理（11 11） 数字滤波（数字滤波（2 2）“函数函数”选板选板 - “信信号处理号处理” - - “波形调波形调理理” -“滤波器滤波器”例：数字滤波例：数字滤波.vifh是高截止频率，以是高截止频率，以Hz为单位。默认值为为单位。默认值为0.45 Hz。 如滤波器类型如滤波器类型为为0（lowpass）或）或1（highpass），），VI将忽略该参数。滤波器类型将忽略该参数

24、。滤波器类型为为2 (Bandpass)或或3 (Bandstop)时，高截止频率：时，高截止频率： fh必须大于低截必须大于低截止频率止频率 fl例题：提取正弦信号例题：提取正弦信号.VI例题：信号谱分析例题：信号谱分析.VI，中值滤波及脉冲测量，中值滤波及脉冲测量.VI左秩是用于计左秩是用于计算左侧中值滤算左侧中值滤波器的元素数。波器的元素数。左秩必须大于左秩必须大于等于等于0。默认值。默认值为为2。 右秩是右秩是用于计算右侧用于计算右侧中值滤波器的中值滤波器的元素数。如右元素数。如右秩小于秩小于0，该，该VI将假定右秩等将假定右秩等于左秩。右阶于左秩。右阶数必须小于数必须小于X。默认值为

25、默认值为-1。 4、时域分析、时域分析例例11-5 脉冲测量脉冲测量.VI峰值检测和显示峰值检测和显示信号的缩放和移位信号的缩放和移位 时域的方法分析工具中有时域的方法分析工具中有6个简单实用的函数。这里个简单实用的函数。这里介绍其中的一个介绍其中的一个 “脉冲测量脉冲测量”函数。它接受一个周期性函数。它接受一个周期性的波形或数组，计算其周期、脉冲宽度、占空比等。的波形或数组，计算其周期、脉冲宽度、占空比等。测量信号的分析与处理（测量信号的分析与处理（1212） 相关分析（相关分析（a a） 为实现对被测信号的有效检测、识别和提取，对被测信号进为实现对被测信号的有效检测、识别和提取，对被测信号

26、进行分析处理中，经常要研究两个信号的相似性，或一信号经过行分析处理中，经常要研究两个信号的相似性，或一信号经过一段延迟后其自身的相似性。这就要对信号进行一段延迟后其自身的相似性。这就要对信号进行“相关相关”分析。分析。相关函数是描述随机函数的重要统计量。相关分析又分相关函数是描述随机函数的重要统计量。相关分析又分“互相互相关关”和和“自相关自相关”。例如通过对一个信号进行自相关分析，可。例如通过对一个信号进行自相关分析，可从噪声很强的信号中检测出是否含有周期成分。从噪声很强的信号中检测出是否含有周期成分。 举例举例：选用选用“卷积和相关卷积和相关”Express VIExpress VI（路径

27、是路径是“函数函数”选板选板“信号处理信号处理”“信号运算信号运算”“卷积和相关卷积和相关”）中的自相）中的自相关函数，来检测由关函数，来检测由“仿真信号仿真信号” Express VI发生的信号（正弦发生的信号（正弦与均有噪声）中与均有噪声）中是否含有周期成分（是否含有周期成分（路径为：路径为：“函数函数”选板选板“信号处理信号处理”“波形生成波形生成”“仿真信号仿真信号”）。）。 测量信号的分析与处理（测量信号的分析与处理（1313） 相关分析（相关分析（b b）例例11-6 相关及谱分析相关及谱分析.VI图图11-12 本例中发生本例中发生3个信号，分别是：个信号，分别是：频率为频率为1

28、0Hz、幅值为、幅值为1的正弦波；的正弦波；幅值为幅值为1的白噪声；的白噪声；前两个信号的叠加。前两个信号的叠加。 另外，设置另外，设置Fs1000，采样数，采样数1000。将这。将这3个信号分别做自个信号分别做自相关，然后再做相关，然后再做FFT，显示自相关后的时域和频域波形。程序，显示自相关后的时域和频域波形。程序框图和前面板如图框图和前面板如图1111和图和图1112所示。所示。 图图11-12中，第一行是正弦波及自相关的波形中，第一行是正弦波及自相关的波形(左左)和功率谱和功率谱(右右)；第二行是白噪声及自相关的波形和功率谱；第三行是正；第二行是白噪声及自相关的波形和功率谱；第三行是正

29、弦波与等幅的白噪声叠加后的波形及自相关的波形和功率谱。弦波与等幅的白噪声叠加后的波形及自相关的波形和功率谱。 由图由图11-12可以看出，周期信号经过自相关后仍然呈现周期可以看出，周期信号经过自相关后仍然呈现周期特征，而白噪声被较大地衰减。因此在它们叠加后形成的波形特征，而白噪声被较大地衰减。因此在它们叠加后形成的波形中可以很明显地发现周期性成分，且从对应的中可以很明显地发现周期性成分，且从对应的3个功率谱上也可个功率谱上也可以看到这一点。以看到这一点。测量信号的分析与处理（测量信号的分析与处理（1414） 曲线拟合（曲线拟合（a a） 在基于计算机实现的测量过程中，为充分利用测在基于计算机实

30、现的测量过程中，为充分利用测得数据，并减小可能的测量误差，经常要对测得结果得数据，并减小可能的测量误差，经常要对测得结果进行所谓曲线拟合。进行所谓曲线拟合。曲线拟合（曲线拟合（curve fittingcurve fitting）是）是指选择适当的曲线类型来拟合观测数据，并用拟合的指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。曲线方程分析两变量间的关系。 曲线拟合的主要作用有：曲线拟合的主要作用有： 1 1）消除测量噪声；）消除测量噪声； 2 2）填充丢失的采样点（例如，如果一个或多个采）填充丢失的采样点（例如，如果一个或多个采样点有丢失，或记录不正确）；样点有丢失

31、，或记录不正确）； 3 3）插值（对采样点之间被测对象可能的取值做出）插值（对采样点之间被测对象可能的取值做出估计）；估计）； 4 4）外推（对采样范围之外被测对象可能的取值进）外推（对采样范围之外被测对象可能的取值进行估计）；行估计）； 5 5）求解某个基于离散数据的被测对象的速度轨迹）求解某个基于离散数据的被测对象的速度轨迹（一阶导数）和加速度轨迹（二阶导数）等。（一阶导数）和加速度轨迹（二阶导数）等。 测量信号的分析与处理（测量信号的分析与处理（1515） 曲线拟合（曲线拟合（b b） 找到找到曲线拟合曲线拟合Express VI的路径：的路径： “函数函数”选板选板“数学数学”“拟合拟合”“曲曲线拟合线拟合