CN110309545A - 基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 - Google Patents
基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 Download PDFInfo
- Publication number
- CN110309545A CN110309545A CN201910480723.9A CN201910480723A CN110309545A CN 110309545 A CN110309545 A CN 110309545A CN 201910480723 A CN201910480723 A CN 201910480723A CN 110309545 A CN110309545 A CN 110309545A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- data
- matrix
- indicate
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computational Linguistics (AREA)
- Algebra (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Health & Medical Sciences (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Databases & Information Systems (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,根据拉依达(Pauta)准则,对公共连接点(PCC)处谐波电压与电流的测量数据进行数据清洗,并随后再依次进行去均值与白化处理;将谐波电压与电流的测量数据作为观测信号,系统侧与用户侧谐波电流作为源信号构建独立成分分析(ICA)模型,并以数据负熵最大为目标函数,采用改进粒子群算法求解出分离矩阵,从测量数据中分离出独立分量信号;以测量数据作为因变量,解混得到的独立分量作为自变量,构建回归方程,利用最小二乘法(OLS)得到混合系数矩阵;根据混合系数与系统谐波阻抗实部与虚部间对应的线性关系得出系统谐波阻抗。与现有技术相比,本发明具有谐波补偿精确度高,计算迅速的优点。
Description
技术领域
本发明涉及谐波补偿技术领域,尤其是涉及一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法。
背景技术
伴随电力电子设备与非线性负荷的大量接入,配电网谐波污染问题愈发严重。谐波责任划分对配电网谐波治理有着极为重要的作用,而准确实现谐波责任量化区分的关键在于如何精确得到系统谐波阻抗值。当前,系统谐波阻抗的估计方法分为“干预试”与“非干预式”两大类,其中由于“非干预式”方法不会影响系统正常运行而受到广泛的研究与应用。
目前“非干预式”方法主要有:波动量法,回归法,协方差特征法,快速独立分析法(FastICA)。波动量法在假设系统谐波波动基本为零的前提下根据PCC处电压与电流波动的比值特征来估计谐波阻抗,然而在实际系统中系统背景谐波一般存在较大波动,假设条件显然不成立,随着波动的增大,该方法的估计误差显著增大。协方差特征法基于PCC处谐波电流与背景谐波电压仅存在弱相关性这一特征进行阻抗估计,该统计特征的依据是PCC处谐波电流主要由用户侧贡献且系统侧谐波电流的方差要小于用户侧,但在实际中用户侧谐波电流波动客观存在,忽略系统侧从而判定弱相关性的做法必然会随着系统侧谐波贡献的增大造成较大的误差。回归法根据PCC处谐波电流与谐波电压之间的线性关系构建方程并通过线性回归估算出阻抗值,然而当系统侧与用户侧谐波电流均存在较大波动时,PCC处谐波电压与谐波电流的线性相关度受到干扰,通过回归算法得到的估计结果无法准确反映出系统谐波阻抗。文献提出采用FastICA对PCC处测量的电压电流数据进行盲源分离,以解混出系统侧与用户侧的谐波电流,再通过混合系数间的线性关系求解得到系统谐波阻抗,该方法不再基于系统侧谐波波动较小这一假设,从而在实际工程运用中可降低估计误差。FastICA的本质是一种基于拉格朗日乘数法与牛顿迭代法相结合的寻优算法,然而牛顿迭代法仅能搜索到局部最优解,且对初值要求严苛、鲁棒性差,在复杂性高、非线性强的实际电力系统环境下,无法准确求解出系统谐波阻抗。文献提出一种改进的FastICA估计方法,该方法通过修正牛顿迭代法以提升FastICA的收敛速度,但并未提升算法的分离精度。粒子群算法(PSO)属于进化算法的一种,它不仅能够搜索到全局最优解,且拥有鲁棒性强、适应性广、易融入其它算法等优点,相较于牛顿法拥有更强大的寻优潜力。然而传统PSO存在易陷入局部最优的缺点,影响算法的寻优精度。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法。
本发明的目的可以通过以下技术方案来实现:
一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,该方法包括以下分步骤:
步骤1:对公共连接点处谐波电压与电流的测量数据进行数据清洗,并对清洗后的数据进行去均值与白化处理;
步骤2:利用经去均值与白化处理的数据,以数据负熵最大作为寻优目标,采用改进粒子群算法进行寻优,从而得到分离矩阵;
步骤3:根据分离矩阵构建回归方程并通过最小二乘法得到混合系数矩阵;
步骤4:根据混合系数矩阵与系统谐波阻抗实部与虚部之间的各自对应关系获取系统谐波阻抗并进一步对系统谐波阻抗对应谐波通过相关设备进行谐波补偿。
进一步地,所述步骤1包括以下分步骤:
步骤11:利用离散傅里叶分析方法采集公共连接点处谐波电压与电流的测量数据;
步骤12:根据拉依达准则清洗测量数据并进行初步数据判断筛选;
步骤13:对经过初步数据判断筛选的数据采用滑差法获取谐波电压与谐波电流的波动数据;
步骤14:根据谐波电压与谐波电流的波动数据构建观测矩阵;
步骤15:对观测矩阵进行去均值处理得到经去均值处理的观测矩阵;
步骤16:对经去均值处理的观测矩阵进行白化处理得到经去均值与白化处理的数据,即经去均值与白化处理的观测矩阵。
进一步地,所述步骤16中的经去均值与白化处理的观测矩阵,其描述公式为:
式中,表示经去均值与白化处理的观测矩阵,V表示白化矩阵且V=ED-0.5ET,其中E表示由数据协方差矩阵的特征向量组成的正交矩阵,D表示由数据协方差矩阵的特征值组成的对角矩阵,表示经去均值处理的观测矩阵。
进一步地,所述步骤2包括以下分步骤:
步骤21:初始化群体规模并随机初始化每个粒子位置及速度;
步骤22:利用经去均值与白化处理的数据,以负熵作为适应度函数开始寻优;
步骤23:比较粒子的适应度函数值以得到当前个体最优位置及全局最优位置;
步骤24:归一化当前个体最优位置对应适应度函数值以得到归一化处理后的当前个体最优位置对应权重;
步骤25:利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新;
步骤26:根据粒子适应度值更新当前个体最优位置及全局最优位置;
步骤27:对得到的更新后的全局最优位置进行归一化处理并根据经过归一化处理后的全局最优位置得到源信号;
步骤28:如果达到结束条件则根据源信号确定分离矩阵,否则返回步骤22按顺序再次循环迭代。
进一步地,所述步骤23中的全局最优位置的计算公式为:
式中,pg表示全局最优位置,pi表示当前个体最优位置,表示经去均值与白化处理的数据,Ng()表示负熵函数。
进一步地,所述步骤24中的归一化处理后的当前个体最优位置对应权重的计算公式为:
ωf=fitness(ωi)/||ωi||
式中,ωf表示归一化处理后的当前个体最优位置对应权重,fitness()表示适应度函数,ωi表示当前个体最优位置对应权重值。
进一步地,所述步骤25中利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新,其描述公式为:
式中,t和t+1分别表示前后对应迭代次数,vi()表示对应粒子速度,wi()表示对应粒子位置,r1()和r2()均为区间[0,1]内均匀随机数,ωp表示惯性权重因子,常数η为区间(0,2.5)内随机数,i为自然数,pg()表示对应全局最优位置,pi()表示对应当前个体最优位置,k表示当前迭代次数,MAX表示最大迭代次数,常数s1和s2均为学习因子。
进一步地,所述步骤3中的混合系数矩阵,其描述公式为:
Ki=(ITI)-1ITXi
式中,Ki表示混合系数矩阵,Xi表示观测数据矩阵,I表示根据分离矩阵得到的独立分量。
进一步地,所述步骤4中的系统谐波阻抗,其描述公式为:
式中,ZS表示系统谐波阻抗,kij表示混合系数矩阵中第i行第j列元素,i和j均为自然数。
本发明的主要原理概括为:首先,根据拉依达(Pauta)准则,对在公共连接点(PCC)处谐波电压与电流的测量数据进行数据清洗,并对清洗后的数据依次进行去均值与白化处理;其次,将谐波电压与电流的测量数据作为观测信号,系统侧与用户侧谐波电流作为源信号构建独立成分分析(ICA)模型,并以数据负熵最大为目标函数,采用改进粒子群算法求解出分离矩阵,从测量数据中分离出独立分量信号;然后,以测量数据作为因变量,解混得到的独立分量作为自变量,构建回归方程,利用最小二乘法(OLS)得到混合系数矩阵;最后,根据混合系数与系统谐波阻抗实部与虚部间对应的线性关系得出系统谐波阻抗并进一步对系统谐波阻抗对应谐波通过相关设备进行谐波补偿。
与现有技术相比,本发明具有以下优点:
(1)本发明方法中的PSO粒子群算法流程步骤包括:初始化群体规模并随机初始化每个粒子位置及速度;利用经去均值与白化处理的数据,以负熵作为适应度函数开始寻优;比较粒子的适应度函数值以得到当前个体最优位置及全局最优位置;归一化当前个体最优位置对应适应度函数值以得到归一化处理后的当前个体最优位置对应权重;利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新;根据粒子适应度值更新当前个体最优位置及全局最优位置;对得到的更新后的全局最优位置进行归一化处理并根据经过归一化处理后的全局最优位置得到源信号;如果达到结束条件则根据源信号确定分离矩阵,否则返回步骤22按顺序再次循环迭代,拥有较好的精度,在背景谐波波动干扰较大的情况下也能够较准确地得出系统谐波阻抗,进一步致使谐波补偿精确度高。
(2)本发明方法中对公共连接点处谐波电压与电流的测量数据进行数据清洗,并对清洗后的数据进行去均值与白化处理;利用经去均值与白化处理的数据,以数据负熵最大作为寻优目标,采用改进粒子群算法进行寻优,从而得到分离矩阵;根据分离矩阵构建回归方程并通过最小二乘法得到混合系数矩阵;根据混合系数矩阵与系统谐波阻抗实部与虚部之间的各自对应关系获取系统谐波阻抗并进一步对系统谐波阻抗对应谐波通过相关设备进行谐波补偿,整个方法的鲁棒性佳,稳定性强,不易崩溃出错。
附图说明
图1为本发明的方法流程图;
图2为本发明方法谐波处理的诺顿等效模型;
图3为本发明方法与其它三种方法的系统谐波阻抗幅值相对误差对比图;
图4为本发明方法与其它三种方法的系统谐波阻抗相角相对误差对比图;
图5为本发明方法与其它三种方法的系统谐波阻抗幅值结果对比图;
图6为本发明方法与其它三种方法的系统谐波阻抗相角结果对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都应属于本发明保护的范围。
实施例
如图1为本发明方法的流程示意图,具体包括以下步骤:
1)利用离散傅里叶分析采集PCC处谐波电压UPCC与谐波电流数据IPCC;
2)根据拉依达(Pauta)准则清洗测量数据,如果测量数据Xii满足则认为Xii为离群值,予以剔除。式中,是数据均值,σ是数据的标准差;
3)如图2所示为本发明方法谐波处理的诺顿等效模型,对其采用滑差法得到谐波电压与谐波电流的波动数据和
式中,d表示设定的固定时间间隔;
4)根据步骤3)得到的谐波电压与谐波电流波动量构建观测矩阵:
式中,与分别是的实、虚部,与分别是的实、虚部;
5)对步骤4)获得的观测矩阵,进行去均值处理:
式中,表示去均值后的数据,n表示采集数据点数;
6)对步骤5)去均值后的数据进行白化处理:
式中,表示经去均值与白化处理的观测矩阵,V表示白化矩阵且V=ED-0.5ET,其中E表示由数据协方差矩阵的特征向量组成的正交矩阵,D表示由数据协方差矩阵的特征值组成的对角矩阵,表示经去均值处理的观测矩阵;
7)以数据负熵最大作为寻优目标,采用改进粒子群算法进行寻优,从而得到分离矩阵W;
71):初始化群体规模NPSO并随机初始化每个粒子位置及速度;
72):利用经去均值与白化处理的数据,以负熵作为适应度函数开始寻优;
Ng(y)≈[E{g(y)}-E{g(ygauss)}]2
根据当前粒子位置计算每个粒子的适应度函数值,其中w表示位置;
73):比较粒子的适应度函数值以得到当前个体最优位置pi=wi及全局最优位置;
式中,pg表示全局最优位置,pi表示当前个体最优位置,表示经去均值与白化处理的数据,Ng()表示负熵函数。
74):归一化当前个体最优位置对应适应度函数值以得到归一化处理后的当前个体最优位置对应权重;
ωf=fitness(ωi)/||ωi||
式中,ωf表示归一化处理后的当前个体最优位置对应权重,fitness()表示适应度函数,ωi表示当前个体最优位置对应权重值。
75):利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新;
式中,t和t+1分别表示前后对应迭代次数,vi()表示对应粒子速度,wi()表示对应粒子位置,r1()和r2()均为区间[0,1]内均匀随机数,ωp表示惯性权重因子,常数η为区间(0,2.5)内随机数,i为自然数,pg()表示对应全局最优位置,pi()表示对应当前个体最优位置,k表示当前迭代次数,MAX表示最大迭代次数,常数s1和s2均为学习因子。
76):根据粒子适应度值更新当前个体最优位置及全局最优位置;
77):对得到的更新后的全局最优位置进行归一化处理并根据经过归一化处理后的全局最优位置得到源信号;
78):如果达到结束条件则根据源信号确定分离矩阵,否则返回步骤22按顺序再次循环迭代。
8)根据步骤7)得到的分离矩阵W算出独立分量I:
I=W·X
9)将观测数据矩阵Xi作为因变量,根据分离矩阵得到的独立分量I作为自变量,构建回归方程,并通过最小二乘法(OLS)得到混合系数矩阵:
Ki=(ITI)-1ITXi
式中,Ki表示混合系数矩阵,Xi表示观测数据矩阵,I表示根据分离矩阵得到的独立分量。
10)通过混合系数矩阵中的各元素计算出系统谐波阻抗的实部与虚部:
式中,ZS表示系统谐波阻抗,kij表示混合系数矩阵中第i行第j列元素,i和j均为自然数。
11)最后,系统谐波阻抗幅值与相位的结果为:
式中,RS和XS分别为系统谐波阻抗的实部和虚部,M表示总估计点数。
图3和图4为基于MATLAB R2016a仿真平台的本发明方法与其它三种方法的系统谐波阻抗幅值与相角估计相对误差的对比结果,其中,方法一为本发明方法,方法二为FastICA法,方法三为基于独立随机矢量协方差特性的阻抗估计方法,方法四为波动量法。从图3中看出,伴随k值的增大(即系统侧背景谐波干扰的增大),方法三与方法四的估计误差有着较大幅度的上升,当k=1.0时,两者的误差分别达到17.63%与36.10%。方法二的幅值估计误差始终保持在10%以内,但相比之下,本发明方法均保持在5%以内,估计误差更小,且无明显上升趋势,估计结果更稳定。图4中,随着系统侧背景谐波干扰的增大,方法三与方法四依旧存在着明显的上升趋势,而本发明方法与方法二均维持在较低水平,且本发明方法估计精度更高。图5与四6为基于本发明方法和另外三种方法得到的某低压配电网PCC处3次系统谐波阻抗幅值与相角估计结果对比图,图中发现,方法四相角与幅值的估计结果均波动剧烈;方法三相较于方法四,估计结果较平稳,但相角与幅值的估算值均在参考值以下,偏离理论值;本发明方法与方法二均在参考值附近上下波动,但本发明方法的波动幅度更小。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (9)
1.一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,该方法包括以下分步骤:
步骤1:对公共连接点处谐波电压与电流的测量数据进行数据清洗,并对清洗后的数据进行去均值与白化处理;
步骤2:利用经去均值与白化处理的数据,以数据负熵最大作为寻优目标,采用改进粒子群算法进行寻优,从而得到分离矩阵;
步骤3:根据分离矩阵构建回归方程并通过最小二乘法得到混合系数矩阵;
步骤4:根据混合系数矩阵与系统谐波阻抗实部与虚部之间的各自对应关系获取系统谐波阻抗并进一步对系统谐波阻抗对应谐波通过相关设备进行谐波补偿。
2.根据权利要求1所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤1包括以下分步骤:
步骤11:利用离散傅里叶分析方法采集公共连接点处谐波电压与电流的测量数据;
步骤12:根据拉依达准则清洗测量数据并进行初步数据判断筛选;
步骤13:对经过初步数据判断筛选的数据采用滑差法获取谐波电压与谐波电流的波动数据;
步骤14:根据谐波电压与谐波电流的波动数据构建观测矩阵;
步骤15:对观测矩阵进行去均值处理得到经去均值处理的观测矩阵;
步骤16:对经去均值处理的观测矩阵进行白化处理得到经去均值与白化处理的数据,即经去均值与白化处理的观测矩阵。
3.根据权利要求2所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤16中的经去均值与白化处理的观测矩阵,其描述公式为:
式中,表示经去均值与白化处理的观测矩阵,V表示白化矩阵且V=ED-0.5ET,其中E表示由数据协方差矩阵的特征向量组成的正交矩阵,D表示由数据协方差矩阵的特征值组成的对角矩阵,表示经去均值处理的观测矩阵。
4.根据权利要求1所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤2包括以下分步骤:
步骤21:初始化群体规模并随机初始化每个粒子位置及速度;
步骤22:利用经去均值与白化处理的数据,以负熵作为适应度函数开始寻优;
步骤23:比较粒子的适应度函数值以得到当前个体最优位置及全局最优位置;
步骤24:归一化当前个体最优位置对应适应度函数值以得到归一化处理后的当前个体最优位置对应权重;
步骤25:利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新;
步骤26:根据粒子适应度值更新当前个体最优位置及全局最优位置;
步骤27:对得到的更新后的全局最优位置进行归一化处理并根据经过归一化处理后的全局最优位置得到源信号;
步骤28:如果达到结束条件则根据源信号确定分离矩阵,否则返回步骤22按顺序再次循环迭代。
5.根据权利要求4所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤23中的全局最优位置的计算公式为:
式中,pg表示全局最优位置,pi表示当前个体最优位置,表示经去均值与白化处理的数据,Ng()表示负熵函数。
6.根据权利要求4所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤24中的归一化处理后的当前个体最优位置对应权重的计算公式为:
ωf=fitness(ωi)/||ωi||
式中,ωf表示归一化处理后的当前个体最优位置对应权重,fitness()表示适应度函数,ωi表示当前个体最优位置对应权重值。
7.根据权利要求4所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤25中利用归一化处理后的当前个体最优位置对应权重对粒子位置与速度进行更新,其描述公式为:
式中,t和t+1分别表示前后对应迭代次数,vi()表示对应粒子速度,wi()表示对应粒子位置,r1()和r2()均为区间[0,1]内均匀随机数,ωp表示惯性权重因子,常数η为区间(0,2.5)内随机数,i为自然数,pg()表示对应全局最优位置,pi()表示对应当前个体最优位置,k表示当前迭代次数,MAX表示最大迭代次数,常数s1和s2均为学习因子。
8.根据权利要求1所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤3中的混合系数矩阵,其描述公式为:
Ki=(ITI)-1ITXi
式中,Ki表示混合系数矩阵,Xi表示观测数据矩阵,I表示根据分离矩阵得到的独立分量。
9.根据权利要求1所述的一种基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法,其特征在于,所述步骤4中的系统谐波阻抗,其描述公式为:
式中,ZS表示系统谐波阻抗,kij表示混合系数矩阵中第i行第j列元素,i和j均为自然数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480723.9A CN110309545A (zh) | 2019-06-04 | 2019-06-04 | 基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910480723.9A CN110309545A (zh) | 2019-06-04 | 2019-06-04 | 基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110309545A true CN110309545A (zh) | 2019-10-08 |
Family
ID=68075027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910480723.9A Pending CN110309545A (zh) | 2019-06-04 | 2019-06-04 | 基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110309545A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441601A (zh) * | 2019-10-10 | 2019-11-12 | 四川大学 | 基于独立分量分析的多谐波源责任量化方法 |
CN111175608A (zh) * | 2020-01-19 | 2020-05-19 | 上海电力大学 | 基于加速独立分量分析的配电网谐波责任定量划分方法 |
CN112001256A (zh) * | 2020-07-24 | 2020-11-27 | 东风汽车集团有限公司 | 混合信号去工频干扰方法及系统 |
CN112464601A (zh) * | 2020-12-03 | 2021-03-09 | 电子科技大学 | 一种多谐振点电容的高频spice模型建立方法 |
CN112858782A (zh) * | 2021-01-07 | 2021-05-28 | 国网河南省电力公司电力科学研究院 | 电力系统多元用户交互影响下的谐波责任量化方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108089059A (zh) * | 2017-12-19 | 2018-05-29 | 西南交通大学 | 一种基于稀疏独立分量分析的谐波阻抗估算方法 |
CN108615068A (zh) * | 2018-03-24 | 2018-10-02 | 西安电子科技大学 | 一种混沌扰动与自适应惯性权重的粒子群优化方法 |
-
2019
- 2019-06-04 CN CN201910480723.9A patent/CN110309545A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108089059A (zh) * | 2017-12-19 | 2018-05-29 | 西南交通大学 | 一种基于稀疏独立分量分析的谐波阻抗估算方法 |
CN108615068A (zh) * | 2018-03-24 | 2018-10-02 | 西安电子科技大学 | 一种混沌扰动与自适应惯性权重的粒子群优化方法 |
Non-Patent Citations (4)
Title |
---|
ZHAO X AND YANG HG: "A new method to calculate the utility harmonic impedance based on FastICA", 《IEEE TRANS ON POWER DELIVERY》 * |
林顺富,李杨,汤波等: "基于改进FastICA及偏最小二乘法的系统谐波阻抗估计", 《电网技术》 * |
谭萍,李伟娟,李明: "基于改进粒子群算法的独立分量分析", 《兰州大学学报》 * |
陈漠: "用于最优化问题的改进粒子群优化算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441601A (zh) * | 2019-10-10 | 2019-11-12 | 四川大学 | 基于独立分量分析的多谐波源责任量化方法 |
CN110441601B (zh) * | 2019-10-10 | 2020-01-07 | 四川大学 | 基于独立分量分析的多谐波源责任量化方法 |
CN111175608A (zh) * | 2020-01-19 | 2020-05-19 | 上海电力大学 | 基于加速独立分量分析的配电网谐波责任定量划分方法 |
CN112001256A (zh) * | 2020-07-24 | 2020-11-27 | 东风汽车集团有限公司 | 混合信号去工频干扰方法及系统 |
CN112001256B (zh) * | 2020-07-24 | 2023-05-26 | 东风汽车集团有限公司 | 混合信号去工频干扰方法及系统 |
CN112464601A (zh) * | 2020-12-03 | 2021-03-09 | 电子科技大学 | 一种多谐振点电容的高频spice模型建立方法 |
CN112464601B (zh) * | 2020-12-03 | 2022-10-11 | 电子科技大学 | 一种多谐振点电容的高频spice模型建立方法 |
CN112858782A (zh) * | 2021-01-07 | 2021-05-28 | 国网河南省电力公司电力科学研究院 | 电力系统多元用户交互影响下的谐波责任量化方法和系统 |
CN112858782B (zh) * | 2021-01-07 | 2022-04-26 | 国网河南省电力公司电力科学研究院 | 电力系统多元用户交互影响下的谐波责任量化方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110309545A (zh) | 基于改进粒子群算法获取系统谐波阻抗的谐波补偿方法 | |
CN105403834B (zh) | 一种发电机动态状态估计方法 | |
CN105719023A (zh) | 一种基于混合高斯分布的风电功率实时预测误差分析方法 | |
CN107658881A (zh) | 基于戴维南等值方法的电压稳定临界点判断方法 | |
CN105974777B (zh) | 一种利用Algos和Kalman组合产生原子时标的方法 | |
CN106980044B (zh) | 一种适应风电接入的电力系统谐波电流估计方法 | |
CN107038292A (zh) | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 | |
CN104112062B (zh) | 基于插值方法的风资源分布的获取方法 | |
CN108462184B (zh) | 一种电力系统线路串补优化配置方法 | |
CN104979829B (zh) | 一种基于改进的FastICA算法的谐波电流估计方法 | |
CN108170885A (zh) | 一种配电网中多谐波源识别方法 | |
CN107025609A (zh) | 基于奇异值分解cdkf的电力系统动态状态估计方法 | |
CN109921426A (zh) | 基于cv-kde的风电并网系统概率潮流计算方法 | |
CN113468792A (zh) | 电磁暂态仿真模型的参数校正方法、装置和电子设备 | |
CN110880757B (zh) | 基于最优分区的主动配电网分布式区间状态估计方法 | |
CN115081299A (zh) | 一种基于upf的电力系统鲁棒辅助预测状态估计方法 | |
CN111709350A (zh) | 基于fcm聚类的低频振荡模态参数识别方法及系统 | |
CN113743452B (zh) | 电能质量评估方法及装置 | |
CN108054781B (zh) | 一种基于db6小波的风电汇聚效应时频分析方法 | |
CN110795841A (zh) | 一种间歇性能源出力不确定性的数学建模方法 | |
CN113919430A (zh) | 电压暂降监测判别方法及装置 | |
CN109284936A (zh) | 一种基于云图形可视化输出的电能质量评估方法 | |
CN106599541B (zh) | 一种动态电力负荷模型的结构和参数在线辨识方法 | |
CN114839586B (zh) | 基于em算法的低压台区计量装置失准计算方法 | |
CN110021931B (zh) | 一种计及模型不确定性的电力系统辅助预测状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |