CN103809023B - 基于二分搜索的电网同步谐波相量测量方法 - Google Patents
基于二分搜索的电网同步谐波相量测量方法 Download PDFInfo
- Publication number
- CN103809023B CN103809023B CN201410036038.4A CN201410036038A CN103809023B CN 103809023 B CN103809023 B CN 103809023B CN 201410036038 A CN201410036038 A CN 201410036038A CN 103809023 B CN103809023 B CN 103809023B
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- signal
- harmonic
- value
- 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.)
- Active
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种基于二分搜索的电网同步谐波相量测量方法,它首先对信号进行采集和预处理,然后利用对二分法对谐波真实频率进行搜索迭代运算,并通过比较和更新谐波相量值保证迭代过程正确,以预先设置好的误差比精度为准进行反复迭代计算,逼近谐波的真实频率,最终迭代算出更为精确的谐波相量估计值。从而在利用传统傅立叶变换传递函数幅频响应特点的同时,得到相对精确的谐波相量估计值。该方法能更逼近谐波的真实频率及相量幅值,应用于电网谐波测量中,有助于提高谐波源定位的准确性,为抑制电网谐波提供更多精确的信息,利于谐波抑制与管理。
Description
技术领域
本发明涉及一种电力行业中同步谐波相量测量方法。
背景技术
近年来,随着工业技术的不断发展,电力系统谐波污染问题日益严重。由于电弧炉、大量电力电子装置以及电力机车等非线性负荷的广泛使用,产生了大量高次谐波,造成电力系统电压、电流畸变,严重危害了电力系统的安全稳定运行,同时影响了用户的用电品质。因此,谐波问题的解决显得十分迫切,而谐波检测、谐波定位是谐波治理的关键前提。
目前,离散傅里叶变换(Discrete Fourier Transform,DFT)算法在静态条件下具有较好应用价值,已被广泛应用到现有谐波测量中。
现有的DFT算法,在非同步采样情况下存在幅值和频率误差,且随着其非同步性的增强,误差急剧增大,往往达不到实际应用的要求。当系统频率是随着系统配置参数以及状态变化而变化时,不能保证采样系统对实测信号一直保持同步采样。对于频率为fs的正弦序列,它的频谱()应该只是在fs处有离散谱。但是,在利用DFT求它的频谱做了截短,结果使信号的频谱不只是在fs处有离散谱,而是在以fs为中心的频带范围内都有谱线出现,它们可以理解为是从fs频率上“泄露”出去的;同时,对一函数实行采样,即是抽取采样点上的对应的函数值。其效果如同透过栅栏的缝隙观看外景一样,只有落在缝隙前的少数景象被看到,其余景象均被栅栏挡住而视为零,这种现象称为栅栏效应。因此应用DFT计算谐波时所引起的频谱泄露和栅栏现象可能会使得DFT算法产生较大的误差,甚至得到一个不可用的结果。DFT算法能够在一定条件和范围下对谐波进行检测,但存在一定的局限性,而谐波测量的准确度直接影响谐波定位。谐波相量测量不准确最终将导致谐波源定位不准,谐波责任难以划分,给谐波治理带来极大困难。
发明内容
本发明的目的是提供一种基于二分搜索的电网同步谐波相量测量方法,该方法能减少离散傅里叶变换算法工作在非同步采样条件的影响,保持其对稳定信号的高精度处理功能,有效减少谐波信号频率误差及提高谐波相量测量精度。
本发明为解决其技术问题,所采用的技术方案为:一种基于二分搜索的电网同步谐波相量测量方法,其步骤为:
A、数据采集与预处理
在一个测量周期内,采样器从安装于电网上的电压或电流互感器采集得到电网的电压信号或电流信号,数值信号处理器对采集得到电网的电压信号或电流信号作为输入信号,对输入信号进行离散傅里叶变换得到输入信号的频谱,并对频谱进行最小二乘拟合得到输入信号的频域信号;
B、数值信号处理器根据A步得到的输入信号的频域信号,分别计算出t1时刻输入信号的相位和t2时刻的输入信号的相角并利用公式计算出电流的基波频率估计值及其K次谐波的频率估计值其中t1、t2分别为间隔小于预计的基波周期的两个时刻,K表示谐波的次数、为大于等于2且小于等于20的整数;
C、对估计频率修正值提取
数值信号处理器根据K次谐波的频率估计值fk和A步得到的输入信号的频域信号,计算得到频率估计值fk对应的幅值Mk;同时,根据A步中的输入信号的频域信号分别计算得到fk-Δf=f1频率和fk+Δf=f3频率对应的幅值M1和M3;其中Δf为迭代频率间隙;
D、二分法迭代计算谐波相量
将频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3,以幅值最大为目标,通过A步中的频域信号进行二分法迭代,得到迭代后的频率估计值fk’及其对应的幅值M’k,即测量出K次谐波的相量值。
与现有技术相比,本发明的有益效果是:
对测量周期内通过傅里叶变换得到的频谱信号进行最小二乘拟合得到连续的频域信号,使得每一个频率点均有对应的幅值。然后通过预计的一个基波周期内的两个时刻点的相位差计算出基波频率估计值,根据该基波频率估计值得出K次谐波频率估计值,并根据上述的频域信号得出对应的幅值及其相隔迭代频率间隙的两个邻近频率的幅值。通过频域信号进行二分迭代搜索出K次谐波频率估计值附近的最大幅值对应的频率值。由于谐波在谐波频率处的幅值最大,因此该最大幅值及其对应的频率值即为测出的谐波幅值及其频率。
总之,本发明通过对频谱信号的拟合,弥补了离散傅里叶算法在非同步采样下得到的频谱信号在谐波频率点的缺失,并通过二分迭代法在拟合后的无缺失频率点的信号上搜索出最大幅值的谐波频率及其幅值。从而能够有效减少谐波信号幅值与频率误差,应用于电网谐波测量,将有助于提高谐波源定位的准确性,为抑制电网谐波提供更多精确的信息,利于谐波抑制与管理。
所述的步骤D中以幅值最大为目标,通过A步中的频域信号进行二分法迭代的具体做法是:
将计算得到的频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3进行比较:
D1、如果幅值Mk最大:
并且幅值M1为最小,则令f1=fk,M1=Mk,并将fk与f3之和的二分之一赋值给fk,f3保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
或者幅值M3为最小,则令f3=fk,M3=Mk,并将f1与fk之和的二分之一赋值给fk,f1保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
D2、如果M3幅值最大:
若f1与fk未经过迭代运算,则直接将f3与f1之和的二分之一赋值给f3;否则,f1与fk分别返回上次迭代运算前的值,并将f3与f1之和的二分之一赋值给f3;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D3、如果M1幅值最大:
若fk与f3未经过迭代运算,则直接将f1与f3之和的二分之一赋值给f1;否则,fk与f3分别返回上次迭代运算前的值,并将f1与f3之和的二分之一赋值给f1;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D4、计算误差比精度λ,其中max{M1,Mk,M3}表示M1、Mk、M3中的最大值;min{M1,Mk,M3}表示M1、Mk、M3中的最小值;
当误差比精度λ≥λ0时,λ0为设定的误差比精度的阈值,重复D1~D3的操作;否则,结束操作,得到的频率估计值fk及其对应的幅值Mk即为测出的K次谐波的相量值。
D2、D3具体迭代算法中,能够使整个迭代过程始终向着估计谐波估计频率值附近的最大幅值逼近,保证算法迅速收敛。同时二分迭代也使得算法简单,收敛迅速。从而使得本发明能够快速准确地得出符合精度要求的谐波相量幅值及其频率。
下面结合具体实施方式对本发明做进一步的详细说明。
具体实施方式
实施例
一种基于二分搜索的电网同步谐波相量测量方法,其步骤为:
A、数据采集与预处理
在一个测量周期内(通过为5~30秒),采样器从安装于电网上的电压或电流互感器采集得到电网的电压信号或电流信号,数值信号处理器对采集得到电网的电压信号或电流信号作为输入信号,对输入信号进行离散傅里叶变换得到输入信号的频谱,并对频谱进行最小二乘拟合得到输入信号的频域信号;
B、数值信号处理器根据A步得到的输入信号的频域信号,分别计算出t1时刻输入信号的相位和t2时刻的输入信号的相角并利用公式计算出电流的基波频率估计值及其K次谐波的频率估计值其中t1、t2分别为间隔小于预计的基波周期的两个时刻,K表示谐波的次数、为大于等于2且小于等于20的整数。预计的基波周期通常为19-21毫秒,t1、t2的间隔一般取15-18毫秒。
C、对估计频率修正值提取
数值信号处理器根据K次谐波的频率估计值fk和A步得到的输入信号的频域信号,计算得到频率估计值fk对应的幅值Mk;同时,根据A步中的输入信号的频域信号分别计算得到fk-Δf=f1频率和fk+Δf=f3频率对应的幅值M1和M3;其中Δf为迭代频率间隙;Δf通常取0.2~0.8Hz;
D、二分法迭代计算谐波相量
将频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3,以幅值最大为目标,通过A步中的频域信号进行二分法迭代,得到迭代后的频率估计值fk’及其对应的幅值M’k,即测量出K次谐波的相量值。
本例的步骤D中以幅值最大为目标,通过A步中的频域信号进行二分法迭代的具体做法是:
将计算得到的频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3进行比较:
D1、如果幅值Mk最大:
并且幅值M1为最小,则令f1=fk,M1=Mk,并将fk与f3之和的二分之一赋值给fk,f3保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
或者幅值M3为最小,则令f3=fk,M3=Mk,并将f1与fk之和的二分之一赋值给fk,f1保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
D2、如果M3幅值最大:
若f1与fk未经过迭代运算,则直接将f3与f1之和的二分之一赋值给f3;否则,f1与fk分别返回上次迭代运算前的值,并将f3与f1之和的二分之一赋值给f3;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D3、如果M1幅值最大:
若fk与f3未经过迭代运算,则直接将f1与f3之和的二分之一赋值给f1;否则,fk与f3分别返回上次迭代运算前的值,并将f1与f3之和的二分之一赋值给f1;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D4、计算误差比精度λ,其中max{M1,Mk,M3}表示M1、Mk、M3中的最大值;min{M1,Mk,M3}表示M1、Mk、M3中的最小值;
当误差比精度λ≥λ0时,λ0为设定的误差比精度的阈值,重复D1~D3的操作;否则,结束操作,得到的频率估计值fk及其对应的幅值Mk即为测出的K次谐波的相量值。
本发明中的谐波次数K,通常为3、5、7、9、11次,每次计算时取其中一个数字;需要测量两个以上谐波相量值时则同步进行两个以上的算法操作;
仿真实验
为研究基于二分搜索谐波相量测量算法在频率小幅波动情况下的运算精度,仿真实验设置在采样频率2400Hz,系统频率49.95Hz的情况下,观察该发明算法通过迭代运算对基波及各次谐波的幅值和频率真实值的逼近情况,并与DFT算法比较。
信号模型为式中:A1=10,A2=0.5,A3=0.1,A4=0.08,分别表示基波和3、5、7次谐波的幅值,同时加入信噪比为30dB的高斯白噪声。
实验结果:2种算法在系统频率为49.95Hz时测得的幅值估计值与真实值的比较见表1。
表1
从表1可以看出,传统DFT算法由于受到系统频率偏差的影响,测量得到的幅值估计值与真实值存在大约-0.04的偏差,偏差较大;本发明算法测得的幅值估计值与真实值之差约为-0.01。3次谐波情况虽有所波动,但本发明算法得出的测量结果基本稳定在0附近,而传统DFT算法存在大约-0.016的偏差;同理,5次谐波和7次谐波情况与基波和3次谐波类似。本发明算法与传统DFT算法比较结果如表1所示。可见本发明算法在频率小幅波动情况下,较现有的DFT算法有更高的准确度。
Claims (2)
1.一种基于二分搜索的电网同步谐波相量测量方法,其步骤为:
A、数据采集与预处理
在一个测量周期内,采样器从安装于电网上的电压或电流互感器采集得到电网的电压信号或电流信号,数值信号处理器对采集得到电网的电压信号或电流信号作为输入信号,对输入信号进行离散傅里叶变换得到输入信号的频谱,并对频谱进行最小二乘拟合得到输入信号的频域信号;
B、数值信号处理器根据A步得到的输入信号的频域信号,分别计算出t1时刻输入信号的相位和t2时刻的输入信号的相角并利用公式计算出电流的基波频率估计值及其K次谐波的频率估计值其中t1、t2分别为间隔小于预计的基波周期的两个时刻,K表示谐波的次数、为大于等于2且小于等于20的整数;
C、对估计频率修正值提取
数值信号处理器根据K次谐波的频率估计值fk和A步得到的输入信号的频域信号,计算得到频率估计值fk对应的幅值Mk;同时,根据A步中的输入信号的频域信号分别计算得到fk-Δf=f1频率和fk+Δf=f3频率对应的幅值M1和M3;其中Δf为迭代频率间隙;
D、二分法迭代计算谐波相量
将频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3,以幅值最大为目标,通过A步中的频域信号进行二分法迭代,得到迭代后的频率估计值fk’及其对应的幅值M’k,即测量出K次谐波的相量值。
2.如权利要求1所述的基于二分搜索的电网同步谐波相量测量方法,其特征在于:所述的步骤D中以幅值最大为目标,通过A步中的频域信号进行二分法迭代的具体做法是:
将计算得到的频率估计值fk对应的幅值Mk、f1频率对应的幅值M1和f3频率对应的幅值M3进行比较:
D1、如果幅值Mk最大:
并且幅值M1为最小,则令f1=fk,M1=Mk,并将fk与f3之和的二分之一赋值给fk,f3保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
或者幅值M3为最小,则令f3=fk,M3=Mk,并将f1与fk之和的二分之一赋值给fk,f1保持不变;根据更新后的频率估计值fk和A步中的频域信号重新计算其幅值,得到更新的幅值Mk;
D2、如果M3幅值最大:
若f1与fk未经过迭代运算,则直接将f3与f1之和的二分之一赋值给f3;否则,f1与fk分别返回上次迭代运算前的值,并将f3与f1之和的二分之一赋值给f3;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D3、如果M1幅值最大:
若fk与f3未经过迭代运算,则直接将f1与f3之和的二分之一赋值给f1;否则,fk与f3分别返回上次迭代运算前的值,并将f1与f3之和的二分之一赋值给f1;
将更新后的频率估计值f1、fk与f3根据A步中的频域信号重新计算其幅值,得到更新的幅值M1、Mk与M3;
D4、计算误差比精度λ,其中max{M1,Mk,M3}表示M1、Mk、M3中的最大值;min{M1,Mk,M3}表示M1、Mk、M3中的最小值;
当误差比精度λ≥λ0时,λ0为设定的误差比精度的阈值,重复D1~D3的操作;否则,结束操作,得到的频率估计值fk及其对应的幅值Mk即为测出的K次谐波的相量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036038.4A CN103809023B (zh) | 2014-01-26 | 2014-01-26 | 基于二分搜索的电网同步谐波相量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410036038.4A CN103809023B (zh) | 2014-01-26 | 2014-01-26 | 基于二分搜索的电网同步谐波相量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103809023A CN103809023A (zh) | 2014-05-21 |
CN103809023B true CN103809023B (zh) | 2016-08-24 |
Family
ID=50706089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410036038.4A Active CN103809023B (zh) | 2014-01-26 | 2014-01-26 | 基于二分搜索的电网同步谐波相量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103809023B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983851B (zh) * | 2014-05-22 | 2016-08-31 | 河海大学 | 主动配电网谐波源定位分析方法 |
CN103983847B (zh) * | 2014-06-12 | 2016-10-05 | 福州大学 | 一种同步相量测量中基于rls的自适应频率跟踪测量方法 |
CN104849548B (zh) * | 2015-06-01 | 2018-05-25 | 海南大学 | 一种电力系统瞬时频率监测方法及系统 |
CN105004920B (zh) * | 2015-07-10 | 2017-11-17 | 国网天津市电力公司 | 傅里叶修正系数频率测量方法 |
CN106501602B (zh) * | 2016-09-28 | 2018-10-30 | 西南交通大学 | 一种基于滑窗频谱分离的基波参数测量方法 |
CN108519512B (zh) * | 2018-03-23 | 2020-05-22 | 深圳市计量质量检测研究院 | 用于高精度谐波测量的方法和系统 |
CN112433088A (zh) * | 2019-08-26 | 2021-03-02 | 中车大连电力牵引研发中心有限公司 | 谐波识别方法、装置、电子设备及存储介质 |
CN112557781B (zh) * | 2020-10-21 | 2024-07-26 | 广东电网有限责任公司广州供电局 | 一种适用于校准器的宽频域信号测量方法 |
CN112557751B (zh) * | 2020-12-03 | 2023-07-18 | 东南大学 | 一种基于dft迭代法的谐波参数估计方法 |
CN113533846B (zh) * | 2021-07-13 | 2024-08-06 | 国网江苏省电力有限公司营销服务中心 | 一种用电负荷数据高速同步采集计算方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236048A (zh) * | 2011-06-16 | 2011-11-09 | 国网电力科学研究院 | 一种电力系统相量频率测量方法 |
CN102401866A (zh) * | 2011-12-13 | 2012-04-04 | 广东电网公司佛山供电局 | 三调谐直流滤波器失谐故障元件检测方法 |
CN102749521A (zh) * | 2012-07-18 | 2012-10-24 | 华北电力大学(保定) | 一种电力系统谐波阻抗计算方法 |
CN102998535A (zh) * | 2012-12-10 | 2013-03-27 | 华北电力大学(保定) | 一种基于极大似然估计理论的系统谐波阻抗计算方法 |
CN103091554A (zh) * | 2013-03-01 | 2013-05-08 | 山东电力集团公司电力科学研究院 | 线路末端同步测量信号多母线微电网稳态谐波分析方法 |
-
2014
- 2014-01-26 CN CN201410036038.4A patent/CN103809023B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102236048A (zh) * | 2011-06-16 | 2011-11-09 | 国网电力科学研究院 | 一种电力系统相量频率测量方法 |
CN102401866A (zh) * | 2011-12-13 | 2012-04-04 | 广东电网公司佛山供电局 | 三调谐直流滤波器失谐故障元件检测方法 |
CN102749521A (zh) * | 2012-07-18 | 2012-10-24 | 华北电力大学(保定) | 一种电力系统谐波阻抗计算方法 |
CN102998535A (zh) * | 2012-12-10 | 2013-03-27 | 华北电力大学(保定) | 一种基于极大似然估计理论的系统谐波阻抗计算方法 |
CN103091554A (zh) * | 2013-03-01 | 2013-05-08 | 山东电力集团公司电力科学研究院 | 线路末端同步测量信号多母线微电网稳态谐波分析方法 |
Non-Patent Citations (5)
Title |
---|
MUSIC法的有限域单向快速频率搜索;蒋毅等;《电子科技大学学报》;20080330;第37卷(第02期);第241-243页 * |
同步电机谐波相量间的基本关系及其应用;张小平;《中国电机工程学报》;19950330;第15卷(第02期);第110-117页 * |
基于Rife-Vincent窗的高准确度电力谐波相量计算方法;曾博等;《电工技术学报》;20090826;第24卷(第08期);第154-158页 * |
改进的数据采集系统性能测试的正弦拟合法;朱仕银等;《北京理工大学学报》;20001230;第20卷(第06期);第757-761页 * |
特高压母线故障暂态高频电流特征和母线保护应对措施;吴崇昊等;《电力自动化设备》;20080410;第28卷(第04期);第10-15页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103809023A (zh) | 2014-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103809023B (zh) | 基于二分搜索的电网同步谐波相量测量方法 | |
Su et al. | Power harmonic and interharmonic detection method in renewable power based on Nuttall double‐window all‐phase FFT algorithm | |
Xia et al. | A complex least squares enhanced smart DFT technique for power system frequency estimation | |
CN103399213B (zh) | 一种基于基波频率修正的高压设备介质损耗因数计算方法 | |
Pazoki | A new DC-offset removal method for distance-relaying application using intrinsic time-scale decomposition | |
CN110333389A (zh) | 基于插值dft的正弦信号频率估计方法 | |
CN106018958B (zh) | 孤网变频系统中二次侧电流电压频率跟踪方法 | |
CN104049144A (zh) | 一种滤除衰减直流分量的同步相量测量的实现方法 | |
CN102818930B (zh) | 一种高精度快速计算电力谐波参数的方法 | |
Wang et al. | Adaptive robust cubature Kalman filter for power system dynamic state estimation against outliers | |
CN104316768B (zh) | 一种三相不平衡扰动源定位的负序阻抗参数估算方法 | |
CN104897960A (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及系统 | |
WO2018227485A1 (en) | Method, system and apparatus for fault detection | |
CN106771586B (zh) | 一种直流控制保护板卡的回路信号分析方法及装置 | |
Kumar et al. | A novel sub-cycle-based method for estimation of decaying DC component and fundamental phasor | |
CN104849552A (zh) | 一种基于自适应陷波滤波器的谐波提取方法 | |
Godse et al. | Real‐time digital filtering algorithm for elimination of the decaying DC component using mathematical morphology | |
CN104407197B (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
Tao et al. | A robust parametric method for power harmonic estimation based on M-estimators | |
Bertocco et al. | Frequency tracking for efficient phasor measurement based on a CSTFM model | |
Dutra et al. | High-precision frequency estimation of real sinusoids with reduced computational complexity using a model-based matched-spectrum approach | |
Serov et al. | Reseach of influency of mains frequency deviation on voltage spectrum measurement error by DFT method | |
CN105486991A (zh) | 一种局部放电脉冲提取方法 | |
Barros et al. | A Fast Gridless algorithm for harmonics and interharmonics estimation | |
Amin et al. | Kalman filter estimation of impedance parameters for medium transmission line |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |