CN105738867A - Illegal radio station rapid positioning and monitoring method adopting high-precision clock - Google Patents

Illegal radio station rapid positioning and monitoring method adopting high-precision clock Download PDF

Info

Publication number
CN105738867A
CN105738867A CN201610065675.3A CN201610065675A CN105738867A CN 105738867 A CN105738867 A CN 105738867A CN 201610065675 A CN201610065675 A CN 201610065675A CN 105738867 A CN105738867 A CN 105738867A
Authority
CN
China
Prior art keywords
receiving front
coordinate
receiving
radio
monitoring
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
Application number
CN201610065675.3A
Other languages
Chinese (zh)
Inventor
焦杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin Institute Of Radio And Television (scientific And Technological Information Center Jilin Bureau Of Press Publication Radio Film And Television)
Original Assignee
Jilin Institute Of Radio And Television (scientific And Technological Information Center Jilin Bureau Of Press Publication Radio Film And Television)
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin Institute Of Radio And Television (scientific And Technological Information Center Jilin Bureau Of Press Publication Radio Film And Television) filed Critical Jilin Institute Of Radio And Television (scientific And Technological Information Center Jilin Bureau Of Press Publication Radio Film And Television)
Priority to CN201610065675.3A priority Critical patent/CN105738867A/en
Publication of CN105738867A publication Critical patent/CN105738867A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/06Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-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/14Determining absolute distances from a plurality of spaced points of known location

Abstract

The invention relates to an illegal radio station rapid positioning and monitoring method adopting a high-precision clock, and belongs to the radio and television and radio monitoring technology field. By adopting the characteristic of the direct ratio between the distance between the illegal radio station and the monitoring front end and the electromagnetic wave transmission time, the position coordinate of the illegal radio station can be calculated, and a monitoring system built by adopting the above mentioned method is constituted by four monitoring front ends and one monitoring center. The method provided by the invention is advantageous in that compared to the monitoring vehicle positioning method, the production cost is low, the operation cost is low, the maintenance is easy, and the positioning speed is fast by adopting the fixed unmanned monitoring point method; the precision of the positioning method is mainly relied on the clock precision, when the clock precision is high, and the error is in a range from 10-10~10-12, the measurement precision can reach the meter level, and the position of the illegal radio station can be located to the concrete building, and even the concrete floor.

Description

The method using high precision clock position monitor quick to pirate radio
Technical field
The invention belongs to radio and television and radio monitoring technical field, be mainly used in radio supervision and believe with radio and television Number monitoring industry, be a kind of method for determining pirate radio position.
Background technology
In recent years, occur in that a lot of unauthorized illegal broadcasting station throughout the country, based on f-m broadcast station, Rare have amplitude modulation broadcasting and television broadcasting station.Original radio and television are the instruments that Party and government carry out propaganda work, and The appearance of pirate radio, has upset normal broadcast and has televised the order with radio control.It pretends to be some country matchmakers Body, is engaged in false advertisement of drugs publicity mostly, and content is vulgar, has had a strong impact on the public credibility of mass media of Party and government. More seriously, some pirate radios of border area also have the figure of hostile forces behind, they carry out various instead Dynamic publicity, destroys stablizing of country, serious harm society.Expressly provide: any without Wen Guangxin office above county level The radio and television station that approval is set up, broadly falls into " the black appliances platform " illegally set up, once finds, all should give in accordance with the law To ban.
Pirate radio is the most all erected in the skyscraper in people district, illegal divided and has usually pretended equipment hidden Cover.The pirate radio also having can be taked to broadcast with frequency with legal radio station, or is all put on automobile mobile broadcast, adds The difficulty of work is banned in supervision location.
Substantially first set out monitoring car carries out bigness scale to the means taked pirate radio location at present, tentatively finds illegal The residential block at place, radio station.In order to improve the efficiency of location, typically can set out two monitoring cars, respectively from one simultaneously The both direction in city starts simultaneously at directional monitoring, needs, during monitoring, exchange monitoring information of communicating with each other.Utilize The information that both direction crosses gradually approaches the region of pirate radio.As the residential neighborhood being positioned near pirate radio Time, the location tasks of monitoring car is just complete, and investigates with artificial handheld portable monitoring machine walking the most again.
Monitoring car self cost is the most at a relatively high, and flat is that maintenance cost is the highest, and traditional monitoring car location mode relates to Department and personnel many, therefore use cost is the highest.Monitoring car needs road to travel, location work be one gradually Near the process of pirate radio, therefore locating speed is slow, and whole position fixing process needs several hours at the soonest, the most then need Take one day.Monitoring car uses DF and location method, and positioning precision is relatively low, can only navigate near community, it is impossible to Navigate to pirate radio place residential building, more can not navigate to the floor at pirate radio place.
Summary of the invention
The present invention provides a kind of method using high precision clock position monitor quick to pirate radio, to solve monitoring car Location cost is high, and locating speed is slow, the problem that positioning precision is low.
The present invention adopts the technical scheme that, comprises the following steps:
(1) in the city of monitoring management, a set of unmanned monitoring front end that is respectively mounted of four different locations is selected to set Standby as four unmanned monitoring points, respectively receiving front-end 1, receiving front-end 2, receiving front-end 3 and receiving front-end 4, And the longitude and latitude of the good each receiving front-end reception antenna position of surveying record and height above sea level.In order to improve at elevation On direction, the precision of location, needs to allow four front ends the most at grade;In each unmanned monitoring front end, installation accuracy is 10-10~10-12High precision clock, and after close adjustment, synchronization, can be carried out further work;
(2) first allow these four receiving front-ends receive the electromagnetic wave signal that same pirate radio P is sent, first simultaneously Determining that this electromagnetic wave signal has been received normally by each front end, then these four receiving front-ends are at high precision clock Synchronous regime under, start simultaneously at the demodulation acoustical signal to this electromagnetic wave signal and sample;
(3) after sampling terminates, first data are carried out lossless compress, then compression data are sent back to center;
(4) after center harvests the data sent of four front ends, respectively decompress after obtain receiving front-end 1, 2, original sampling data sequence U of 3,41n、U2n、U3n、U4n, n in subscript represents the sequence number of data sequence;
To data U1nAnd U2nIt is analyzed calculating its dependency:
γ 12 = Σ n = 1 N U 1 n U 2 n
From U1nSecond data start, again with U2nCalculate dependency;The most again from U1nThe 3rd data open Begin to calculate.So obtain the sequence of the dependency of the sampled data of summary front end 1,2, found in this sequence Maximum, counter can release the U that maximum is corresponding1nAnd U2nThe subscript sequence number of two groups of data differences, before receiving Corresponding clock times during end 1,2 sampling, can obtaining front end 1 and front end 2 in conjunction with sample frequency, to receive this non- The time difference Δ of method radio signals12
In the same way, calculating receiving front-end 1 and receiving front-end 3 receive the time difference Δ of pirate radio signal13, Receiving front-end 1 and receiving front-end 4 is also had to receive the time difference Δ of pirate radio signal14
(5) pirate radio P coordinates are (X, Y, Z), and the radio signal launched at moment t will be at ti It is (X that moment arrives coordinatei,Yi,Zi) No. i-th receiving front-end, i=1,2,3,4, pirate radio and No. i-th The distance of receiving front-end is:
( X - X i ) 2 + ( Y - Y i ) 2 + ( Z - Z i ) 2 = C ( t i - t )
Wherein C is the light velocity;
The range equation of any two front end is subtracted each other, and disappear moment t, obtain receiving front-end 1 and receiving front-end 2 away from Deviation equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 2 ) 2 + ( Y - Y 2 ) 2 + ( Z - Z 2 ) 2 = CΔ 12
Front end 1 and front end 3 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 3 ) 2 + ( Y - Y 3 ) 2 + ( Z - Z 3 ) 2 = CΔ 13
Front end 1 and front end 4 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 4 ) 2 + ( Y - Y 4 ) 2 + ( Z - Z 4 ) 2 = CΔ 14
Wherein, (X1,Y1,Z1)、(X2,Y2,Z2)、(X3,Y3,Z3) and (X4,Y4,Z4) be respectively code name be 1,2, 3,4 these four receiving front-ends are at three-dimensional coordinate, and (X, Y, Z) is the coordinate position of pirate radio;
(6) calculate the rough position of pirate radio in the plane with approximatioss, appoint with kilometer integral multiple and take a point coordinates (X, Y), substitutes into receiving front-end 1 and receiving front-end 2 range difference equation, is calculated a time difference Δ, the most again Positive direction along X changes a coordinate again, again calculates a time difference Δ of getting back ', if Δ ' than Δ closer to Δ12, So illustrate that the direction approached is correct, be exactly otherwise to have walked anti-direction;Until finally giving closest to Δ12That Individual coordinate points, then changes a Y at this closest to coordinate points, again obtains immediate point with above-mentioned approximatioss, The most repeatedly, institute in the receiving front-end 1 in units of 1 kilometer and receiving front-end 2 range difference equation is obtained the thickest Solution sequence S;Coordinate datas all in S are substituted into successively the range difference equation of receiving front-end 1 and receiving front-end 3, To range difference Δ " and Δ13Compare, find that coordinate immediate, it is simply that the rough position of pirate radio in the plane Put;
(7) near pirate radio rough position, appoint in units of the integer-bit of 50 meters and take a point coordinates (X, Y), Use the mode of last step, again substitute into receiving front-end 1 and receiving front-end 2 range difference equation, and at this 1 kilometer Within change Y, obtain the receiving front-end 1 in units of 50 meters and receiving front-end 2 range difference equation likely solve Sequence E, afterwards with each coordinate figure in the range difference equation cycle tests E of receiving front-end 1 and receiving front-end 3, Find and Δ13That coordinate that error is minimum, it is simply that the pirate radio plane coordinates in units of 50 meters;
(8) it is calculated near pirate radio plan-position in 300 meters in previous step, adds the coordinate of Z-direction, In units of 25 meters, divide solid space become the mesh space of 12 × 12 × 12, all seats in this space Punctuate substitutes into receiving front-end 1 and receiving front-end 2, receiving front-end 1 and receiving front-end 3, receiving front-end 1 one by one and connects Receive the range difference equation of front end 4, obtain all time differences and measurement result that coordinate the most immediate, it is simply that be illegal The coordinate of place, radio station three-dimensional space position.
In these four receiving front-ends of the present invention, need the height above sea level of one of them compared with other receiving front-end, Difference in height is 50~300 meters.
One embodiment of the present invention is: in step (five), when a City Terrain rises and falls the least, Z1=Z2=Z3=Z4
It is an advantage of the invention that the spy utilizing pirate radio to be directly proportional to the electromagnetic transmission time to the distance of monitoring front end Property, the method calculating pirate radio position coordinates, because the method that have employed fixing unmanned monitoring point, so with Monitoring car localization method is compared, and has cost much lower, and operating cost is low, safeguards easily, the advantage that locating speed is fast; The precision of this localization method depends primarily on the precision of clock, and when clock accuracy is higher, error is 10-10~10-12Time, The precision measured, up to meter level, is possible not only to position pirate radio position reaches which building, building concrete, it might even be possible to navigate to Concrete floor.
Accompanying drawing explanation
Fig. 1 is system primary structure of the present invention, and system is utilized network by central server and unmanned headend equipment Coupling composition, radio supervision person on duty manipulates this set monitoring system on central server, positions in three dimensions Needing the unmanned monitoring headend equipment of four diverse locations, on two dimensional surface, location is that three dimensions depth displacement is equal to The special case of zero, at this moment has only to the unmanned monitoring headend equipment of three diverse locations;
Fig. 2 describes two receiving front-ends 1,2 and receives what the situation of same signal, curve 1 and curve 2 received Same signal is all initial at arrow from figure, but has difference on a timeline, and curve 1 is moved on a timeline, When moving to dotted line position, the numerical value γ of two groups of data calculating dependencys reaches maximum, now curve 1 data subscript The numerical value that sequence number moves is multiplied by the time difference Δ that the cycle of data sampling is exactly two diverse location front ends of signal arrival;
Fig. 3 describes the pirate radio P that coordinate is (X, Y, Z), at the aerogram that moment t launches Number will be at tiIt is (X that moment arrives coordinatei,Yi,Zi) i-th receiving front-end, the distance etc. of pirate radio and front end It is multiplied by the time in the light velocity.
Detailed description of the invention
Comprise the following steps:
(1) in the city of monitoring management, a set of unmanned monitoring front end that is respectively mounted of four different locations is selected to set Standby as four unmanned monitoring points, respectively receiving front-end 1, receiving front-end 2, receiving front-end 3 and receiving front-end 4, And the longitude and latitude of the good each receiving front-end reception antenna position of surveying record and height above sea level.In order to improve at elevation On direction, the precision of location, needs to allow four front ends the most at grade;In these four receiving front-ends, need it In the height above sea level of compared with other receiving front-end, difference in height is 50~300 meters;In supervision center machine room, peace Filling a set of monitoring center system, utilize network to interconnect with each unmanned monitoring point, accompanying drawing 1 describes the basic of this system Structure, in each unmanned monitoring front end, installation accuracy is 10-10~10-12High precision clock, and through close adjustment, with Further work is can be carried out after step;
(2) first allow these four receiving front-ends receive the electromagnetic wave signal that same pirate radio P is sent, first simultaneously Determine that this electromagnetic wave signal has been received normally by each front end, before then recycling network advertisement these four receives End, under the synchronous regime of high precision clock, starts simultaneously at the demodulation acoustical signal to this electromagnetic wave signal and adopts Sample;Different according to signal receiving quality, sampling time length was from several seconds to tens seconds.During signal to noise ratio height, several seconds just Enough;During poor signal to noise, the sampling time longer precision being favorably improved location measurement;
(3) after sampling terminates, first data are carried out lossless compress, then compression data are sent back to center, utilize The program solution harmonic signal frequency that data sampling frequency is sent far above pirate radio P, the number of adjacent two sampled points According to the feature that difference is less;First carry out Run-Length Coding with to data, be compressed with huffman coding the most again, afterwards Recycling network transmission;The data that so can in minutes realize that each front-end sampling is obtained are sent to prison On the server of tube hub machine room;
(4) after center harvests the data sent of four front ends, respectively decompress after obtain receiving front-end 1, 2, original sampling data sequence U of 3,41n、U2n、U3n、U4n, n in subscript represents the sequence number of data sequence;
To data U1nAnd U2nIt is analyzed calculating its dependency:
γ 12 = Σ n = 1 N U 1 n U 2 n
From U1nSecond data start, again with U2nCalculate dependency;The most again from U1nThe 3rd data open Begin to calculate.So obtain the sequence of the dependency of the sampled data of summary front end 1,2, found in this sequence Maximum, counter can release the U that maximum is corresponding1nAnd U2nThe subscript sequence number of two groups of data differences, before receiving Corresponding clock times during end 1,2 sampling, can obtaining front end 1 and front end 2 in conjunction with sample frequency, to receive this non- The time difference Δ of method radio signals12, accompanying drawing 2 describes this process;
In the same way, calculating receiving front-end 1 and receiving front-end 3 receive the time difference Δ of pirate radio signal13, Receiving front-end 1 and receiving front-end 4 is also had to receive the time difference Δ of pirate radio signal14
(5) accompanying drawing 3 describes a pirate radio P coordinate is (X, Y, Z), in the nothing that moment t launches The line signal of telecommunication will be at tiIt is (X that moment arrives coordinatei,Yi,Zi) No. i-th receiving front-end, i=1,2,3,4, non- Method radio station with the distance of No. i-th receiving front-end is:
( X - X i ) 2 + ( Y - Y i ) 2 + ( Z - Z i ) 2 = C ( t i - t )
Wherein C is the light velocity;
The range equation of any two front end is subtracted each other, and disappear moment t, obtain receiving front-end 1 and receiving front-end 2 away from Deviation equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 2 ) 2 + ( Y - Y 2 ) 2 + ( Z - Z 2 ) 2 = CΔ 12
Front end 1 and front end 3 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 3 ) 2 + ( Y - Y 3 ) 2 + ( Z - Z 3 ) 2 = CΔ 13
Front end 1 and front end 4 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 4 ) 2 + ( Y - Y 4 ) 2 + ( Z - Z 4 ) 2 = CΔ 14
Wherein, (X1,Y1,Z1)、(X2,Y2,Z2)、(X3,Y3,Z3) and (X4,Y4,Z4) be respectively code name be 1,2, 3,4 these four receiving front-ends are at three-dimensional coordinate, and (X, Y, Z) is the coordinate position of pirate radio;
If a City Terrain rises and falls the least, then may be considered and positions on a two dimensional surface, namely exists Z1=Z2=Z3=Z4Situation.
(6) rough position of pirate radio is calculated in the plane with approximatioss, in a city, with kilometer integer Appoint again and take a point coordinates (X, Y), substitute into receiving front-end 1 and receiving front-end 2 range difference equation, when being calculated one Between difference Δ, the positive direction along X changes a coordinate more again, again calculates a time difference Δ of getting back ', if Δ ' Than Δ closer to Δ12, then illustrate that the direction approached is correct, be exactly otherwise to have walked anti-direction;Until final To closest to Δ12That coordinate points, then change a Y at this closest to coordinate points, again obtain with above-mentioned approximatioss To immediate point, the most repeatedly, the receiving front-end 1 in units of 1 kilometer and receiving front-end 2 range difference side are obtained Institute likely thick solution sequence S in journey;Coordinate datas all in S are substituted into receiving front-end 1 and receiving front-end 3 successively Range difference equation, obtains range difference Δ " and Δ13Compare, find that coordinate immediate, it is simply that in the plane Pirate radio rough position;
(7) near pirate radio rough position, appoint in units of the integer-bit of 50 meters and take a point coordinates (X, Y), Use the mode of last step, again substitute into receiving front-end 1 and receiving front-end 2 range difference equation, and at this 1 kilometer Within change Y, obtain the receiving front-end 1 in units of 50 meters and receiving front-end 2 range difference equation likely solve Sequence E, afterwards with each coordinate figure in the range difference equation cycle tests E of receiving front-end 1 and receiving front-end 3, Find and Δ13That coordinate that error is minimum, it is simply that the pirate radio plane coordinates in units of 50 meters;
(8) it is calculated near pirate radio plan-position in 300 meters in previous step, adds the coordinate of Z-direction, In units of 25 meters, divide solid space become the mesh space of 12 × 12 × 12, all seats in this space Punctuate substitutes into receiving front-end 1 and receiving front-end 2, receiving front-end 1 and receiving front-end 3, receiving front-end 1 one by one and connects Receive the range difference equation of front end 4, obtain all time differences and measurement result that coordinate the most immediate, it is simply that be illegal The coordinate of place, radio station three-dimensional space position.
In general, the pirate radio coordinate precision so obtained has met actually used needs, if in system Clock accuracy sufficiently high (10-12), can be near the calculated coordinate position of upper step, with the side as previous step Method, is that unit divides space again with rice, can be calculated more accurate pirate radio position;It addition, in plane On the location that carries out be the special case of three-dimensional space location, when positioning in the plane, it is only necessary to have two range equations The most much of that, say, that plane positioning only need to have three monitoring points just can realize positioning function.

Claims (3)

1. the method using high precision clock position monitor quick to pirate radio, it is characterised in that include following Step:
(1) in the city of monitoring management, a set of unmanned monitoring front end that is respectively mounted of four different locations is selected to set Standby as four unmanned monitoring points, respectively receiving front-end 1, receiving front-end 2, receiving front-end 3 and receiving front-end 4, And the longitude and latitude of the good each receiving front-end reception antenna position of surveying record and height above sea level.In order to improve at elevation On direction, the precision of location, needs to allow four front ends the most at grade;In each unmanned monitoring front end, installation accuracy is 10-10~10-12High precision clock, and after close adjustment, synchronization, can be carried out further work;
(2) first allow these four receiving front-ends receive the electromagnetic wave signal that same pirate radio P is sent, first simultaneously Determining that this electromagnetic wave signal has been received normally by each front end, then these four receiving front-ends are at high precision clock Synchronous regime under, start simultaneously at the demodulation acoustical signal to this electromagnetic wave signal and sample;
(3) after sampling terminates, first data are carried out lossless compress, then compression data are sent back to center;
(4) after center harvests the data sent of four front ends, respectively decompress after obtain receiving front-end 1, 2, original sampling data sequence U of 3,41n、U2n、U3n、U4n, n in subscript represents the sequence number of data sequence;
To data U1nAnd U2nIt is analyzed calculating its dependency:
γ 12 = Σ n = 1 N U 1 n U 2 n
From U1nSecond data start, again with U2nCalculate dependency;The most again from U1nThe 3rd data open Begin to calculate.So obtain the sequence of the dependency of the sampled data of summary front end 1,2, found in this sequence Maximum, counter can release the U that maximum is corresponding1nAnd U2nThe subscript sequence number of two groups of data differences, before receiving Corresponding clock times during end 1,2 sampling, can obtaining front end 1 and front end 2 in conjunction with sample frequency, to receive this non- Time difference △ of method radio signals12
In the same way, calculating receiving front-end 1 and receiving front-end 3 receive time difference △ of pirate radio signal13, Receiving front-end 1 and receiving front-end 4 is also had to receive time difference △ of pirate radio signal14
(5) pirate radio P coordinates are (X, Y, Z), and the radio signal launched at moment t will be at ti It is (X that moment arrives coordinatei,Yi,Zi) No. i-th receiving front-end, i=1,2,3,4, pirate radio and No. i-th The distance of receiving front-end is:
( X - X i ) 2 + ( Y - Y i ) 2 + ( Z - Z i ) 2 = C ( t i - t )
Wherein C is the light velocity;
The range equation of any two front end is subtracted each other, and disappear moment t, obtain receiving front-end 1 and receiving front-end 2 away from Deviation equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 2 ) 2 + ( Y - Y 2 ) 2 + ( Z - Z 2 ) 2 = CΔ 12
Front end 1 and front end 3 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 3 ) 2 + ( Y - Y 3 ) 2 + ( Z - Z 3 ) 2 = CΔ 13
Front end 1 and front end 4 range difference equation:
( X - X 1 ) 2 + ( Y - Y 1 ) 2 + ( Z - Z 1 ) 2 - ( X - X 4 ) 2 + ( Y - Y 4 ) 2 + ( Z - Z 4 ) 2 = CΔ 14
Wherein, (X1,Y1,Z1)、(X2,Y2,Z2)、(X3,Y3,Z3) and (X4,Y4,Z4) be respectively code name be 1,2, 3,4 these four receiving front-ends are at three-dimensional coordinate, and (X, Y, Z) is the coordinate position of pirate radio;
(6) calculate the rough position of pirate radio in the plane with approximatioss, appoint with kilometer integral multiple and take a point coordinates (X, Y), substitutes into receiving front-end 1 and receiving front-end 2 range difference equation, is calculated time difference △, the most again Positive direction along X changes a coordinate again, again calculate time difference △ of getting back ', if △ ' than △ closer to △12, So illustrate that the direction approached is correct, be exactly otherwise to have walked anti-direction;Until finally giving closest to △12That Individual coordinate points, then changes a Y at this closest to coordinate points, again obtains immediate point with above-mentioned approximatioss, The most repeatedly, institute in the receiving front-end 1 in units of 1 kilometer and receiving front-end 2 range difference equation is obtained the thickest Solution sequence S;Coordinate datas all in S are substituted into successively the range difference equation of receiving front-end 1 and receiving front-end 3, To range difference △ " and △13Compare, find that coordinate immediate, it is simply that the rough position of pirate radio in the plane Put;
(7) near pirate radio rough position, appoint in units of the integer-bit of 50 meters and take a point coordinates (X, Y), Use the mode of last step, again substitute into receiving front-end 1 and receiving front-end 2 range difference equation, and at this 1 kilometer Within change Y, obtain the receiving front-end 1 in units of 50 meters and receiving front-end 2 range difference equation likely solve Sequence E, afterwards with each coordinate figure in the range difference equation cycle tests E of receiving front-end 1 and receiving front-end 3, Find and △13That coordinate that error is minimum, it is simply that the pirate radio plane coordinates in units of 50 meters;
(8) it is calculated near pirate radio plan-position in 300 meters in previous step, adds the coordinate of Z-direction, In units of 25 meters, divide solid space become the mesh space of 12 × 12 × 12, all seats in this space Punctuate substitutes into receiving front-end 1 and receiving front-end 2, receiving front-end 1 and receiving front-end 3, receiving front-end 1 one by one and connects Receive the range difference equation of front end 4, obtain all time differences and measurement result that coordinate the most immediate, it is simply that be illegal The coordinate of place, radio station three-dimensional space position.
The method of use high precision clock the most according to claim 1 position monitor quick to pirate radio, its Be characterised by: in these four receiving front-ends described, need the height above sea level of one of them compared with other receiving front-end, Difference in height is 50~300 meters.
The method of use high precision clock the most according to claim 1 position monitor quick to pirate radio, its It is characterised by step (five), when a City Terrain rises and falls the least, Z1=Z2=Z3=Z4
CN201610065675.3A 2016-01-31 2016-01-31 Illegal radio station rapid positioning and monitoring method adopting high-precision clock Pending CN105738867A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610065675.3A CN105738867A (en) 2016-01-31 2016-01-31 Illegal radio station rapid positioning and monitoring method adopting high-precision clock

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610065675.3A CN105738867A (en) 2016-01-31 2016-01-31 Illegal radio station rapid positioning and monitoring method adopting high-precision clock

Publications (1)

Publication Number Publication Date
CN105738867A true CN105738867A (en) 2016-07-06

Family

ID=56247075

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610065675.3A Pending CN105738867A (en) 2016-01-31 2016-01-31 Illegal radio station rapid positioning and monitoring method adopting high-precision clock

Country Status (1)

Country Link
CN (1) CN105738867A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112883895A (en) * 2021-03-08 2021-06-01 山东大学 Illegal electromagnetic signal detection method based on self-adaptive weighted PCA and realization system thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201656990U (en) * 2010-04-08 2010-11-24 成都西锐科技有限公司 Intelligent radio interference search system
CN201773183U (en) * 2010-05-07 2011-03-23 北京邮电大学 Portable radio direction finding device
CN102694607A (en) * 2011-03-25 2012-09-26 北京海曼无限信息技术有限公司 Networked software radio signal monitoring method aiming at mobile interference source and system thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201656990U (en) * 2010-04-08 2010-11-24 成都西锐科技有限公司 Intelligent radio interference search system
CN201773183U (en) * 2010-05-07 2011-03-23 北京邮电大学 Portable radio direction finding device
CN102694607A (en) * 2011-03-25 2012-09-26 北京海曼无限信息技术有限公司 Networked software radio signal monitoring method aiming at mobile interference source and system thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高远谋: "基于射频传感器网络的无线信号监测与定位技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112883895A (en) * 2021-03-08 2021-06-01 山东大学 Illegal electromagnetic signal detection method based on self-adaptive weighted PCA and realization system thereof
CN112883895B (en) * 2021-03-08 2022-03-25 山东大学 Illegal electromagnetic signal detection method based on self-adaptive weighted PCA and realization system thereof

Similar Documents

Publication Publication Date Title
JP6426222B2 (en) Ephemeris Extension System and how to use in GNSS
US7212160B2 (en) Locating an object of interest using back tracked-rays derived from multipath signals
Guzman-Acevedo et al. GPS, accelerometer, and smartphone fused smart sensor for SHM on real-scale bridges
US20090121940A1 (en) System for determining position over a network
CN101536596B (en) Method and device for determination of the position of a terminal in a mobile communication network
Ogle et al. Accuracy of global positioning system for determining driver performance parameters
CN1383492A (en) Time of arrival estimation for positioning systems
CN101836080A (en) A method of processing positioning data
CN101742609A (en) Positioning method based on AOA (Activity on Arrows)and terminal moving track under single-cell environment
CN102253379B (en) System for monitoring landslide based on radio interference technology
Sharp et al. Positional accuracy measurement and error modeling for mobile tracking
CN105334495A (en) Non-line-of-sight stable positioning method based on signal arrival time in wireless network
CN101860384A (en) Ionosphere hybrid modeling method in short-wave ray tracing technique
CN102749637A (en) Method for realizing accurate positioning of vehicle-mounted GPS (Globe Positioning System)
CN102226700B (en) Method for matching electronic map of flyover road network
CN104217588A (en) Method, server and system for acquiring real-time traffic information
CN103810855A (en) Multi-characteristic quantity comprehensive traffic detection system oriented to multi-user group
CN1310036C (en) Apparatus and method for providing a telematics service having and agps function
WO2006088390A1 (en) Method for on-line monitoring and controlling moving objects
CN105738867A (en) Illegal radio station rapid positioning and monitoring method adopting high-precision clock
CN1417557A (en) Intelligent vehicle navigator
CN100516927C (en) Calibration and correction system for satellite position location systems
CN115698625A (en) Method for estimating bridge displacement and electronic device for estimating bridge displacement
CN105182393A (en) Realization method based on Beidou satellite positioning technology and internet-of-vehicles technology
Yokota Localization algorithm based on altitude time series in GNSS-denied 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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160706