CN1946069A - Detecting and analyzing method for multi system frequency shift key control signal - Google Patents
Detecting and analyzing method for multi system frequency shift key control signal Download PDFInfo
- Publication number
- CN1946069A CN1946069A CNA2006101140851A CN200610114085A CN1946069A CN 1946069 A CN1946069 A CN 1946069A CN A2006101140851 A CNA2006101140851 A CN A2006101140851A CN 200610114085 A CN200610114085 A CN 200610114085A CN 1946069 A CN1946069 A CN 1946069A
- Authority
- CN
- China
- Prior art keywords
- signal
- centerdot
- frequency
- sequence
- delta
- 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.)
- Granted
Links
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Complex Calculations (AREA)
Abstract
This invention relates to automatic identification to communication signals and analysis technology aiming at testing and analyzing multi-bit frequency shift key controlled signals characterizing in comparing the mean value of the absolute value of the variance rate sequence of the normalized instant frequencies with the set threshold value to judge MFSK signal, when analyzing, Smooth filtration is used in reducing shake of the MFSK signal transient frequency then the entropy principle is used in analyzing the histogram for the distribution of transient frequencies to estimate the state number of code elements and frequency sphere with an equal probability distribution principle then to use the equal frequency distance to get more accurate code element state number, frequency sphere and central frequency.
Description
Technical field
The invention belongs to the signal of communication analysis field, the particularly automatic identification and the analysis of multi system frequency shift key control (MFSK) signal in the digital modulation mode.
Background technology
Modulation system identification automatically is the important step that electronic warfare (signal reconnaissance), radio monitoring, adaptive communications etc. are used, and is an important topic in non-collaboration communication field, and the various theories of signal of communication analysis have obtained extensive use in this direction.Frequency shift keying modulation system (FSK) is being played the part of important role in the shortwave/ultra short wave communication of army, government and amateur radio communication, being a kind of important digital signal modulation mode, is the important content of radio signal monitoring and military communication scouting to the detection of such signal and analysis.The present invention relates generally to the detection and the analysis of frequency shift keyed signals, the i.e. analysis of the identification of frequency shift keyed signals and order of modulation.
In recent years, at multi system frequency shift key control (MFSK) input and analysis, following several method has been proposed:
1. the method for analyzing based on signal envelope analysis and instantaneous frequency
People such as E.E.Azzouz in Automatic modulation recognition (Journal Of The FranklinInstitute-Engineering And Applied Mathematics, Mar 1997), proposed to use signal envelope spectrum peak γ
MaxPermanent envelope (frequency modulated signal) signal and non-constant envelope signal are classified, use the standard deviation of normalization instantaneous frequency absolute value
AfTo BFSK signal and the classification of QFSK signal, step is as follows:
1). the signal that receives is converted into analytic signal, removes the weak signal section;
2). calculate zero center normalization instantaneous amplitude a
Cn
3). the spectrum peak γ of signal calculated envelope
Max, should be worth with predefine thresholding τ
γCompare, if γ
Max<τ
γThen received signal is judged to be frequency modulated signal;
4). calculate its zero center normalization instantaneous frequency f by the instantaneous phase of analytic signal
Cn
5). calculate the standard deviation of normalization instantaneous frequency absolute value
Af, should be worth with predefine thresholding τ
σ afCompare, if σ
α f<τ
σThen be judged to be the BFSK signal, otherwise be judged to be the QPSK signal.
The problem of this method is:
● because factor affecting such as frequency selective fadings, the MFSK sample of signal of many actual reception is not a constant envelope signal;
● this method can only be discerned BFSK and QFSK signal.
Cao Zhi has just waited the people " to need not the digital signal automatic Modulation Recognition method commonly used of priori " and also related in (patent No. 02123627.5) detection and the analysis of MFSK signal in relevant patent.Relevant step is as follows:
1). the R parameter of signal calculated sample, if this parameter less than the predefine thresholding, then is classified as constant envelope signal with this signal;
2). to constant envelope signal, the instantaneous phase of signal calculated s (n)=x (n)+iy (n):
3). go volume folded to instantaneous phase, add the phasing sequence for the phase sequence of mould 2 π:
θ(n)=(n)+C
k(n)
4). the instantaneous frequency f of signal calculated
1(n):
5). to f
1(n) carry out the statistics of absolute value, obtain statistics number C greater than pi/2
p, satisfy C
p<T
CpSignal determining be the signal that does not have phase hit, i.e. MFSK signal;
6). to f
1(n) carry out interpolation processing, obtain f
2(n):
7). to f
2(n) carry out normalized and obtain normalized frequency f (n):
8). the distribution to f (n) is added up, and obtains the peak value number N of f (n)
f, if N
f=2, then be judged to be the BFSK signal, otherwise be judged to be the FM signal.
The problem of this method is:
● because factor affecting such as frequency selective fadings, the MFSK sample of signal of many actual reception is not a constant envelope signal;
● this patent does not provide the peak value number N of the Distribution Statistics of calculating f (n)
fMethod, and can only discern the BFSK signal.
2. based on the method for zero passage sampling
People such as S.-Z.Hsue are at paper Automatic Modulation Classification Using Zero Crossing (Radar andSignal Processing, IEE Proceedings For, vol.37, No.6, pp.459-464, Dec.1990.K.) proposed in to use the zero passage sampler to extract the inverse of the instantaneous frequency of signal, detect the MFSK signal by analyzing this variance reciprocal then, analyze the method for MFSK signal by analyzing the instantaneous frequency distribution histogram, step is as follows:
1). received signal, the time point sequence x (i) of the zero passage of usefulness zero passage sampler tracer signal;
2). by time x (i) calculated value sequences y (i)=x (the i+1)-x (i) and value sequence z (i)=y (the i+1)-y (i) of signal zero passage;
3). set thresholding τ
z, will satisfy z (i)>τ
zSample point be defined as " intersymbol saltus step (IST) sample ", remove corresponding y (i), obtain revised value sequence y
a(i);
4). calculate y
a(i) variance G, and set thresholding T, the sample of signal that satisfies G>T is judged as the MFSK signal;
5). with the y between the IST sample point
a(i) get average, calculate its distribution histogram;
6). calculate the order of modulation of MFSK signal by analyzing this histogrammic peak value number.
The problem of this method is:
● the accuracy of zero passage Frequency Estimation is very responsive to SNR, and in relevant report, the author has only provided the simulation result of CNR 〉=15dB;
● zero passage method is measured instantaneous frequency needs very high sample rate;
● under the low signal-to-noise ratio condition, calculate accurately relatively difficulty of order of modulation by the peak value number of analyzing the resulting frequency distribution block diagram of zero passage Frequency Estimation.
3. based on the method for signal discrete Fourier transform
People such as Zaihe Yu are at paper A practical classification algorithm for M-ary frequency shift keyingsignals (Military Communications Conference, 2004.MILCOM 2004.IEEE Volume 2,31 Oct.-3Nov.2004 Page (s): proposed a kind of method that obtains MFSK signal order of modulation by direct analysis MFSK power spectrum signal 1123-1128 Vol.2).Do not have the detection method of MFSK signal in this report, but suppose that signal to be analyzed is the MFSK signal.This method step is as follows:
1). the discrete Fourier transform (DFT) of signal calculated sample (DFT);
2). according to thresholding rule, local extremum rule, equidistant rule, perfect information rule, rule of classification, adaptive threshold rule and confidence level rule analysis sample of signal DFT peak value number, obtain order of modulation M.
The problem of this method is:
● this method supposes that the distribution of each code element on power equates, still, in the test environment of reality, because factors such as frequency selective fadings, each symbol power distributes and often do not equate;
● when calculating DFT, the author states the FFT length that use is long as far as possible, and this is unfavorable for the realization of real system, and under the signal bandwidth condition of unknown, estimates that a rational FFT length also is very difficult.
Summary of the invention
Main purpose of the present invention is identification and analyzes multi system frequency shift key control (MFSK) signal, comprises a kind of detection method of MFSK signal and a kind of analytical method of MFSK signal.
The invention is characterized in that described method realizes successively according to the following steps in computer:
Step (12) receives pending data, obtains real signal x (n);
Step (13) converts real signal x (n) to analytic signal s (n) according to the following steps:
Step (13.1) is calculated FFT after real signal x (n) zero padding, obtains the discrete Fourier transform (DFT) F of signal
x(n), the length N of zero padding
zObtain by following formula:
N
z=N
FFT-N
x, N
xBe the length of real signal x (n),
Expression is greater than the smallest positive integral of " ", N
FFTBe the length of carrying out the FFT computing;
Step (13.2) is got described F
x(n) first half, n=1 ..., N
FFT/ 2, obtain the discrete Fourier transform (DFT) F of signal s (n)
s(n);
Step (13.3) is F
s(n) first half, n=1 ..., N
FFT/ 4 and latter half n=N
FFT/ 4+1 ..., N
FFT/ 2 transposings;
The F that step (13.4) calculation procedure (2.3) obtains
s(n) IFFT obtains analytic signal s (n), s (n)=IFFT (F
s(n));
The analytic signal s's (n) that step (13.5) obtains step (2.4)
Part is clipped, and remainder is the analytic signal that subsequent analysis is used, and its effective length is
Step (14) is according to the predefine thresholding, and segmentation detects the amplitude of s (n), removes the weak signal section, and its steps in sequence is as follows:
The signal s (n) that step (14.1) obtains step (2.5)=x (n)+jy (n) is divided into isometric N
SegSection is with a sequence S
i, i=1,2 ..., N
SegExpression, N
Seg=5~20;
Described each the segment signal amplitude of step (14.2) calculation procedure (3.1) and m
i:
Step (14.3) is set thresholding τ
m:
Step (14.4) is with m
i<τ
mSignal segment be judged to be the weak signal section, remove all weak signal sections, remaining signal segment is connected into new analytic signal sequence s (n), length is N;
The normalization instantaneous frequency sequence f of the described analytic signal s of step (15) calculation procedure (3.4) (n)
n(n), realize according to the following steps:
Step (15.1) is calculated as follows analytic signal sequence s (n)=x (n)+jy (n), n=1, and 2 ..., the instantaneous phase sequence (n) of N:
The instantaneous phase sequence (n) that step (15.2) obtains step (4.1) does and gets the surplus instantaneous frequency sequence f (n) that obtains signal behind the calculus of differences:
f(n)=mod((n+1)-(n)+π,2π),n=1,2,……,N-1;
Step (15.3) is normalized to zero-mean, unit variance sequence f to this instantaneous frequency sequence f (n)
n(n):
Wherein,
The f that step (16) calculation procedure (4.3) obtains
nThe rate of change sequence Δ f of normalization instantaneous frequency (n)
n(n) and the average m of absolute value
Δ fn, wherein:
Δf
n(n)=|f
n(n+2)-f
n(n)|,n=1,2,……,N-3,
Step (17) is described m
Δ fnWith predefine thresholding τ
mRelatively:
If
Then described signal x (n) is judged to be the MFSK signal, otherwise, just non-MFSK signal;
Step (18) is carried out following steps successively and is reduced the shake of instantaneous frequency to utilize mean filter:
Step (18.1) utilizes mean filter that the instantaneous frequency value sequence is carried out smothing filtering, the instantaneous frequency sequence f that obtains upgrading
n(n), filter length l
f, by the possible maximum symbol rate f of analyzed signal
DmaxSetting sample frequency f with signal
sObtain by following formula:
The instantaneous frequency sequence of the renewal that step (18.2) obtains according to step (7.1) recomputates the normalization instantaneous frequency rate of change sequence Δ f of signal
n(n): Δ f
n(n)=| f
n(n+2)-f
n(n) |;
The variance τ of the normalization instantaneous frequency rate of change sequence that step (18.3) obtains according to following formula calculation procedure (7.2)
Δ f:
Step (18.4) search Δ f
n(n) sequence, if
F then
n(n+1) belong to symbol stable region Г
Si, i=1,2 ..., N
s(N
sNumber for the symbol stable region), if
F then
n(n+1) belong to the interval Г of symbol saltus step
Ti, i=1,2 ..., N
t(N
tNumber for symbol saltus step interval);
Step (18.5) changes the instantaneous frequency value of each symbol stable region in the step (7.4) average of all instantaneous frequency values of this interval into, gives up the instantaneous frequency value in all symbol saltus step intervals, constitutes new instantaneous frequency value sequence f
n(n);
Step (19) basis " entropy principle " is the histogram of the described instantaneous frequency sequential value of analytical procedure (7.5) according to the following steps:
Step (19.1) is according to normalization instantaneous frequency sequence f
n(n) set up distribution histogram
With this as f
nThe estimation of probability density function (n):
Wherein, N
FnBe f
n(n) length,
Be the length of distribution histogram,
Count (x, y} represent the to satisfy condition number of x value of y
N=10~100:
Step (19.3) in the hunting zone each, value carries out step:
Step (19.3.1) search is satisfied
F
i, continuous f
iBe classified as one group;
Step (19.3.2) is calculated the probability density function of each group and is estimated
Sum obtains each code element of MFSK signal
Step (19.3.3) is calculated should τ
pEntropy e:
Step (19.4) is got the τ of corresponding entropy maximum
p, with corresponding f
iThe number of group is as code element state number
initially estimate
Meter, every group minimum instantaneous frequency
With the maximum instantaneous frequency
Frequency range as each code element state
Step (20) basis " equiprobability distribution principle ", accurate according to the following steps estimating code element state number
Exact value:
Step (20.1) is set thresholding
Step (20.2) is removed
State;
Step (20.3) recomputates code element state number
The initial estimation and the frequency range of each code element state
Step (21) is according to accurate estimating code element state number such as " frequency difference principle "
With each state centre frequency, its step is as follows:
Step (21.3) is calculated the weight w of each frequency difference
Kl:
Weight equals to obtain the smaller value that the probability distribution of two code element states of this frequency difference is estimated, the probability distribution that deducts all the code element states between these two code element states is estimated sum, if difference less than 0, then weight is 0:
Step (21.4) obtains the estimation of frequency difference to each frequency difference weighted average:
Step (21.5) is set thresholding τ
Min, τ
Max:
Step (21.6) for
Adjacent states, merge code element state k and k+1;
Step (21.7) for
Adjacent states, between code element state k and k+1, insert
Individual code element state (
Expression is greater than the smallest positive integral of " ");
Step (22) subsequent treatment:
Wherein [] computing is represented to get from " " nearest integer.
The present invention has following advantage:
● do not rely on the hypothesized model that the MFSK signal is a constant envelope signal, allow the decline of signal frequency of occurrences selectivity, the scope of application is greater than other existing algorithms;
● caught the MFSK signal to have the substantive characteristics of instantaneous frequency saltus step, discerned and analyzed the MFSK signal as essential characteristic, had robustness preferably with instantaneous frequency;
● at the test of computer mould analog signal with at the test shows of actual signal, the algorithm that the present invention proposes has good recognition effect and stronger practicality.
The effect that the present invention reaches:
Multi system frequency shift key control (MFSK) is the important digital signal modulation mode of a class, and the present invention relates generally to the detection and the analysis of the signal of communication of such modulation system.Computer simulation experiment shows, to order of modulation is 2~32 MFSK signal, and under the condition of signal to noise ratio greater than 5dB, the detection accuracy of institute of the present invention extracting method has reached more than 99%, under the condition of signal to noise ratio greater than 7.5dB, analyze accuracy and reached more than 95%.To the test shows of the sample of signal of actual reception, the detection accuracy of institute of the present invention extracting method has reached more than 95%, analyzes accuracy and has reached more than 90%.
Description of drawings
MFSK input and analysis process figure that Fig. 1 proposes for the present invention.
Fig. 2 is the hardware elementary diagram that adopts when actual signal is tested among the embodiment.
Fig. 3 is the identification test result of emulation MFSK signal, and wherein dotted line is represented the m of fsk signal
Δ fnValue, solid line is represented the m of non-fsk signal
Δ fnValue.
Fig. 4 is the analytical test result of emulation MFSK signal, and every curve is represented the analysis accuracy of a kind of modulation system under different signal to noise ratio conditions.
Fig. 5 for the identification test result of actual reception MFSK signal wherein dotted line represent the m of fsk signal
Δ fnValue, solid line is represented the m of non-fsk signal
Δ fnValue.
Embodiment
MFSK input proposed by the invention and analytical method overall procedure may further comprise the steps as shown in Figure 1:
1) receives pending data, obtain real signal x (n);
2) real signal x (n) is become analytic signal s (n);
3) according to the predefine thresholding, segmentation detects s (n) amplitude, removes the weak signal section;
4) the normalization instantaneous frequency sequence f of calculating s (n)
n(n);
5) the rate of change sequence Δ f of calculating normalization instantaneous frequency
n(n) and the average m of absolute value
Δ fn
6) with m
Δ fnWith predefine thresholding τ
mCompare, if
Be the MFSK signal then, otherwise judge that this signal is not the MFSK signal this signal determining;
7) utilize mean filter to reduce the shake of instantaneous frequency;
8) analyze the instantaneous frequency distribution histogram according to " entropy principle ", obtain code element state number
Initial estimation;
9) according to " equiprobability distribution principle " accurate estimating code element state number
10) according to accurate estimating code element state number such as " frequency difference principle "
With each state centre frequency;
11) subsequent treatment.
Here provide at the test of Computer Simulation signal with at two embodiment of test of actual reception signal.
1. at the test of Computer Simulation signal
A) identification of MFSK signal
Under different signal to noise ratio conditions, generate two groups of signals and be used for test.One group is fsk signal, and one group is non-fsk signal:
● the fsk signal group comprises 2FSK, 4FSK, 8FSK, 16FSK, 32FSK signal, 256 code elements of signal length, and 100 sampled points of each code element, frequency difference 100Hz between the code element, sample rate is 4 times of maximum modulating frequency;
● non-fsk signal group comprises BPSK, QPSK, 8PSK, 16QAM, 64QAM, OQPSK, CW signal, and 256 code elements of signal length, sample frequency are 8 times of signal baud rate, adopts the raised cosine window shaping pulse of rolloff-factor 0.8.
The signal that receives in the actual application environment, signal to noise ratio are unknown, also are unsettled.For the situation of more realistic application,, added the white noise of 5dB to 30dB at random to two groups of test signals that generate.Obtain two groups of simulate signals at last, respectively comprise 500 segment signal samples, calculate the m of every segment signal respectively
Δ fnValue is estimated rational thresholding τ
m, obtain test result as accompanying drawing 3.Transverse axis is the sample of signal numbering in the accompanying drawing 3, and the longitudinal axis is m
Δ fnValue.As can be seen, the fsk signal group is obviously separated with non-fsk signal group, and reasonably the threshold value scope is τ
m=0.8~1.2.Work as τ
m=0.84 o'clock, the accuracy of identification reached 100%.
B) analysis of MFSK signal
Under 0~20dB signal to noise ratio condition, generation 2FSK, 4FSK, 8FSK, 16FSK, 32FSK signal are used for test.Signal length is 1024 code elements, and 100 sampled points of each code element, frequency difference 200Hz between the code element, signal sampling rate are 4 times of maximum modulating frequency.With the algorithm that the present invention proposes the signal that generates is analyzed, every kind of order of modulation/signal to noise ratio combination is repeated 100 tests, obtained the test result as accompanying drawing 4.Accompanying drawing 4 is that this algorithm of utilization is analyzed the curve of the analysis accuracy that obtains to the test signal of different modulating exponent number under different signal to noise ratio conditions, and transverse axis is the signal to noise ratio of test signal, and the longitudinal axis is to analyze accuracy, and every curve is represented a kind of modulation system.As can be seen, when SNR 〉=7.5dB, the analysis accuracy of 2FSK, 4FSK, 8FSK, 16FSK, 32FSK modulation signal has all been reached more than 95%.
2. at the test of actual reception signal (test is seen accompanying drawing 2 with hardware elementary diagram)
A) identification of MFSK signal
With seemingly, test data is divided into two groups at the test class of simulate signal:
● the MFSK signal, 2FSK (58), 4FSK (16), 6FSK (3), 8FSK (17), 12FSK (3), 13FSK (2), 16FSK (1), 20FSK (2), 32FSK (5) signal (being the sample number of this order of modulation signal in the bracket) have been comprised, the sample rate scope is 6000~8000Hz, chip rate scope 13Bd~1000Bd.
● non-MFSK signal, comprised BPSK (8), QPSK (2), 8PSK (9) signal (being the sample number of this modulation system signal in the bracket), the sample rate scope is 6000~8000Hz, the chip rate scope is 30Bd~2400Bd.
From two groups of samples, the random extraction time span is the signal segment of 1s (6000~8000 sampled point), calculates its m
Δ fnValue is estimated rational thresholding τ
m, obtain test result as accompanying drawing 5.In the accompanying drawing 5, transverse axis is the sample signal segment number, and the longitudinal axis is m
Δ fnValue.As can be seen from the figure, two groups of samples are well separated basically, and reasonably the threshold value scope is τ
m=0.8~1.2, conform to simulation result.Get τ
mIn the time of=1.17, the accuracy of separation has reached 95.6%.
B) analysis of MFSK signal
We utilize the actual reception signal that the parser of the MFSK signal of the present invention's proposition is tested, because sufficiently long 16FSK signal and 32FSK signal are difficult to obtain, test has included only 2FSK, 4FSK, 8FSK signal with set of signals.In the test every segment signal sample is divided into the data segment that length is 1s (6000~8000 sample points), analyzes with the MFSK signal analysis algorithm that the present invention proposes, the test result that obtains is as follows:
Signal type | The data hop count | Correct number of results | Accuracy (%) |
2FSK | 1078 | 988 | 91.7 |
4FSK | 127 | 115 | 90.6 |
8FSK | 169 | 155 | 91.8 |
Amount to | 1374 | 1258 | 91.56 |
Test result from last table as can be seen, this algorithm has reached than higher accuracy the analysis of actual signal.
Claims (1)
1. the detection of a multi system frequency shift key control signal and analytical method is characterized in that, described method realizes in computer successively according to the following steps:
Step (1) receives pending data, obtains real signal x (n);
Step (2) converts real signal x (n) to analytic signal s (n) according to the following steps:
Step (2.1) is calculated FFT after real signal x (n) zero padding, obtains the discrete Fourier transform (DFT) F of signal
x(n), the length N of zero padding
zObtain by following formula:
N
z=N
FFT-N
x, N
xBe the length of real signal x (n),
Expression is greater than the smallest positive integral of " ", N
FFTBe the length of carrying out the FFT computing;
Step (2.2) is got described F
x(n) first half, n=1 ..., N
FFT/ 2, obtain the discrete Fourier transform (DFT) F of signal s (n)
s(n);
Step (2.3) is F
s(n) first half, n=1 ..., N
FFT/ 4 and latter half n=N
FFT/ 4+1 ..., N
FFT/ 2 transposings;
The F that step (2.4) calculation procedure (2.3) obtains
s(n) IFFT obtains analytic signal s (n), s (n)=IFFT (F
s(n));
The analytic signal s's (n) that step (2.5) obtains step (2.4)
Part is clipped, and remainder is the analytic signal that subsequent analysis is used, and its effective length is
Step (3) is according to the predefine thresholding, and segmentation detects the amplitude of s (n), removes the weak signal section, and its steps in sequence is as follows:
The signal s (n) that step (3.1) obtains step (2.5)=x (n)+jy (n) is divided into isometric N
SegSection is with a sequence S
i, i=1,2 ..., N
SegExpression, N
Seg=5~20;
Described each the segment signal amplitude of step (3.2) calculation procedure (3.1) and m
i:
Step (3.3) is set thresholding τ
m:
Step (3.4) is with m
i<τ
mSignal segment be judged to be the weak signal section, remove all weak signal sections, remaining signal segment is connected into new analytic signal sequence s (n), length is N;
The normalization instantaneous frequency sequence f of the described analytic signal s of step (4) calculation procedure (3.4) (n)
n(n), realize according to the following steps:
Step (4.1) is calculated as follows analytic signal sequence s (n)=x (n)+jy (n), n=1, and 2 ..., the instantaneous phase sequence (n) of N:
The instantaneous phase sequence (n) that step (4.2) obtains step (4.1) does and gets the surplus instantaneous frequency sequence f (n) that obtains signal behind the calculus of differences:
f(n)=mod((n+1)-(n)+π,2π),n=1,2,……,N-1;
Step (4.3) is normalized to zero-mean, unit variance sequence f to this instantaneous frequency sequence f (n)
n(n):
Wherein,
The f that step (5) calculation procedure (4.3) obtains
nThe rate of change sequence Δ f of normalization instantaneous frequency (n)
n(n) and the average m of absolute value
Δ fn, wherein:
Step (6) is described m
Δ fnWith predefine thresholding τ
mRelatively:
If
Then described signal x (n) is judged to be the MFSK signal, otherwise, just non-MFSK signal;
Step (7) is carried out following steps successively and is reduced the shake of instantaneous frequency to utilize mean filter:
Step (7.1) utilizes mean filter that the instantaneous frequency value sequence is carried out smothing filtering, the instantaneous frequency sequence f that obtains upgrading
n(n), filter length l
fBy the possible maximum symbol rate f of analyzed signal
DmaxSetting sample frequency f with signal
sObtain by following formula:
The instantaneous frequency sequence of the renewal that step (7.2) obtains according to step (7.1) recomputates the normalization instantaneous frequency rate of change sequence Δ f of signal
n(n): Δ f
n(n)=| f
n(n+2)-f
n(n) |;
The variance τ of the normalization instantaneous frequency rate of change sequence that step (7.3) obtains according to following formula calculation procedure (7.2)
Δ f:
Step (7.4) search Δ f
n(n) sequence, if
F then
n(n+1) belong to symbol stable region Γ
Si, i=1,2 ..., N
s(N
sNumber for the symbol stable region), if
F then
n(n+1) belong to the interval Γ of symbol saltus step
Ti, i=1,2 ..., N
t(N
tNumber for symbol saltus step interval);
Step (7.5) changes the instantaneous frequency value of each symbol stable region in the step (7.4) average of all instantaneous frequency values of this interval into, gives up the instantaneous frequency value in all symbol saltus step intervals, constitutes new instantaneous frequency value sequence f
n(n);
Step (8) basis " entropy principle " is the histogram of the described instantaneous frequency sequential value of analytical procedure (7.5) according to the following steps:
Step (8.1) is according to normalization instantaneous frequency sequence f
n(n) set up distribution histogram
With this as f
nThe estimation of probability density function (n):
Wherein, N
FnBe f
n(n) length,
Be the length of distribution histogram,
Step (8.2) is determined thresholding τ
pThe hunting zone be
Step-size in search
N=10~100;
Step (8.3) is to each τ in the hunting zone
pValue is carried out step:
Step (8.3.1) search is satisfied
F
i, continuous f
iBe classified as one group;
Step (8.3.2) is calculated the probability density function of each group and is estimated
Sum obtains the estimated value of the distribution probability of each code element of MFSK signal
Step (8.3.3) is calculated should τ
pEntropy e:
Step (8.4) is got the τ of corresponding entropy maximum
p, with corresponding f
iThe number of group is as code element state number
Initial estimation, every group minimum instantaneous frequency
With the maximum instantaneous frequency
Frequency range as each code element state
Step (9) basis " equiprobability distribution principle ", accurate according to the following steps estimating code element state number
Exact value:
Step (9.1) is set thresholding
Step (9.2) is removed
State;
Step (9.3) recomputates code element state number
The initial estimation and the frequency range of each code element state
Step (10) is according to accurate estimating code element state number such as " frequency difference principle "
With each state centre frequency, its step is as follows:
Step (10.2) is calculated the frequency difference between each state
Step (10.3) is calculated the weight w of each frequency difference
Kl:
Weight equals to obtain the smaller value that the probability distribution of two code element states of this frequency difference is estimated, the probability distribution that deducts all the code element states between these two code element states is estimated sum, if difference less than 0, then weight is 0:
Step (10.4) obtains the estimation of frequency difference to each frequency difference weighted average:
Step (10.5) is set thresholding τ
Min, τ
Max:
Step (10.6) for
Adjacent states, merge code element state k and k+1;
Step (10.7) for
Adjacent states, between code element state k and k+1, insert
Individual code element state (
Expression is greater than the smallest positive integral of " ");
Step (11) subsequent treatment:
Be calculated as follows the order of modulation that obtains the MFSK signal
Wherein [] computing is represented to get from " " nearest integer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101140851A CN100521670C (en) | 2006-10-27 | 2006-10-27 | Detecting and analyzing method for multi system frequency shift key control signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006101140851A CN100521670C (en) | 2006-10-27 | 2006-10-27 | Detecting and analyzing method for multi system frequency shift key control signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1946069A true CN1946069A (en) | 2007-04-11 |
CN100521670C CN100521670C (en) | 2009-07-29 |
Family
ID=38045286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006101140851A Expired - Fee Related CN100521670C (en) | 2006-10-27 | 2006-10-27 | Detecting and analyzing method for multi system frequency shift key control signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100521670C (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102113283A (en) * | 2008-06-03 | 2011-06-29 | 汤姆森特许公司 | Apparatus and method for determination of signal format |
CN102763339A (en) * | 2010-02-16 | 2012-10-31 | 飞思卡尔半导体公司 | Detector and method for detecting an oscillatory signal among noise |
CN101605001B (en) * | 2009-07-13 | 2012-11-21 | 中国船舶重工集团公司第七一五研究所 | Doppler measurement and correction method for MFSK underwater acoustic communication |
CN109543643A (en) * | 2018-11-30 | 2019-03-29 | 电子科技大学 | Carrier signal detection method based on one-dimensional full convolutional neural networks |
CN109581429A (en) * | 2018-12-18 | 2019-04-05 | 中国电子科技集团公司第五十四研究所 | A kind of GNSS signal acquisition performance analysis method |
CN111756038A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
CN111800360A (en) * | 2020-07-03 | 2020-10-20 | 国仪量子(无锡)技术有限公司 | FSK software decoding method based on frequency identification |
CN111800359A (en) * | 2020-09-07 | 2020-10-20 | 中国人民解放军国防科技大学 | Method, device, equipment and medium for identifying communication signal modulation mode |
CN112737992A (en) * | 2020-09-23 | 2021-04-30 | 青岛科技大学 | Underwater sound signal modulation mode self-adaptive in-class identification method |
-
2006
- 2006-10-27 CN CNB2006101140851A patent/CN100521670C/en not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102113283B (en) * | 2008-06-03 | 2014-09-03 | 汤姆森特许公司 | Apparatus and method for determination of signal format |
CN102113283A (en) * | 2008-06-03 | 2011-06-29 | 汤姆森特许公司 | Apparatus and method for determination of signal format |
CN101605001B (en) * | 2009-07-13 | 2012-11-21 | 中国船舶重工集团公司第七一五研究所 | Doppler measurement and correction method for MFSK underwater acoustic communication |
CN102763339A (en) * | 2010-02-16 | 2012-10-31 | 飞思卡尔半导体公司 | Detector and method for detecting an oscillatory signal among noise |
CN102763339B (en) * | 2010-02-16 | 2015-04-29 | 飞思卡尔半导体公司 | Detector and method for detecting an oscillatory signal among noise |
US9094084B2 (en) | 2010-02-16 | 2015-07-28 | Freescale Semiconductor, Inc. | Detector and method for detecting an oscillatory signal among noise |
CN109543643B (en) * | 2018-11-30 | 2022-07-01 | 电子科技大学 | Carrier signal detection method based on one-dimensional full convolution neural network |
CN109543643A (en) * | 2018-11-30 | 2019-03-29 | 电子科技大学 | Carrier signal detection method based on one-dimensional full convolutional neural networks |
CN109581429A (en) * | 2018-12-18 | 2019-04-05 | 中国电子科技集团公司第五十四研究所 | A kind of GNSS signal acquisition performance analysis method |
CN111800360A (en) * | 2020-07-03 | 2020-10-20 | 国仪量子(无锡)技术有限公司 | FSK software decoding method based on frequency identification |
CN111800360B (en) * | 2020-07-03 | 2022-07-26 | 国仪量子(无锡)技术有限公司 | FSK software decoding method based on frequency identification |
CN111756038A (en) * | 2020-07-09 | 2020-10-09 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
CN111756038B (en) * | 2020-07-09 | 2021-08-13 | 西安交通大学 | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics |
CN111800359B (en) * | 2020-09-07 | 2020-12-04 | 中国人民解放军国防科技大学 | Method, device, equipment and medium for identifying communication signal modulation mode |
CN111800359A (en) * | 2020-09-07 | 2020-10-20 | 中国人民解放军国防科技大学 | Method, device, equipment and medium for identifying communication signal modulation mode |
CN112737992A (en) * | 2020-09-23 | 2021-04-30 | 青岛科技大学 | Underwater sound signal modulation mode self-adaptive in-class identification method |
CN112737992B (en) * | 2020-09-23 | 2023-02-03 | 青岛科技大学 | Underwater sound signal modulation mode self-adaptive in-class identification method |
Also Published As
Publication number | Publication date |
---|---|
CN100521670C (en) | 2009-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1946069A (en) | Detecting and analyzing method for multi system frequency shift key control signal | |
CN1746973A (en) | Distributed speech recognition system and method | |
CN1083183C (en) | Method and apparatus for reducing noise in speech signal | |
CN1914839A (en) | Reception device | |
CN1223109C (en) | Enhancement of near-end voice signals in an echo suppression system | |
CN104468001B (en) | Signal identification method and system based on radio signal frequency spectrum feature template | |
CN101079266A (en) | Method for realizing background noise suppressing based on multiple statistics model and minimum mean square error | |
CN101064571A (en) | Apparatus for enhancing channel evaluation in OFDM receiver and its method | |
CN1754308A (en) | Coherent am demodulator using a weighted lsb/usb sum for interference mitigation | |
CN101047995A (en) | Channel switchover method and adaptive method of interference detection threshold | |
CN101080911A (en) | Calibrating amplitude and phase imbalance and DC offset of an analog i /q modulator in a high-frequency transmitter | |
CN1794757A (en) | Telephone and method for processing audio single in the telephone | |
CN1624767A (en) | Noise reduction apparatus and noise reducing method | |
CN100340122C (en) | Channel estimation method for a mobile communication system | |
CN1764094A (en) | Decreasing computational complexity of TD-SCDMA measurement process | |
CN101075845A (en) | Method and apparatus for realizing down synchronization in first search of area | |
CN1623186A (en) | Voice activity detector and validator for noisy environments | |
CN100338969C (en) | Microphone-loudspeaker device | |
CN1174592C (en) | Data communication apparatus and method based on orthogonal FDMA | |
CN101053192A (en) | Method and system for determining a frequency offset | |
CN109076038B (en) | Method for estimating parameters of a signal contained in a frequency band | |
CN1611963A (en) | Signal analysis | |
CN1905383A (en) | Shared frequency cell channel estimating apparatus and method | |
CN1949754A (en) | Method for estimating OFDM integer frequency shift based on virtual subcarrier and frequency domain differential sequence | |
CN1736039A (en) | Echo suppression device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090729 Termination date: 20091127 |