The AAPS is comprised of a master AIS base station (BS), some slave base stations, a shipborn AIS equipment and an ASF correction system. The system configuration is shown in
Figure 1. The AIS BSs require time synchronization. The master AIS BS and other slave BSs takes turns sending signals. The shipborn AIS equipment receives the VHF AIS signal from the BSs. Furthermore, an ASF correction system is set up in order to improve the positioning accuracy. In this section, the principle of AAPS is investigated in detail.
Figure 1.
AIS autonomous positioning system configuration.
2.1. Position Algorithm
The ranging—mode position technique in AIS is based on a signal propagation delay technique. Signal propagation time position techniques can be categorized as time of arrival (TOA) [
12] and time difference of arrival (TDOA) [
13].
In the TOA technique, the fundamental objective is to estimate the vessel position by measuring the signal transmission delay from BSs to the vessel. Then the distance can be obtained by multiplying by the speed of light in the free space
c. However, in practice deviations exist in the calculated distance values due to the clock bias between the vessel and the AIS BSs. Thus the positioning equation used for the TOA technique should be written as:
where
is the measured distance from the
ith BS to the vessel; (
x,
y) is the unknown true position of the vessel; (
xi,
yi) is the known precise position of the
ith BS; Δ
t is the clock bias between the vessel and the BSs. The clock bias augments the two-dimensional location vector forming a three dimensional state vector, so there are three unknowns (
x,
y and Δ
t) in the above equation. Thus, to estimate the vessel position, more than three independent TOA measurements are required at the same time from the different BSs. Otherwise, there are only two channels in the shipborn AIS equipment, so it is very difficult to get more than three measurements at the same time.
In the TDOA technique, the vessel position is estimated by measuring the relative time delay between the master and the slave BSs. This does not require knowledge of the precise absolute time that the signal left the BS. Thus, the vessel doesn’t need to synchronize with the BSs, so it is more feasible for users in practice. The positioning equation used for the TDOA technique is given by:
where (
xi,
yi) and (
xm,
ym) are the known precise coordinates of the
ith slave BS and the master BS, respectively. As there are two unknowns in the position equation, at least two BS couples are required for the distance difference measurements at the same time.
Yi [
14,
15] analysed the mean square error for potential position algorithms based on the ranging-mode AIS used in AAPS. The performance of the TDOA technique is always superior to the TOA technique without time synchronization. Thus, the positioning algorithm for AAPS uses the TDOA technique. The details of the position estimation algorithm based TDOA technique in AAPS is presented in [
16].
2.2. Signal Measurement Technique
In the existing AIS, the maritime information is transmitted in the VHF mobile band. Specifically, two channels located at 161.975 and 162.025 MHz have been designated for this purpose [
17]. The two channels organized into time slots are shared by all BSs using time division multiple access (TDMA). Each channel has 2250 slots per minute. The data signaling rate is 9600 bits per second. The structure of the AIS signal is illustrated in
Figure 2.
Figure 2.
Structure of AIS signals.
Figure 2.
Structure of AIS signals.
The non-coherent demodulation technique is widely used in traditional shipborne AIS receivers for communication. There is no strict requirement for the decision time of each bit, thus the measurement accuracy of the bit transition is theoretically one bit period. This cannot meet the accuracy requirements of the positioning system, so the carrier frequency and the phase of the AIS signals from BSs should also be tracked in the receiver for AAPS. However, the problem of the cycle ambiguity is difficult to resolve, as the carriers of the AIS signals are not the continuous wave signals. The carrier-phase measurement alone is not feasible. The receivers used in AAPS estimate the carrier frequency, the carrier phase and the time of bit transition. Then they can provide precise times for the bit transitions with assistance of carrier tracking. The corresponding block diagram is shown in
Figure 3. The method is called GMSK demodulation based on carrier phase tracking. The key technology used in this demodulation method is the carrier extraction and tracking, which has been investigated in [
18].
Figure 3.
Block diagram of GMSK demodulation in an AIS receiver.
Figure 3.
Block diagram of GMSK demodulation in an AIS receiver.
GMSK is a noncausal system [
19]. The step response function of Gauss filter is given as:
where
ω denotes the bandwidth. The simulation curves with different
ω are shown in
Figure 4.
Figure 4.
Step response of a Gauss filter.
Figure 4.
Step response of a Gauss filter.
It can be seen from
Figure 4 that the location of each zero crossing of the baseband signal corresponds to the location of the bit transition of serial data signals. Each zero crossing of the baseband signal can be detected in GMSK demodulation based on carrier phase tracking. Therefore, the precise transition time from one bit to the next bit in the serial data signals can be measured in GMSK demodulation based on a FPGA circuit.
In the AIS receiver based on FPGA, a special counter is used to record the bit transition time of serial data signals. It starts from zero at the beginning of a slot. However, the measurement error of the single measurement according to the counter cannot be ignored due to the existence of random errors and quantization errors. In order to reduce the measurement error of the single sampled value, a number of sampled values in one slot are measured. The first sampled value of time delay measurement corresponds to the last bit of the start flag (7E). After that, each bit is sampled in one time slot. There are 193 sample datapoints from the last bit of the start flag to the end flag in one time slot. The measurement of time delay
Ti corresponding to the
ith bit is given by:
where
B is the time of one bit, which is calculated according to the bit synchronization signal of the AIS receiver; Δ
T is the time delay from the BS to the vessel. 8, 24 and 7 are the lengths of the ramp up, the training sequence and the first seven bits of the start flag shown in
Figure 2. The unit of these parameters in Equation (3) is the clock cycle count with a frequency of 174.72 MHz. The use of a clock frequency of 174.72 MHz is a compromise. The maximum frequency of the designed AIS receiver based on FPGA is less than 200 MHz. As the intermediate frequency of AIS signals is 455 KHz, the clock frequency should be the integer times of the intermediate frequency of AIS signal, so to obtain higher time resolution, 174.72 MHz is chosen as the clock frequency.
To evaluate the accuracy of the measurement, a total of 193 sample data in one slot is analyzed during 700 time slots. The curve of the standard deviation of time delay measurement according to every sample data is drawn in
Figure 5. It can be seen that the standard deviation of the measurement corresponding to every sample bit is very large. The average reaches 179.9156 clock cycles count, corresponding to 1.04 µs, so the accuracy of the distance measurement is 312 m. It cannot meet the requirements of positioning accuracy.
Figure 5.
Standard deviation of measurement for sample data.
Figure 5.
Standard deviation of measurement for sample data.
Under the ideal condition, all of the time data on the transition from one bit to the next bit in a slot should be on a straight line, with a slope of the baud rate of replication synchronization bit. Thus, the least squares method is used for the linear fit of all of the time data on the transition from one bit to the next bit in a slot. The time of the bit transitions can be measured more accurately by reducing the random turbulence. According to the Equation (3), the first order coefficient after the linear fitting is the corrected value of the bit synchronization signal, and the constant is the estimation of the time delay. Experimental results after the linear fitting based on the least square method are shown in
Figure 6. In
Figure 6, the horizontal axis represents the number of the slot and the vertical axis represents the constant calculated by the linear fitting corresponding to each slot, which indicates the estimation of the time delay. The standard deviation of measurement corresponding to the time delay after the linear fitting is 42.70 clock cycle counts, corresponding to 0.24 µs, so the accuracy of the distance measurement is reduced to 73 m by the linear fitting based on the least squares method.
Figure 6.
Experimental results after the linear fitting based on the least squares method.
Figure 6.
Experimental results after the linear fitting based on the least squares method.
2.3. Geometric Dilution of Precision
Accuracy is a crucial parameter for evaluating a positioning system. The accuracy depends on the measurement accuracy and the geometry of the BSs, which is described by the GDOP [
20]. In general, positioning errors can be expressed as the product of the GDOP at the vessel location and the root mean square deviation of the measurement errors [
21]. The measurement error has been discussed above in
Section 2.2. The GDOP of AAPS will be investigated below. The existing AIS BSs have already been established for a communication system, but the layout of the location of AIS BSs is not suitable for positioning, without considering the effects of the geometry, so we need to evaluate the GDOP of AIS BSs in AAPS.
The location equation (Equation (2)) is nonlinear. We obtain the linear location equation using a Taylor-series, keeping only terms below second order:
where
and
denote the measurement value and its estimation for the
ith pair of BSs, respectively. Stacking all the measurements from different BSs, the above Equation (4) can be written in matrix form as follows:
where:
where
γi is the estimated azimuth angle directed from the estimated position of the vessel to the
ith BS. As:
According to Equation (8), then the GDOP for AAPS can be calculated by:
We take the established verification system of AAPS as an example in this paper. The location information of BSs in the verification system is shown in
Table 1.
Table 1.
Location information of base stations.
Table 1.
Location information of base stations.
BS Name | Latitude | Longitude |
---|
Lingjing | 38°50'21.309"N | 121°30'46.005"E |
Heishijiao | 38°51'58.689"N | 121°33'05.779"E |
Fujiazhuang | 38°51'52.660"N | 121°36'49.543"E |
The GDOP of the established verification system of AAPS is evaluated as shown in
Figure 7. GDOP is calculated based on the vessel’s bearings to each of the BSs. Since the GDOP can be interpreted as a multiplier on the measurement error, lower numbers are better. According to the conclusions of
Section 2.2, the accuracy of the distance measurement is 73 m by the linear fitting based on the least squares method, therefore, the positioning precision for the verification system of AAPS is about 100 m (2
σ) in the area where GDOP is less than 1.5.
Figure 7.
GDOP for the verification system of AAPS.
Figure 7.
GDOP for the verification system of AAPS.
2.4. Time Synchronization Technique
The most prominent advantage of the TDOA technique is that the vessel doesn’t need to synchronize with the BSs. This is more feasible for users in reality. Even so, it still requires synchronization between all the BSs. As part of their definition as a TDMA communications link, AIS transmissions are synchronized. The BSs are synchronized with the coordinated universal time (UTC) given by the GNSS receivers. According to the time synchronization specification in IEC 62320-1 [
22], the AIS BS has a time source synchronization to UTC better than 50 µs. This obviously cannot meet the accuracy requirement for the position system.
Time synchronization in AAPS could be achieved using the GNSS disciplined rubidium clock technology [
23]. The principle of the GNSS disciplined rubidium clock technology is described briefly below. The phase difference between pulse per seconds (PPS) given by the rubidium clock and the GNSS time receiver is compared in real-time. According to the phase difference, the phase locked loop adjusts the rubidium clock frequency, which is used for time synchronization of the BSs. The technology of GNSS disciplined rubidium clock accomplishes a couple of things: One is making the rubidium clock output frequency of 10 MHz more accurate. The other is making the PPS output by the rubidium clock more accurate and stable, which is aligned with the PPS output from the GNSS time receiver.
When the phase of the PPS output by the rubidium clock is aligned to the PPS output from the GNSS receiver, it’s called synchronization mode. It could be used to provide the time synchronization up until the point in time that the GNSS signal is lost. When the GNSS signal is not present, it is called holdover mode. After this point, the signal would remain in tolerance for 24 h for a rubidium clock in the holdover mode. When the GNSS signal is received again, it can quickly resume the high precision synchronization.
In the experiment, the rubidium clock had been locked to the GNSS time receiver for the first 2 h, but lost it later.
Figure 8 illustrates the experimental results in the holdover mode for 24 h, in which each division along the horizontal axis corresponds to one minute, and each division along the vertical axis corresponds to 5 ns.
Figure 8.
Experimental results of GNSS disciplined rubidium clock technology.
Figure 8.
Experimental results of GNSS disciplined rubidium clock technology.
When the rubidium clock is not locked to any reference, such as the GNSS time receiver, the phase difference between the PPS outputted by the rubidium clock and the GNSS receiver is changing with time as denoted by the blue dots in
Figure 8. The frequency stability of the rubidium clock is influenced by environmental factors, especially the temperature. Thus, the laboratory must meet strict work environment requirements. The phase deviation is less than 300 ns within 24 h when the rubidium clock is running individually. This is consistent with performance levels given by [
24]. The phase differences are represented by the green and the red dots in the synchronization mode and the holdover mode of GNSS disciplined rubidium clock technology, respectively. In the synchronization mode, the PPS given by the rubidium clock is always aligned with the PPS output from the GNSS time receiver. Therefore, the phase difference of a GNSS disciplined rubidium clock is very small. In the holdover mode, the rubidium clock has already been disciplined by GNSS. Although without the PPS output from the GNSS time receiver, the accurate frequency drift correction model can be used, whose parameters are obtained according to long-term statistics of the character of the rubidium clock and short-term statistics derived from synchronization mode data. The experimental results in
Figure 8 show that the phase difference is less than 30 ns within 24 h in the holdover mode using the GNSS disciplined rubidium clock technology. However, the period of validity of short-term statistical data in the correction model is 24 h. The phase error will increase quickly after 24 h.
2.5. Additional Secondary Factor Correction
Due to the signals propagation over paths of varying conductivity, topography and weather, there is a significant factor limiting the accuracy of the positioning system. Thus, in the ground-based positioning system, ASF correction is generally employed. In the investigation of ASF corrections, the simplest approach is a theoretical model construction, but in that case the calculation error is large and it cannot meet the accuracy requirements of AAPS. Thus, the ASF correction system has been set up [
25]. In the ASF correction system, the ASF transmitter is time synchronized with the ASF receiver. The ASF transmitter sends its precise position and ranging code to the ASF receiver. Therefore, the theoretical distance
RASF0 from the ASF transmitter to the receiver can be calculated according to the precise position information of the ASF transmitter and the ASF receiver. The actual measurement distance
RASF can be obtained by measuring the propagation time of signals from the ASF transmitter to the ASF receiver. The ASF correction value is calculated by comparing the distance error between the theoretical and actual measurement distance from the ASF transmitter to the receiver. Then the real-time correction value is broadcasted to the vessels in the area by the AIS base station to improve the positioning accuracy. The ASF correction system is connected with AAPS as shown in
Figure 9.
Figure 9.
Connection diagram between the ASF correction system and AAPS.
Figure 9.
Connection diagram between the ASF correction system and AAPS.
The measurement variances of ASF and AIS receiver are denoted by
and
, respectively. The theoretical and actual measurement distances from the AIS BS to the vessel are denoted by
RAIS0 and
RAIS, respectively. The relationship between them is given by:
where
RASF,
RASF0 and
σASF obtained by the ASF correction system are all sent to the AIS master BS through a serial port according to the interface protocol. The AIS BS sends these parameters to the vessels in the same area. The AAPS receiver corrects the measured distance
RAIS according to Equation (11) using these ASF parameters and its own
σAIS: