CN109506773A - 基于数字解码的高精度外差激光干涉振动校准方法 - Google Patents
基于数字解码的高精度外差激光干涉振动校准方法 Download PDFInfo
- Publication number
- CN109506773A CN109506773A CN201811633851.4A CN201811633851A CN109506773A CN 109506773 A CN109506773 A CN 109506773A CN 201811633851 A CN201811633851 A CN 201811633851A CN 109506773 A CN109506773 A CN 109506773A
- Authority
- CN
- China
- Prior art keywords
- signal
- vibration
- output voltage
- pdif
- voltage signal
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000005284 excitation Effects 0.000 claims abstract description 46
- 238000005070 sampling Methods 0.000 claims abstract description 25
- 230000001133 acceleration Effects 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 238000002955 isolation Methods 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了基于数字解码的高精度外差激光干涉振动校准方法,该方法包括:以固定采样率同步采集外差激光干涉仪输出的S/PDIF信号与被校振动传感器及测量仪的输出电压信号;通过对采集的S/PDIF信号进行解码,并利用混合插值算法剔除解码信号的噪声点;再对输出电压信号与解码的振动激励信号对齐处理,保证被校振动传感器及测量仪的校准精度;最后基于正弦逼近法(SAM)拟合解码的振动激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励加速度峰值与初相及其输出电压信号的峰值与初相,以实现振动校准。相比于现有的外差激光干涉振动校准方法,本方法在保证校准精度的前提下,实现稳定、可靠、快速的校准。
Description
技术领域
本发明涉及一种基于数字解码的高精度外差激光干涉振动校准方法,属于激光干涉技术领域。
背景技术
振动传感器及测量仪被广泛用于多个领域的实时监控与振动参数测量,如桥梁建筑、地震、汽车、航空航天等。为了保证振动传感器及测量仪所得测量结果的可靠性与准确性,需要定期对振动传感器及测量仪进行校准。
常规的外差激光干涉振动校准方法有基于Nyquist采样定理的NS采集方式校准方法、使用模拟混频器与低通滤波器的MLPFS采集方式校准方法、基于带通采样的BPS采集方式校准方法。NS采集方式的校准方法具有所需采样频率高、采样数据量大、采集信号测量精度高等特点;基于MLPFS采集方式的校准方法具有所需采样率低,采集数据量小,但引入外部模拟器件导致校准精度有限。基于BPS采集方式的校准方法具有所需采样率低,采集数据量小、采集信号测量精度高,不需要外部模拟器件而不存在相位延时,但是需要计算最佳采样频率,增加了校准过程的复杂性。
因此,针对目前外差激光干涉振动校准方法存在的所需采样率高且处理数据量大、存在相位延时、灵敏度相位校准精度有限、校准过程复杂等缺点,本发明提出一种所需采样率低、不需要外部模拟器件、测量精度高、校准过程简单的外差激光干涉振动校准方法。
发明内容
本发明的目的在于针对目前的外差激光干涉振动校准方法存在所需采样率高且处理数据量大、需要外部模拟器件而引入相位延时导致校准精度有限、校准过程复杂等不足,提供了一种基于数字解码的高精度外差激光干涉振动校准方法。
为实现上述目的,本发明采用的技术方案为基于数字解码的高精度外差激光干涉振动校准方法,该方法包括:以固定采样率同步采集外差激光干涉仪输出的S/PDIF信号与被校振动传感器及测量仪的输出电压信号,对采集到S/PDIF信号进行解码处理,对解码过程中出现的噪声点利用混合插值算法处理,得到振动激励信号。
引入振动激励信号与输出电压信号的对齐处理过程,保证两信号的对齐,提高校准精度。
基于SAM拟合对振动激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励峰值与初相及其输出电压信号的峰值与初相,以实现振动校准。
基于数字解码的高精度外差激光干涉仪振动校准方法包括以下步骤,
S1:以固定采样率同步采集外差激光干涉仪的S/PDIF信号与被校振动传感器及测量仪输出电压信号;
S2:对S/PDIF信号进行解码,获取被校振动传感器及测量仪的输入激励信号,对解码的激励信号使用混合插值算法剔除其噪声点;
S3:对采集到的输出电压信号与解码的激励信号进行对齐处理,以保证被校振动传感器及测量仪的校准精度;
S4:基于SAM拟合激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励峰值与初相及其输出电压信号的峰值与初相。
利用混合插值算法剔除解码激励信号的噪声点,具体包括,
(1)线性插值;
当振动频率小于1kHz时,根据S/PDIF信号编码规则,在固定采样率下12.288MHz下,每个振动周期的采样点数大于96个,连续3个点之间的关系可近似为直线,使用线性插值算法剔除噪声点,计算如下:
其中,Lj(x)为连接(xj,xj+1)的线段,i、j分别为采样点序号。
(2)样条插值;
当振动频率大于等于1kHz时,利用趋于正弦分布的样条插值算法剔除噪声点,以提高校准精度,在区间[xi,xi+1],插值y计算如下:
y=Ayi+1+Byi+1+Cy″i+Dy″i+1 (2)
其中,
B=1-A
对采集到的输出电压信号与解码的振动激励信号进行对齐处理的方法:
通过判断S/PDIF信号中第一个左声道头码的位置,计算出第一个左声道头码之前的无效位位数,相应地剔除采集到的输出电压信号相同的位数,即实现输出电压信号与振动激励信号的对齐。
基于SAM拟合振动激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励加速度峰值与初相及输出电压信号的峰值与初相。
本发明外差激光干涉振动校准方法具有如下优势:
(1)本发明方法具有所需采样率低、处理数据量小、实时性好、校准精度高的优点;
(2)本发明方法只需以固定采样率采集,无需其他参数设置与计算,校准过程简单;
(3)本发明方法不需要使用任何外部模拟器件,固不存在外部模拟器件引入的相位延时;
(4)本发明方法基于SAM对振动激励信号与传感器及测量仪输出电压信号拟合,保证振动校准的准确性。
附图说明
图1为本发明方法具体实施实例安装装置示意图;
图2为基于数字解码的高精度外差激光干涉振动校准方法流程图;
图3为外差激光干涉仪的S/PDIF信号解码与利用混合插值方法剔除解码激励信号噪声点流程图;
图4为振动激励信号与振动传感器及测量仪输出电压信号对齐方法流程图;
图5为本发明方法具体实施实例的相位灵敏度校准结果;
图6为本发明方法具体实施实例的幅值灵敏度校准结果;
具体实施方式
针对目前外差激光干涉振动校准方法存在的所需采样率高且处理数据量大、存在相位延时、灵敏度相位校准精度有限、校准过程复杂等缺点,本发明提出一种所需采样率低、不需要外部模拟器件、测量精度高、校准过程简单的外差激光干涉振动校准方法,下面结合附图和具体的实施实例对本发明做出详细描述。
参考图1位本发明方法的实施实例装置示意图,该装置主要包括振动台1、振动台工作台面2、被校传感器及测量仪3、反光纸4、外差激光干涉仪5、隔振台6、数据传输线7、信号采集设备8、信号处理及显示设备9;隔振台6用于固定外差激光干涉仪5并使激光束垂直于振动台1的工作台面2;被校振动传感器及测量仪3固定在振动台工作台面2;反光纸4用于反射外差激光干涉仪5发出的激光光束;被校振动传感器及测量仪3、外差激光干涉仪5、信号采集设备8、信号处理及显示设备9通过数据传输线7连接;数据传输线7用于实时传输被校传感器2、外差激光干涉仪2、信号采集设备8信号;信号处理及显示设备9用于处理S/PDIF信号与被校振动传感器及测量仪的输出电压信号,并显示波形与校准结果;
参考图2为基于数字解码的高精度外差激光干涉振动校准方法流程图。本发明方法主要包括以下步骤:
S20:以固定采样率同步采集外差激光干涉仪输出的S/PDIF信号与被校振动传感器及测量仪的输出电压信号,对采集的S/PDIF信号进行解码并使用插值算法剔除解码的激励信号噪声点;
S40:对采集的输出电压信号与解码的振动激励信号进行对齐处理,以提高校准精度;
S60:基于SAM拟合振动激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励加速度峰值与初相及其输出电压信号的峰值与初相;
S80:利用不同频率下的激励信号峰值与初相及输出电压信号计算响应灵敏度,获得振动传感器及测量仪的频率特性。
参考图3为S/PDIF信号解码与利用混合插值方法剔除噪声点流程图。本发明外差激光干涉仪的S/PDIF信号解码与使用混合插值方法剔除噪声点包括如下步骤:
步骤S21:对采集到的S/PDIF信号向上取整,得到规范的二进制数S/PDIF信号;
步骤S22:依据头码的编码规则设置模板,通过模板获取S/PDIF信号的第一个左声道头码;
步骤S23:以第一个左声道头码为起始位置,获取S/PDIF信号中所有的完整左声道;
步骤S24:判断每个左声道的有效位元是否为0,若为0,此左声道有效,跳至步骤S25,否则此左声道无效,跳至步骤S27;
步骤S25:判断每个左声道是否符合S/PDIF编码规则,符合则跳至步骤S26,否则跳至步骤S27;
步骤S27:无效的左声道与不符合编码规则的左声道为噪声点,使用十进制数0替换;
步骤S28:振动频率小于1kHz,跳至步骤S29,振动频率大于等于1kHz,跳至步骤S30;
步骤S29:利用线性插值算法计算噪声点位置的实际值;
步骤S30:使用样条插值算法计算噪声点位置的实际值;
步骤S31:将插值算法获取的值替换原噪声点位置的0值;
步骤S32:输出被校振动传感器及测量仪的激励信号。
参考图4为输出电压信号与振动激励信号对齐处理方法流程图。本发明对采集到输出电压信号与振动激励信号对齐处理包括如下步骤:
步骤S41:以获取S/PDIF信号的第一个左声道头码为起始位置;
步骤S42:计算第一个左声道头码起始位置之前无效的S/PDIF编码位数;
步骤S43:剔除与无效的S/PDIF编码位数相同的输出电压信号位数。
参考图5-6为本发明方法的具体实施实例的传感器振动校准结果图。本实施实例装置的具体参数为:德国保利泰克OFV-5000外差式激光干涉仪,振动频率为5Hz-20kHz的PCB 396C11振动台,凌华9846数据采集卡,被校传感器为PCB 396C11内置传感器,隔振台为一级光学隔振平台。本次校准使用图1对应的安装方式,在正弦激励输入下的振动台工作台面实现传感器校准,以每个频点测量10次的平均值得出校准结果。
上述详细描述为本发明方法的一个具体实施实例,其并非用以限定本发明的应用范围。本领域相关技术人员可以在本发明的基础上可做出一系列的优化与改进、等同修改等。因此本发明的保护范围应由所附权利要求来限定。
Claims (8)
1.基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:所述方法包括以下步骤,
S1:以固定采样率同步采集外差激光干涉仪的S/PDIF信号与被校振动传感器及测量仪输出电压信号;
S2:对S/PDIF信号进行解码,获取被校振动传感器及测量仪的输入激励信号,对解码的激励信号使用混合插值算法剔除其噪声点;
S3:对采集的输出电压信号与解码的激励信号进行对齐处理,以保证被校振动传感器及测量仪的校准精度;
S4:基于SAM拟合激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励峰值与初相及其输出电压信号的峰值与初相。
2.根据权利要求1所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:利用混合插值算法剔除解码激励信号的噪声点,具体包括,
(1)线性插值;
当振动频率小于1kHz时,根据S/PDIF信号编码规则,在固定采样率12.288MHz下,每个振动周期的采样点数大于96个,连续3个点之间的关系近似为直线,使用线性插值算法剔除噪声点,计算如下:
其中,Lj(x)为连接(xj,xj+1)的线段,i、j分别为采样点序号;
(2)样条插值;
当振动频率大于等于1kHz时,利用趋于正弦分布的样条插值算法剔除噪声点,以提高校准精度,在区间[xi,xi+1],插值y计算如下:
y=Ayi+1+Byi+1+Cy”i+Dy”i+1 (2)
其中,
B=1-A
3.根据权利要求1所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:
对采集的输出电压信号与解码的激励信号进行对齐处理:
通过判断S/PDIF信号中第一个左声道头码的位置,计算出第一个左声道头码之前的无效位位数,相应地剔除采集到的输出电压信号相同的位数,实现输出电压信号与振动激励信号的对齐。
4.根据权利要求1所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:基于SAM拟合解码的激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励加速度峰值与初相及输出电压信号的峰值与初相。
5.根据权利要求1所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:
振动校准方法包括以下步骤:
S20:以固定采样率同步采集外差激光干涉仪输出的S/PDIF信号与被校振动传感器及测量仪的输出电压信号,对采集的S/PDIF信号进行解码并使用插值算法剔除解码的激励信号噪声点;
S40:对采集的输出电压信号与解码的振动激励信号进行对齐处理,以提高校准精度;
S60:基于SAM拟合振动激励信号与输出电压信号,获得被校振动传感器及测量仪的输入激励加速度峰值与初相及其输出电压信号的峰值与初相;
S80:利用不同频率下的激励信号峰值与初相及输出电压信号计算响应灵敏度,获得振动传感器及测量仪的频率特性。
6.根据权利要求5所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:
对采集的S/PDIF信号解码与剔除噪声点包括如下步骤:
步骤S21:对采集到的S/PDIF信号向上取整,得到规范的二进制数S/PDIF信号;
步骤S22:依据头码的编码规则设置模板,通过模板获取S/PDIF信号的第一个左声道头码;
步骤S23:以第一个左声道头码为起始位置,获取S/PDIF信号中所有的完整左声道;
步骤S24:判断每个左声道的有效位元是否为0,若为0,此左声道有效,跳至步骤S25,否则此左声道无效,跳至步骤S27;
步骤S25:判断每个左声道是否符合S/PDIF编码规则,符合则跳至步骤S26,否则跳至步骤S27;
步骤S27:无效的左声道与不符合编码规则的左声道为噪声点,使用十进制数0替换;
步骤S28:振动频率小于1kHz,跳至步骤S29,振动频率大于等于1kHz,跳至步骤S30;
步骤S29:利用线性插值算法计算噪声点位置的实际值;
步骤S30:使用样条插值算法计算噪声点位置的实际值;
步骤S31:将插值算法获取的值替换原噪声点位置的0值;
步骤S32:输出被校振动传感器及测量仪的激励信号。
7.根据权利要求6所述的基于数字解码的高精度外差激光干涉振动校准方法,其特征在于:
对采集的输出电压信号与解码的振动激励信号对齐处理包括如下步骤:
步骤S41:以获取S/PDIF信号的第一个左声道头码为起始位置;
步骤S42:计算第一个左声道头码起始位置之前无效的S/PDIF编码位数;
步骤S43:剔除与无效的S/PDIF编码位数相同的输出电压信号位数。
8.利用权利要求1所述方法进行的基于数字解码的高精度外差激光干涉振动校准装置,其特征在于:
该装置包括振动台(1)、振动台工作台面(2)、被校传感器(3)、反光纸(4)、外差激光干涉仪(5)、隔振台(6)、数据传输线(7)、信号采集设备(8)、信号处理及显示设备(9);隔振台(6)用于固定外差激光干涉仪(5)并使激光束垂直于振动台(1)的工作台面(2);被校振动传感器及测量仪(3)固定于振动台的工作台面(2);反光纸(4)用于反射外差激光干涉仪(5)发出的激光光束;被校振动传感器及测量仪(3)、外差激光干涉仪(5)、信号采集设备(8)、信号处理及显示设备(9)通过数据传输线(7)连接;数据传输线(7)用于实时传输被校传感器(2)、外差激光干涉仪(2)、信号采集设备(8)信号;信号处理及显示设备(9)用于处理S/PDIF信号与被校振动传感器及测量仪的输出电压信号,并显示波形与校准结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633851.4A CN109506773B (zh) | 2018-12-29 | 2018-12-29 | 基于数字解码的高精度外差激光干涉振动校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633851.4A CN109506773B (zh) | 2018-12-29 | 2018-12-29 | 基于数字解码的高精度外差激光干涉振动校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109506773A true CN109506773A (zh) | 2019-03-22 |
CN109506773B CN109506773B (zh) | 2020-01-24 |
Family
ID=65756876
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811633851.4A Active CN109506773B (zh) | 2018-12-29 | 2018-12-29 | 基于数字解码的高精度外差激光干涉振动校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109506773B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098026A (zh) * | 2020-09-08 | 2020-12-18 | 杭州亿恒科技有限公司 | 一种噪声及环路检测加速方法及其系统 |
CN112556827A (zh) * | 2020-10-22 | 2021-03-26 | 中国计量科学研究院 | 一种基于激光干涉法的高加速度振动校准方法及装置 |
CN112649087A (zh) * | 2020-11-04 | 2021-04-13 | 北京航天计量测试技术研究所 | 基于声振解耦的噪声传感器振动灵敏度校准方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11351958A (ja) * | 1998-06-11 | 1999-12-24 | Akashi Corp | 振動ピックアップの校正装置 |
CN104501940A (zh) * | 2014-12-17 | 2015-04-08 | 中国计量科学研究院 | 一种外差激光的信号解调方法及其系统 |
US20170122798A1 (en) * | 2014-04-10 | 2017-05-04 | Korea Research Institute Of Standards And Science | Method for measuring vibration displacement using state variation principle |
CN106895904A (zh) * | 2016-12-16 | 2017-06-27 | 中国计量科学研究院 | 一种高精度的外差式激光振动校准灵敏度相位测量方法 |
CN107356321A (zh) * | 2017-08-25 | 2017-11-17 | 浙江省计量科学研究院 | 三维光电振动测量仪电路计量校准方法及其校准装置 |
-
2018
- 2018-12-29 CN CN201811633851.4A patent/CN109506773B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11351958A (ja) * | 1998-06-11 | 1999-12-24 | Akashi Corp | 振動ピックアップの校正装置 |
US20170122798A1 (en) * | 2014-04-10 | 2017-05-04 | Korea Research Institute Of Standards And Science | Method for measuring vibration displacement using state variation principle |
CN104501940A (zh) * | 2014-12-17 | 2015-04-08 | 中国计量科学研究院 | 一种外差激光的信号解调方法及其系统 |
CN106895904A (zh) * | 2016-12-16 | 2017-06-27 | 中国计量科学研究院 | 一种高精度的外差式激光振动校准灵敏度相位测量方法 |
CN107356321A (zh) * | 2017-08-25 | 2017-11-17 | 浙江省计量科学研究院 | 三维光电振动测量仪电路计量校准方法及其校准装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112098026A (zh) * | 2020-09-08 | 2020-12-18 | 杭州亿恒科技有限公司 | 一种噪声及环路检测加速方法及其系统 |
CN112556827A (zh) * | 2020-10-22 | 2021-03-26 | 中国计量科学研究院 | 一种基于激光干涉法的高加速度振动校准方法及装置 |
CN112649087A (zh) * | 2020-11-04 | 2021-04-13 | 北京航天计量测试技术研究所 | 基于声振解耦的噪声传感器振动灵敏度校准方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109506773B (zh) | 2020-01-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10634554B2 (en) | Integrated vibration measurement and analysis system | |
CN109506773B (zh) | 基于数字解码的高精度外差激光干涉振动校准方法 | |
CN106895904B (zh) | 一种高精度的外差式激光振动校准灵敏度相位测量方法 | |
CN211348423U (zh) | 一种高频信号测量装置 | |
CN106767745A (zh) | 一种光电传感器测角系统的信号处理方法 | |
CN102624456B (zh) | 一种光纤干涉型传感信号解调的系统及方法 | |
CN108132448B (zh) | 用于相位发生器相位零点校准的装置及方法 | |
CN204255494U (zh) | 桥梁振动监测装置 | |
CN109490794A (zh) | 一种电源纹波测试电路和测试方法 | |
CN104880597A (zh) | 一种基于可编程逻辑电路的微弱信号检测方法 | |
CN104678170B (zh) | 一种基于谐波分析仪的电力谐波分析方法和谐波分析仪 | |
CN102095940A (zh) | 一种介质损耗角的测量方法和装置 | |
CN104734715A (zh) | 一种提高a/d转换器分辨率的方法 | |
CN109541309B (zh) | 一种频谱分析仪及其信号处理方法 | |
CN101660925A (zh) | 数字编码器解码过程中干扰信号的处理方法 | |
CN103197153A (zh) | 一种基于矢量三角形的电容电感参数测量电路及其测量方法 | |
CN117213619A (zh) | 一种宽频带振动校准外差激光干涉法及其测量装置 | |
CN108414001B (zh) | 非均匀采样正弦波形失真度的确定方法 | |
CN106595841B (zh) | 基于带通采样的外差式激光测振方法 | |
CN104215597B (zh) | 甲烷浓度参考信息获取方法及装置 | |
Li et al. | A design of inductosyn angle measurement system based on AD2S80A and FPGA | |
CN112857349A (zh) | 一种应用于液浮陀螺仪的高精度信号采集系统及方法 | |
JP2009128334A (ja) | インピーダンス測定装置、及び検出方法 | |
CN219349111U (zh) | 一种双通道电压信号源超低误差高精度测量系统 | |
CN110380415A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |