CN104734692A - 基于频差法和dds实现的高采样率等效采样方法和*** - Google Patents
基于频差法和dds实现的高采样率等效采样方法和*** Download PDFInfo
- Publication number
- CN104734692A CN104734692A CN201510043845.3A CN201510043845A CN104734692A CN 104734692 A CN104734692 A CN 104734692A CN 201510043845 A CN201510043845 A CN 201510043845A CN 104734692 A CN104734692 A CN 104734692A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- sampled
- 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.)
- Granted
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
一种等效采样方法和***，该方法基于频差法原理，在同一频率基准上产生的两路信号，经过处理，一路作为采样时钟信号，一路作为待采样信号。可以使两路信号频率差很小，每个待采样信号周期进行一次采样，通过多次顺序采样，获得完整的、时间展宽的待采样信号波形，实现顺序等效采样；还可以使两个信号频率有倍数关系，每个待采样信号周期进行多次采样，再对采样数据进行排列组合，恢复出时间展宽的待采样信号波形，在实现伪随机等效采样的同时，降低等效采样时间。本发明可以广泛应用在采用周期性脉冲信号作为测量载体的穿墙雷达、防撞雷达等，采用周期性连续波信号作为测量载体的探地雷达、时域反射仪等。
Description
技术领域
本发明涉及一种基于频差法和DDS(Direct Digital frequency Synthesis，直接数字频率合成)实现的高采样率等效采样方法和***，具体而言，设计了一种对周期性信号进行等效采样的方法和***。该项方法适用于周期性信号的等效采样。
背景技术
对信号波形的数据采集通常有两种方法：实时采样和等效采样；其中，等效采样又分为三种，顺序等效采样、伪随机等效采样和随机等效采样，其中，伪随机等效采样可以看作是改进的顺序等效采样。
目前的顺序等效采样方法可以分为两种实现方式，一种时域方式，一种频域方式。时域实现方式包含采用延时芯片实现、采用斜波比较法实现和采用特殊工艺法实现等；频域实现方式包含采用频差法实现、采用频率步进法实现等。如专利申请号为201410001964.8的专利所述，采用级联延时芯片实现；如文章名为《时域反射仪高精度步进延迟***设计》所述，采用了斜波比较法产生步进延时；如文章名为《基于DDFS和恒精度测频法的等效采样***设计》所述，采用了频差法实现等效采样；如论文名为《频率步进穿墙雷达***设计及仿真研究》所述，采用频率步进法实现等效采样。
目前的时域方式实现顺序等效采样方法在应用中存在以下问题：
1、步进延时步长长，等效采样率低。采用低采样率的芯片实现等效高采样率，需要精确控制采样时钟信号/待采样信号的步进延时。目前实现步进延时的芯片或者斜波比较法，其等效步进延时步长长，不能够实现高精度的步进延时步长控制，最终实现的等效采样率低；
2、步进延时步长精度低。常用的延时电路采用专用延时芯片来实现，存在延时芯片延时精度不高，不同延时芯片的延时单元量有偏差，这将导致延时精度低，并且线性度差，影响信号的采样质量，进而影响***的测量精度；
3、需要步进延时***，控制复杂。实现步进延时的方法包括延时芯片法、斜波比较法和特殊工艺法等，这些方法实现起来控制相对复杂，尤其是总的步进延时比较大，需要进行级联控制时，控制复杂。
目前频域方式实现顺序等效采样方法在应用中存在以下问题：
1、频差法需要一个高精度的输入信号频率测量模块。传统意义的频差法测量，需要一个测量输入信号频率测量模块，之后产生一个与输入信号频率相近采样时钟信号，实现等效采样；
2、频率步进法需要产生一个连续的步进频率信号，该信号的周期在连续改变，会伴随一定的杂散频率，最终输出时会带来较大的抖动，降低信号质量，导致等效采样率低。
针对上述问题，本专利采用频差法实现方式，提出了独特的思路，将时域的问题转移到频域处理，并且将频域的问题处理简单化。本发明提出的一种基于频差法和DDS实现的高采样率等效采样方法和***，解决了高速采样芯片价格高且难以获取的问题，解决了现有时域法实现等效采样需要延时***且延时精度不高的问题，降低了***复杂度，提高了***稳定性；本发明还解决了频差法需要高精度信号频率测量模块的问题，解决了频率步进法的杂散频率引入的信号抖动，极大地提高了***的等效采样率。
本发明提出的一种基于频差法和DDS实现的高采样率等效采样方法和***，可以实现很高的采样率，适用于对周期性的连续波信号和脉冲信号进行顺序采样，配以高带宽、低抖动的取样头或者取样电路，可以实现高精度时空分辨率的测量。在电子学抖动很小或者不考虑抖动的情况下，可以实现1000GS/s以上的等效采样率。
本发明提出的方法和***，可以广泛应用在各类型雷达中，例如采用周期性脉冲信号作为测量载体的穿墙雷达、防撞雷达等，采用周期性连续波信号作为测量载体的探地雷达、时域反射仪等。
对数据采集时间有要求的场合，可以同时提高采样时钟频率和待采样信号频率，也可以只提高采样时钟频率，在第二种情况下，采样是伪随机等效采样，对采样数据的恢复组合，就相对麻烦一些，具体的数据恢复组合方式可以参考专利申请公开号为CN103634006A的专利所述。
发明内容
本发明克服现有技术存在的问题，提出一种基于频差法和DDS实现的高采样率等效采样方法和***，该方法通过频域的频差法来实现。该方法的采用，解决了现有时域法实现等效采样需要延时***且延时精度不高的问题，降低了***复杂度，提高了***稳定性；本发明解决了频差法需要高精度信号频率测量模块的问题，解决了频率步进法的杂散频率引入的信号抖动，极大地提高了***的等效采样率。
本发明的技术解决方案是：
一种基于频差法和DDS实现的高采样率等效采样方法，其特别之处在于：包括以下步骤：
1)产生一个高稳时钟基准；
2)根据高稳时钟基准产生两路信号，对信号进行滤波、波形转换处理；其中一路作为采样时钟信号，另一路经过待测模块的反射，作为待采样信号；
所述两路信号的频率按照采样方式不同按以下方式确定：
若进行顺序等效采样，则两路信号之间的频率差很小且采样时钟信号频率小于待采样信号频率；
若进行伪随机等效采样，则采样时钟信号频率为待采样信号频率的近似倍数N；
3)根据采样时钟信号的控制，数据采集***对待采样信号进行采样，得到采样数据；
4)若进行顺序等效采样，则对采样数据按顺序排列并输出；
若进行伪随机等效采样，则对采样数据进行排列组合并输出。
还包括调节等效采样率的步骤：通过调节两路信号之间的频率差，实现不同的等效采样率。
还包括降低等效采样时间的步骤：通过调节采样时钟信号频率和待采样信号频率之间的N倍数关系，可以N倍降低***的总的等效采样时间。
还包括调节采样起始点的步骤：通过调节相位差和延时，实现不同的起始点开始采样。
一种基于频差法和DDS实现的高采样率等效采样***，其特别之处在于：包括时钟基准模块U1、信号产生模块U2、输入输出模块U3、数据采集模块U4、数据处理模块U5和控制模块U6；
所述时钟基准模块U1产生一个高稳时钟基准，并送入信号产生模块U2；
所述信号产生模块U2根据高稳时钟基准产生两路信号，并对信号进行滤波、波形转换等处理，其中一路信号送入数据采集模块U4作为采样时钟信号，另一路信号经过输入输出模块U3送入待测模块U7，其信号反射部分经过输入输出模块U3作为数据采集模块U4的待采样信号；
所述控制模块U6控制信号产生模块U2的输出信号，控制数据处理模块U5对采样数据进行排列组合；
所述数据处理模块U5对数据采集模块U4输出的采样数据进行排列组合，恢复出时间展宽的待采样信号波形。
所述信号产生模块U2包括带有双通道输出、高频率控制字、高相位控制字、高分辨率DAC的DDS芯片，或者是两片单道输出、高频率控制字、高相位控制字、高分辨率DAC的DDS芯片；所述控制模块U6通过信号产生模块U2调节两路信号之间的频率差、相位差和延时参数。
本发明所具有的优点：
1、本发明采用频差法实现等效高采样率，相对于传统的时域、频域实现方法，降低了***复杂度，提高了***的稳定性，极大地提高了***的等效采样率；
2、本发明采用频差法实现，相对于时域方法实现，优势明显。时域方法对于延时的控制，小于10ps已经很不易实现，频差法则不同，通过调节频率差，可以很轻易实现10ps的步进延时，在电子学抖动很小或者不考虑抖动的情况下，甚至可以达到优于1ps的步进延时，从而实现很高的顺序等效采样率；
3、本发明采用高稳晶振作为时钟源基准，采用DDS芯片产生两路信号，两路信号之间的频率差、相位差和延时等参数可调，并且在传输过程中保持稳定，其中一路作为采样时钟信号，另一路作为待采样信号；
4、本发明采用的DDS芯片，具有很高的频率控制字、相位控制字，可以实现高的频率分辨率，可以根据需要，更改DDS的输出的两路信号之间的频率差。当两路信号之间的频率差很小，通过调节两路信号之间的频率差，可以实现不同的顺序等效采样率，通过调节相位差和延时可以从不同的起始点开始测量；当两路信号之间的频率有近似倍数N的关系，实现伪随机等效采样，在相同的等效采样率下，伪随机等效采样时间是顺序等效采样所需时间的1/N，通过调节两路信号频率，实现不同的伪随机等效采样率，并缩短伪随机等效采样所需的时间。
5、本发明采用的DDS芯片，可以根据需要，产生各种连续波和脉冲信号，实现不同的测量方式，应用于不同场合。
附图说明
图1是本发明等效采样***原理图；
图2是信号产生模块原理图；
图3是顺序等效采样原理图；
图4是伪随机等效采样数据排列组合原理图。
具体实施方式
本发明用于实现等效高采样率，该等效采样方法可以广泛应用在各类型雷达中，例如采用周期性脉冲信号作为测量载体的穿墙雷达、防撞雷达等，采用周期性连续波信号作为测量载体的探地雷达、时域反射仪等。
参见图1，该***主要包括时钟基准模块U1、信号产生模块U2、输入输出模块U3、数据采集模块U4、数据处理模块U5、控制模块U6、待测模块U7。
控制模块U6控制信号产生模块U2的两路输出信号，控制数据处理模块U5，对采样数据进行排列，恢复出时间展宽的待采样信号波形，根据波形推出待测模块U7的性质。
信号产生模块U2产生两路频率差很小/两路频率近似成倍数N关系的信号，最终两路信号分别作为数据采集模块U4的采样时钟信号、待采样信号，当两路信号之间的频率差很小，采样时钟信号频率小于待采样信号频率，通过数据处理模块U5对数据顺序排列，实现顺序等效采样，通过调整两路信号之间的频率差，实现不同的顺序等效采样率；当两路信号之间的频率有近似倍数N的关系，采样时钟信号频率为待采样信号频率的倍数，通过数据处理模块U5对数据进行组合排列，实现伪随机等效采样；在相同的等效采样率下，伪随机等效采样时间是顺序等效采样所需时间的1/N。
根据恢复出的时间展宽的待采样信号波形，推出待测模块U7的性质。
以下结合本发明，对顺序等效采样方法和伪随机等效采样方法的原理及其实施作具体介绍。
信号产生模块U2产生两路频率差很小的信号，最终两路信号分别作为数据采集模块U4的采样时钟信号、待采样信号，两路信号之间的频率差很小，采样时钟信号频率小于待采样信号频率。两路信号fclk、fdata之间的频率差Δf很小，反映到时域上，两路信号的周期有微小时间差别Δt，Δt就是等效采样的步进延时，其原理如图3所示。
下面举例说明顺序等效采样，采用10MS/s的实时采样率实现10GS/s的顺序等效采样率。将采样时钟信号频率设置为10MHz，fclk＝10MHz＝10000000Hz，周期Tclk＝100000ps。要达到10GS/s的顺序等效采样率，Δt＝100ps，那么，对应待采样信号的周期为Tdata＝Tclk-Δt＝99900ps，fdata≈10010010Hz，Δf＝fdata-fclk＝10010Hz。同样以采用10MS/s的实时采样率实现500GS/s的顺序等效采样率为例，fclk＝10MHz＝10000000Hz，周期Tclk＝100000ps。要达到500GS/s的顺序等效采样率，Δt＝2ps，那么，对应待采样信号的周期为Tdata＝Tclk-Δt＝99998ps，fdata≈10000200Hz，Δf＝fdata-fclk＝200Hz。通过信号产生模块U2产生两路信号，频率分别为10000000Hz、10000200Hz，可以采用10MS/s的实时采样率实现500GS/s的顺序等效采样率，因为每个待采样信号周期99998ps只采样一个点，所以，最终采集完成并恢复出展宽的待采样信号波形，需要较长的总的等效采样时间，如果采样点为65536点，那么总的等效采样时间为65536×99998ps≈6.55ms，这个量级的总的等效采样时间是可以接受的。但是，当待采样信号频率为fdata＝10000Hz，对应周期为100us，那么，在相同的等效采样率和采样点的情况下，采样时钟信号频率为fclk＝9999.9998Hz，Δf＝fdata-fclk＝0.0002Hz，等效采样时间为6.55s。这时候，可以通过伪随机等效采样，即保证了高的等效采样率，又降低了总的等效采样时间。
以频率为fdata＝10000Hz的待采样信号为例，说明伪随机等效采样方法，其数据排列组合原理如图4所示。等效采样率为500GS/s，采样点数65536点，对顺序等效采样来讲，其对应的fclk＝9999.9998Hz，总的等效采样时间约为6.55s。假设伪随机等效采样率500GS/s，采样点为65536点，如果采样时钟信号频率为fclk＝29999.9994Hz，和待采样信号频率大约有3倍关系，N＝3，那么总的等效采样时间约为6.55s/3＝2.18s；如果采样时钟信号频率为fclk＝29999999.4Hz，和待采样信号频率大约有3000倍关系，N＝3000，那么总的等效采样时间约为6.55s/3000＝2.18ms。在伪随机等效采样中，N的倍数越大，总的等效采样时间越小，相应地，其伪随机等效采样数据排列组合越复杂。
Claims (6)
1.一种基于频差法和DDS实现的高采样率等效采样方法，其特征在于：包括以下步骤：
1)产生一个高稳时钟基准；
2)根据高稳时钟基准产生两路信号，对信号进行滤波、波形转换处理；其中一路作为采样时钟信号，另一路经过待测模块的反射，作为待采样信号；
所述两路信号的频率按照采样方式不同按以下方式确定：
若进行顺序等效采样，则两路信号之间的频率差很小且采样时钟信号频率小于待采样信号频率；
若进行伪随机等效采样，则采样时钟信号频率为待采样信号频率的近似倍数N；
3)根据采样时钟信号的控制，数据采集***对待采样信号进行采样，得到采样数据；
4)若进行顺序等效采样，则对采样数据按顺序排列并输出；
若进行伪随机等效采样，则对采样数据进行排列组合并输出。
2.根据权利要求1所述基于频差法和DDS实现的高采样率等效采样方法，其特征在于：
还包括调节等效采样率的步骤：通过调节两路信号之间的频率差，实现不同的等效采样率。
3.根据权利要求1或2所述基于频差法和DDS实现的高采样率等效采样方法，其特征在于：
还包括降低等效采样时间的步骤：通过调节采样时钟信号频率和待采样信号频率之间的N倍数关系，可以N倍降低***的总的等效采样时间。
4.根据权利要求3所述基于频差法和DDS实现的高采样率等效采样方法，其特征在于：
还包括调节采样起始点的步骤：通过调节相位差和延时，实现不同的起始点开始采样。
5.一种基于频差法和DDS实现的高采样率等效采样***，其特征在于：包括时钟基准模块(U1)、信号产生模块(U2)、输入输出模块(U3)、数据采集模块(U4)、数据处理模块(U5)和控制模块(U6)；
所述时钟基准模块(U1)产生一个高稳时钟基准，并送入信号产生模块(U2)；
所述信号产生模块(U2)根据高稳时钟基准产生两路信号，并对信号进行滤波、波形转换等处理，其中一路信号送入数据采集模块(U4)作为采样时钟信号，另一路信号经过输入输出模块(U3)送入待测模块(U7)，其信号反射部分经过输入输出模块(U3)作为数据采集模块(U4)的待采样信号；
所述控制模块(U6)控制信号产生模块(U2)的输出信号，控制数据处理模块(U5)对采样数据进行排列组合；
所述数据处理模块(U5)对数据采集模块(U4)输出的采样数据进行排列组合，恢复出时间展宽的待采样信号波形。
6.根据权利要求5所述的基于频差法和DDS实现的高采样率等效采样***，其特征在于：所述信号产生模块(U2)包括带有双通道输出、高频率控制字、高相位控制字、高分辨率DAC的DDS芯片，或者是两片单道输出、高频率控制字、高相位控制字、高分辨率DAC的DDS芯片；所述控制模块(U6)通过信号产生模块(U2)调节两路信号之间的频率差、相位差和延时参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510043845.3A CN104734692B (zh) | 2015-01-28 | 2015-01-28 | 基于频差法和dds实现的高采样率等效采样方法和*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510043845.3A CN104734692B (zh) | 2015-01-28 | 2015-01-28 | 基于频差法和dds实现的高采样率等效采样方法和*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104734692A true CN104734692A (zh) | 2015-06-24 |
| CN104734692B CN104734692B (zh) | 2017-09-29 |
Family
ID=53458201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510043845.3A Active CN104734692B (zh) | 2015-01-28 | 2015-01-28 | 基于频差法和dds实现的高采样率等效采样方法和*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104734692B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106443630A (zh) * | 2016-11-03 | 2017-02-22 | 长沙太电子科技有限公司 | 一种高效的组合式顺序等效采样方法 |
| CN106772269A (zh) * | 2017-03-03 | 2017-05-31 | 南京邮电大学 | 一种应用探地雷达回波信号采集的等效采样电路 |
| CN109324215A (zh) * | 2018-09-21 | 2019-02-12 | 北京无线电计量测试研究所 | 一种基于dds的标准相位产生方法和装置 |
| CN109324313A (zh) * | 2018-09-25 | 2019-02-12 | 北京无线电计量测试研究所 | 一种用于冲激雷达的等效时间采样***和方法 |
| CN110736877A (zh) * | 2019-09-26 | 2020-01-31 | 山东信通电子股份有限公司 | 一种时域反射信号的高速采集方法及装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6194945B1 (en) * | 1994-07-15 | 2001-02-27 | Unisys Corporation | Dual threshold digital receiver with large noise margin |
| CN101144835A (zh) * | 2007-10-17 | 2008-03-19 | 电子科技大学 | 一种等效采样装置 |
| CN101312357A (zh) * | 2008-07-03 | 2008-11-26 | 电子科技大学 | 一种用于Chirp超宽带通信***的调制方法 |
| CN101666819A (zh) * | 2009-07-21 | 2010-03-10 | 秦轲 | 数字示波器信号等效采样方法 |
| CN101888258A (zh) * | 2010-06-02 | 2010-11-17 | 航天恒星科技有限公司 | 高动态环境下基于3g的geo卫星移动通信时隙同步***及方法 |
-
2015
- 2015-01-28 CN CN201510043845.3A patent/CN104734692B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6194945B1 (en) * | 1994-07-15 | 2001-02-27 | Unisys Corporation | Dual threshold digital receiver with large noise margin |
| CN101144835A (zh) * | 2007-10-17 | 2008-03-19 | 电子科技大学 | 一种等效采样装置 |
| CN101312357A (zh) * | 2008-07-03 | 2008-11-26 | 电子科技大学 | 一种用于Chirp超宽带通信***的调制方法 |
| CN101666819A (zh) * | 2009-07-21 | 2010-03-10 | 秦轲 | 数字示波器信号等效采样方法 |
| CN101888258A (zh) * | 2010-06-02 | 2010-11-17 | 航天恒星科技有限公司 | 高动态环境下基于3g的geo卫星移动通信时隙同步***及方法 |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106443630A (zh) * | 2016-11-03 | 2017-02-22 | 长沙太电子科技有限公司 | 一种高效的组合式顺序等效采样方法 |
| CN106443630B (zh) * | 2016-11-03 | 2017-09-12 | 湖南太康电子信息技术有限公司 | 一种高效的组合式顺序等效采样方法 |
| CN106772269A (zh) * | 2017-03-03 | 2017-05-31 | 南京邮电大学 | 一种应用探地雷达回波信号采集的等效采样电路 |
| CN109324215A (zh) * | 2018-09-21 | 2019-02-12 | 北京无线电计量测试研究所 | 一种基于dds的标准相位产生方法和装置 |
| CN109324313A (zh) * | 2018-09-25 | 2019-02-12 | 北京无线电计量测试研究所 | 一种用于冲激雷达的等效时间采样***和方法 |
| CN110736877A (zh) * | 2019-09-26 | 2020-01-31 | 山东信通电子股份有限公司 | 一种时域反射信号的高速采集方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104734692B (zh) | 2017-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104734692A (zh) | 基于频差法和dds实现的高采样率等效采样方法和*** | |
| US9077513B2 (en) | High precision synchronized measured value acquisition | |
| CN101976037B (zh) | 一种多次同步模拟内插的时间间隔测量方法和装置 | |
| CN102882673B (zh) | 多通道高速dac同步实现方法 | |
| CN105549379A (zh) | 一种基于高精度时间基准触发的同步测量装置及方法 | |
| US10177897B2 (en) | Method and system for synchronizing and interleaving separate sampler groups | |
| CN110350892B (zh) | 一种基于dds时钟移相技术的延时装置及方法 | |
| CN106647926B (zh) | 用于冷原子干涉仪激光时序控制的dds跳频装置 | |
| US9553600B1 (en) | Skew detection and correction in time-interleaved analog-to-digital converters | |
| CN110445493A (zh) | 一种基于fpga tdc的数据采集同步装置及方法 | |
| CN103487649A (zh) | 一种兼容连续波和脉冲调制载波频率测量的方法及装置 | |
| CN112968690B (zh) | 一种高精度低抖动延时脉冲发生器 | |
| CN102723931A (zh) | 一种宽动态高精度边沿时间可调的脉冲波产生方法 | |
| CN102859880A (zh) | 用于信号转换器的源同步测试的装置和方法 | |
| CN103728893B (zh) | 一种探地雷达高精度时序控制电路 | |
| CN202421767U (zh) | 一种时统设备 | |
| CN103595522A (zh) | 一种双通道数模转换器同步方法 | |
| CN114063505A (zh) | 一种同步控制方法、装置、设备、同步***及存储介质 | |
| EP3179225B1 (en) | Apparatus and method for low latency, reconfigurable and picosecond resolution time controller | |
| CN104133409A (zh) | 一种对称性可调的三角波合成装置 | |
| CN102916681B (zh) | 一种脉宽可调的nrz/rz码转换装置 | |
| Sudars et al. | Signal Analog-to-Event-to-Digital converting based on periodic sampling and precise event timing | |
| CN109445272A (zh) | 时钟信号同步方法、调节信号频率的方法及装置 | |
| US20070188363A1 (en) | Equivalent time asynchronous sampling arrangement | |
| CN102916682B (zh) | 一种脉宽可调的不归零码/归一码转换装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |