质量保证协议.docVIP

变压器传递函数的系统辩识 朱明林　朱子述 摘　要：利用最小二乘法，由变压器绕组频率响应辩识其传递函数的解析表达式，分析了影响辩识准确性的因素，提出了改进措施，提高了辩识的准确性。关键词：系统辩识　病态矩阵　变压器　传递函数 System Identification of Transformer′s Transfer Function Abstract：A methodology to identify analytical expression of transformer winding′s transfer function from frequency response is presented based on least squares. Some factors influencing result are discussed. Then improvement steps are put forward to improve the accuracy.Key words：system identification　morbid matrix　transformer　transfer function▲ 0　引　言 　　据上海电网资料统计，1991～1998年期间，因抗短路能力不强，造成15台次220 kV以上主变损坏，发生绕组变形的现象更加普遍。检测变压器绕组变形通常采用频响法(FRA)［1］，即同一台变压器和历史记录的频响比较、同一台变压器不同相之间的频响比较、相同型号变压器之间的频响比较等，都是通过比较来判断变压器绕组的状况，如资料少则很难判断。因此，需要定量地分析变压器的频率响应，以建立数学模型，即根据变压器绕组频响辩识出其传递函数的解析表达式。 1　方法和误差分析 1.1　最小二乘法辩识变压器传递函数［2］　　所谓系统辩识，也叫动态建模，就是利用系统的输入、输出数据，经计算处理后，估计出系统的数学模型。变压器绕组可以当成一个由电阻、电感和电容等参数组成的无源线性双口网络，在频域上可以用传递函数H(jω)来描述。 　(1) 其中于是μ＜v，一般情况下，μ＝v－1。当扫频点频率ω＝ωk时， 　(2) 其中(3) 令式(2)实部和虚部分别相等得： 　(4) 未知系数有a0，a1，…，aν－1和b1，b2，…，bν共n=2ν个，假设频响扫频点共有m个，则求解传递函数H(jω)系数的方程为： Ax=b，　(5) 其中 (6) 　　通常频响扫频观测点数m传递函数H(jω)系数的总个数n，方程(5)存在最小二乘解： x=(ATA）－1ATb，　(7) 1.2　最小二乘辩识变压器传递函数的误差分析　　测量一台变压器C相低压绕组传递函数频响曲线的接线图见图1，测量结果见图2曲线1，扫频范围为0.5～500 kHz，共1 500个扫频点。同时，采用最小二乘法辩识出该变压器的传递函数，画出的频响曲线见图2中的曲线2，其中阶数ν＝20，比较曲线1、2，可见两者有较大误差。 图1　传递函数频响测量接线图 图2　变压器绕组传递函数频响曲线 　　因式(5)中b的元素是测得的，矩阵A的元素是在b的基础上计算得的，测量观测误差，加上计算舍入误差，使矩阵A和b形成微小的扰动δA和δb，这些扰动都会对最终解x产生影响。假设δx是x由δA和δb引起的变化，则式(5)变为： (A+δA）（x＋δx）＝b＋δb，　(8) 进一步得到：(A+δA）δx＝δb－δA.x。设A－1是矩阵A的广义逆，也约等于(A+δA)的广义逆，则： δx＝A－1（δb－δA．x)　(9) 两边取范数并取不等式： 　(10) 可见，矩阵A和b的相对误差在最终解中可能被放大‖A－1‖．‖A‖倍，‖A－1‖．‖A‖被定义为矩阵A的条件数。通常使用的条件数有：　　1) cond(A)∞＝‖A－1‖∞．‖A‖∞，其中：称为A的行范数；　　2) cond(A)2＝‖A‖2．‖A－1‖2＝其中λmax和λmin分别是矩阵ATA的最大特征值和最小特征值。　　当矩阵A的条件数很大时，矩阵A和b存在的微小扰动δA和δb，会引起方程组Ax=b解的很大变化，则称Ax=b为病态方程组，称矩阵A为病态矩阵。矩阵A的条件数愈大，方程组病态程度愈严重，也就愈难得到方程组准确的解［3］。因此，图2中实测频响曲线和最小二乘辩识出传递函数的频响曲线之间的误差，正是由于方程(5)中的矩阵A是病态矩阵的结果。 2　影响因素和改进措施 2.1　影响因素　　由式(6)矩阵A元素的组成可以看出，最小二乘辩识变压器传递函数方法中，矩阵A的条件数主要依赖于以下三个因素：　　1) 传递函数的阶数ν　　因矩阵A的列数n=2ν，所以传递函数的阶数ν将影响矩阵A的大小；同时随着ν的增大，矩阵A中