雷达物位计-软件算法说明

法兰距离计算流程图信号加窗信号加窗用于减小频谱泄露，可选择三种窗函数之一：矩形窗、汉宁窗、布莱克曼窗。假设中频信号电压采样数组为v[N]，采样点数为N（N=1199）；加窗实际上是构建一个N点的数组w[N]，将v[N]和w[N]进行点乘运算；信号加窗后的输出数组x[N]可表示为： xn矩形窗矩形窗的w[N]是一个全1的数组，相当于对v[N]不进行加窗处理。汉宁窗汉宁窗的w[N]表示如下： wn布莱克曼窗布莱克曼窗的w[N]表示如下： wn信号补零将N（N=1199）点的数组x[N]扩展成M（M=4096）点的数组y[M]： ynFFT运算对M点的数组y[M]进行实数快速傅里叶变换并求模，具体算法请参考数字信号处理相关书籍；TMS320F28335有现成的库函数RFFT_f32(RFFT_F32_STRUCT*)、RFFT_f32_mag(RFFT_F32_STRUCT*)供调用，FFT运算的输出为数组Y[M/2]。频谱峰值探测根据MATLAB信号处理工具箱中PeakFinder的算法思路：处于两个波谷间且大于波谷一定阈值的点即为波峰；对应的流程图为： A0回波跟踪如果选择首回波算法，则回波跟踪过程可以用如下状态机表示：如果选择回波搜索算法，则回波跟踪过程可以用如下状态机表示：回波跟踪的输出为液面状态、位置和幅度。滤波处理有三种滤波方法：跳跃滤波、阻尼、最小二乘法滤波；其中阻尼不能和最小二乘法滤波同时使用，但跳跃滤波能和阻尼或最小二乘法滤波同时使用。跳跃滤波跳跃滤波用于减少液面的大幅度快速抖动，跳跃滤波用到两个参数：跳跃滤波阈值JFilt_Thresh、跳跃滤波延迟JFilt_Delay，其算法思路如下：如果当前液位与前一时刻液位输出值相距大于JFilt_Thresh，跳跃滤波计数器累加；否则跳跃滤波计数器清零,用当前液位值作为当前时刻的输出。如果跳跃滤波计数器大于JFilt_Delay，用当前液位值作为当前时刻的输出；否则用前一时刻液位输出值作为当前时刻的输出。阻尼阻尼用距离滤波因子来度量，施加阻尼后系统的输出为： Dn'= 式中，Dn'— Dn-1'— Dn—当前时刻的 DistFiltFactor—距离滤波因子距离滤波因子和阻尼时间常数之间可用如下关系式表示： DistFiltFactor=1-e-∆TT 式中，T—时间常数，单位为s ∆T—液位计算的时间间隔(0.1s)最小二乘法滤波最小二乘法滤波用于快速跟踪液位的变化，其算法步骤如下：利用前N（N=10）个时刻的输出值D'0、D'∆T、D'2∆T D''t=a+b*t利用拟合直线，计算当前时刻的预测值： D''N∆T=a+b*N∆T根据当前时刻的预测值和测量值，计算当前时刻的输出值：D'N∆T=D式中，DN∆T注：最小二乘法线性拟合见附录A。液位计算以用户定义的液位零点为参考点，计算出的液位值为： Level=TankHeight_R-Flange_D-OffsetDist_G+CalibrationLevel+LevelOffset(2-1)式中，TankHeight_R—储罐高度 Flange_D—上一步计算出来的法兰距离 OffsetDist_G—距离偏移量 CalibrationLevel—液位校准值 LevelOffset—液位修正量距离计算以用户定义的距离零点为参考点，计算出的距离值为： Ullage=Flange_D+OffsetDist_G-CalibrationLevel+UllageOffset(3-1)式中，UllageOffset—距离修正量体积计算有4种体积计算方法：插值法、理想球罐法、理想垂直圆柱罐法、理想水平圆柱罐法。插值法使用拉格朗日插值法进行体积计算；假设体积插值表可表示为L0，V0、L1，V1、…、Level_V=TankHeight_R-Flange_D-OffsetDist_G+CalibrationLevel+LevelOffs(4-1)式中，LevelOffs—插值表液位偏移量线性插值寻找点Li，VV'= 如果Level_V<L0，用如果Level_V≥LN-1，用二次曲线插值法寻找点Li，Vi、 V + +Level_V-LiLevel_V-Li+1Li+2-如果Level_V<L0如果Level_V≥LN-2插值法的体积输出值为： Volume=V式中，VolumeOffs—体积修正量理想球罐法Level_V=TankHeight_R-Flange_D-OffsetDist_G+CalibrationLevel+ZeroLevelToBottomVolume=1式中，ZeroLevelToBottom—液位零点到罐底的距离L1—球罐理想垂直圆柱罐法Level_V=TankHeight_R-Flange_D-OffsetDist_G+CalibrationLevel+ZeroLevelToBottomVolume=1式中，L1—圆柱体理想水平圆柱罐法Level_V=TankHeight_R-Flange_D-OffsetDist_G+CalibrationLevel+ZeroLevelToBottomVolume=-+VolumeOffs(4-7) 式中，L1— L2—液位速率计算液位速率由下式计算： LevelRate=Level_V- 式中，Level_V' ∆T—液位计算的时间间隔(0.1s)

附录A最小二乘法线性拟合最小二乘法线性拟合就是将一组符合Y=a+bX关系的测量数据，用计算的方法求出最佳的a和b。假设直线方程的表达式为：要根据测量数据求出最佳的a和b。对满足线性关系的一组等精度测量数据（xi，yi），假定自变量xi的误差可以忽略，则在同一xi下，测量点yi和直线上的点a+bxi的偏差di如下：显然最好测量点都在直线上（即d1=d2=……=dn=0），求出的a和b是最理想的，但测量点不可能都在直线上，这样只有考虑d1、d2、……、dn为最小，也就是考虑d1+d2+……+dn为最小，但因d1、d2、……、dn有正有负，加起来可能相互抵消，因此不可取；而|d1|+|d2|+……+|dn|又不好解方程，因而不可行。现在采取一种等效方法：当d12+d22+……+dn2对a和b为最小时，d1、d2、……、dn也为最小。取（d12+d22+……+dn2）为最小值，求a和b的方法叫最小二乘法。令＝D对a和b分别求一阶偏导数为： 再求二阶偏导数为：；显然：；满足最小值条件，令一阶偏导数为零：引入平均值：； ；； 则：解得： 将a、b值带入线性方程，即得到回归直线方程。会签第一次评