Summary of the invention
In order to overcome above deficiency; Particularly mine positioning system poor anti jamming capability, setting accuracy are low in order to overcome, there are problems such as NLOS propagation in mine positioning equipment limited volume, ultra-broadband signal under mine; The utility model provides ultra broadband positioning system under a kind of accurate mine; The characteristics of this systems incorporate lane space utilize limited reference mode and simple algorithm to realize in real time accurate target localization, have satisfied the needs of in time rescuing after mine production scheduling and the calamity.
To achieve these goals, the utility model adopts following technical scheme:
Ultra broadband positioning system under a kind of coal mine comprises telemonitoring layer and on-site supervision layer, and said telemonitoring layer comprises browser server, data server, ground monitoring terminal, switch; Said on-site supervision layer comprises destination node, reference mode, gateway and base station; Data server and ground monitoring terminal constitute spider lines through switch and base station, and send the real-time positioning data to browser server through the Internet network; Destination node, reference mode, gateway and base station constitute the super wideband wireless fixer network; The base station is articulated on the switch through optical fiber or bus; Destination node is carried by underground operators or mechanical hook-up;
Said data server receives and storage destination node location data;
Said ground monitoring terminal; Record reference mode deployed position information; According to the distance measure that receives, carry out NLOS and differentiate, and distribute according to the deployed position and the tunnel of NLOS identification result and reference mode; Calculate the exact position of destination node, and positional information is redispatched to data server and browser server;
Locating information is responsible for transmitting in said switch and base station;
Said gateway is responsible for accomplishing the forwarding of super wideband wireless locating information and the transmission of configuration information;
Said reference mode is placed in the gateway device wireless signal coverage, provides a packets of information that comprises own position coordinates and reception signal strength indication value (RSSI) to give destination node, accomplishes the calculating of time of arrival of transmission ultra-broadband signal;
Said destination node, corresponding with personnel in the pit's identity, accomplish choosing of reference mode;
The underground equipment that relates in the described system all is an intrinsically safe equipment.
Said gateway comprises: microprocessor, storer, bus module, ultra-wideband antenna, power supply, communication interface, LCD, selection key.
Said reference mode comprises: processor, power supply, storer, ultra-broadband signal generation unit, energy measuring receiver structure, ultra-wideband antenna, on-off circuit.
Said destination node comprises: processor, power supply, storer, ultra-broadband signal generation unit, energy measuring receiver structure, ultra-wideband antenna, on-off circuit.
Said energy measuring receiver structure comprises: ultra-wideband antenna, low noise amplifier LNA, BPF. BPF, squarer, integrator, sampling switch, decision device.
To achieve these goals; The utility model provides a kind of and differentiates the UWB location technology with weighted least-squares method based on nonparametric NLOS: according to colliery underworkings actual features, certain apart from the linear deployment reference mode, during target localization on the edge, both sides of the edge in tunnel; At first confirm to participate in the pinpoint reference mode of destination node according to receiving signal strength values; And estimate at every turn the time of arrival (TOA) from the ultra broadband positioning signal of destination node, calculating M corresponding distance value, utilize nonparametric NLOS discrimination method that the channel between destination node and the datum node is carried out NLOS again and differentiate; According to the NLOS identification result; The weighting factor that distributes the least square method location algorithm is constructed the objective function of weighting factor least square method location algorithm, and calculates the position of destination node.
Said UWB location technology based on nonparametric NLOS discriminating and weighted least-squares method may further comprise the steps:
A. confirm to participate in the pinpoint reference mode of destination node, and calculate at every turn TOA time of arrival, repeat M time from the ultra broadband positioning signal of destination node;
B. M the TOA value that estimates according to each datum node calculates M corresponding distance value, utilizes nonparametric NLOS discrimination method that the channel between destination node and the datum node is carried out NLOS again and differentiates;
C. according to the NLOS identification result, distribute the weighting factor of least square method location algorithm;
D. according to the weighting factor that distributes, construct the objective function of weighting factor least square method location algorithm, and the position of estimating target node.
Said steps A comprises the following steps:
A1. in the both sides of the edge in tunnel along certain apart from the linear deployment reference mode;
A2. after destination node MS got into the tunnel, the interior towards periphery reference mode of MS was periodically broadcasted positioning request signal, sends M time continuously; After reference mode is received positioning request signal; Therefrom propose one group of RSSI value, then this group RSSI value is asked on average, and write down this mean value;
A3. destination node MS broadcasts the RSSI information of once collecting after having sent M Location Request; At this moment reference mode can be replied the RSSI mean value of its record to this destination node around; After destination node has been collected RSSI information; If the number of the RSSI value of collecting less than 3, is then broadcasted positioning request signal more again; Otherwise, the RSSI value of collecting is sorted, and choose maximum 3 nodes as a reference, and to these 3 reference modes positioning signal of redispatching M time;
A4. after reference mode receives positioning signal at every turn, with the signal r (t) that receives through carrying out integration sampling, the energy sampled sequence of picked up signal behind the squarer again:
N=T
f/ T
b, wherein n representes the sequence number of energy block, T
bThe expression integration period, T
fIn the expression frame period, N representes a frame self-energy piece number;
A5. the threshold value θ according to described energy sampled sequence and setting compares, and first is exceeded the sampling block of the energy block of thresholding as direct-path DP place, promptly
Y wherein
nBe the energy sampled sequence of signal,
The energy block at expression DP place, θ representes judgement threshold,
Expression DP arrives the time of target.
Said step B comprises the following steps:
B1. datum node has been after having estimated M TOA, can calculate M corresponding distance value, for i (i=1 ..., M) inferior measurement, r under sighting distance (LOS) situation
i=s+n
Noise, r under the NLOS situation
i=s+n
Noise+ n
NLOS, r wherein
iRepresent the i time distance measure, s representes the actual distance between destination node and the reference mode, n
NoiseNoise, n are measured in expression
NLOSExpression NLOS error;
B2. calculate and measure noise n
NoiseProbability density function (PDF), under the LOS situation, n
NoiseCan be expressed as the Gaussian random variable that 0 average value standard deviation is σ, its probability density function is expressed as:
B3. according to the probability density function P of said measurement noise
Noise(x), calculate distance measure r under the LOS situation
i(i=1 ..., PDF M):
Wherein s representes the actual distance between destination node and the reference mode;
B4. according to described M distance measure r
i(i=1 ..., M), the calculating measuring distance
Wherein Φ (x) representes the window function of Gaussian window or rectangular window, and L representes zooming parameter, r
iRepresent i distance measure;
B5. according to described P
Noise(x-s) and
Calculate P
Noise(x-s) and
Between minor increment:
Wherein
Represent M distance measure r
i(i=1 ..., mean value M);
B6. according to described minor increment
The size that compares this minimum value and threshold delta, when
Judgement is LOS, when
Then judgement is NLOS, and wherein threshold delta is represented as follows: δ '=δ-L
2σ
ω 2/ (2 σ
2) (σ
ω 2=∫ x
2Φ (x) dx), L representes the zooming parameter of selected window function, and σ representes to measure under the LOS situation noise n
NoiseStandard deviation.
Among the said step C, be its weights assigned factor, when j reference mode is judged as los state, for the weighting factor of its distribution is k according to the NLOS identification result of j reference mode
j=1; When this reference mode is judged as the NLOS state, then do for the weighting factor of its distribution
Wherein i (i=1,2 ..., M) expression is measured for the i time, and j (j=1,2,3) representes j reference mode, r
JiThe i time distance measure representing j reference mode,
M the distance measure r that representes j reference mode
JiMean value.
Said step D comprises the following steps:
D1. according to described weighting factor, construct the objective function of weighting factor least square method location algorithm for the pinpoint reference mode distribution of participation destination node:
Wherein (x y) representes position coordinates variable, X to X
j(x
j, y
j) j reference mode position coordinates of (j=1,2,3) expression, k
jJ weighting factor that reference mode is assigned to of (j=1,2,3) expression,
M the distance measure r that representes j reference mode
JiMean value;
D2. according to described objective function F (x), use the position coordinates that least square method is calculated destination node
The beneficial effect of the utility model is:
1. circumstance complication under the mine; The conventional positioning system can't reach pinpoint requirement, through adopting the ultra-broadband signal of strong time resolution, because signal dutyfactor is extremely low; Multipath signal is separable in time, has improved the antijamming capability of localization method and system greatly.The use of super-broadband tech significantly reduces power consumption of terminal equipment simultaneously.And because need not adopt complicated modulation-demodulation technique, make system equipment simple in structure, volume is little, is fit to the use of down-hole particular surroundings.The UWB technology has the incomparable advantage of other wireless communication technologys, and its combining with wireless sensor network technology can be applied in the mine easily, reaches high orientation precision.
2. topmost problem is that the NLOS transmission is serious under the mine; For the bigger NLOS positioning error that reduces to produce; Through adopting nonparametric NLOS discrimination method; The main thought of this method is the degree of closeness of the PDF of PDF and distance measure under the comparison LOS situation, can be used for NLOS easily and effectively and differentiate.
3. in order to reduce the influence that the NLOS error worsens location estimation; Through adopting a kind of new weighting factor least square method location algorithm; The weighting factor of this algorithm is confirmed by the result that NLOS differentiates; Size has reflected NLOS and LOS datum node shared proportion in location estimation, can estimate the position coordinates of destination node more accurately.
Embodiment
For the technical scheme, characteristics and the advantage that make the utility model is clearer, the utility model is further specified below in conjunction with accompanying drawing and embodiment.
Fig. 1 be the utility model under the NLOS situation based on the positioning principle figure of TOA.
As shown in Figure 1, in the UWB positioning system based on TOA, each, τ determined that a center of circle is r for the datum node radius time of arrival
iThe circle of=c * τ (i=1,2,3), wherein c is the speed that ultra-broadband signal is propagated in the tunnel.If there is not the NLOS sum of errors to measure noise error, 3 circles by the TOA decision can intersect at a point fully in two-dimensional space, and this point promptly is the position of destination node; Measure noise error if there is the NLOS sum of errors; Its measuring distance is that actual distance s, NLOS sum of errors are measured the noise error sum; Under the reason condition, the NLOS error can make measuring distance produce a bigger overgauge with respect to actual distance, and this can cause the location ambiguity of destination node; As shown in Figure 1, the position of destination node MS is positioned at the overlapping region of 3 circles.In order to reduce the influence that the NLOS error worsens location estimation; At first the channel between destination node and the reference mode being carried out NLOS differentiates; Then through adopting a kind of new weighting factor least square method location algorithm; The weighting factor of this algorithm is confirmed the exact position of calculating destination node by the result that NLOS differentiates.
Fig. 2 is the positioning system composition frame chart of the utility model.
As shown in Figure 2, the utility model, ultra broadband positioning system under a kind of coal mine comprises telemonitoring layer and on-site supervision layer, the telemonitoring layer comprises data server 1, ground monitoring terminal 2, switch 3, browser server 4; The on-site supervision layer comprises base station 5, gateway 6, reference mode 7, destination node 8; Data server 1 constitutes spider lines with ground monitoring terminal 2 through switch 3 and base station 5, and sends the real-time positioning data to browser server 4 through the Internet network; Destination node 8, reference mode 7, gateway 6 constitutes the super wideband wireless fixer network with base station 5; Base station 5 is articulated on the switch 3 through optical fiber or bus; Destination node 8 is carried by underground operators or mechanical hook-up.
Behind the destination node 8 entering tunnels, interior towards periphery reference mode is periodically broadcasted positioning request signal, according to the RSSI mean value of its record of reference mode answer, confirms to participate in the pinpoint reference mode of destination node.Reference mode 7 is accomplished the calculating of time of arrival of transmission ultra-broadband signal, and the guestimate of distance between completion and the destination node, and sends distance measure to ground monitoring terminal 2 through gateway 6, base station 5 and switch 3.Ground monitoring terminal 2; According to the distance measure that receives, carry out NLOS and differentiate, and distribute according to the deployed position and the tunnel of NLOS identification result and reference mode 7; Calculate the exact position of destination node, and positional information is redispatched to data server 1 and browser server 4.
Fig. 3 is the UWB location technology FB(flow block) based on nonparametric NLOS discriminating and weighted least-squares method of the utility model.
As shown in Figure 3, the UWB location technology based on nonparametric NLOS discriminating and weighted least-squares method may further comprise the steps:
A. confirm to participate in the pinpoint reference mode of destination node, and calculate at every turn TOA time of arrival, repeat M time from the ultra broadband positioning signal of destination node;
B. M the TOA value that estimates according to each datum node calculates M corresponding distance value, utilizes nonparametric NLOS discrimination method that the channel between destination node and the datum node is carried out NLOS again and differentiates;
C. according to the NLOS identification result, distribute the weighting factor of least square method location algorithm;
D. according to the weighting factor that distributes, construct the objective function of weighting factor least square method location algorithm, and the position of estimating target node.
Steps A comprises the following steps:
A1. in the both sides of the edge in tunnel along certain apart from the linear deployment reference mode;
A2. after destination node MS got into the tunnel, the interior towards periphery reference mode of MS was periodically broadcasted positioning request signal, sends M time continuously; After reference mode is received positioning request signal; Therefrom propose one group of RSSI value, then this group RSSI value is asked on average, and write down this mean value;
A3. destination node MS broadcasts the RSSI information of once collecting after having sent M Location Request; At this moment reference mode can be replied the RSSI mean value of its record to this destination node around; After destination node has been collected RSSI information; If the number of the RSSI value of collecting less than 3, is then broadcasted positioning request signal more again; Otherwise, the RSSI value of collecting is sorted, and choose maximum 3 nodes as a reference, and to these 3 reference modes positioning signal of redispatching M time;
A4. after reference mode receives positioning signal at every turn, with the signal r (t) that receives through carrying out integration sampling, the energy sampled sequence of picked up signal behind the squarer again:
N=T
f/ T
b, wherein n representes the sequence number of energy block, T
bThe expression integration period, T
fIn the expression frame period, N representes a frame self-energy piece number;
A5. the threshold value θ according to described energy sampled sequence and setting compares, and first is exceeded the sampling block of the energy block of thresholding as direct-path DP place, promptly
Y wherein
nBe the energy sampled sequence of signal,
The energy block at expression DP place, θ representes judgement threshold,
Expression DP arrives the time of target.
Step B comprises the following steps:
B1. datum node has been after having estimated M TOA, can calculate M corresponding distance value, for i (i=1 ..., M) inferior measurement, r under sighting distance (LOS) situation
i=s+n
Noise, r under the NLOS situation
i=s+n
Noise+ n
NLOS, r wherein
iRepresent the i time distance measure, s representes the actual distance between destination node and the reference mode, n
NoiseNoise, n are measured in expression
NLOSExpression NLOS error;
B2. calculate and measure noise n
NoiseProbability density function (PDF), under the LOS situation, n
NoiseCan be expressed as the Gaussian random variable that 0 average value standard deviation is σ, its probability density function is expressed as:
B3. according to the probability density function P of said measurement noise
Noise(x), calculate distance measure r under the LOS situation
i(i=1 ..., PDF M):
Wherein s representes the actual distance between destination node and the reference mode;
B4. according to described M distance measure r
i(i=1 ..., M), the calculating measuring distance
Wherein Φ (x) representes the window function of Gaussian window or rectangular window, and L representes zooming parameter, r
iRepresent i distance measure;
B5. according to described P
Noise(x-s) and
Calculate P
Noise(x-s) and
Between minor increment:
Wherein
Represent M distance measure r
i(i=1 ..., mean value M);
B6. according to described minor increment
The size that compares this minimum value and threshold delta, when
Judgement is LOS, when
Then judgement is NLOS, and wherein threshold delta is represented as follows: δ '=δ-L
2σ
ω 2/ (2 σ
2) (σ
ω 2=∫ x
2Φ (x) dx), L representes the zooming parameter of selected window function, and σ representes to measure under the LOS situation noise n
NoiseStandard deviation.
Among the said step C, be its weights assigned factor, when j reference mode is judged as los state, for the weighting factor of its distribution is k according to the NLOS identification result of j reference mode
j=1; When this reference mode is judged as the NLOS state, then do for the weighting factor of its distribution
Wherein i (i=1,2 ..., M) expression is measured for the i time, and j (j=1,2,3) representes j reference mode, r
JiThe i time distance measure representing j reference mode,
M the distance measure r that representes j reference mode
JiMean value.
Step D comprises the following steps:
D1. according to described weighting factor, construct the objective function of weighting factor least square method location algorithm for the pinpoint reference mode distribution of participation destination node:
Wherein (x y) representes position coordinates variable, X to X
j(x
j, y
j) j reference mode position coordinates of (j=1,2,3) expression, k
jJ weighting factor that reference mode is assigned to of (j=1,2,3) expression,
M the distance measure r that representes j reference mode
JiMean value;
D2. according to described objective function F (x), use the position coordinates that least square method is calculated destination node
Fig. 4 is the circuit block diagram of the gateway of the utility model.
As shown in Figure 4; Gateway is responsible for accomplishing the forwarding of super wideband wireless locating information and the transmission of configuration information, comprises microprocessor 9, storer 10, bus module 11, ultra-wideband antenna 12, power supply 13, communication interface 14, LCD 15, selection key 16.Gateway is by the mainboard 17 that comprises microprocessor 9, storer 10, bus module 11 and comprise that the base plate 18 of power supply 13, communication interface 14, LCD 15, selection key 16 and ultra-wideband antenna 12 form.
Fig. 5 is the circuit block diagram of the energy measuring receiver structure of the utility model.
As shown in Figure 5; The energy measuring receiver structure comprises ultra-wideband antenna 19, low noise amplifier LNA20, BPF. BPF21, squarer 22, integrator 23, sampling switch 24, decision device 25, is used to estimate to receive TOA time of arrival of ultra-broadband signal.
Fig. 6 is the reference mode circuit block diagram identical with destination node of the utility model.
As shown in Figure 6, the circuit block diagram of reference mode and destination node is identical, comprises processor 26, power supply 27, storer 28, ultra-broadband signal generation unit 29, energy measuring receiver structure 30, ultra-wideband antenna 31, on-off circuit 32.Power supply 27, storer 28, ultra-broadband signal generation unit 29, energy measuring receiver structure 30 all link to each other with processor 26.Processor 26 is with storer 28, and the data that responsible butt joint sending and receiving are sent are stored and handled, and control the operation of other unit.Ultra-broadband signal generation unit 29 is used for armed data-switching is become impulse ultra-wideband signal.On-off circuit 32 links to each other with ultra-broadband signal generation unit 29 ultra-wideband antenna 31 with energy measuring receiver structure 30, be used for the conversion between ultra-broadband signal and the reception.
Fig. 7 is the experiment scene arrangement plan of the utility model.
As shown in Figure 7, for the performance based on nonparametric NLOS discriminating and weighted least-squares method UWB location algorithm that is proposed is assessed, carried out experiment test.The length in experiment tunnel is respectively 100m, 3m and 3m, and with reference to the both sides that are deployed in the tunnel of the mutual dot interlace of knot, the horizontal range of adjacent datum node is 7.5m, and promptly the distance of homonymy datum node is 15m.If the actual position of destination node is (x
0, y
0), the estimated position does
Then positioning error is defined as E:
MS physical location (x
0, y
0) get (0.5,0.5), (3,1), (5,1.5), (7.5,2), (10,2.5), (15,2.75) respectively, estimate the position with this location algorithm
And calculate corresponding positioning error E.
Fig. 8 is the test data of experiment of the utility model.
Fig. 8 is the one group of test data of experiment that under experiment scene configuring condition shown in Figure 7, obtains, comprising the number of LOS datum node, and the actual position (x of destination node
0, y
0), the estimated position
Positioning error definition E.
Fig. 9 is the simulation result figure of the utility model.
As shown in Figure 9; In order to react the performance of this location algorithm more intuitively; Drawn out corresponding data and curves according to the data among Fig. 8 with Matlab software, the last figure that obtains is estimated position under actual position and the various LOS situation, has reflected the relation between positioning error and the LOS datum node number; Have the visible many more positioning errors of LOS datum node of figure more little, promptly the estimated position is more near physical location; The figure below that obtains is the cumulative distribution function (cdf) under the various LOS situation, has figure visible, when selected 3 datum nodes all are LOS, 90% positioning error is arranged less than 0.25m
2, when selected datum node all is NLOS, 50% positioning error is arranged less than 0.25m
2, 90% positioning error is less than 0.30m
2This shows that no matter under which kind of situation in Fig. 8, this system can reach high orientation precision.