1.一种基于快速独立分量分析和互信息的谐波源定位方法，其步骤为：A、量测所有母线的谐波电压Vn(t)，得到母线的谐波电压矩阵V(t)，V(t)＝[V1(t)，TV2(t)，…，Vn(t)，…，VN(t)]，其中n为母线的序号，t为时间，T为矩阵的转置；

B、用移动平滑滤波器对谐波电压Vn(t)进行滤波，得到谐波电压Vn(t)中的慢速电压分量Vn，s(t)，将谐波电压Vn(t)与其相减，得到谐波电压Vn(t)中的快速电压分量Vn，f(t)；

C、将快速电压分量Vn，f(t)组成的快速电压分量矩阵Vf(t)，Vf(t)＝[V1，f(t)，V2，Tf(t)，…，Vn，f(t)，…，Vn，f(t)] 去均值并白化，得到具有零均值和单位协方差的快速电压分量白化矩阵V′f(t)，V′f(t)＝[V′1，f(t)，V′2，f(t)，…，V′n，f(t)，…，V′N，f(t)]T；

D、对快速电压分量白化矩阵V′f(t)应用快速独立分量分析算法得到分离矩阵W，即所有母线的导纳构成的导纳矩阵W；然后将导纳矩阵W与谐波电压矩阵V(t)相乘，得到与谐波电压矩阵V(t)对应的谐波电流矩阵I(t)，I(t)＝[I1(t)，I2(t)，…，In(t)，…，TIN(t)]，提取谐波电流矩阵中大于零的谐波电流In(t)，得到非零谐波电流矩阵I′(t)＝T[I′1(t)，I′2(t)，…，I′h(t)，…，I′H(t)]，其中h为大于零的谐波电流In(t)的序号，H为大于零的谐波电流In(t)的个数；

E、计算非零谐波电流矩阵I′(t)和谐波电压矩阵V(t)的互信息矩阵M(I′(t)，找出互信息矩阵中每一列h的最大互信息值 这些最大互信息值 对应的谐波电压Vn′(t)所在的母线n′，即为谐波源所在母线。

2.根据权利要求1所述的基于快速独立分量分析和互信息的谐波源定位方法，其特征在于：所述的D步中，对快速电压分量白化矩阵V′f(t)应用快速独立分量分析算法得到导纳矩阵W的具体做法是：选择一个随机的初始矩阵W，按以下迭代公式迭代计算：+ +

W＝W/||W||

+ + + T

式中，W 是W的迭代值，||W||为W 的2范数；E[·]为均值运算；g(WV′f(t))为T T TWV′f(t)的双曲正切函数， 为g(WV′f(t))对WV′f(t)的偏导数。

3.根据权利要求1所述的基于快速独立分量分析和互信息的谐波源定位方法，其特征在于：所述的E步中计算非零谐波电流矩阵I′(t)和谐波电压矩阵V(t)的互信息矩阵M(I′(t)， 的公式是：M(I′(t)，V(t))＝H(I′(t))-H(I′(t)|V(t))式中，H(I′(t))为所有非零谐波电流I′h(t)的概率密度函数p(I′h(t))的微分熵的矩阵，H(I′(t)|V(t))为所有非零谐波电流I′h(t)与所有量测谐波电压Vn(t)之间的微分熵的矩阵， 其中p(Vn(t))为谐波电压Vn(t)的概率密度函数，p(I′h(t)|Vn(t)为谐波电压Vn(t)条件下非零谐波电流I′h(t)的概率密度函数。