第13章labview+数学分析与信号处理

1、1第第13章章 数学分析与信号处理数学分析与信号处理2引子nLabVIEW作为自动化测试、测量领域的专业软件，其内部集成了600多个分析函数，用于信号生成、频率分析、概率、统计、数学运算、曲线拟合、插值、数字信号处理等等各种数据分析应用。n此外，LabVIEW还提供了附加工具软件专业应用于某些信号处理应用中，如声音与振动、机器视觉、RF/通信测量、瞬态/短时持续信号分析等等。3本章内容n13.1 数学分析数学分析n13.1.1 图形化编程与数学分析图形化编程与数学分析n13.1.2 基本数学函数基本数学函数n13.1.3 线性代数线性代数n13.1.4 曲线拟合曲线拟合n13.1.7 概率与统

2、计概率与统计n13.1.9 常微分方程常微分方程n13.1.12 MathScript4本章内容n13.2 数字信号处理数字信号处理n13.2.1 信号发生信号发生n13.2.3 波形测量波形测量n13.2.5 频域分析频域分析n13.2.7 数字滤波器数字滤波器n13.2.8 逐点分析库逐点分析库513.1.1 图形化编程与数学分析图形化编程与数学分析nLabVIEW具有强大的数学分析能力数学分析数学分析VI函数面板函数面板 613.1.1 图形化编程与数学分析图形化编程与数学分析子面板名称子面板名称描述描述Numeric最基本的数学操作，例如加减乘除、类型转换和数据操作等。Elementa

3、ry & Special Functions一些常用的数学函数，例如正余弦函数、指数函数、双曲线函数、离散函数和贝塞尔函数等。Linear Algebra线性代数，主要是矩阵操作的相关函数Fitting曲线拟合和回归分析Interpolation & Extrapolation一维和二维的插值函数，包括分段插值、多项式插值和傅立叶插值Integration & Differentiation积分与微分函数Probability & Statistics概率与统计Optimization最优化Differential Equations解常微分方程Geometry

4、几何Polynomial多项式计算和分析Scripts & Formulas脚本节点、公式节点以及公式解析的相关函数713.1.2 基本数学函数基本数学函数n基本数学函数分为12类：三角函数、指数函数、双曲线函数、门函数、离散数学函数、贝塞尔函数、函数、超几何分布函数、椭圆积分、指数函数、误差函数和椭圆抛物函数。基本数学函数面板 813.1.3 线性代数线性代数n强大的矩阵运算能力线性代数函数面板913.1.3 线性代数线性代数n例例13.2 解线性方程组解线性方程组Ax=b,其中其中2 . 05 . 32185 . 07328 . 032A=b= 1013.1.4 曲线拟合曲线拟合n

5、二维曲线拟合就是根据输入数据的坐标（xi,yi），即X数组和Y数组，找出yi和xi的函数关系y=f(x)。对于不同的对象，有不同的拟合方法：n线性拟合（Linear Fit）、指数拟合(Exp fit)、幂拟合(Power Fit)、高斯拟合（Gauss Peak Fit）、对数拟合（Logarithm Fit）、多项式拟合（Polynomial Fit）、最小二乘法拟合(Gen. LS Lin. Fit)和非线性拟合（Nonlinear Curve Fit）等等。1113.1.4 曲线拟合曲线拟合n例例13.3 最小二乘法曲线拟合举例最小二乘法曲线拟合举例利用最小二乘法拟合曲线，将因变量y与

6、自变量x的关系表达为本例中：假设猜测函数为：ninniixfaxfaxfaxfaxafy01100)()()()(),(Noisexxxxy14)cos(3)sin(2)()()()()(4433221100 xfaxfaxfaxfaxfay1213.1.4 曲线拟合曲线拟合其中：下面我们通过最小二乘法拟合函数 General LS Linear Fit.vi 来求解回归系数4432210)(1)()cos(3)()sin()(1)(xxfxxxfxxfxxfxf131413.1.7 概率与统计概率与统计概率与统计函数面板概率与统计函数面板 1513.1.7 概率与统计概率与统计n例例13.5

7、 概率与统计函数举例概率与统计函数举例n该例中首先通过该例中首先通过Gaussian White Noise.vi产生一产生一个满足高斯分布的随机数序列，然后通过个满足高斯分布的随机数序列，然后通过Create Histogram和和Statistic两个两个Express VI对该随机对该随机序列进行分析。序列进行分析。 1613.1.7 概率与统计概率与统计1713.1.9 常微分方程常微分方程常微分方程函数面板常微分方程函数面板 1813.1.9 常微分方程常微分方程函数名称函数名称功能功能ODE Solver.vi解带初值的常微分方程：X=F(X,t)ODE Runge Kutta 4

8、th Order.vi用龙格库塔方法解带初值的常微分方程ODE Cash Karp 5th Order.vi用Cash Karp方法解带初值的常微分方程ODE Euler Method.vi用欧拉方法解带初值的常微分方程ODE Linear nth Order Numeric.vi用数值解法解n阶线性齐次常微分方程ODE Linear nth Order Symbolic.vi用符号解法解n阶线性齐次常微分方程ODE Linear System Numeric.vi解一个带有常系数微分方程的n维齐次线性系统，结果为数值解ODE Linear System Symbolic.vi解一个带有常系数

9、微分方程的n维齐次线性系统，结果为符号解1913.1.9 常微分方程常微分方程n例例13.7 常微分方程数值解举例常微分方程数值解举例设河边点O的正对岸为点A，河宽OAh，两岸为平行直线，水流速度为a，有一鸭子从点A游向点O，设鸭子（在静水中）的游速为b（ba）,且鸭子游动方向始终朝着点O.求鸭子游过的迹线方程。鸭子游过的迹线鸭子游过的迹线2013.1.9 常微分方程常微分方程通过分析得到迹线微分方程：使用ODE Runge Kutta 4th Order.vi函数求解2222yxbydtdyyxbxadtdx2113.1.9 常微分方程常微分方程222313.1.12 MathScriptn

10、MathScript是LabVIEW 8以后版本推出的面向数学的文本编程语言，它带有交互式的窗口和可编程的接口。通过MathScript，喜欢文本编程的用户可以在LabVIEW中编写并执行MATLAB式的文本代码（.m文件）并能与图形化编程无缝结合。n新的MathScript包含了600多个数学分析与信号处理函数，并增加和增强了丰富的图形功能。 2413.1.12 MathScriptn使用MathScript的方法有两种：n使用MathScript交互式窗口；通过交互式窗口，你可以像使用MATLAB一样执行命令、编译运行.m脚本文件、查看运行结果等。n在程序框图中使用MathScript节点

11、。 2513.1.12 MathScriptn使用MathScript交互式窗口2613.1.12 MathScriptn运行结果2713.1.12 MathScriptn在程序框图中使用MathScript节点2813.1.12 MathScript2913.2 数字信号处理数字信号处理n高效、灵活、强大的数字信号处理能力数字信号处理函数面板数字信号处理函数面板3013.2 数字信号处理数字信号处理子面板子面板描述描述Waveform Generation通过该VI函数面板可以产生各种不同类型的波形信号Waveform Conditioning用于波形信号的数字滤波和窗函数等信号调理Wave

12、form Measurements波形信号测量面板，用来实现常见的时域和频域的测量，譬如直流交流成分分析、振幅测量、傅立叶变换、功率谱计算、谐波畸变分析、频率响应和信号提取等Signal Generation按照具体的波形模式产生一维实数数组表示的信号。Signal Operation对信号进行各种操作，例如卷积、自相关分析等。Windows窗函数分析Filters实现IIR、FIR和非线性滤波Spectral Analysis实现基于数组的谱分析Transforms信号处理中各种常见的变化函数Point By Point逐点分析函数库3113.2.1 信号发生信号发生n波形发生函数可以用来模

13、拟产生你需要的各种波形。nLabVIEW有两个信号发生函数面板，其中Waveform Generation用于产生波形数据类型表示的波形信号，Signal Generation用于产生一维数组表示的波形信号。 3213.2.1 信号发生信号发生nWaveform Generation3313.2.1 信号发生信号发生nSignal Generation3413.2.1 信号发生信号发生3513.2.1 信号发生信号发生3613.2.3 波形测量波形测量n波形测量面板提供的VI函数用于对波形的各种信息进行测量，譬如直流交流分析、振幅测量、脉冲测量、傅立叶变换、功率谱测量、谐波畸变分析、过渡分析、

14、频率响应等。 3713.2.3 波形测量波形测量n例例13.9 测量波形的直流分量和有效值测量波形的直流分量和有效值3813.2.3 波形测量波形测量3913.2.5 频域分析频域分析n频域分析函数被划分为两个面板：nTransforms面板实现的函数功能主要有傅立叶变换、Hilbert变换、小波变换、拉普拉斯变换等。nSpectral Analysis面板包含的函数主要包括功率谱分析、联合时频分析等。 4013.2.5 频域分析频域分析4113.2.5 频域分析频域分析n例例13.11 信号的傅立叶变换信号的傅立叶变换4213.2.5 频域分析频域分析4313.2.5 频域分析频域分析n转变

15、为单边傅立叶变换4413.2.5 频域分析频域分析4513.2.5 频域分析频域分析n例例13.12 非均匀采样数据的功率谱计算非均匀采样数据的功率谱计算4613.2.5 频域分析频域分析4713.2.7 数字滤波器数字滤波器n滤波器分为模拟滤波器和数字滤波器。传统模拟滤波器的输入与输出都是连续的。而数字滤波器的输入与输出都是离散时间信号。n数字滤波器具有如下好处：n软件可编程，因此易于搭建和测试n只需要加减乘三种基本数学操作n不随外界环境条件变化而漂移，也不会老化n有非常高的性价比4813.2.7 数字滤波器数字滤波器nLabVIEW提供的IIR滤波器类型有Butterworth、Cheby

16、shev、Inverse Chebyshev、Elliptic和Bessel滤波器。它们都有各自的特点，用途也不尽相同。4913.2.7 数字滤波器数字滤波器nLabVIEW还提供了高级IIR和FIR滤波器子面板。在高级面板中，滤波器的设计部分和执行部分是分开的。由于滤波器的设计很费时间，而滤波过程则很快。在含有循环结构的程序中，可以将滤波器的设计放在循环外，将设计好的滤波器参数传递到循环内，在循环内进行滤波，从而提高程序的运行效率。 5013.2.7 数字滤波器数字滤波器n滤波器选择向导滤波器选择向导 5113.2.7 数字滤波器数字滤波器n滤波器函数面板滤波器函数面板 5213.2.7 数

17、字滤波器数字滤波器n例例13.15 低通滤波举例低通滤波举例5313.2.7 数字滤波器数字滤波器5413.2.8 逐点分析库逐点分析库n逐点分析中，数据分析是针对每个数据点的，对采集到的每一点数据都可以立即进行分析，而且分析可以是连续进行的。因此通过实时分析，用户可以实时的观察到当前采集数据的分析结果，从而使用户能够跟踪和处理实时事件。 5513.2.8 逐点分析库逐点分析库5613.2.8 逐点分析库逐点分析库n逐点分析函数的基本结构逐点分析函数的基本结构 5713.2.8 逐点分析库逐点分析库n例例13.16 基于逐基于逐点分析的实时滤波点分析的实时滤波5813.2.8 逐点分析库逐点分析库59小试身手小试身手n1. LabVIEW在数学分析与信号处理方面具有哪些优势？n2. 利用基本数学函数面板提供的函数画出如下公式在-2pi,2pi之间的曲线。(利用Signal Generation下的函数产生-2pi,2pi之间的均匀采样点)xexyxcos560小试身手小试身手n3. 求解线性方程组求解