CN110389327A - The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error - Google Patents
The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error Download PDFInfo
- Publication number
- CN110389327A CN110389327A CN201910689238.2A CN201910689238A CN110389327A CN 110389327 A CN110389327 A CN 110389327A CN 201910689238 A CN201910689238 A CN 201910689238A CN 110389327 A CN110389327 A CN 110389327A
- Authority
- CN
- China
- Prior art keywords
- formula
- receiving station
- away
- target
- biradical
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses the more external illuminators-based radars of multistation under a kind of receiving station's location error are biradical away from localization method.The metric data of the more external radiation source radar systems of multistation is grouped by the present invention according to receiving station, and the distance of selection target to receiving station is auxiliary variable, biradical away from the linearisation of non-linear measurement equation puppet by every group.Receiving station position is measured into noise and biradical incorporate in target location algorithm away from measurement noise statistics designs weight, establishes the cost function based on weighted least-squares.On this basis, being associated with as constraint condition, building constraint weighted least-squares location model using auxiliary variable and target position, and optimized by method of Lagrange multipliers.Various estimated result Weighted Fusions are finally obtained target position finally to estimate.The present invention reduces the complexity of positioning under the premise of guaranteeing to estimate performance, reduces influence of the error to target positioning performance.
Description
Technical field
The invention belongs to radar data process fields, and in particular to the more external sort algorithms of multistation under a kind of receiving station's location error
Radar is biradical away from localization method.
Background technique
External illuminators-based radar utilizes third party non-co-operation signal source (such as mobile communication signal, television broadcasting signal and enemy
Radar information etc.) detection target, signal and direct wave from target scattering third party's radiation emission are received by receiving station
Signal obtains the time delay and Doppler frequency difference of signal using Coherent processing.Time delay observation corresponds directly to signal from external sort algorithm
After target reflects reach receiving station propagation distance and, it is also known as biradical away from (Bistatic Range, BR).In multiple-input multiple-output body
It makes in lower external illuminators-based radar positioning system, biradical under multiple external sort algorithms-reception source (T/R) can be obtained away from fusion treatment
To the positioning realized to static target.
Currently, being widely studied using time delay/Doppler frequency difference target radiation source orientation problem.However, utilizing
The research of the non-cooperative location problem of BR/BRR is then relatively fewer.Godrich, H. team propose a kind of based on maximum likelihood
The direct localization method of estimation, the likelihood function of direct construction target position and speed obtain mesh by multi-dimensional grid search
Cursor position parameter Estimation.Noroozi etc. proposes one kind for the moving target orientation problem in distributed MIMO radar system
Two step weighted least square algorithms.Zhao Yongsheng etc. proposes a kind of three step weighted least square algorithms in succession.Most based on multistep weighting
Small two multiply estimation method when measurement noise is smaller, and positioning accuracy can reach a carat Metro lower bound, but when noise increases to one
When determining degree, position error be increased dramatically, and threshold effect occurs.
Above-mentioned estimation method assumes that the accurately known of the position of receiving station.And receiving station is normal in practical non-cooperative location
It is often installed on the motion platforms such as satellite, aircraft, naval vessels or surface car, position often can not accurately obtain, receiving station
There are certain errors for position.Receiving station's location error is very big on target location accuracy influence, and therefore, external illuminators-based radar positioning is asked
Receiving station's location error must be considered when topic research.Patent 201910048330.0 proposes a kind of receiving station's location error and places an order station
More external illuminators-based radars are biradical away from orientation problem, and it is biradical away from positioning that the present invention designs the more external illuminators-based radars of multistation under a kind of error
Algorithm.
Summary of the invention
The present invention is directed to the more external illuminators-based radars of multistation, considers influence of receiving station's location error to positioning performance, proposes
The more external illuminators-based radars of multistation are biradical away from positioning under receiving station's location error of the one kind based on constraint weighted least-squares (CWLS)
Method.The distance of the algorithms selection target to receiving station is auxiliary variable, is linearized biradical away from non-linear measurement equation puppet, will
Receiving station's position noise and biradical incorporate in target location algorithm away from measurement noise statistics design weight, establish based on weighting
The cost function of least square.On this basis, being associated with as constraint condition, building using auxiliary variable and target position
CWLS location model, and optimized by method of Lagrange multipliers.
Specific steps of the method for the invention are:
Each receiving station n, which is received, in the more external illuminators-based radar nets of step 1. multistation (M external sort algorithm and N number of receiving station) comes
From target scattering third side emitter m emit signal, obtain target away from external sort algorithm m and away from radar receiving station n it is biradical away from
Measurement information dm,n。
Measurement is divided into N group according to different receiving stations by step 2., and corresponding one group of each receiving station n (n=1 ..., N) double
Cardinal distance measures dm,n(m=1,2 ..., M).The distance of selection target to receiving station n are auxiliary variable Rn, will be non-linear biradical away from amount
Pseudo- linearisation is surveyed, estimation equation is obtained
Step 3. selects least-square residuals quadratic sum for cost function, will measure noise statistics and incorporates location algorithm
Middle design weight constructs cost functionWeightConsider auxiliary
Variable RnWith the correlation of target position, structure constraint condition θn TΩθn=0, the constraint weighted least-squares for constructing quadratic form are estimated
Meter problem.
Step 4. utilizes Lagrange relaxation, converts unconfined optimization problem for constrained optimization problemOptimization Solution
Step 5. utilizes minimum mean square error criterion, and each group of receiving station is measured the target position estimated value obtained and is melted
It closes, obtains globally optimal solution.
Beneficial effects of the present invention:
1. by introduce auxiliary variable, rationally by it is non-linear it is biradical be converted into pseudo square-free function model away from measurement model, obtain
Obtain the enclosed analytic solutions of target position estimation.
2. according to receiving station's location error and it is biradical away from measure Noise Design optimizing index weight, to reduce error to mesh
The influence of positioning performance is marked, target position estimated accuracy is improved.
3. consider the relevance between auxiliary variable and unknown variable, structure constraint condition, design constraint weighting minimum two
Multiplication algorithm further decreases evaluated error.
Specific embodiment
It is biradical away from localization method, the party that the present invention devises the more external illuminators-based radars of multistation under a kind of receiving station's location error
Method the following steps are included:
Step 1: in the more external illuminators-based radar nets of multistation, it is assumed that the transmitting station You Ge M, m emit station locationN number of receiving station, the actual position of n-th of receiving stationActually due to depositing
It is disturbing, the actual position of receiving station is unable to get, and be can only obtain the position comprising noise and isAnd For the position error vector of receiving station n, and assume to be independent Gauss zero-mean white noise,
Covariance isTarget position to be estimated is Starget=[x, y, z]T, then target away from external sort algorithm m and away from
The biradical of receiving station n be away from measurement information
In formula, dm,nFor the biradical measured value away from measurement;| | | | it is Euclidean distance, the distance of target to receiving station n areDistance of the target to external sort algorithm mem,nTo be biradical away from error in measurement, obey
Zero-mean gaussian distribution;
Measurement is divided into N group according to different receiving stations by step 2., and corresponding one group of each receiving station n biradical away from measurement dm,n;
Formula (1) is subjected to pseudo- linearisation
In formula, Measurement distance for target away from receiving station n, For corresponding gradient,
Write formula (2) as matrix form
In formula,en=[e1,n e2,n … eM,n]T
Step 3. establishes constraint weighted least square model;
Step 3.1. selects least-square residuals quadratic sum for cost function J (θn)
Noise statistics will be measured and incorporate design weight W in location algorithmn
In formula, n-th group is biradical away from error in measurement covarianceReceiving station n location error covariance
Step 3.2. considers auxiliary variableWith the correlation of target position, constraint condition between the two is established
Then have
θn TΩθn=0 (7)
In formula,Ω=diag [1 1 1-1];
The constraint weighted least square problem that step 3.3. constructs quadratic form is as follows
Step 4. utilizes Lagrange relaxation, and Optimization Solution constrains weighted least square problem;
Step 4.1. introduces Lagrange multiplier λn, unconstrained optimization problem is converted by constrained optimization problem, that is, formula (8),
Establish Lagrangian
Step 4.2. solves parameter θ to be estimated using Lagrangian Relaxation Algorithmn;
To LagrangianL (θn,λn) local derviation is sought, enable partial derivativeIt is zero, can obtains
Due to λnIt is unknown, formula (10) are substituted into formula (7), can be obtained
It, will be in formula (11) using Eigenvalues Decomposition methodDiagonalization
In formula, Λn=diag { γ1,…,γ4, and γi(i=1 ..., 4) is characteristic value, UnFor corresponding eigenvalue composition
Feature vector;
Formula (12) are substituted into formula (11), can be obtained
In formula,
By convolution and polynomial rooting operation, λ is acquirednMultiple;Choose the λ of real rootnSubstitution formula (10), can be in the hope of
Several θ outn, by θnValue substitutes into cost function formula (4), selects cost function J (θn) the smallest θnValue
Step 4.3. is worth available target locator value according to optimal estimation are as follows:
In formula,For based on n-th group target position estimated value,For estimated valueFirst three items;
Step 5. utilizes minimum mean square error criterion, and the target position estimated value obtained is measured to each group of receiving stationWeighted Fusion obtains the globally optimal solution of target position;
Step 5.1. calculates n-th group and measures lower estimated valueCovariance;
By HnH is decomposed by columnn=[Hn,1:3,Hn,4], wherein Hn,1:3And Hn,4Respectively indicate Hn1~3 column of matrix and the 4th
Column;Therefore cost function formula (4) indicates again are as follows:
In formula, gn=Hn,1:3vn+Hn,4(vn Tvn)1/2-Zn,
To gnThe small noise disturbance analysis of single order is carried out, can be obtained
In formula, Gn=Hn,1:3vn+Hn,4vn T(vn Tvn)-1/2;
N-th group can be obtained as a result, measures lower estimated valueCovariance be
Step 5.2. hypothesis is biradical uncorrelated away from the estimation for measuring lower target position according to N group, according to unbiased lowest mean square
Poor criterion obtains final global optimization value to N group target state estimator value Weighted Fusion
Claims (1)
1. the more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error, it is characterised in that: this method includes
Following steps:
Step 1: in the more external illuminators-based radar nets of multistation, it is assumed that the transmitting station You Ge M, m-th of transmitting station locationN number of receiving station, the actual position of n-th of receiving stationActually due to depositing
It is disturbing, the actual position of receiving station is unable to get, and be can only obtain the position comprising noise and isAnd For the position error vector of receiving station n, and assume to be independent Gauss zero-mean white noise,
Covariance isTarget position to be estimated is Starget=[x, y, z]T, then target away from external sort algorithm m and away from
The biradical of receiving station n be away from measurement information
In formula, dm,nFor the biradical measured value away from measurement;| | | | it is Euclidean distance, the distance of target to receiving station n areDistance of the target to external sort algorithm mem,nTo be biradical away from error in measurement, obey
Zero-mean gaussian distribution;
Measurement is divided into N group according to different receiving stations by step 2., and corresponding one group of each receiving station n biradical away from measurement dm,n;By formula
(1) pseudo- linearisation is carried out
In formula, Measurement distance for target away from receiving station n, For corresponding gradient,
Write formula (2) as matrix form
In formula,en=[e1,n e2,n … eM,n]T
Step 3. establishes constraint weighted least square model;
Step 3.1. selects least-square residuals quadratic sum for cost function J (θn)
Noise statistics will be measured and incorporate design weight W in location algorithmn
In formula, n-th group is biradical away from error in measurement covarianceReceiving station n location error covariance
Step 3.2. considers auxiliary variableWith the correlation of target position, constraint condition between the two is established
Then have
θn TΩθn=0 (7)
In formula,Ω=diag [1 1 1-1];
The constraint weighted least square problem that step 3.3. constructs quadratic form is as follows
Step 4. utilizes Lagrange relaxation, and Optimization Solution constrains weighted least square problem;
Step 4.1. introduces Lagrange multiplier λn, unconstrained optimization problem is converted by constrained optimization problem, that is, formula (8), is established
Lagrangian
L(θn,λn)=(Zn-Hnθn)TWn -1(Zn-Hnθn)+λnθn TΩθn (9)
Step 4.2. solves parameter θ to be estimated using Lagrangian Relaxation Algorithmn;
To LagrangianL (θn,λn) local derviation is sought, enable partial derivativeIt is zero, can obtains
Due to λnIt is unknown, formula (10) are substituted into formula (7), can be obtained
It, will be in formula (11) using Eigenvalues Decomposition methodDiagonalization
In formula, Λn=diag { γ1,…,γ4, and γi(i=1 ..., 4) is characteristic value, UnFor the spy of corresponding eigenvalue composition
Levy vector;
Formula (12) are substituted into formula (11), can be obtained
In formula,
By convolution and polynomial rooting operation, λ is acquirednMultiple;Choose the λ of real rootnSubstitution formula (10), if can find out
Dry θn, by θnValue substitutes into cost function formula (4), selects cost function J (θn) the smallest θnValue
Step 4.3. is worth available target locator value according to optimal estimation are as follows:
In formula,For based on n-th group target position estimated value,For estimated valueFirst three items;
Step 5. utilizes minimum mean square error criterion, and the target position estimated value obtained is measured to each group of receiving stationWeighting
Fusion, obtains the globally optimal solution of target position;
Step 5.1. calculates n-th group and measures lower estimated valueCovariance;
By HnH is decomposed by columnn=[Hn,1:3,Hn,4], wherein Hn,1:3And Hn,4Respectively indicate Hn1~3 column and the 4th column of matrix;Cause
This cost function formula (4) indicates again are as follows:
In formula, gn=Hn,1:3vn+Hn,4(vn Tvn)1/2-Zn,
To gnThe small noise disturbance analysis of single order is carried out, can be obtained
In formula, Gn=Hn,1:3vn+Hn,4vn T(vn Tvn)-1/2;
N-th group can be obtained as a result, measures lower estimated valueCovariance be
Step 5.2. hypothesis is biradical uncorrelated away from the estimation for measuring lower target position according to N group, quasi- according to unbiased Minimum Mean Square Error
Then, final global optimization value is obtained to N group target state estimator value Weighted Fusion
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910689238.2A CN110389327A (en) | 2019-07-29 | 2019-07-29 | The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910689238.2A CN110389327A (en) | 2019-07-29 | 2019-07-29 | The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110389327A true CN110389327A (en) | 2019-10-29 |
Family
ID=68287834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910689238.2A Pending CN110389327A (en) | 2019-07-29 | 2019-07-29 | The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110389327A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111123197A (en) * | 2019-12-21 | 2020-05-08 | 杭州电子科技大学 | TDOA-based target radiation source positioning method |
CN111308530A (en) * | 2020-02-17 | 2020-06-19 | 中国人民解放军战略支援部队信息工程大学 | Short wave multi-station and single-satellite cooperative direct positioning method based on two-dimensional direction of arrival |
CN111371432A (en) * | 2020-03-24 | 2020-07-03 | 宁波飞拓电器有限公司 | Noise-related nonlinear observability degree analysis method |
CN111983561A (en) * | 2020-06-30 | 2020-11-24 | 江西锐迪航空科技发展有限公司 | TDOA (time difference of arrival) positioning method for multiple unmanned aerial vehicle targets under receiver position error |
CN112526449A (en) * | 2020-11-27 | 2021-03-19 | 中国人民解放军海军工程大学 | Method for calibrating position information of receiving station by utilizing moving target |
CN112683265A (en) * | 2021-01-20 | 2021-04-20 | 中国人民解放***箭军工程大学 | MIMU/GPS integrated navigation method based on rapid ISS collective filtering |
CN115508774A (en) * | 2022-10-12 | 2022-12-23 | 中国电子科技集团公司信息科学研究院 | Time difference positioning method and device based on two-step weighted least square and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105091884A (en) * | 2014-05-08 | 2015-11-25 | 东北大学 | Indoor moving robot route planning method based on sensor network dynamic environment monitoring |
CN106405533A (en) * | 2016-08-30 | 2017-02-15 | 西安电子科技大学 | Radar target combined synchronization and positioning method based on constraint weighted least square |
US20170286730A1 (en) * | 2016-04-04 | 2017-10-05 | Mojix, Inc. | Location Estimation and Tracking for Passive RFID and Wireless Sensor Networks Using MIMO Systems |
CN109633592A (en) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | The external illuminators-based radar time difference and frequency difference co-located method under movement observations station error |
CN109633591A (en) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | External illuminators-based radar is biradical away from localization method under a kind of observation station location error |
-
2019
- 2019-07-29 CN CN201910689238.2A patent/CN110389327A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105091884A (en) * | 2014-05-08 | 2015-11-25 | 东北大学 | Indoor moving robot route planning method based on sensor network dynamic environment monitoring |
US20170286730A1 (en) * | 2016-04-04 | 2017-10-05 | Mojix, Inc. | Location Estimation and Tracking for Passive RFID and Wireless Sensor Networks Using MIMO Systems |
CN106405533A (en) * | 2016-08-30 | 2017-02-15 | 西安电子科技大学 | Radar target combined synchronization and positioning method based on constraint weighted least square |
CN109633592A (en) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | The external illuminators-based radar time difference and frequency difference co-located method under movement observations station error |
CN109633591A (en) * | 2019-01-18 | 2019-04-16 | 杭州电子科技大学 | External illuminators-based radar is biradical away from localization method under a kind of observation station location error |
Non-Patent Citations (1)
Title |
---|
周成等: "一种利用MIMO雷达的CWLS目标定位算法", 《西安电子科技大学学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111123197A (en) * | 2019-12-21 | 2020-05-08 | 杭州电子科技大学 | TDOA-based target radiation source positioning method |
CN111308530A (en) * | 2020-02-17 | 2020-06-19 | 中国人民解放军战略支援部队信息工程大学 | Short wave multi-station and single-satellite cooperative direct positioning method based on two-dimensional direction of arrival |
CN111371432A (en) * | 2020-03-24 | 2020-07-03 | 宁波飞拓电器有限公司 | Noise-related nonlinear observability degree analysis method |
CN111371432B (en) * | 2020-03-24 | 2024-05-31 | 宁波飞拓电器有限公司 | Nonlinear observability analysis method with noise correlation |
CN111983561A (en) * | 2020-06-30 | 2020-11-24 | 江西锐迪航空科技发展有限公司 | TDOA (time difference of arrival) positioning method for multiple unmanned aerial vehicle targets under receiver position error |
CN112526449A (en) * | 2020-11-27 | 2021-03-19 | 中国人民解放军海军工程大学 | Method for calibrating position information of receiving station by utilizing moving target |
CN112526449B (en) * | 2020-11-27 | 2022-12-27 | 中国人民解放军海军工程大学 | Method for calibrating position information of receiving station by utilizing moving target |
CN112683265A (en) * | 2021-01-20 | 2021-04-20 | 中国人民解放***箭军工程大学 | MIMU/GPS integrated navigation method based on rapid ISS collective filtering |
CN115508774A (en) * | 2022-10-12 | 2022-12-23 | 中国电子科技集团公司信息科学研究院 | Time difference positioning method and device based on two-step weighted least square and storage medium |
CN115508774B (en) * | 2022-10-12 | 2023-07-28 | 中国电子科技集团公司信息科学研究院 | Time difference positioning method and device based on two-step weighted least square and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110389327A (en) | The more external illuminators-based radars of multistation are biradical away from localization method under receiving station's location error | |
CN109633592B (en) | External radiation source radar time difference and frequency difference co-location method under motion observation station error | |
CN106610483B (en) | MIMO radar angle estimation algorithm based on tensor subspace and spectrum peak search | |
Noroozi et al. | Weighted least squares target location estimation in multi‐transmitter multi‐receiver passive radar using bistatic range measurements | |
CN107015213B (en) | Bistatic MIMO radar angle evaluation method based on MUSIC algorithm | |
CN109633591A (en) | External illuminators-based radar is biradical away from localization method under a kind of observation station location error | |
CN108387876B (en) | External radiation source radar network double-base-distance error registration method based on CTLS | |
CN110389326A (en) | The more external illuminators-based radar moving target localization methods of multistation under a kind of reception station error | |
CN107991659B (en) | Method for measuring height of low-elevation target of meter wave radar based on dictionary learning | |
CN109283522A (en) | A kind of co-located MIMO radar method for tracking target of joint spatial-temporal resource management | |
CN108957387A (en) | A kind of satellite-signal two-dimentional angle estimation method and system | |
CN104950297A (en) | Array element error estimation method based on matrix 1-norm fitting | |
CN109683151A (en) | Tenth of the twelve Earthly Branches rooting MUSIC angle estimating method under non-uniform noise environment based on matrix completion | |
CN102213766A (en) | Method and device for avoiding multi-path errors in satellite navigation receiver | |
CN103852749A (en) | Robust waveform optimization method for improving MIMO-STAP detection performance | |
CN111107626A (en) | DOA positioning method based on time reversal | |
Liu et al. | Target localization in high-coherence multipath environment based on low-rank decomposition and sparse representation | |
CN105445732A (en) | Object track initialization method using multipath observation under dense clutter condition | |
CN108761384B (en) | Target positioning method for robust sensor network | |
CN114690175B (en) | Target direct detection and tracking method based on passive external radiation source radar | |
CN112986975B (en) | Distance weighting-based passive radar network centralized detection method | |
CN105204018B (en) | A kind of two-dimentional DOA trackings using many frame informations | |
Li et al. | Two-step Bayesian detection for MIMO radar in compound-Gaussian clutter with Gamma texture | |
Tang et al. | Bistatic MIMO radar height estimation method based on adaptive beam-space RML data fusion | |
Kilpatrick et al. | Sidelobe suppression and super resolution for MIMO imaging radar |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191029 |