基于数字滤波技术的虚拟频率补偿仪设计.docVIP

基于数字滤波技术的虚拟频率补偿仪设计 1．设计原理 动态特性补偿目的 对测量系统(传感器)来说，大量被测物理量是随时间变化的动态信号，即被测量是时间t的函数x(t)，而不是常量。系统的动态特性反映它测量动态信号的能力。理想的传感器系统，其输出量y(t)与输入量x(t)的时间函数表达式相同。但实际上，两者只能在一定频率范围内，在允许的动态误差条件下保持所谓的一致。 常用的传感器多是一阶、二阶系统，频率特性都有一定的频带限制，频谱幅度随着频率的增加总体呈下降的趋势。这样就导致了实际测量中存在动态误差。两类系统的动态误差可用下式描述: 上式中， 为测量系统的幅频特性，为理想幅频特性，对于一阶、二阶系统为直流放大倍数。 不同角频率的输入信号通过系统都将产生不同程度的动态误差。两类系统的动态误差值列入表5－1中。 注：二阶系统的一组值是在阻尼系数为0的情况下得到的。表中ωc为一阶系统转折角频率， ωn为二阶系统固有角频率。 从表5—1中看到，如果允许动态误差为1％，则一阶系统输入信号频率必须小于转折频率的1／7，二阶系统输入信号频率必须小于固有频率的1／10。所以，当被测信号的频率高，而测量系统的工作频带不能满足测量允许误差的要求时，则希望扩展系统的频带以改善系统的动态性能，称为系统动态特性补偿。 对系统的频带拓宽目前采用滤波法或频域校正法等方法进行补偿，这些方法需要事先知道系统的频率特性。所以，用硬件实现补偿虽然未尝不可，但缺乏灵活性，无法针对不同系统设计成任意需要的扩展倍数；而用软件实现，可任意设置补偿目标和补偿参数，具有硬件方法不可比拟的优越性。所以，这里尝试用虚拟仪器的形式实现系统动态特性补偿仪。 (2)补偿原理 此处采用了数字滤波法，其补偿原理框图如下图所示。 Y(s)U(s) Y(s) U(s) H’(s) H’(s) H(s) 设传感器系统的传递函数为H(s)，在系统后串接一个补偿环节H’(s)，使得最后的输出Y(s)比传感器输出U(s)更接近被测量X(s)。 一阶系统的补偿原理 当系统是一阶系统时，。若补偿环节，则补偿的结果是Y(s)=X(s)。但是这个补偿环节是不稳定的，因为从该系统的幅频特性： （1） 可以看出，当时系统频率响应趋向∞，所以完全的补偿是无法实现的。 通常的做法是先确定系统频带欲扩展的倍数A，由A确定目的传递函数，即 （2） 再由W(s)=H(s)×H’(s)，得到补偿环节的传递函数： （3） = 2 \* GB3 ②二阶系统的补偿原理 当系统是二阶系统时，传递函数为，和一阶系统一样，对二阶系统也无法实现稳定的完全补偿。通常也是先确定目的函数，然后推出补偿环节的传递函数。二阶系统比一阶系统稍微复杂一点，因为它的幅频特性并不是单调的下降趋势，而是一个峰值，动态误差有正有负，且同时和两个频域参数有关——固有频率ωn和阻尼系数ξ。 从二阶系统动态误差表达式可以知道，在工作频带内，系统动态误差γ和ξ、ωn成反比，所以，要改善系统的动态误差，需要增大系统的固有频率ωn和阻尼系数ξ。 假设补偿后的系统的固有频率为ω’n、阻尼比系数ξ’，即有 （4） 则补偿环节的传递函数 （5） 但是，在补偿时，如果同时考虑两个可变参数，情况过于复杂，补偿环节也不好实现和控制，因而考虑对其中一个参数取一个合适的固定值。考虑到固有频率和系统截止频率有着更密切的关系，而阻尼系数则更多与某一频段的系统响应幅值有关，所以可将阻尼系数取一个固定值。 当固有频率取值固定时，阻尼系数ξ和系统截止频率ωh的关系如图： ξ－ωh曲线图 （ωn＝1） 从图看到，系统频带随阻尼系数的增大而减小，当ωh＝ωn’时，两者有一对一的关系，且此时系统处于临界阻尼状态，超调量小，很近似一阶系统， 所以取补偿后系统阻尼系数为ξ＝0.707最合适。 = 3 \* GB3 ③数字滤波器的离散表达式 上面的补偿环节可以用模拟滤波器实现，也可以用软件设计数字滤波器实现。采用后向差分法时，通过等效变换可以得到数字滤波器的表达式。 由式（3）得到一阶系统的补偿环节等效数字滤波器表达式为： , ,b=Ac 式中τ为系统原时间常数，A为系统频带扩展倍数。软件实现的差分方程表达式为： （6） 由式(5)得到二阶系统的补偿环节等效数字滤波器表达式为： （7） 其中，ωn’是补偿后的系统固有频率，有ωn’＝Aωn，ξ’是补偿后的系统阻尼系数，有ξ’=0.707。 软件实现的差分方程表达式为： （8） 其中：K＝ b0＝ b1＝－ b2＝