本发明属于电气工程和电磁兼容领域，具体指一种宽带载波通讯设备的辐射噪声诊断方法。
背景技术：
低压电力线宽带载波通信是利用低压电力配电线作为信息传输媒介进行语音或数据传输的一种特殊通信方式。该技术是把载有信息的高频信号加载于电流，然后利用各种等级的电力线传输，接受信息的调制解调器再把高频信号从电流中分离出来，并传送到电力线宽带用户终端。该技术在不需要重新布线的基础上，在现有电线上实现数据、语音和视频等多业务的承载。
然而，宽带载波通讯设备的通讯频段为2-12mhz，通讯速率可达10mbps，其通讯频段内和通讯频段外会产生大量的辐射噪声，严重影响周围电子设备和电网的安全可靠性。现有辐射源诊断与测试方法，是利用emi接收机、电波暗室和远场天线能够确定非同频辐射噪声源，但无法诊断出同频辐射噪声源。此外，中国发明专利(zl200810242655.4)公开了一种《诊断辐射emi机理的实验台及辐射emi机理的简易诊断方法》，能够通过测量近场波阻抗，判定辐射噪声源的机理，包括共模噪声和差模噪声，但该专利技术也无法诊断出同频辐射噪声源。
技术实现要素：
为了克服上述不足，本发明提供一种宽带载波通讯设备的辐射噪声诊断方法，包括以下步骤：
步骤1：采用多通道高数数字示波器和m个近场磁场探头同时测量载波通讯设备内部电路板不同位置产生的辐射磁场，测量结果分别为xij(t)，构造辐射磁场时域信号矩阵x(t)：
其中，j为1至m，xij(t)为时间t的函数，xij(t)为测试场点j测量得到的第i个辐射磁场时域信号，每组xij(t)的采样点个数为n，i为1至n；
步骤2：第二步：求解x(t)第j列的数学期望值ej：
步骤3：对x(t)按式(3)进行信号处理，得到y(t)：
步骤4：求解y(t)的协方差矩阵c(t)：
式中，yt(t)的是y(t)的转置矩阵；
步骤5：求解c(t)的特征值λj及其对应的特征向量uj，得到λ和u矩阵，
u＝[u1u2…un]\*mergeformat(6)
步骤6：求解白化矩阵z，
z＝vy(t)＝λuty(t)\*mergeformat(7)
式中，v＝λut；
步骤7：设任意矩阵m是n×n的矩阵，
式中，qij是四阶累积量矩阵q第i行第j列对应的值，kijkl(z)是z中第i，j，k，l四个量的四阶张量，mkl是矩阵m第k行第l列对应的值；
步骤8：重复迭代求解w，使得式(9)中d最小，
d＝wtqw\*mergeformat(9)
步骤9：辐射噪声源的分量s：
s＝wty(t)\*mergeformat(10)
进一步，所述步骤1中多通道高数数字示波器为tektronixdpo5204b。
进一步，所述的宽带载波通讯设备的辐射噪声诊断方法，其特征在于：所述步骤1中近场磁场探头为rohde&schwarz公司的hz-11。
与现有技术相比，本发明的有益效果是：先利用多通道高速数字示波器同时测量被测电路在不同场点产生的辐射噪声时域信号，再利用信号分析方法得到引起辐射噪声超标的噪声源特征，从而实现宽带载波设备的辐射噪声源的有效诊断，有助于提高电网的安全可靠性；本发明主要应用于电力线宽带载波通信设备，亦可拓展应用于医疗器械设备、电力设备、大功率变流器、汽车电子、轨道交通装备、网络通讯设备等电子设备中，能够有效诊断出相关电子设备辐射噪声超标的原因，为后续辐射噪声抑制提供理论和实践依据，从而使相关电子设备通过gb9254、gb4824、yy0505等电磁兼容标准检测。
附图说明
图1为实验测试系统的示意图
具体实施方式
本发明先利用多通道高速数字示波器同时测量被测电路在不同场点产生的辐射噪声时域信号，再利用信号分析方法得到辐射噪声的时频二重域信号，然后根据超标频点的大小提取辐射噪声超标频率的时域信号，最后利用非线性信号分离算法确定辐射噪声源的分量，从而实现同频辐射噪声源的有效诊断，有助于降低电子设备产生的辐射噪声。
下面结合附图和具体实施例做进一步说明，本发明包括以下步骤：
步骤1：采用多通道高数数字示波器tektronixdpo5204b和m个近场磁场探头rohde&schwarz公司的hz-11同时测量载波通讯设备内部电路板3个不同位置产生的辐射磁场，测量结果分别为xij(t)，j为1至3，xij(t)为时间t的函数，xij(t)为测试场点j测量得到的第i个辐射磁场时域信号，每组xij(t)的采样点个数n为1024，i为1至1024，构造辐射磁场时域信号矩阵x(t)，
步骤2：求解x(t)第j列的数学期望值ej
步骤3：对x(t)按式(3)进行信号处理，得到y(t)
步骤4：求解y(t)的协方差矩阵c(t)
式中，yt(t)的是y(t)的转置矩阵；
步骤5：求解c(t)的特征值λj及其对应的特征向量uj，得到λ和u矩阵，
u＝[u1u2…un]\*mergeformat(16)
步骤6：求解白化矩阵z，
z＝vy(t)＝λuty(t)\*mergeformat(17)
式中，v＝λut；
步骤7：设任意矩阵m是n×n的矩阵，
式中，qij是四阶累积量矩阵q第i行第j列对应的值，kijkl(z)是z中第i，j，k，l四个量的四阶张量，mkl是矩阵m第k行第l列对应的值；
步骤8：重复迭代求解w，使得式(9)中d最小，
d＝wtqw\*mergeformat(19)
步骤9：辐射噪声源的分量s为
s＝wty(t)\*mergeformat(20)
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。