A METHOD FOR EVALUATING QUALITY OF RADIO FREQUENCY SIGNALS FOR STELLITE NAVIGATION SYSTEM
20170359132 · 2017-12-14
Assignee
Inventors
Cpc classification
G01S19/396
PHYSICS
B60W2556/50
PERFORMING OPERATIONS; TRANSPORTING
International classification
Abstract
A method for evaluating the quality of a radio frequency signal for a satellite navigation system, the method comprising: sampling the payload radio frequency signal of a satellite to obtain an intermediate frequency signal, and filtering the signal; downconverting the filtered intermediate signal to obtain the corresponding actual baseband signal; generating the signal components of the ideal baseband signal branches on the basis of the obtained actual baseband signal and the signal system thereof; reproducing the ideal baseband signal, which is used to evaluating the actual baseband signal; establishing a correlation function between the actual baseband signal to be evaluated and the ideal baseband signal, and obtaining through corresponding calculations a series of quality evaluation indexes such as spurious transmission in the band and related loss, thereby enabling the evaluation of the quality of the radio frequency signal. The present invention clarifies the relation between signal quality indexes and the capturing, tracking, and demodulation performance of a signal, and can clearly and regularly evaluate the effect of the signal quality on navigation performance.
Claims
1. A method for assessing a quality of a radio frequency signal for a satellite navigation system, the method comprising: 1) sampling a radio frequency signal of satellite payload to obtain an intermediate frequency signal, then filtering the intermediate frequency signal; 2) downconverting a filtered intermediate frequency signal to obtain a corresponding actual baseband signal; 3) generating respective branch signal components of an ideal baseband signal having a same sampling rate and one-code-cycle length on the basis of the actual baseband signal obtained in 2) and a signal framework thereof; and by utilizing the actual baseband signal and the respective branch signal components of the ideal baseband signal, reproducing an ideal baseband signal to be used to assess the actual baseband signal; 4) establishing a correlation function between an actual baseband signal to-be-assessed and the ideal baseband signal as follows:
2. The method of claim 1, wherein the in-band spurious and the related losses are calculated according to the following formulas:
r.sub.ST[dBc]=10×log.sub.10(1−max(|CCF(τ)|.sup.2)) wherein r.sub.ST[dB] denotes a dB value of in-band spurious,
L.sub.CCF[dB]=20×log.sub.10(max(|CCF(τ)|.sup.2)) wherein L.sub.CCF[dB] denotes a dB value of correlation loss.
3. The method of claim 1, further comprising a process for quantitative assessment of pseudo-code ranging error and code tracking performance in 5), and the process is specifically as follows: for pseudo-code-ranging error, as one metric, it may be obtained as follows: firstly, according to correlation peaks of the ideal baseband signal and the actual baseband signal, calculating a corresponding discriminant function, then performing linear fitting to zero-crossing point offset of the discriminant function; the correspondingly obtained data can be used to reflect the pseudo-code ranging error; and for code tracking performance, as another metric, it may be obtained as follows: firstly, according to correlation peaks of the ideal baseband signal and the actual baseband signal, calculating corresponding discriminant functions respectively, then performing linear fitting to zero-crossing point slopes of the two discriminant functions respectively; the correspondingly obtained ratio of the slops can be used to reflect the code tracking performance.
4. The method of claim 1, further comprising a process for calculations of a signal-component phase difference and a signal-component magnitude difference in 5), and the calculations are executed respectively according to the following formulas:
Δθ.sub.i=θ.sub.i.sup.t−θ.sub.i wherein Δθ.sub.i denotes a phase difference of an i-th signal component of the actual baseband signal, θ.sub.i.sup.t denotes a phase angle of the maximum-modulus point in the cross-correlation sequence between the i-th signal components of the actual baseband signal and the ideal baseband signal, and θ.sub.i denotes a design value of the carrier phase of the i-th signal component of the ideal baseband signal;
5. The method of any of claim 1, further comprising a process for calculations of signal-component code phase consistency and frequency-point code phase consistency in 5), and the calculations are executed respectively according to the following formulas:
Δε.sub.i-th component=ε.sub.i-th component−ε.sub.0 wherein Δε.sub.i-th component denotes a code phase deviation of an i-th signal component of the actual baseband signal, ε.sub.i-th component denotes a zero-crossing point offset value of the discriminant function of the i-th signal component, and ε.sub.0 denotes a zero-crossing point offset value of the discriminant function of the reference signal component;
Δε.sub.i-th frequency-point=ε.sub.i-th frequency-point−ε.sub.1 wherein Δε.sub.i-th frequency-point denotes a code phase deviation of an i-th frequency point of the actual baseband signal, ε.sub.i-th frequency-point denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the i-th frequency point, and ε.sub.1 denotes a zero-crossing point offset value of the discriminant function of aggregate signal at the reference frequency point.
6. The method of claim 5, further comprising processes of calculations of a series of metrics comprising: signal power spectrum, spectrum distortion, constellation diagram, error vector magnitude, and carrier leakage in 5).
7. The method of claim 2, further comprising a process for quantitative assessment of pseudo-code ranging error and code tracking performance in 5), and the process is specifically as follows: for pseudo-code-ranging error, as one metric, it may be obtained as follows: firstly, according to correlation peaks of the ideal baseband signal and the actual baseband signal, calculating a corresponding discriminant function, then performing linear fitting to zero-crossing point offset of the discriminant function; the correspondingly obtained data can be used to reflect the pseudo-code ranging error; and for code tracking performance, as another metric, it may be obtained as follows: firstly, according to correlation peaks of the ideal baseband signal and the actual baseband signal, calculating corresponding discriminant functions respectively, then performing linear fitting to zero-crossing point slopes of the two discriminant functions respectively; the correspondingly obtained ratio of the slops can be used to reflect the code tracking performance
8. The method of claim 2, further comprising a process for calculations of a signal-component phase difference and a signal-component magnitude difference in 5), and the calculations are executed respectively according to the following formulas:
Δθ.sub.i=θ.sub.i.sup.t−θ.sub.i wherein Δθ.sub.i denotes a phase difference of an i-th signal component of the actual baseband signal, θ.sub.i.sup.t denotes a phase angle of the maximum-modulus point in the cross-correlation sequence between the i-th signal components of the actual baseband signal and the ideal baseband signal, and θ.sub.i denotes a design value of the carrier phase of the i-th signal component of the ideal baseband signal;
9. The method of claim 3, further comprising a process for calculations of a signal-component phase difference and a signal-component magnitude difference in 5), and the calculations are executed respectively according to the following formulas:
Δθ.sub.i=θ.sub.i.sup.t−θ.sub.i wherein Δθ.sub.i denotes a phase difference of an i-th signal component of the actual baseband signal, θ.sub.i.sup.t denotes a phase angle of the maximum-modulus point in the cross-correlation sequence between the i-th signal components of the actual baseband signal and the ideal baseband signal, and θ.sub.i denotes a design value of the carrier phase of the i-th signal component of the ideal baseband signal;
10. The method of claim 2, further comprising a process for calculations of signal-component code phase consistency and frequency-point code phase consistency in 5), and the calculations are executed respectively according to the following formulas:
Δε.sub.i-th component=ε.sub.i-th component−ε.sub.0 wherein Δε.sub.i-th component denotes a code phase deviation of an i-th signal component of the actual baseband signal, ε.sub.i-th component denotes a zero-crossing point offset value of the discriminant function of the i-th signal component, and denotes a zero-crossing point offset value of the discriminant function of the reference signal component;
Δε.sub.i-th frequency-point=ε.sub.i-th frequency-point−ε.sub.1 wherein Δε.sub.i-th frequency-point denotes a code phase deviation of an i-th frequency point of the actual baseband signal, ε.sub.i-th frequency-point denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the i-th frequency point, and denotes a zero-crossing point offset value of the discriminant function of aggregate signal at the reference frequency point.
11. The method of claim 3, further comprising a process for calculations of signal-component code phase consistency and frequency-point code phase consistency in 5), and the calculations are executed respectively according to the following formulas: Δε.sub.i-th component=ε.sub.i-th component−ε.sub.0 wherein Δε.sub.i-th component denotes a code phase deviation of an i-th signal component of the actual baseband signal, ε.sub.i-th component denotes a zero-crossing point offset value of the discriminant function of the i-th signal component, and ε.sub.0 denotes a zero-crossing point offset value of the discriminant function of the reference signal component;
Δε.sub.i-th frequency-point=ε.sub.i-th frequency-point−ε.sub.1 wherein Δε.sub.i-th frequency-point denotes a code phase deviation of an i-th frequency point of the actual baseband signal, ε.sub.i-th frequency-point denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the i-th frequency point, and denotes a zero-crossing point offset value of the discriminant function of aggregate signal at the reference frequency point.
12. The method of claim 4, further comprising a process for calculations of signal-component code phase consistency and frequency-point code phase consistency in 5), and the calculations are executed respectively according to the following formulas:
Δε.sub.i-th component=ε.sub.i-th component−ε.sub.0 wherein Δε.sub.i-th component denotes a code phase deviation of an i-th signal component of the actual baseband signal, ε.sub.i-th component denotes a zero-crossing point offset value of the discriminant function of the i-th signal component, and ε.sub.0 denotes a zero-crossing point offset value of the discriminant function of the reference signal component;
Δε.sub.i-th frequency-point=ε.sub.i-th frequency-point−ε.sub.1 wherein Δε.sub.i-th frequency-point denotes a code phase deviation of an i-th frequency point of the actual baseband signal, ε.sub.i-th frequency-point denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the i-th frequency point, and ε.sub.1 denotes a zero-crossing point offset value of the discriminant function of aggregate signal at the reference frequency point.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0028]
DETAILED DESCRIPTION OF THE EMBODIMENTS
[0029] For further illustrating the invention, experiments detailing a method for assessing a quality of a radio frequency signal for a satellite navigation system are described below. It should be noted that the following examples are intended to describe and not to limit the invention.
[0030]
[0031] Step 1
[0032] Sampling a RF signal of satellite payload to obtain an IF signal, then filtering the IF signal.
[0033] Specifically, firstly, acquiring a satellite signal to-be-assessed by sampling a RF signal of satellite payload, for example, by intercepting a tow-code-cycle signal; according to the sampling rate and the center frequency of the RF signal, calculating the IF frequency of the sampling signal, and filtering the RF-sampling signal, for example, through an ideal block filter, to obtain the IF signal to-be-assessed. Then, according to the signal framework of the signal to be assessed, generating a one-code-cycle ideal baseband aggregate signal and respective branch signal components thereof.
[0034] Step 2
[0035] Downconverting the filtered IF signal to obtain a corresponding actual baseband signal.
[0036] Step 3
[0037] On the basis of the actual baseband signal to-be-assessed obtained in step 2 and the signal framework thereof, generating respective branch signal components of an ideal baseband signal having a same sampling rate and one-code-cycle length; then, by utilizing the actual baseband signal and the respective branch signal components of the ideal baseband signal, reproducing an ideal baseband signal to be used to assess the actual baseband signal.
[0038] Step 4
[0039] Establishing a correlation function between an actual baseband signal to-be-assessed and the ideal baseband signal, to be used in subsequent processes for calculation of a plurality of assessment metrics, the correlation function being shown as below:
where, s.sub.rec(t) denotes a variable that represents the actual baseband signal s.sub.rec and is determined by a time variable t, s.sub.0(t) denotes a variable that represents the ideal baseband signal s.sub.0 and is determined by the time variable t, t denotes the time variable, τ denotes a correlation time delay, Tp denotes a length of time of the actual baseband signal s.sub.rec, s.sub.0*(t−τ) denotes a conjugate operation executed on a variable that represents the ideal baseband signal s.sub.0 and is determined by both the time variable t and the correlation time delay τ, CCF(τ) denotes a variable that represents a cross-correlation sequence CCF between the actual baseband signal s.sup.rec and the ideal baseband signal s.sub.0 and is determined by the correlation time delay τ.
[0040] Step 5
[0041] By utilizing the correlation function established in above step, calculating quality-assessment metrics related to in-band spurious and correlation loss, thereby accomplishing quality assessment procedure of the RF signal.
[0042] According to a preferred embodiment of the present invention, for the calculation of in-band spurious, it is preferably executed according to the following formula:
r.sub.ST[dBc]=10×log.sub.10(1−max(|CCF(τ)|.sup.2))
where, r.sub.ST[dBc] denotes a dB value of in-band spurious, CCF(τ) denotes a variable that represents a cross-correlation sequence between an actual baseband signal and an ideal baseband signal and is determined by correlation time delay τ;
[0043] Likewise, for the calculation of correlation loss, it is preferably executed according to the following formula:
L.sub.CCF[dBc]=20×log.sub.10(max(|CCF(τ)|.sup.2))
where, L.sub.CCF[dBc] denotes a dB value of correlation loss, CCF(τ) denotes a variable that represents a cross-correlation sequence between an actual baseband signal and an ideal baseband signal and is determined by correlation time delay τ.
[0044] In order to perform a more comprehensive and accurate assessment of a RF signal of satellite navigation system, apart from the above basic metrics such as in-band spurious and correlation loss, the present invention also designs several other relevant quality-assessment metrics, of which the calculation or processing processes are explained respectively and specifically as follows:
[0045] Pseudo-Code Ranging Error
[0046] According to a preferred embodiment of the present invention, this metric preferably may be obtained by a process as follows: firstly, according to a correlation peak of the ideal baseband signal, calculating a corresponding discriminant function, then using fitting, for example linear fitting, to obtain zero-crossing point offset of the discriminant function; thus the pseudo-code ranging error of the navigation signal is:
e.sub.DB=ε.sub.b(σ)×c
where, e.sub.DB denotes ranging deviation caused by signal distortion, with unit m, ε.sub.b(σ) denotes zero-crossing point offset of the discriminator, with unit s, c denotes electromagnetic wave propagation speed.
[0047] Code Tracking Performance
[0048] According to a preferred embodiment of the present invention, this metric preferably may be obtained by a process as follows: firstly, according to correlation peaks of both the ideal baseband signal and the actual baseband signal, calculating corresponding discriminant functions respectively, then using fitting, for example linear fitting, to obtain zero-crossing point slopes of the two discriminant functions respectively; thus, the ratio of the two slopes can be used to reflect the code tracking performance.
[0049] Signal-Component Phase Difference
[0050] According to a preferred embodiment of the present invention, this assessment metric may be calculated according to the following formula:
Δθ.sub.i=θ.sub.i.sup.t−θ.sub.i
where, Δθ.sub.i denotes a phase difference of an i-th signal component of the actual baseband signal, θ.sub.i.sup.t denotes a phase angle of the maximum-modulus point in the cross-correlation sequence between the i-th signal components of the actual baseband signal and the ideal baseband signal, θ.sub.i denotes a design value of the carrier phase of the i-th signal component of the ideal baseband signal;
[0051] Signal-Component Magnitude Difference
[0052] According to a preferred embodiment of the present invention, this assessment metric may be calculated according to the following formula:
where, Δp.sub.i[dB] denotes the dB value of the magnitude difference of the i-th signal component of the actual baseband signal, p.sub.i denotes a design value of the proportion of the i-th signal component of the actual baseband signal to the total power, p.sub.i.sup.t denotes the actual value of the proportion of the i-th signal component of the actual baseband signal to the total power.
[0053] Signal-Component Code Phase Consistency
[0054] According to a preferred embodiment of the present invention, this assessment metric may be calculated according to the following formula:
Δε.sub.i-th component=ε.sub.i-th component−ε.sub.0
where, Δε.sub.i-th component denotes a code phase deviation of an i-th signal component of the actual baseband signal ε.sub.i-th component denotes a zero-crossing point offset value of the discriminant function of the i-th signal component, ε.sub.0 denotes a zero-crossing point offset value of the discriminant function of the reference signal component.
[0055] Frequency-Point Code Phase Consistency
[0056] According to a preferred embodiment of the present invention, this assessment metric may be calculated according to the following formula:
Δε.sub.i-th frequency-point=ε.sub.i-th frequency-point−ε.sub.1
where, Δε.sub.i-th frequency-point denotes a code phase deviation of an i-th frequency point of the actual baseband signal, ε.sub.i-th frequency-point denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the i-th frequency point, ε.sub.1 denotes a zero-crossing point offset value of the discriminant function of the aggregate signal at the reference frequency point.
[0057] In addition, the quality of the RF signal assessment method according to the present invention preferably further comprises processes for involving a series of metrics—including signal power spectrum, spectrum distortion, constellation diagram, error vector magnitude (EVM), carrier leakage, etc.—in the quality-assessment procedure; the calculation and processing processes of these metrics are explained respectively and specifically as follows:
[0058] Plotting Signal Power Spectrum
[0059] Firstly, an FFT algorithm, for example, may be used to analyze the Fourier transform S.sub.rec(f) of the actual baseband signal s.sub.rec(t), thus the signal power spectrum is
and accordingly, the power spectrum can be plotted, where, T denotes the length of time of the sampling signal, here it is a one-code-cycle length, Fs denotes the signal sampling rate.
[0060] In addition, a power spectrum envelope may be extracted. Based on inverse Fourier transform of the product of the signals in frequency domain, time-domain correlation is calculated, then the time-domain correlation is time-windowed, with the time-window width generally being taken as 5 code-chip widths (the size of the time window determines the fineness degree of the power spectrum envelope; the wider the time window, the finer the power spectrum envelope curve), and subsequently, through fast Fourier transform of the time-windowed time-domain correlation, the power spectrum envelope is obtained.
[0061] Signal Spectrum Distortion
[0062] Firstly, the actual baseband signal spectrum S.sub.rec(f) and the ideal baseband signal spectrum S.sub.s(f) are calculated respectively, then the spectrum distortion may be calculated according to the following formula, and accordingly, a distortion diagram can be plotted.
[0063] Plotting Constellation Diagram
[0064] Firstly, both the ideal baseband signal and the actual baseband signal undergo magnitude normalization, that is, the actual baseband signal is divided by its average magnitude, thereby obtaining the normalized actual signal, and the normalized ideal signal can be obtained likewise.
[0065] Then, in a coordinate, taking the I branch of the normalized actual signal as horizontal axis, the Q branch as vertical axis, a scatter diagram of the actual signal can be plotted, with the range of the horizontal and vertical axes being [X.sub.min, X.sub.max, Y.sub.min, Y.sub.max], where, X.sub.min denotes the minimum abscissa of the signal points, X.sub.max denotes the maximum abscissa of the signal points, Y.sub.min denotes the minimum ordinate of the signal points, Y.sub.max denotes the maximum ordinate of the signal points. The plane of the axes is divided into N*N equal-area square grids, for example, taking N=100. The number of signal points in each square grid is counted, and the color depth of each square grid is set according to the number of points; Likewise, the ideal signal points are plotted in the same diagram for comparison.
[0066] Error Vector Magnitude
[0067] Firstly, both the actual baseband signal to-be-assessed and the ideal baseband signal may be normalized, so that the average magnitude of the normalized signal is 1;
[0068] Next, the ideal baseband signal points are numbered according to their distribution positions on the coordinate and denoted as V.sub.k=1.sub.k+jQ.sub.k, then according to the distances from the actual baseband signal points to these K ideal baseband signal points, all of the actual baseband signal points and the ideal baseband signal points that pertain to a closest distance are sorted into one class and denoted as below:
{(I.sub.rec(m), Q.sub.rec(m))|(I.sub.ideal(m), Q.sub.ideal(m))=V.sub.k, m=1,2, . . . , M.sub.k}
where, M.sub.k denotes, in a sample, the number of the sampling points corresponding to the signal vector V.sub.k, (I.sub.rec(m), Q.sub.rec(m)) denotes the actual signal points, (I.sub.ideal(m), Q.sub.ideal(m)), denotes the ideal signal points;
[0069] Consequently, the vector error magnitude of the to-be-assessed signal corresponding to the signal vector Vie is calculated according to the following formula:
EVM.sub.k=√{square root over ((I.sub.rec.sup.k−I.sub.k).sup.2+(Q.sub.rec.sup.k−Q.sub.k).sup.2)}
where, (I.sub.rec.sup.k, Q.sub.rec.sup.k) denotes the k-class of actual signal points, i.e., the actual signal points corresponding to the ideal signal vector V.sub.k.
[0070] Carrier Leakage
[0071] The carrier leakage is calculated according to the following formula, where, G(f) denotes the calculated signal power spectrum in step 6), f.sub.c denotes the carrier frequency, with unit Hz, [f.sub.l, f.sub.h] denotes the frequency range of the specified frequency band:
[0072] For the above-mentioned series of metrics, including signal power spectrum, spectrum distortion, constellation diagram, error vector magnitude, carrier leakage, etc., their basic calculation and processing formulas are common knowledge well known to those skilled in the art, and thus will not be described here.
[0073] Unless otherwise indicated, the numerical ranges involved in the invention include the end values. While particular embodiments of the invention have been shown and described, it will be obvious to those skilled in the art that changes and modifications may be made without departing from the invention in its broader aspects, and therefore, the aim in the appended claims is to cover all such changes and modifications as fall within the true spirit and scope of the invention.