10章信号处理-频域分析解读.ppt

六、曲线拟合 曲线拟合（Curve Fitting）技术是用于科学、工程、商业、医学和其他数据分析领域中的常见技术。该技术主要是用于从一组数据中提取曲线参数或者系数，以得到这组数据的函数表达式。通过曲线拟合，可以用连续模型表示离散的数据。曲线拟合的实际应用很广泛。 曲线拟合过程中的主要算法就是最小二乘法。 LabVIEW的分析软件库提供了多种线性和非线性的曲线拟合算法，例如线性拟合、指数拟合、通用多项式拟合、非线性Levenberg-Marquardt拟合等。 LabVIEW的曲线拟合模块位于All Functions→Analyze→Mathematics→Curve Fitting子模板 。 曲线拟合子模板 [例] 给定数据进行线性拟合。 这个例子收集了10对实验数据x和y，如表格中所示，采用线性关系进行拟合。 X[i] 0 1 2 3 4 5 6 7 8 9 Y[i] 0.9 1.4 3.3 4 4.92 7 7.4 8.1 9 10 线性拟合的前面板设置和框图程序 (1)在主框图程序中，使用了如下的子程序: Linear Fit.vi子程序（在Analyze→Mathematics→Curve Fitting子模板）。在本例中，Linear Fit.vi把实验数据拟合为一条直线，求出系数，以满足；还得出了实验结果和拟合结果之间的误差的均方根值 (2)注意：输入数据是一个二维数组，是采用DAQ卡采集数据时的通用格式。我们用Index Array子程序可以分解得到两个一维数组和。 另外，本例求出了“MSE”，它表示误差均方根值。误差越小，拟合结果越好。 (3)运行该程序。曲线将显示实验数据和拟合结果。 [例]给定数据进行多项式拟合。 实验数据和上例中数据一样，采用多项式进行拟合。 (1) 本例中假设实验数据之间为多项式关系，即 (2) 当多项式为一阶时，有两个参数，结果是一个线性拟合，与前一个例子相同。当阶数为二阶时，结果是一个二阶多项式带三个参数，多项式参数存放于数组。 （3） 在框图程序中使用了如下的子程序： General Polynomial Fit.vi子程序（在Analyze→Mathematics→Curve Fitting子模板）。在本例中,该函数的功能是把实验数据拟合为二阶多项式曲线，返回结果为拟合数据[i]，多项式系数a[i]以及拟合误差的均方根值MSE。 （4）本例使用多项式拟合，多项式阶数为2，求出了多项式系数。一般情况下，我们总是尽可能使用最低阶的多项式。 （5）运行程序，将显示实验数据和拟合结果的数据。 多项式拟合的前面板设置和框图程序 信号的频域分析 傅里叶变换是信号处理和数据处理中的一个重要的分析工具，其意义在于将时域与频域信号联系起来。 LabVIEW高级分析程序库中的频域分析子模板提供了丰富的时域和频域转换函数，包括傅里叶变换、Hilbert变换、Hartley变换、功率谱分析、谐波分析等。频域分析子模板位于All Functions模板下Analyze→Signal Processing→Frequency Domain子模板。 频域分析子模板 一、 FFT变换 FFT的输出都是双边的，它同时显示了正负频率的信息。通过只使用一半FFT输出采样点可转换成单边FFT。FFT的采样点之间的频率间隔是fs/N，这里fs是采样频率。 Analyze库中有两个可以进行FFT的VI，分别是Real FFT VI和Complex FFT VI 。 例 双边和单边FFT变换 双边FFT变换的前面板设置和框图程序 (1) 框图程序中的各函数作用如下： Sine Waveform.vi（All Functions→Analyze→Waveform Generation）：产生时域正弦波形。 Real FFT.vi（实数快速傅里叶变换）：用于计算输入数据的FFT，将时域信号转换为频域信号。输入为实数数组，输出为复数数组。 Array Size 函数（All Functions→Array模板）：用于根据采样点数N对FFT输出的结果进行处理。将FFT输出除以N，可获得正确的频率幅度信息。 Complex to Polar.vi（All Functions→Numeric→Complex子模板）：将输入数据从复数坐标系转换到极坐标系。此例将FFT输出分解为实部和虚部（幅值和相位），相位的单位是弧度，这里只需显示FFT的幅值 (2) 把该VI保存为LabVIEW\lianxi中的FFT_2sided.vi。 (3)