CN106597095B - 一种双cpu等间隔采样的频率跟踪实现方法 - Google Patents
一种双cpu等间隔采样的频率跟踪实现方法 Download PDFInfo
- Publication number
- CN106597095B CN106597095B CN201611199964.9A CN201611199964A CN106597095B CN 106597095 B CN106597095 B CN 106597095B CN 201611199964 A CN201611199964 A CN 201611199964A CN 106597095 B CN106597095 B CN 106597095B
- Authority
- CN
- China
- Prior art keywords
- frequency
- cpu
- resampling
- sampling
- module
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种双采样CPU在等间隔采样时实现频率跟踪的方法,通过配置信息中测频通道接入的位置确定两个CPU中哪一个作为测频CPU,哪一个作为跟随CPU;测频CPU在测频模块中计算出模拟量的频率,通过数据交互模块将计算出的频率传输给跟随CPU,在重采样模块中使用上次算出的频率得到重采样间隔,对原始采样值进行重采样;跟随CPU从数据交互模块中读取频率并计算重采样间隔,在重采样模块中对原始采样值进行重采样。本发明解决在数字化装置中采用双采样CPU等间隔采样时,两个CPU都需要接入测频通道才能实现频率跟踪的问题。
Description
技术领域
本发明属于电力系统自动化技术领域,具体涉及高压继电保护装置和智能测控装置在数据采集中的频率跟踪方法。
背景技术
随着国内电力需求的日益增大,电力技术水平不断提高,电力保障设备的使用显得尤为重要。而高压继电保护和测控装置数据量大、要求采样精度高,尤其是在一些数据量较大的智能化装置上,通常采用多CPU分布采集。在电网频率波动时,要保证采集数据的同步以及准确性。
传统保护测控装置由本身进行数据采集,一般是根据测频计算出采样间隔,调整采样中断,从而实现的频率跟踪,是不等间隔采样。数字化保护、测控装置的采样数据一般来自合并单元,采样频率高,采样间隔恒定。合并单元本身不做频率跟踪,这就要求使用采样数据的后端设备自行做频率跟踪。目前,现有保护测控装置的等间隔采样频率跟踪方法都为单CPU实现,如果有2个采样CPU的话要将测频通道同时接入2个CPU,然后各自实现频率跟踪;尚无只接入一个测频通道实现双CPU等间隔采样的频率跟踪的方法。
发明内容
针对现有技术中存在的问题,本发明的目的是提供一种双采样CPU在等间隔采样时实现频率跟踪的方法,用以解决在数字化装置中采用双采样CPU时,两个CPU都需要接入测频通道的问题。
为实现上述目的,本发明采用的技术方案为:一种双CPU等间隔采样的频率跟踪实现方法,包括:通过配置信息中测频通道接入的位置确定两个CPU中哪一个作为测频CPU,哪一个作为跟随CPU;测频CPU在测频模块中计算出本次频率和重采样间隔,在重采样模块中使用上次算出的重采样间隔对原始采样值进行重采样,通过数据交互模块将计算出的频率和重采样间隔传输给跟随CPU;跟随CPU从数据交互模块中读取重采样间隔,在重采样模块中对原始采样值进行重采样。以上所述的所有工作在一个中断中完成。
重采样的值具有采样编号,用于后端同步判断。初次采样由上端控制,开始前对两个CPU的重采样编号清零,确保两个采样编号同步。
所述的测频模块,具体方法如下:①根据近一周波内原始采样点数据计算出傅里叶实部值,②计算傅里叶实部值过零点的时间差,③计算本周波频率。
所述步骤①中的傅里叶实部值得计算方法公式为:
其中N为原始采样点每周波采样点数,i为当前采样点序号,h=i-(N-1)为一周波前的采样点序号,Smpk为第k点的采样值;
所述步骤②中计算实部过零点时间差方法,公式为计算频率条件,当满足该条件时,记录下Rei、Rei-1、Rei-2;并记录一个周波内不满足该公式条件的采样点个数M。设Rei-2时刻Ti-2=0,则Rei-1时刻Ti-1=TN、Rei时刻Ti=2*TN(TN为原始采样点采样间隔时间),利用拉格朗日二阶插值公式计算出实部值为零时对应的采样时刻:
则有:Tpos=2*TN-Tzero;Tnag=Tzero-TN;
傅里叶实部过零点的时间差:Tdelta=M*TN+T’pos+Tnag
其中,T’pos为上次计算得到的Tpos值。
所述步骤②完成过零点时间差计算后将M赋为0。
所述步骤③中的本周波的频率f’F=1s/Tdelta。
所述的测频模块在每收到一个原始采样点时调用一次。
所述的数据交互模块:两个CPU之间的交互数据要求实时性非常高,每次采样中断都需要进行数据交换,本发明采用双口RAM,对两个CPU之间的数据传输提供支撑,且外部中断信号源确保两个CPU同时进入采样中断。
当刚进入采样中断时,两个CPU都从双口RAM的指定地址读取本次中断中应使用的频率值,该值由上次中断结束时测频CPU写入,初次上电时该频率值默认50Hz。为确保两个CPU对双口RAM的读写操作不产生冲突,本发明采用了分时读写和读写控制两种措施:
数据交互模块中所述的分时读写措施:测频CPU在中断结束前执行写操作,程序内部保证中断执行时间不小于50us(该时间大于写RAM的执行时间即可);两颗CPU在刚进入中断时执行读操作。这样通过时间控制避免双口RAM的读写冲突;
数据交互模块中所述的读写控制措施:写数据之前,先将写操作标志WR_BUSY置1,写结束最后将WR_BUSY清0。当读操作时判断WR_BUSY,若为0,直接读取;若为1,则等待10us后再读,若等待时间超过50us,认为写操作超时,放弃读取,采用上次中断中使用的值进行运算并告警。
所述的重采样模块:根据原始采样点,通过拉格朗日二阶插值法得到新的采样值。①记录并存储原始采样点的采样值Sn和采样时刻Tn;②在数据交换模块中得到新的频率fF,算出重采样间隔TF,由上一个重采样点的时刻T’Re得到新的重采样时刻TRe;③如果TRe≤Tn,在原始采样点缓冲区找出TRe时刻之前的两个采样点(采样值S0、S1,采样时刻T0、T1)和之后的一个采样点(采样值S2,采样时刻T2),根据拉格朗日二阶插值法得到TRe时刻的采样值SRe;④判断下一个重采样点的时刻是否在Tn之前,如果是,TRe加上TN后重复上述步骤③。⑤将重采样后的数据和本次中断重采样点数目发送给后端;重采样后的点具有采样编号,便于后端对两个CPU的重采样数据进行同步识别。
所述步骤②中重采样间隔计算公式:TF=F/fF,其中F为每周波重采样点数;新的重采样时刻计算公式为:TRe=T′Re+TF。
所述步骤③中采样值SRe的计算公式为:
所述步骤④在每次进入重采样模块后只执行一次,确保每个中断最多只做两次重采样插值,保证中断执行时间不超限。
所述步骤⑤中的采样编号为0~(50*F-1)循环,不受频率的变化影响。
本发明的有益效果:当保护测控装置接收采样的数据量较大时,对CPU的处理能力要求很高,尤其是数字化装置,还要对网络接口做安全保护机制。这样,对处理器的性能要求就很高,如果使用多CPU进行采样的方案,又面临着频率跟踪中测频通道需要分别接入两个CPU的问题。本发明使用两个CPU,可以降低大量的处理数据对CPU的性能要求,而又顺利解决了数字化装置中两个CPU间实现频率跟综功能,测频通道不好分两接入的问题。
附图说明
图1是本发明实施例的工作流程图。
图2是本发明实施例的测频模块的过零点时差计算示意图。
图3是本发明实施例的重采样示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步描述。
本实施例的双CPU等间隔采样的频率跟踪实现方法,两个CPU分为测频CPU和跟随CPU,其步骤如下:在程序初始化时根据配置信息中测频通道接入的位置确定哪一个作为测频CPU,哪一个作为跟随CPU。两个CPU之间的交互数据要求实时性非常高,每次采样中断都需要进行数据交换,本发明采用双口RAM,对两个CPU之间的数据传输提供支撑,且外部中断信号源确保两个CPU同时进入采样中断。进入中断时,两个CPU先从双口RAM中读取本次应使用重采样频率fF,该频率由测频CPU在上次中断中计算得出,并在上次中断结束前写入双口RAM;然后读取原始采样点。测频CPU根据原始采样率进入测频模块计算当前模拟量频率;然后进入重采样模块对原始数据线进行重采样,跟随CPU读取原始采样点后直接进入重采样模块。重采样前判断当前中断可以处理处几个点的重采样值,为确保中断执行时间不超限,每次最多只做两次重采样插值,多出的点由下次中断再处理,之后将本次中断计算出的重采样值和点数进行处理,发送给后端;重采样结束后测频CPU将计算出的当前频率写入双口RAM,待下次中断中使用,退出中断结束。本发明解决在数字化装置中采用双采样CPU等间隔采样时,两个CPU都需要接入测频通道才能实现频率跟踪的问题。
图1是本实施例的工作流程,所有工作内容在一个中断中完成,该中断为间隔等于标准频率下重采样间隔的等间隔中断。进入中断时,两个CPU先从双口RAM中读取本次应使用重采样频率fF,该频率由测频CPU在上次中断中计算得出,并在上次中断结束前写入双口RAM;然后读取原始采样点,一般原始采样点采样率高于重采样率,每次读取的原始采样点数多于1个,记录每个原始采样点的采样值Sn和采样时刻Tn。
测频CPU根据原始采样率进入测频模块计算当前模拟量频率;然后进入重采样模块对原始数据线进行重采样,此时使用的重采样频率为刚进入中断时从双口RAM中读到的频率,重采样前判断当前中断可以处理处几个点的重采样值,为确保中断执行时间不超限,每次最多只做两次重采样插值,多出的点由下次中断再处理,之后将本次中断计算出的重采样值和点数进行处理,发送给后端;重采样结束后测频CPU将计算出的当前频率写入双口RAM,待下次中断中使用,退出中断结束。跟随CPU在读取完原始采样点后直接进入重采样模块,进入模块后的处理方式与测频CPU相同,处理结束退出中断。
图2给出了测频模块计算过零点时间差的示意,步骤如下:
①每收到一个原始采样点时,根据公式计算出近一周波内原始采样点的傅里叶实部值。
其中N为原始采样点每周波采样点数,i为当前采样点序号,h=i-(N-1)为一周波前的采样点序号,Smpk为第k点的采样值。
②计算傅里叶实部值不满足条件的点数M,当满足该条件时,记录下Rei、Rei-1、Rei-2;设Rei-2时刻Ti-2=0,则Rei-1时刻Ti-1=TN、Rei时刻Ti=2*TN(TN为原始采样点采样间隔时间),利用拉格朗日二阶插值公式计算出实部值为零时对应的采样时刻:
则有:Tpos=2*TN-Tzero;Tnag=Tzero-TN;Tpos为过零点采样间隔的正轴时间,Tnag为过零点采样间隔的负轴时间。
傅里叶实部过零点的时间差:Tdelta=M*TN+T’pos+Tnag
其中,T’pos为上次计算得到的Tpos值,本次算出的Tpos留至下次计算频率时用,计算结束后将M赋为0。
③计算出本周波频率f’F=1s/Tdelta。
图3是重采样超值找到的差值点,重采样步骤如下:
①根据进入中断时读到的fF和重采样率F算出重采样间隔TF=F/fF,由上一个重采样点的时刻T’Re得到新的重采样时刻TRe=T′Re+TF;
②如果TRe≤Tn,在原始采样点缓冲区找出TRe时刻之前的两个采样点(采样值S0、S1,采样时刻T0、T1)和之后的一个采样点(采样值S2,采样时刻T2),根据拉格朗日二阶插值法得到TRe时刻的采样值
③判断下一个重采样点的时刻是否在Tn之前,如果是,TRe加上TN后重复上述步骤②(只重复执行一次,确保中断执行时间不超);将重采样后的数据和本次中断重采样点数目发送给后端,重采样时对采样值进行编号,编号范围为0~(50*F-1)循环,不受频率的变化影响。用于后端对两个CPU的重采样数据进行同步识别。
两个CPU对双口RAM的读取操作和测频CPU对双口RAM的写入操作构成了本发明的数据交互模块。当刚进入采样中断时,两个CPU都从双口RAM的指定地址读取本次中断中应使用的频率值,该值由上次中断结束时测频CPU写入,初次上电时该频率值默认50Hz。为确保两个CPU对双口RAM的读写操作不产生冲突,采用了分时读写和读写控制两种措施。
分时读写措施:测频CPU在中断结束前执行写操作,程序内部保证中断执行时间不小于50us(该时间大于写RAM的执行时间即可);两颗CPU在刚进入中断时执行读操作。这样通过时间控制避免双口RAM的读写冲突;
读写控制措施:写数据之前,先将写操作标志WR_BUSY置1,写结束最后将WR_BUSY清0。当读操作时判断WR_BUSY,若为0,直接读取;若为1,则等待10us后再读,若等待时间超过50us,认为写操作超时,放弃读取,采用上次中断中使用的值进行运算并告警。
虽然本发明已以较佳实施例公开如上,但它们并不是用来限定本发明,任何熟习此技艺者,在不脱离本发明之精神和范围内,自当可作各种变化或润饰,但同样在本发明的保护范围之内。因此本发明的保护范围应当以本申请的权利要求保护范围所界定的为准。
Claims (7)
1.一种双CPU等间隔采样的频率跟踪实现方法,其特征在于:通过配置信息中测频通道接入的位置确定两个CPU中哪一个作为测频CPU,哪一个作为跟随CPU;测频CPU在测频模块中计算出本次频率,在重采样模块中使用上次算出的频率计算出重采样间隔,对原始采样值进行重采样,通过数据交互模块将计算出的频率和重采样间隔传输给跟随CPU;跟随CPU从数据交互模块中读取重采样间隔,在重采样模块中对原始采样值进行重采样;
所述测频模块,计算出原始采样点的变换值,通过找出变换值的过零点时间差,得到所测频率;
具体方法如下:①根据近一周波内原始采样点数据计算出傅里叶实部值,②计算傅里叶实部值过零点的时间差,③计算本周波频率;
所述步骤①中的傅里叶实部值的计算公式为:
其中N为原始采样点每周波采样点数,i为当前采样点序号,h=i-(N-1)为一周波前的采样点序号,Smpk为第k点的采样值;
所述步骤②中计算实部过零点时间差方法,公式为计算频率条件,当满足该条件时,记录下Rei、Rei-1、Rei-2;并记录一个周波内不满足该公式条件的采样点个数M;设Rei-2时刻Ti-2=0,则Rei-1时刻Ti-1=TN、Rei时刻Ti=2*TN,TN为原始采样点采样间隔时间,利用拉格朗日二阶插值公式计算出实部值为零时对应的采样时刻:
则有:Tpos=2*TN-Tzero;Tnag=Tzero-TN;
傅里叶实部过零点的时间差:Tdelta=M*TN+T’pos+Tnag;
其中,T’pos为上次计算得到的Tpos值;
所述步骤②完成过零点时间差计算后将M赋为0;
所述步骤③中的本周波的频率f’F=1s/Tdelta;
所述数据交互模块,两个CPU通过双口RAM进行数据交换,实现数据迅速及时的传输;通过对数据操作的分时处理和互斥机制确保了数据交换的准确性;
所述的重采样模块,实现在一次中断中得到最多两个重采样结果,每个重采样结果都有重采样编号。
2.根据权利要求1所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:数据采集、测频、数据交互、重采样、数据传输在一个等间隔中断中完成。
3.根据权利要求2所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:两个CPU的中断进入时刻和中断间隔时间相同。
4.根据权利要求1所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:所述测频模块,每收到一个原始采样点时调用一次。
5.根据权利要求1所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:所述数据交互模块的分时处理,测频CPU在中断结束前执行写操作,程序内部保证中断执行时间不小于50us;两个CPU在刚进入中断时执行读操作;这样通过时间控制避免双口RAM的读写冲突。
6.根据权利要求1所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:所述数据交互模块的互斥机制,写数据之前,先将写操作标志WR_BUSY置1,写结束最后将WR_BUSY清0;当读操作时判断WR_BUSY,若为0,直接读取;若为1,则等待10us后再读,若等待时间超过50us,认为写操作超时,放弃读取,采用上次中断中使用的值进行运算并告警。
7.根据权利要求1所述的双CPU等间隔采样的频率跟踪实现方法,其特征在于:所述重采样模块,重采样编号为0~(50*F-1)循环,不受频率的变化影响,其中F为每个周波的重采样率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611199964.9A CN106597095B (zh) | 2016-12-22 | 2016-12-22 | 一种双cpu等间隔采样的频率跟踪实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611199964.9A CN106597095B (zh) | 2016-12-22 | 2016-12-22 | 一种双cpu等间隔采样的频率跟踪实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106597095A CN106597095A (zh) | 2017-04-26 |
CN106597095B true CN106597095B (zh) | 2019-09-10 |
Family
ID=58601001
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611199964.9A Active CN106597095B (zh) | 2016-12-22 | 2016-12-22 | 一种双cpu等间隔采样的频率跟踪实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106597095B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108037357B (zh) * | 2017-11-29 | 2019-11-01 | 西电通用电气自动化有限公司 | 一种电力系统自动化装置的频率跟踪方法 |
CN111736014B (zh) * | 2020-07-20 | 2020-11-24 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种频率信号采集和数据处理方法及装置 |
CN113466552B (zh) * | 2021-07-14 | 2024-02-02 | 南京海兴电网技术有限公司 | 一种定间隔采样下频率跟踪方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2385341Y (zh) * | 1999-03-24 | 2000-06-28 | 山东艾诺仪器有限公司 | 智能电力谐波分析仪 |
CN101350518A (zh) * | 2008-09-05 | 2009-01-21 | 国电南瑞科技股份有限公司 | 数字化保护测控装置采样值处理方法 |
CN102135567A (zh) * | 2011-02-25 | 2011-07-27 | 深圳市业通达实业有限公司 | 一种电力系统交流采样实时频率跟踪及谐波测量方法 |
WO2014013116A1 (es) * | 2012-07-18 | 2014-01-23 | Universitat Rovira I Virgili | Método y dispositivo para el seguimiento del punto de potencia máxima en sistemas fotovoltaicos |
CN103575980A (zh) * | 2012-07-26 | 2014-02-12 | 施耐德电器工业公司 | 系统频率测量方法、同步相量测量方法及设备 |
CN103675522A (zh) * | 2013-11-12 | 2014-03-26 | 国电南瑞科技股份有限公司 | 一种面向间隔的智能变电站多功能二次装置及采样平台 |
CN106018957A (zh) * | 2016-05-18 | 2016-10-12 | 南京大全自动化科技有限公司 | 一种适用于馈线终端装置的硬件测频及频率跟踪方法 |
-
2016
- 2016-12-22 CN CN201611199964.9A patent/CN106597095B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2385341Y (zh) * | 1999-03-24 | 2000-06-28 | 山东艾诺仪器有限公司 | 智能电力谐波分析仪 |
CN101350518A (zh) * | 2008-09-05 | 2009-01-21 | 国电南瑞科技股份有限公司 | 数字化保护测控装置采样值处理方法 |
CN102135567A (zh) * | 2011-02-25 | 2011-07-27 | 深圳市业通达实业有限公司 | 一种电力系统交流采样实时频率跟踪及谐波测量方法 |
WO2014013116A1 (es) * | 2012-07-18 | 2014-01-23 | Universitat Rovira I Virgili | Método y dispositivo para el seguimiento del punto de potencia máxima en sistemas fotovoltaicos |
CN103575980A (zh) * | 2012-07-26 | 2014-02-12 | 施耐德电器工业公司 | 系统频率测量方法、同步相量测量方法及设备 |
CN103675522A (zh) * | 2013-11-12 | 2014-03-26 | 国电南瑞科技股份有限公司 | 一种面向间隔的智能变电站多功能二次装置及采样平台 |
CN106018957A (zh) * | 2016-05-18 | 2016-10-12 | 南京大全自动化科技有限公司 | 一种适用于馈线终端装置的硬件测频及频率跟踪方法 |
Non-Patent Citations (2)
Title |
---|
"基于插值和时标变换的智能变电站保护数据重采样算法";蔡超 等;《电力系统自动化》;20131010;第80-85页 |
"电力有源滤波器中的几种谐波实时检测方法研究";张俊敏;《CNKI硕士学位论文数据库》;20030301;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106597095A (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597095B (zh) | 一种双cpu等间隔采样的频率跟踪实现方法 | |
CN106338706B (zh) | 一种电能计量装置整体误差检测的方法、装置和系统 | |
CN108595301B (zh) | 一种基于机器学习的服务器能耗预测方法和系统 | |
CN106154816B (zh) | 一种自动装置高精度守时方法 | |
CN108874535B (zh) | 一种任务调节方法、计算机可读存储介质及终端设备 | |
CN105243068A (zh) | 数据库系统的查询方法、服务器和能耗测试系统 | |
CN104619009B (zh) | 一种定位数据采样周期的调整方法、装置及移动终端 | |
CN111998918A (zh) | 一种误差校正方法、误差校正装置及流量传感系统 | |
CN103033675A (zh) | 一种交流同步采集系统及其实现方法 | |
CN103365965A (zh) | 一种数据的汇总处理方法和装置 | |
CN107656851A (zh) | 一种基于部件能耗模型的云服务器能耗测算方法及系统 | |
CN110990235A (zh) | 异构存储设备的性能数据管理方法、装置、设备及介质 | |
CN111948454B (zh) | 一种抗衰减直流分量干扰的同步相量测量方法及系统 | |
CN103018547B (zh) | 归一化多处理器电功率计量方法 | |
CN109885599B (zh) | 基于多设备乱序数据的实时计算方法及系统 | |
CN107220166A (zh) | 一种cpu占用率的统计方法及装置 | |
CN105160446A (zh) | 一种获取借款额度的方法及装置 | |
CN104316893A (zh) | 互感器负荷箱的全自动校验方法 | |
CN104317704A (zh) | 一种生产线自动测试系统的测量软件集成方法 | |
CN104698929B (zh) | 费控智能电能表的实时控制方法 | |
CN104573136B (zh) | 电网数据处理方法及装置 | |
CN203893955U (zh) | 恒温槽性能测试系统 | |
CN110780114A (zh) | 一种台区理论线损分析方法及可读存储介质 | |
CN111062004A (zh) | 一种预测模型调整方法、装置、预测模拟器以及存储介质 | |
CN108804640A (zh) | 基于最大化iv的数据分组方法、装置、储存介质及设备 |
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 |