CN110412504A - It is associated with based on angle with the passive track-corelation of time difference information and localization method - Google Patents
It is associated with based on angle with the passive track-corelation of time difference information and localization method Download PDFInfo
- Publication number
- CN110412504A CN110412504A CN201910738048.5A CN201910738048A CN110412504A CN 110412504 A CN110412504 A CN 110412504A CN 201910738048 A CN201910738048 A CN 201910738048A CN 110412504 A CN110412504 A CN 110412504A
- Authority
- CN
- China
- Prior art keywords
- target
- observation station
- angle
- error
- cost matrix
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/0257—Hybrid positioning
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/10—Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements, e.g. omega or decca systems
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/14—Determining absolute distances from a plurality of spaced points of known location
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention belongs to electronic countermeasure technology fields, and in particular to a kind of to be associated with based on angle with the passive track-corelation of time difference information and localization method.Method of the invention mainly uses multi radiation sources data association algorithm, is associated to obtain cost matrix C by the angular observation data of multiple observation stations and time difference observation data simultaneously respectively1(p, q) and C2(p, q), then it is utilized respectively angle measurement and the standard deviation of range difference measurement value error obtains total cost matrix C (p, q), C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.Beneficial effects of the present invention are that the present invention accurately can complete to be associated with to multistation data, and finally accurately estimate the position of target, and method is simple, works well.
Description
Technical field
The invention belongs to electronic countermeasure technology fields, and in particular to a kind of passive multistation based on angle and time difference information is more
Target association and localization method.
Background technique
In electronic countermeasure field, due to the development of perturbation technique, active location technology poor anti jamming capability lacks hidden
Property weakness start gradually to become obvious, and passive location technology has very strong superiority relative to active location technology.
It has many advantages, such as remote operating distance, lobe-on-receive, strong antijamming capability and is not easy to be realized by other side, for improving system in electricity
Survival ability and fighting efficiency under son war environment have a very important role.Passive location mainly passes through radiation emission
Electromagnetic wave signal reaches the time difference of observation station obtaining the azimuth of radiation source, pitch angle or electromagnetic wave and realizes.Therefore
Passive location generally requires multiple base station collaborations and completes.But when there are multiple targets, passive multistation location is easy to appear void
Decoy point causes to position mistake.So carrying out data correlation for Multiple targets observation data, false target point is rejected for standard
Determine that position is of great significance.
In passive location system, direction cross positioning is presented as that a plurality of straight line crosses in space, be with most extensively,
A kind of longest localization method of search time.But meanwhile it there is also some problems, such as in a practical situation it measurement miss
Difference is larger, can have biggish position error particularly with long distance positioning;If it is in the case where carrying out Multi-target position,
Multiple observation stations can generate a large amount of false crosspoint.And positioning using TDOA is presented as that multiple hyperboloids cross in space, with it
Anti-interference, the advantages that Yi Shixian and precision are high, occupies more and more important position in passive location.Positioning using TDOA is to pass through
Measurement emitter Signals reach the time difference of different observation stations to determine target position, need multiple base stations to cooperate with, therefore also hold
Easily there is false target point.So angle information and time difference information joint are carried out data correlation, for realizing that target is accurately fixed
Position has great research significance.
Summary of the invention
In view of the above-mentioned problems, being closed the invention proposes a kind of based on the passive track-corelation of angle information and time difference information
Connection and localization method.
The technical solution adopted by the present invention is that:
It is a kind of that the positioning of track-corelation is realized based on being associated with localization method for angle information and time difference information, it is described
Be in the positioning that crossed using angle and the straight line of the progress of time difference information and hyperbola, to false target reject to
Realize the data correlation of multi-source data.The track-corelation location model of this method is as shown in Figure 1.
It has been known that there is 3 observation stations, and wherein observation station 1 is main website, and observation station 2 and observation station 3 are slave station, their coordinate point
Not are as follows: (x1,y1)、(x2,y2)、(x3,y3), the coordinate of target 1 and target 2 is unknown.Assuming that 3 observation stations observe target respectively
1 azimuth is θ1,1、θ2,1、θ3,1If the coordinate of radiation source 1 is (xp,yq), then by formula:
When only existing single target, the coordinate of radiation source 1 can be found out by azimuth information, wherein niIt is i-th
The angle measurement error of observation station.
If being positioned using time difference information to single target, it may be assumed that when known to 1 position of target, target 1 reaches
The distance between main website 1 and slave station 2,3 difference are respectively △ d21、△d31, according to hyp definition as can be seen that target 1 is located at
It is fixed point with main website and slave station, is that fixed length is formed by hyperbola with the distance between the fixed point of target 1 to two difference.When target position
When setting unknown, it is assumed that the coordinate of unknown object 2 is (x in Fig. 1p,yq), with main website 1, the distance of slave station 2 and slave station 3 is distinguished
For Indicate target 2 reach the distance between main website and slave station k it is poor, then just like
Lower equation:
When being related to single goal positioning, above equation can directly calculate the position (x of unknown object 2p,yq), target position
Set be two-by-two observation station the intersection point on hyperbola is formed by using range difference as fixed length.
Above method can only be positioned with single target, when there are multiple targets, each observation
The problem of azimuth information or time difference information of multiple targets can all be observed by standing, this relates to data correlation.For example sentence
Which angle information is the corresponding same target in disconnected each observation station out, could more accurately be positioned in this way, otherwise
The influence of false target may be will receive, to cause positioning mistake.Since angle measurement error is opposite in actual measurement
It is larger, if will cause biggish position error in the farther away situation in target range observation station, therefore, in order to ensure that target positions
Accuracy, the present invention use angle information and the united method of time difference information, association algorithm flow chart as shown in Fig. 2, tool
Body method is as follows:
Assuming that observation station and target are all located on X/Y plane, if main website position coordinates (x1,y1), the coordinate of slave station is (xk,
yk), wherein k=2,3 ..., K, K are observation station sum.Then there is j-th of target to the distance of k-th observation station and arrives main website
Apart from its difference are as follows:Wherein k=2,3 ..., K, j=1,2 ..., J, J be target sum,Table
Show measurement error.It enablesFor the azimuth information of j-th of target of k-th of observation station measurement, wherein k=1,
2 ..., K, j=1,2 ..., J, θk,jIndicate k-th of observation station to j-th of target azimuth information true value, nk,jTable
Show measurement noise.
Objective plane is divided into the grid of P × Q range, each mesh point represents a position coordinates in objective plane
(xp,yq), wherein p=1,2 ..., P, q=1,2 ... Q traverses each mesh point in grid plan:
1, point (x is calculatedp,yq) azimuth relative to each observation station:Meanwhile
Calculate this for the distance of main website and slave station and they the distance between it is poor:Its
Middle k=2,3 ..., K.
2, each Searching point (x is calculatedp,yq) azimuth relative to observation stationThe side observed with observation station
Parallactic angleBetween errorWherein k=1,2 ..., K, j=1,2 ..., J.
Range difference of this relative to observation station is calculated simultaneouslyThe range difference observed with observation stationBetween
ErrorWherein k=2,3 ..., K, j=1,2 ..., J.
3, the cost matrix being made of the angle overall error searched for every time is calculatedWherein
K=1,2 ..., K, p=1,2 ..., P, q=1,2 ..., Q and the cost matrix being made of range difference overall errorWherein k=2,3 ..., K, p=1,2 ..., P, q=1,2 ..., Q.
4, combine the cost matrix C based on angle1With the cost matrix C based on range difference2Obtain total cost matrix C (p,
Q)=1/ σ1*C1+1/σ2*C2, wherein σ1And σ2Respectively represent the standard of observation station angle measurement and range difference measurement value error
Difference.
5, C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.
Beneficial effects of the present invention are that can accurately be associated to multistation data, and finally accurately estimate target
Position, method is simple, works well.
Detailed description of the invention
Fig. 1 is multistation angle and time difference Model for joint positioning figure;
Fig. 2 is multi radiation sources data association algorithm flow chart;
Fig. 3 is the Multi-target position puppet spectral peak figure based on angle;
Fig. 4 is the Multi-target position puppet spectral peak figure based on angle and the time difference;
Fig. 5 is positioning result figure;
Fig. 6 is correct localization with threshold variation curve graph.
Specific embodiment
The present invention is described in detail below with reference to specific example:
It is fixed that this example intersects the above-mentioned passive track-corelation direction finding based on angle information and time difference information using matlab
Position method is verified, and for simplicity, is made the following assumptions to algorithm model:
1. main website, slave station and target are all in X/Y plane;
2. main website and slave station range accuracy having the same;
3. all Engineering Errors are all added in range error;
4. assuming that target is static or movement velocity is extremely low.
Assuming that target area is the rectangular region of 100km × 80km, 500 Monte Carlo emulation, 4 mesh are carried out
Target coordinate is respectively (15,70), (35,75), (58,75), (85,70), unit km.Use 1 main website and 3 slave stations pair
Above-mentioned target is positioned, and the coordinate of main website and slave station is respectively (5,30), (30,15), (70,15), (95,30), and unit is
Km, the range error at all stations obey the Gaussian Profile that mean value is zero, and mutually indepedent between the range error at all stations.500
Multistation data correlation effect and locating effect are obtained after secondary Monte Carlo emulation.
Track-corelation data correlation effect:
To compare direction cross positioning method based on angle information and being closed based on the target of angle information and time difference information
Connection localization method, in addition sets (10,50) for target, (13,50), (50,30), (85,60), wherein the first two target position
Relatively.It is to be obtained after the angle information of each target observed to multiple observation stations is associated as shown in Figure 3
Target puppet spectral peak figure, as can be seen from the figure since two of them target range closely causes very much corresponding two spectrums cutting edge of a knife or a sword to become one
Spectrum cutting edge of a knife or a sword to can not be accurately positioned, Fig. 4 be in situation identical with the target position Fig. 3 to the angle information of each target and when
What poor information obtained after being associated has the pseudo- spectral peak figure of 4 peak values, Fig. 3 is compared, in Fig. 4 corresponding to 4 target positions
It is high-visible to compose cutting edge of a knife or a sword, it follows that the object localization method based on angle information and time difference information is excellent for the resolution ratio of distance
In the direction cross positioning method based on angle information.
For the locating effect for more intuitively showing algorithm, Fig. 5 gives the 100th locating effect.In figure, border circular areas
Correct localization region is represented, estimated location is located in border circular areas then identification and correct positioning has been carried out, can from Fig. 5
To find out, estimated location and target actual position very close to, it is believed that accurate target position estimation may be implemented in the algorithm.
Error minimum value, max value of error and average error representated by three curves respectively indicate multiple in Fig. 6
The minimum value of the distance between target state estimator position and actual position, maximum value and to multiple target state estimator positions and true position
The distance between set the average value that summation is averaged.Thresholding is the judgment criteria of target positioning correctness, when target is estimated
It indicates that positioning is correct when counting position and actual position distance lower than thresholding, and then indicates to position mistake beyond thresholding.Accuracy is
Refer to and positions the percentage that correct number accounts for total degree in 500 positioning carried out in the case where some determines thresholding.It can be with by Fig. 6
Find out, correct localization is gradually promoted with the increase of thresholding, the energy when actual position and Estimated Position Error are in 3km or so
Reach very high correct localization, it can be seen that the passive track-corelation direction finding based on angle information and time difference information is handed over
Fork localization method works well.
Claims (1)
1. being associated with based on angle with the passive track-corelation of time difference information and localization method, observation station and target is enabled to be all located at XY
In plane, main website position coordinates (x1,y1), the coordinate of slave station is (xk,yk), wherein k=2,3 ..., K, K are observation station sum;
There is distance of j-th of target to k-th of observation stationWith the distance for arriving main websiteDifference are as follows:Wherein
J=1,2 ..., J, J are target sums,It indicates measurement error, enablesFor the jth of k-th of observation station measurement
The azimuth information of a target, wherein k=1,2 ..., K, j=1,2 ..., J, θk,jIndicate k-th of observation station to j-th of mesh
The true value of target azimuth information, nk,jIndicate measurement noise.Objective plane is divided into the grid of P × Q range, each net
Lattice point represents a position coordinates (x in objective planep,yq), wherein p=1,2 ..., P, q=1,2 ... Q, traversal grid are flat
Each mesh point in face;It is characterized in that, the described method comprises the following steps:
S1, Searching point (x is calculatedp,yq) azimuth relative to each observation station:Meanwhile it calculating
This for the distance of main website and slave station and they the distance between it is poor:
S2, each Searching point (x is calculatedp,yq) azimuth relative to observation stationThe azimuth observed with observation stationBetween errorRange difference of this relative to observation station is calculated simultaneouslyWith
The range difference that observation station is observedBetween error
The cost matrix that S3, calculating are made of the angle overall error searched for every timeAnd by away from
The cost matrix of deviation overall error composition
S4, combine the cost matrix C based on angle1With the cost matrix C based on range difference2Obtain total cost matrix C (p, q)
=1/ σ1*C1+1/σ2*C2, wherein σ1And σ2Respectively represent the standard deviation of observation station angle measurement and range difference measurement value error;
S5, C'=1./C is enabled, the position by solving the local maximum that C' corresponds to number of targets positions target.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910738048.5A CN110412504A (en) | 2019-08-12 | 2019-08-12 | It is associated with based on angle with the passive track-corelation of time difference information and localization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910738048.5A CN110412504A (en) | 2019-08-12 | 2019-08-12 | It is associated with based on angle with the passive track-corelation of time difference information and localization method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110412504A true CN110412504A (en) | 2019-11-05 |
Family
ID=68367061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910738048.5A Pending CN110412504A (en) | 2019-08-12 | 2019-08-12 | It is associated with based on angle with the passive track-corelation of time difference information and localization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110412504A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111308418A (en) * | 2020-03-10 | 2020-06-19 | 慧众行知科技(北京)有限公司 | Steady method for two-dimensional positioning of target with unknown height |
CN111337875A (en) * | 2020-04-13 | 2020-06-26 | 北京航天长征飞行器研究所 | Space multi-station multi-system dimension reduction positioning method |
CN113390406A (en) * | 2021-06-16 | 2021-09-14 | 电子科技大学 | Multi-target data association and positioning method based on passive multi-sensor system |
CN115524662A (en) * | 2022-10-27 | 2022-12-27 | 中国电子科技集团公司信息科学研究院 | Direction finding time difference combined positioning method and system, electronic equipment and storage medium |
RU2790347C1 (en) * | 2022-04-14 | 2023-02-16 | Общество с ограниченной ответственностью "Специальный Технологический Центр" | Difference-range method for determining the location of objects |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108061877A (en) * | 2017-12-14 | 2018-05-22 | 电子科技大学 | A kind of passive track-corelation direction cross positioning method based on angle information |
CN108089148A (en) * | 2017-12-14 | 2018-05-29 | 电子科技大学 | A kind of passive track-corelation direction cross positioning method based on time difference information |
CN108152789A (en) * | 2018-01-03 | 2018-06-12 | 电子科技大学 | Utilize the passive track-corelation data correlation and localization method of RSS information |
CN108802674A (en) * | 2018-07-19 | 2018-11-13 | 中国人民解放军战略支援部队信息工程大学 | It is a kind of for the combined method for searching and device that directly position |
CN109001670A (en) * | 2018-06-01 | 2018-12-14 | 中国人民解放军战略支援部队信息工程大学 | A kind of distributed passive location method and device for combining the time difference and angle |
CN109917333A (en) * | 2019-01-31 | 2019-06-21 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Merge the passive location method of AOA observed quantity and TDOA observed quantity |
-
2019
- 2019-08-12 CN CN201910738048.5A patent/CN110412504A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108061877A (en) * | 2017-12-14 | 2018-05-22 | 电子科技大学 | A kind of passive track-corelation direction cross positioning method based on angle information |
CN108089148A (en) * | 2017-12-14 | 2018-05-29 | 电子科技大学 | A kind of passive track-corelation direction cross positioning method based on time difference information |
CN108152789A (en) * | 2018-01-03 | 2018-06-12 | 电子科技大学 | Utilize the passive track-corelation data correlation and localization method of RSS information |
CN109001670A (en) * | 2018-06-01 | 2018-12-14 | 中国人民解放军战略支援部队信息工程大学 | A kind of distributed passive location method and device for combining the time difference and angle |
CN108802674A (en) * | 2018-07-19 | 2018-11-13 | 中国人民解放军战略支援部队信息工程大学 | It is a kind of for the combined method for searching and device that directly position |
CN109917333A (en) * | 2019-01-31 | 2019-06-21 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Merge the passive location method of AOA observed quantity and TDOA observed quantity |
Non-Patent Citations (2)
Title |
---|
JIHAO YIN: "A simple and accurate TDOA- AOA localization method using two stations", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING LETTERS》 * |
熊健: "一种改进的多站AOA/TDOA联合无源定位闭合形式算法", 《电讯技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111308418A (en) * | 2020-03-10 | 2020-06-19 | 慧众行知科技(北京)有限公司 | Steady method for two-dimensional positioning of target with unknown height |
CN111308418B (en) * | 2020-03-10 | 2021-11-23 | 慧众行知科技(北京)有限公司 | Steady method for two-dimensional positioning of target with unknown height |
CN111337875A (en) * | 2020-04-13 | 2020-06-26 | 北京航天长征飞行器研究所 | Space multi-station multi-system dimension reduction positioning method |
CN113390406A (en) * | 2021-06-16 | 2021-09-14 | 电子科技大学 | Multi-target data association and positioning method based on passive multi-sensor system |
CN113390406B (en) * | 2021-06-16 | 2022-05-24 | 电子科技大学 | Multi-target data association and positioning method based on passive multi-sensor system |
RU2790347C1 (en) * | 2022-04-14 | 2023-02-16 | Общество с ограниченной ответственностью "Специальный Технологический Центр" | Difference-range method for determining the location of objects |
CN115524662A (en) * | 2022-10-27 | 2022-12-27 | 中国电子科技集团公司信息科学研究院 | Direction finding time difference combined positioning method and system, electronic equipment and storage medium |
CN115524662B (en) * | 2022-10-27 | 2023-09-19 | 中国电子科技集团公司信息科学研究院 | Direction finding time difference joint positioning method, system, electronic equipment and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110412504A (en) | It is associated with based on angle with the passive track-corelation of time difference information and localization method | |
CN108089148B (en) | A kind of passive track-corelation direction cross positioning method based on time difference information | |
CN108061877B (en) | A kind of passive track-corelation direction cross positioning method based on angle information | |
CN110988923B (en) | Satellite navigation interference source passive positioning method based on multi-platform cooperation | |
CN109581317A (en) | One kind being based on the matched corner object localization method of echo-peak | |
CN106526549A (en) | False target identification method with combination of two-coordinate radar and three-coordinate radar statistics | |
CN108152789B (en) | Utilize the passive track-corelation data correlation and localization method of RSS information | |
CN107229033A (en) | Multiple target reaching time-difference localization method based on height dimension sectioning search | |
CN106872942A (en) | For the positioning precision calculation method of Distributed Multi positioning monitoring system | |
CN105188082A (en) | Evaluation method for RSS (Received Signal Strength)/AOA (Angle of Arrival)/TDOA (Time Difference of Arrival) positioning performance under indoor WLAN (Wireless Local Area Network) environment | |
CN108614268A (en) | The acoustics tracking of low altitude high speed airbound target | |
CN103376447A (en) | Method for achieving three-dimension positioning of non-cooperative double multistatic radar target | |
CN104330788A (en) | Radio location method based on reach path reverse tracking | |
Lei et al. | Optimal distribution for four-station TDOA location system | |
CN104021292A (en) | Dim target detection and tracking method based on formation active networking | |
CN106680783B (en) | Method based on station site error's blending algorithm confrontation false targets | |
CN104793177A (en) | Microphone array direction finding method based on least square methods | |
CN105787081A (en) | Radiation platform correlation method based on radiation source space positions | |
Giacometti et al. | Accuracy studies for TDOA-AOA localization of emitters with a single sensor | |
Kazemi et al. | Efficient closed-form solution for 3-D hybrid localization in multistatic radars | |
CN104569915A (en) | Positioning method used in multiple photoelectric detection systems and based on target movement model | |
CN108562893A (en) | A kind of external illuminators-based radar multistation combined tracking method | |
CN103728608A (en) | Antenna arrangement method for improving MIMO-OTH radar detecting performance in ionized layer double-Gaussian model | |
Ding et al. | Multistatic ultra-wideband localization for NLOS environments | |
Wang et al. | The research of indoor three-dimensional positioning algorithm based on ultra-wideband technology |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191105 |
|
RJ01 | Rejection of invention patent application after publication |