CN105301114A - Acoustic coating layer insertion loss measurement method based on multi-channel space-time inverse filtering technology - Google Patents
Acoustic coating layer insertion loss measurement method based on multi-channel space-time inverse filtering technology Download PDFInfo
- Publication number
- CN105301114A CN105301114A CN201510683611.5A CN201510683611A CN105301114A CN 105301114 A CN105301114 A CN 105301114A CN 201510683611 A CN201510683611 A CN 201510683611A CN 105301114 A CN105301114 A CN 105301114A
- Authority
- CN
- China
- Prior art keywords
- signal
- insertion loss
- liftering
- acoustic
- channel
- 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
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Provided is an acoustic coating layer insertion loss measurement method based on a multi-channel space-time inverse filtering technology. The method comprises the following steps that 1, a multi-channel space-time inverse filtering emission signal is generated; 2, a direct signal is collected on the condition that a test sample does not exist; 3, a transmission signal is collected on the condition that the test sample exists; 4, an insertion loss measured value is calculated. According to the acoustic coating layer insertion loss measurement method based on the multi-channel space-time inverse filtering technology, a channel impulse response function containing information of a circuit channel and an underwater acoustic channel is obtained through constant interval emission of an emission array transducer and sequential receiving of a hydrophone, and the purpose of receiving a space acoustic energy focus and a time domain impulse waveform focus of the hydrophone is achieved through the simple and effective multi-channel space-time inverse filtering technology, so that the purposes of multi-path signal separation and reverberation suppression in acoustic coating layer insertion loss measurement are achieved, the space resolution of measurement parameters is improved, and the method is particularly suitable for acoustic performance tests of low-mid frequency materials. Meanwhile, it is verified that the method is effective in acoustic coating layer insertion loss measurement.
Description
Technical field
The present invention relates to a kind of measuring method of underwater sound overlayer full-page proof insertion loss.Underwater sound overlayer uses extensive and vital parts under water in Underwater Acoustics Engineering, and be the insertion loss performance of assessment acoustic stimulation sample, under laboratory condition, the insertion loss measurement of large part sample is main means.
Background technology
Acoustic stimulation refers to dedicated functionality acoustical material on the parts under water of being laid in and structure, is that parts realize an important means of Sound stealth under water.Along with the fast development of viscoelastic material and going deep into of Sound stealth technical research, acoustic stimulation product is comparatively complete, but can divide for by its different acoustic function watt to be shaken watt by sound eliminating tile, sound insulation, decoupling watt, array exhaust box etc.The energy of active probe sound wave can significantly absorb by acoustic stimulation on the one hand, and active sonar sound wave reflection is reduced, and reduces the acoustic target intensity of submerged structure, can suppress structure self radiated noise on the other hand.Thus, acoustic stimulation uniquely to reduce at present the effective technology that target strength can suppress again radiated noise, also just becomes our problem in the urgent need to address to its acoustical behavior research in underwater sound environment for use and parameters,acoustic measurement.
Insertion loss is one of important indicator weighing overlayer acoustical behavior.Existing acoustic stimulation insertion loss is measured must meet free-field condition.Large sample in free field is measured limited for the measurement capability of frequency acoustic performance parameter, along with the decline of frequency, because the acoustic absorbant of measurement environment surrounding cannot absorb sound wave completely, border all-bottom sound wave reflection can not be ignored, and measuring-signal is by reverberation serious interference.Meanwhile, because beam main lobe width can increase gradually with the reduction of frequency, namely main lobe resolving power can reduce gradually, and sample to be tested size is limited, when main lobe width is greater than sample to be tested size, can cause comparatively big error to test result.In addition, for the model of multilayer shell, the transmission performance usually need assessing every one deck sample respectively, to realize optimum comprehensive sound absorbing capabilities, thus has higher requirement to measurement.
Summary of the invention
In order to overcome existing acoustic stimulation insertion loss measuring technique in the test of medium and low frequency section, the deficiency that measuring result error is comparatively large, precision is lower, the invention provides one to be applicable to medium and low frequency section, effectively to reduce measuring error, to put forward high-precision acoustic stimulation full-page proof insertion loss measuring method, liftering technology when the method is empty by hyperchannel, the dimensional energy realizing transmitting at sample place focuses on and time domain impulse focuses on.
The present invention is as follows for solving the problems of the technologies described above adopted technical scheme.
An acoustic stimulation insertion loss measuring method for liftering technology time empty based on hyperchannel, described measuring method comprises the steps:
1) when generation hyperchannel is empty, liftering transmits: without in sample situation, each transducer of transducer array launches initialize signal successively, and nautical receiving set is Received signal strength successively, utilizes single channel liftering technology to obtain single channel liftering and transmits; According to the time delay spacing that transducer transmits, the liftering corresponding to each passage transmits and carries out time delay and even up, generate hyperchannel empty time liftering transmit;
2) without the collection of incoming signal in sample situation: liftering signal when transducer transmitting hyperchannel is empty, then according to liftering principle during sky, transmit and produce space-focusing and time domain is waveform of sharp pulses, the incoming signal p of hydrophone recordings sample at nautical receiving set place
i;
3) there is the collection of transmission signal in sample situation: the sample laying acoustic stimulation is put into test environment, repeat
2) step in, records transmission focusing signal p during test sample
t.
4) calculating of insertion loss measured value: utilize the computing formula of insertion loss to obtain
Wherein I
lrepresent insertion loss value.
Technical conceive of the present invention is: launched by initial steering signal and nautical receiving set reception, obtain the Received signal strength comprising circuit channel and underwater acoustic channel information, by hyperchannel liftering technology, obtain transmitting of each array element of transducer array, penetrate through paroxysm, obtain receiving hydrophone place space-focusing and time domain is the signal of impulse form, thus realize space-focusing and the Time-Domain Pulse Compression of measuring-signal, and suppress reverberation.
Compared with existing insertion loss measuring method, advantage of the present invention is:
1) utilize liftering technology to carry out inversion calculation acquisition to the whole signal transmission channel comprising circuit channel and underwater acoustic channel to transmit, realize the Time-Domain Pulse Compression of incoming signal, be conducive to the separation of multipath signal.
2) synchronized transmissions of multi channel signals is utilized, improve incident wave main lobe resolution, reduce reverberation to the impact of measurement result, also can increase the spatial resolution of measurement parameter simultaneously, reduce the requirement to specimen size and the test space, be specially adapted to the test of medium and low frequency section Material acoustic performance.
3) test environment water body internal mechanism can because liftering technology during sky is set up equally under nonhomogeneous media environment, therefore, be eliminated on the impact of measuring-signal.
Accompanying drawing explanation
Fig. 1 is the measuring system overall schematic of the inventive method.
The space-focusing oscillogram of liftering technology when Fig. 2 is hyperchannel sky.
When Fig. 3 is hyperchannel sky, the time domain impulse of liftering technology focuses on oscillogram.
Fig. 4 is 5 millimeters thick steel plate sample insertion loss theories and measures comparison of test results figure.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.
With reference to Fig. 1 ~ Fig. 4, a kind of acoustic stimulation insertion loss measuring method based on liftering technology during hyperchannel sky, for the measurement of finite space acoustic stimulation full-page proof insertion loss.The technical scheme of a whole set of measuring method is as follows:
1) when generation hyperchannel is empty, liftering transmits;
Without in sample situation, each transducer of transducer array launches initialize signal successively, and nautical receiving set is Received signal strength successively, if the frequency domain that transmits can be expressed as s (f), wherein f represents signal frequency, and Received signal strength is expressed as x
n(f), wherein n=1,2 ..., N represents that transducer is numbered, and N represents transducer number;
According to formula
Estimate the channel response function H between each transmission channel and receiving hydrophone
n(f); For avoiding the singular problem occurred in inversion calculation, utilize formula
Estimate the inverse function obtaining channel response
wherein ε is a little value constant, unusual to avoid inverting; Inverse Fourier transform is done to this inverse function and obtains time domain Secondary Emission signal e
n(t), wherein t represents the time;
Utilize the signal transmission time interval of each transmitting transducer to liftering signal e
nt () is carried out time delay and is evened up, the e thus the liftering of acquisition hyperchannel transmits
n(t-τ
n), wherein τ
nfor the delay volume that each transmission channel is corresponding;
2) without the collection of direct signal during sample;
As shown in Figure 1, instrument is exported by computer system control hyperchannel individual control signal, through multichannel power amplifier, utilize transducer array synchronized transmissions by step 1) generate empty time liftering signal, meanwhile, the direct signal p when hydrophone recordings of computer system control focal position does not test sample
i.
3) collection of transmission signal during sample is had;
The sample laying acoustic stimulation is put into test environment, repeats 2) middle step, now record transmission focusing signal p during test sample
t.
4) calculating of insertion loss measured value;
The measurement data calculated needed for insertion loss can be obtained by above three steps.Utilize 2) and 3) the incoming signal p that gathers
iand transmission signal p
t, can be calculated insertion loss measured value is
Example illustrates: be focusing effect during hyperchannel liftering technology empty in checking the present invention, carried out simulation analysis.Simulation pond measurement environment, if pond degree of depth 1.2m, length 15m, the water body velocity of sound is c
w=1480m/s, the density of water is ρ
w=1g/cm
3, be the thick sand bed of 20cm below water layer, its velocity of sound is respectively c
p=1700m/s, density is ρ
s=1.8g/cm
3, attenuation coefficient is α
s=0.5dB/ λ is the semi-infinite space of the hard end below sedimentary deposit.4 sound sources lay respectively at degree of depth 0.3m, 0.5m, 0.7m and 0.9m, receiving hydrophone range transmission battle array 10m, degree of depth 0.55m, and Received signal strength average signal-to-noise ratio is 15dB.In liftering algorithm, ε gets 0.05.When what Fig. 2 showed is hyperchannel sky, liftering is launched, the dimensional energy distribution figure of signal, Fig. 3 are the time domain waveform that receiving hydrophone receives, and wherein signal sampling rate is 50kHz.As seen from the figure, when hyperchannel is empty, liftering technology achieves space-focusing and the time domain impulse focusing of signal preferably.
For the validity of liftering technology in acoustic stimulation insertion loss is measured during checking hyperchannel sky, carry out basin test checking.Laboratory water tank is about 8m, wide 4m, depth of water 1.6m, and acoustic absorbant is posted at waveguide four sides, and bottom is covered with the fine sand of 0.25m.In experiment, ternary launches battle array with test sample at a distance of 7m.Test sample is steel plate, and its physical dimension is 1.1m × 1.0m × 5mm, and density is 7.84 × 10
3kg/m
3, velocity of sound 5470m/s.The insertion loss data acquisition and processing (DAP) of 0.5kHz ~ 20kHz frequency has been carried out in this test.What Fig. 4 showed is steel plate sample insertion loss theory and test findings comparison diagram, as seen from the figure, even if under the pond guided-wave conditions that multipath phenomenon is very serious, measurement result of the present invention is basic consistent with calculated value, and each relative error measuring frequency corresponding is all less than 10%.
Claims (2)
1., based on an acoustic stimulation insertion loss measuring method for liftering technology during hyperchannel sky, described measuring method comprises the steps:
1) when generation hyperchannel is empty, liftering transmits;
2) without the collection of direct signal during sample: export instrument by computer system control hyperchannel individual control signal, through multichannel power amplifier, utilize transducer array synchronized transmissions by step 1) generate empty time liftering signal, meanwhile, the direct signal p when hydrophone recordings of computer system control focal position does not test sample
i;
3) have the collection of transmission signal during sample: the sample laying acoustic stimulation is put into test environment, repeat 2) middle step, now record transmission focusing signal p during test sample
t;
4) calculating of insertion loss measured value: utilize the incoming signal p gathered
iand transmission signal p
t, calculate insertion loss measured value
2. as claimed in claim 1 a kind of empty based on hyperchannel time liftering technology acoustic stimulation insertion loss measuring method, it is characterized in that: described step 1) in, the process that when generation hyperchannel is empty, liftering transmits is: without in sample situation, each transducer of transducer array launches initialize signal successively, nautical receiving set is Received signal strength successively, if transmit, frequency domain can be expressed as s (f), and wherein f represents signal frequency, and Received signal strength is expressed as x
n(f), wherein n=1,2 ..., N represents that transducer is numbered, and N represents transducer number; According to formula
estimate the channel response function H between each transmission channel and receiving hydrophone
n(f); For avoiding the singular problem occurred in inversion calculation, utilize formula
estimate the inverse function obtaining channel response
wherein ε is a little value constant, unusual to avoid inverting; Inverse Fourier transform is done to this inverse function and obtains time domain Secondary Emission signal e
n(t), wherein t represents the time; Utilize the signal transmission time interval of each transmitting transducer to liftering signal e
nt () is carried out time delay and is evened up, the e thus the liftering of acquisition hyperchannel transmits
n(t-τ
n), wherein τ
nfor the delay volume that each transmission channel is corresponding.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510683611.5A CN105301114B (en) | 2015-10-20 | 2015-10-20 | A kind of acoustic stimulation insertion loss measuring method based on multichannel space-time liftering technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510683611.5A CN105301114B (en) | 2015-10-20 | 2015-10-20 | A kind of acoustic stimulation insertion loss measuring method based on multichannel space-time liftering technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105301114A true CN105301114A (en) | 2016-02-03 |
CN105301114B CN105301114B (en) | 2018-01-12 |
Family
ID=55198622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510683611.5A Active CN105301114B (en) | 2015-10-20 | 2015-10-20 | A kind of acoustic stimulation insertion loss measuring method based on multichannel space-time liftering technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105301114B (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562649A (en) * | 2018-03-21 | 2018-09-21 | 哈尔滨工程大学 | A kind of device and method measuring approximately level cylindrical shell backing coating low frequency insertion loss |
CN109186735A (en) * | 2018-09-21 | 2019-01-11 | 中国船舶工业系统工程研究院 | The measurement method of structural radiation acoustical power insertion loss based on vector hydrophone |
CN109211392A (en) * | 2018-09-21 | 2019-01-15 | 中国船舶工业系统工程研究院 | The measurement method of radiated insertion loss automation based on sound intensity probe |
CN109490425A (en) * | 2018-12-18 | 2019-03-19 | 浙江大学 | A kind of passive material bass reflex coefficient measuring method based on Green's function reconfiguration technique |
CN109696239A (en) * | 2018-12-03 | 2019-04-30 | 浙江大学 | Underwater sound passive material echo applied to pressure noise elimination water pot reduces method for automatic measurement |
CN109696238A (en) * | 2018-12-03 | 2019-04-30 | 浙江大学 | Underwater sound passive material insertion loss method for automatic measurement applied to pressure noise elimination water pot |
CN110399696A (en) * | 2019-08-01 | 2019-11-01 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | The double-current domain of the deposited acoustic stimulation in part couples three dimensional sound elasticity test verification method |
CN110988490A (en) * | 2019-12-20 | 2020-04-10 | 北京无线电计量测试研究所 | Power filter differential loss time domain measurement system and method |
CN111007152A (en) * | 2019-12-24 | 2020-04-14 | 无锡吉兴汽车声学部件科技有限公司 | Acoustic performance comprehensive evaluation method and system |
CN113091877A (en) * | 2021-04-07 | 2021-07-09 | 上海海洋大学 | Device and method for measuring target strength of pressure-controlled underwater acoustic scatterer |
CN115266942A (en) * | 2022-07-20 | 2022-11-01 | 哈尔滨工程大学 | Insertion loss measuring device and method for large-size sound insulation structure based on reverberation method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183150A (en) * | 2007-11-30 | 2008-05-21 | 浙江大学 | Time reversal transmitting focusing inverse beamforming method based on high-frequency time reversal array |
CN102539530A (en) * | 2012-01-12 | 2012-07-04 | 浙江大学 | Method for measuring echo reduction/reflection coefficient of underwater sound passive material based on time reversal focusing |
CN102590349A (en) * | 2012-01-12 | 2012-07-18 | 中国船舶工业集团公司船舶系统工程部 | Method for measuring insertion loss/ transmission coefficient of underwater sound passive material based on time reversal focusing |
CN104180891A (en) * | 2014-08-26 | 2014-12-03 | 中国船舶重工集团公司第七一五研究所 | Method for measuring ocean sound transmission loss in real time based on sound matrix |
CN104792877A (en) * | 2015-04-03 | 2015-07-22 | 浙江大学 | Method for measuring decoupling performance of underwater decoupling and vibration reducing material |
-
2015
- 2015-10-20 CN CN201510683611.5A patent/CN105301114B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101183150A (en) * | 2007-11-30 | 2008-05-21 | 浙江大学 | Time reversal transmitting focusing inverse beamforming method based on high-frequency time reversal array |
CN102539530A (en) * | 2012-01-12 | 2012-07-04 | 浙江大学 | Method for measuring echo reduction/reflection coefficient of underwater sound passive material based on time reversal focusing |
CN102590349A (en) * | 2012-01-12 | 2012-07-18 | 中国船舶工业集团公司船舶系统工程部 | Method for measuring insertion loss/ transmission coefficient of underwater sound passive material based on time reversal focusing |
CN104180891A (en) * | 2014-08-26 | 2014-12-03 | 中国船舶重工集团公司第七一五研究所 | Method for measuring ocean sound transmission loss in real time based on sound matrix |
CN104792877A (en) * | 2015-04-03 | 2015-07-22 | 浙江大学 | Method for measuring decoupling performance of underwater decoupling and vibration reducing material |
Non-Patent Citations (5)
Title |
---|
J.-F. AUBRY, M.ET AL: "Optimal focusing by spatio-temporal inverse filter.II. Experiments. Application to focusing through absorbing and reverberating media", 《THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERIA》 * |
M. TANTER.ET AL: "Optimal focusing by spatio-temporal inverse filter.I. Basic principles", 《THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》 * |
李建龙等: "基于时反算子分解的时反高分辨率定位技术研究", 《声学学报》 * |
李建龙等: "水声无源材料回声降低测量时反聚焦方法", 《声学学报》 * |
潘翔等: "基于时间反转处理的空—时逆滤波聚焦研究", 《传感技术学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108562649A (en) * | 2018-03-21 | 2018-09-21 | 哈尔滨工程大学 | A kind of device and method measuring approximately level cylindrical shell backing coating low frequency insertion loss |
CN108562649B (en) * | 2018-03-21 | 2020-09-11 | 哈尔滨工程大学 | Device and method for measuring low-frequency insertion loss of near-water-surface cylindrical shell backing covering layer |
CN109186735A (en) * | 2018-09-21 | 2019-01-11 | 中国船舶工业系统工程研究院 | The measurement method of structural radiation acoustical power insertion loss based on vector hydrophone |
CN109211392A (en) * | 2018-09-21 | 2019-01-15 | 中国船舶工业系统工程研究院 | The measurement method of radiated insertion loss automation based on sound intensity probe |
CN109696239A (en) * | 2018-12-03 | 2019-04-30 | 浙江大学 | Underwater sound passive material echo applied to pressure noise elimination water pot reduces method for automatic measurement |
CN109696238A (en) * | 2018-12-03 | 2019-04-30 | 浙江大学 | Underwater sound passive material insertion loss method for automatic measurement applied to pressure noise elimination water pot |
CN109490425B (en) * | 2018-12-18 | 2020-05-19 | 浙江大学 | Passive material low-frequency reflection coefficient measuring method based on Green function reconstruction technology |
CN109490425A (en) * | 2018-12-18 | 2019-03-19 | 浙江大学 | A kind of passive material bass reflex coefficient measuring method based on Green's function reconfiguration technique |
CN110399696A (en) * | 2019-08-01 | 2019-11-01 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | The double-current domain of the deposited acoustic stimulation in part couples three dimensional sound elasticity test verification method |
CN110399696B (en) * | 2019-08-01 | 2021-11-12 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | Double-current-domain coupling three-dimensional acoustoelastic test verification method for local acoustic covering layer |
CN110988490A (en) * | 2019-12-20 | 2020-04-10 | 北京无线电计量测试研究所 | Power filter differential loss time domain measurement system and method |
CN111007152A (en) * | 2019-12-24 | 2020-04-14 | 无锡吉兴汽车声学部件科技有限公司 | Acoustic performance comprehensive evaluation method and system |
CN113091877A (en) * | 2021-04-07 | 2021-07-09 | 上海海洋大学 | Device and method for measuring target strength of pressure-controlled underwater acoustic scatterer |
CN115266942A (en) * | 2022-07-20 | 2022-11-01 | 哈尔滨工程大学 | Insertion loss measuring device and method for large-size sound insulation structure based on reverberation method |
Also Published As
Publication number | Publication date |
---|---|
CN105301114B (en) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105301114B (en) | A kind of acoustic stimulation insertion loss measuring method based on multichannel space-time liftering technology | |
CN102539530B (en) | Method for measuring echo reduction/reflection coefficient of underwater sound passive material based on time reversal focusing | |
Sabra et al. | Arrival-time structure of the time-averaged ambient noise cross-correlation function in an oceanic waveguide | |
CN102590349B (en) | Method for measuring insertion loss/ transmission coefficient of underwater sound passive material based on time reversal focusing | |
Siderius et al. | Modeling broadband ocean acoustic transmissions with time-varying sea surfaces | |
CN102183435B (en) | Method for measuring submarine density and sound velocity based on multi-path reflection theory | |
CN104808208B (en) | Measurement system and method for detecting orientation and size of underwater target on basis of laser sound source | |
CN105181800B (en) | Acoustic stimulation echo based on multichannel space-time liftering technology reduces measuring method | |
Rouseff | Effect of shallow water internal waves on ocean acousticstriation patterns | |
CN101620004B (en) | Method for measuring direct wave signal sound pressure of sound wave transmitted in limited regional medium | |
CN105223271A (en) | Anti-during underwater sound passive material oise insulation factor under a kind of limited space conditions focus on measuring method | |
CN109001297B (en) | Method for measuring acoustic reflection coefficient of large-sample underwater acoustic material based on single-vector hydrophone | |
CN108680234A (en) | A kind of water-depth measurement method of quarice layer medium | |
Kargl et al. | Synthetic aperture sonar imaging of simple finite targets | |
CN103941250A (en) | Ocean channel simulation method in passive sonar array signal simulation | |
Jiang et al. | Geoacoustic inversion of broadband data by matched beam processing | |
CN106769711A (en) | The focusing distribution of sediment on-line measurement instrument of frequency sweep type ultrasonic | |
CN104280737B (en) | Weighted broadband time reversal operator resolution acoustic imaging method | |
CN115236592B (en) | Ice sound positioning method based on single-array-element time-frequency curve matching | |
Choi et al. | Measurement and simulation of the channel intensity impulse response for a site in the East China Sea | |
CN112834623B (en) | Acoustic covering layer echo reduction measurement method based on compressed sensing technology | |
CN109238440A (en) | The acoustic stimulation insertion loss measurement method of focusing technology when based on optimal sky | |
Li et al. | Measurement of acoustic properties for passive-material samples using multichannel inverse filter | |
CN109238441A (en) | It is a kind of based on it is optimal empty when focusing technology acoustic stimulation echo reduce measurement method | |
CN113126029A (en) | Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment |
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 |