CN103051367A - Clustering-based blind source separation method for synchronous orthogonal frequency hopping signals - Google Patents
Clustering-based blind source separation method for synchronous orthogonal frequency hopping signals Download PDFInfo
- Publication number
- CN103051367A CN103051367A CN2012105250644A CN201210525064A CN103051367A CN 103051367 A CN103051367 A CN 103051367A CN 2012105250644 A CN2012105250644 A CN 2012105250644A CN 201210525064 A CN201210525064 A CN 201210525064A CN 103051367 A CN103051367 A CN 103051367A
- Authority
- CN
- China
- Prior art keywords
- frequency
- time
- signal
- frequency hopping
- domain
- 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
- Mobile Radio Communication Systems (AREA)
Abstract
The invention discloses a clustering-based blind source separation method for synchronous orthogonal frequency hopping signals. The method comprises the following steps of: acquiring M sampled paths of discrete time-domain mixed signals; obtaining M time-frequency domain matrixes of the mixed signals; preprocessing the time-frequency domain matrixes of the frequency hopping mixed signals; estimating frequency hopping moments, normalized mixed matrix column vectors and frequency hopping frequency; estimating time-frequency domain frequency hopping source signals by utilizing the estimated normalized mixed matrix column vectors; splicing the time-frequency domain frequency hopping source signals between different frequency hopping points; and recovering time-domain source signals according to time-frequency domain estimate values of the source signals. According to the method, the frequency hopping source signals are estimated only according to the received mixed signals of a plurality of frequency hopping signals under the condition of unknown channel information, and the frequency hopping signals can be subjected to blind estimation under the condition that the number of receiving antennae is smaller than that of the source signals; short-time Fourier transform is utilized, so that the method is low in computation amount; and the frequency hopping signals are subjected to blind separation, and meanwhile, a part of parameters can also be estimated, so that the method is high in practicability.
Description
Technical field
The invention belongs to communication and signal processing technology field, relate in particular to a kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster.
Background technology
The separation of blind source refers to is not knowing under the condition of any channel information, only estimates the process of source signal according to the mixed signal that observes.When observation signal number during greater than the source signal number, be called the blind separation of overdetermination; When the observation signal number equals the source signal number, be called suitable fixed blind separation; When observation signal number during less than the source signal number, become and owe fixed blind separation.
Blind source separate technology obtains using more and more widely in the signal of communication process field, the at present research of blind source separate technology more is more than or equal to the overdetermination of source signal number or suitable fixed blind separation for the observation signal number, EASI (Equivariant adaptive source separation such as classics, Deng the variation self adaptation) algorithm, FastICA (Fast Independent Component Analysis, Fast Independent Component Analysis) algorithm etc.Owing to satisfy statistical independence between the different Frequency Hopping Signals, and satisfy certain correlation between the different frequency hopping points of identical Frequency Hopping Signal, the thought that therefore has document to propose to utilize blind source to separate is carried out blind estimation to a plurality of Frequency Hopping Signals that receive.
The blind separation algorithm of Frequency Hopping Signal that provides in the existing document requires the observation signal number more than or equal to the source signal number.But in the practical communication process, because the restriction of reception antenna number, when frequency hopping synthesizer signal number was more, the observation signal number was often less than the source signal number, therefore existing algorithm can't in the less situation of reception antenna number, carry out blind estimation to a plurality of Frequency Hopping Signals.
Summary of the invention
The invention provides a kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster, be intended to solve the blind separation algorithm requirement of existing Frequency Hopping Signal observation signal number more than or equal to the source signal number, but in the practical communication process, because the restriction of reception antenna number, when frequency hopping synthesizer signal number is more, the observation signal number is often less than the source signal number, and therefore existing algorithm can't be in the less situation of reception antenna number, the problem of a plurality of Frequency Hopping Signals being carried out blind estimation.
The object of the present invention is to provide a kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster, the method may further comprise the steps:
Step 3 is to the frequency-hopping mixing signal time-frequency domain matrix that obtains in the step 2
Carry out preliminary treatment;
Step 4 utilizes clustering algorithm to estimate jumping moment and corresponding normalized hybrid matrix column vector, the frequency hopping frequency of each jumping of each jumping;
Step 6 is spliced the time-frequency domain frequency hopping synthesizer signal between the different frequency hopping points;
Step 7, root signal time-frequency domain estimated value is recovered the time domain source signal.
Further, in step 2, (p, q) expression time-frequency index, concrete time-frequency value is
Here N
FftThe length of expression FFT conversion, P represents the windowing number of times, C is integer, the sampling number at expression Short Time Fourier Transform windowing interval, general C<N
Fft, and K
c=N
Fft/ C is integer, and what that is to say employing is the Short Time Fourier Transform of overlapping windowing.
Further, in step 3, to frequency-hopping mixing signal time-frequency domain matrix
Carry out preliminary treatment, specifically comprise following two steps:
The first step is right
Go low-yield preliminary treatment, namely at each sampling instant p, will
Amplitude sets to 0 less than the value of thresholding ε, obtains
The setting of thresholding ε can be determined according to the average energy that receives signal;
Second step is found out the constantly time-frequency domain data of (p=0,1,2, L P-1) non-zero of p, uses
Expression, wherein
Constantly time-frequency response of expression p
Corresponding frequency indices when non-zero to these non-zero normalization preliminary treatment, obtains pretreated vectorial b (p, q)=[b
1(p, q), b
2(p, q), L, b
M(p, q)]
T, wherein
Further, in step 4, when utilizing clustering algorithm to estimate that the jumping moment of each jumping and each are jumped corresponding normalized hybrid matrix column vector, frequency hopping frequency, may further comprise the steps:
The first step, p (p=0,1,2 ... P-1) constantly, right
The frequency values of expression carries out cluster, the cluster centre number that obtains
The carrier frequency number that expression p exists constantly,
Individual cluster centre then represents the size of carrier frequency, uses respectively
Expression;
Second step, to each sampling instant p (p=0,1,2 ... P-1), utilize clustering algorithm pair
Carry out cluster, can obtain equally
Individual cluster centre is used
Expression;
In the 4th step, find out
The moment, use p
hExpression is to the p of the continuous value of each section
hAsk intermediate value, use
Represent the l section p that links to each other
hIntermediate value, then
Represent the estimation constantly of l frequency hopping;
The 5th step is according to what estimate in the second step to obtain
P ≠ p
hAnd the 4th estimate that the frequency hopping obtain constantly estimates each and jumps corresponding in the step
Individual hybrid matrix column vector
Concrete formula is:
Here
It is corresponding to represent that l jumps
Individual hybrid matrix column vector estimated value;
The 6th step, estimate the carrier frequency frequency that each jumping is corresponding, use
It is corresponding to represent that l jumps
Individual frequency estimation, computing formula is as follows:
Further, in step 5, estimate time-frequency domain frequency hopping synthesizer signal according to the normalization hybrid matrix column vector of estimating in the step 4 to obtain, concrete steps are as follows:
The first step judges that to all sampling instant index p which this moment index belongs to and jump, and concrete grammar is: if
Represent that then p belongs to the l jumping constantly; If
Represent that then p belongs to the 1st jumping constantly;
Second step is to l (l=1,2, all moment p that L) jump
l, estimate that this jumps the time-frequency domain data of each frequency hopping synthesizer signal, computing formula is as follows:
Further, in step 6, the time-frequency domain frequency hopping synthesizer signal between the different frequency hopping points is spliced,
Concrete steps are:
The first step, it is corresponding to estimate that l jumps
Individual incident angle is used
Represent that l jumps n the incident angle that source signal is corresponding,
Computing formula as follows:
Second step, judge the 1st (1=2,3 ...) jump the source signal and first of estimating and jump corresponding relation between the source signal of estimation, judgment formula is as follows:
M wherein
n (l)Represent that l jumps the m that estimates
n (l)Individual signal and first is jumped n the signal of estimating and is belonged to same source signal;
The 3rd step, with different frequency hopping point estimation to the signal that belongs to same source signal be stitched together, estimate as final time-frequency domain source signal, use Y
nN source signal of (p, q) expression the time time-frequency domain estimated value on frequency (p, q), p=0,1,2 ...., P, q=0,1,2 ..., N
Fft-1, namely
Further, in step 7, when root signal time-frequency domain estimated value was recovered the time domain source signal, concrete steps were as follows:
The first step, to each sampling instant p (p=0,1,2 ...) and frequency domain data Y
n(p, q), q=0,1,2, L, N
Fft-1 is N
FftThe IFFT conversion of point obtains time domain source signal corresponding to p sampling instant, uses y
n(p, q
t) (q
t=0,1,2, L, N
Fft-1) expression;
Second step is to above-mentioned all time domain source signal y that constantly obtain
n(p, q
t) merge processing, obtain final time domain source signal and estimate, concrete formula is as follows:
Here K
c=N
Fft/ C, the sampling number at C windowing interval.
Synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster provided by the invention, do not knowing under the condition of any channel information, only according to the mixed signal of a plurality of Frequency Hopping Signals of receiving, estimate the frequency hopping synthesizer signal, can be under the condition of reception antenna number less than the source signal number, a plurality of Frequency Hopping Signals are carried out blind estimation, only utilized Short Time Fourier Transform, amount of calculation is little, realize that easily the method can also be estimated partial parameters when Frequency Hopping Signal is carried out blind separation, practical, have stronger propagation and employment and be worth.
Description of drawings
Fig. 1 is the realization flow figure based on the synchronized orthogonal Frequency Hopping Signal blind source separation method of cluster that the embodiment of the invention provides;
Fig. 2 is the 1st source signal and the estimated signal time domain waveform figure thereof that the embodiment of the invention provides;
Fig. 3 is the 2nd source signal and the estimated signal time domain waveform figure thereof that the embodiment of the invention provides;
Fig. 4 is the 3rd source signal and the estimated signal time domain waveform figure thereof that the embodiment of the invention provides;
Fig. 5 is the 4th source signal and the estimated signal time domain waveform figure thereof that the embodiment of the invention provides;
Fig. 6 is the embodiment of the invention original Frequency Hopping Signal that provides and Frequency Hopping Signal coefficient correlation and the Between Signal To Noise Ratio curve chart estimated.
Embodiment
In order to make purpose of the present invention, technical scheme and advantage clearer, below in conjunction with drawings and Examples, the present invention is further described in detail.Should be appreciated that specific embodiment described herein only in order to explaining the present invention, and be not used in and limit invention.
The realization flow based on the synchronized orthogonal Frequency Hopping Signal blind source separation method of cluster that Fig. 1 shows that the embodiment of the invention provides.
The method may further comprise the steps:
Step S101 utilizes and contains the array antenna received of M array element from the Frequency Hopping Signal of a plurality of synchronized orthogonal frequency hopping radio sets, respectively each road frequency-hopping mixing signal is sampled, and establishing sample frequency is f
s, the sampling interval is T
s=1/f
s, the discrete data of field frequency-hopping mixing signal is used during m road after the sampling
Expression then has
M=1,2, L, M;
Step S102 is to the discrete time domain mixed signal in M road
Carry out overlapping windowing Short Time Fourier Transform, obtain the time-frequency domain matrix of M frequency-hopping mixing signal
P=0,1, L P-1, q=0,1, L N
Fft-1;
Step S103 is to the time-frequency domain matrix of frequency-hopping mixing signal
Carry out preliminary treatment, obtain vectorial a (p, q);
Step S104 utilizes clustering algorithm to estimate jumping moment and corresponding normalized hybrid matrix column vector, the frequency hopping frequency of each jumping of each jumping;
Step S105, the normalization hybrid matrix column vector that estimation obtains according to step S104 is estimated the time-frequency domain source signal;
Step S106 splices the time-frequency domain frequency hopping synthesizer signal between the different frequency hopping points;
Step S107 according to frequency hopping synthesizer signal time-frequency domain estimated value, recovers time domain frequency hopping synthesizer signal.
In embodiments of the present invention, in step S102, (p, q) expression time-frequency index, concrete time-frequency value is
Here N
FftThe length of expression FFT conversion, P represents the windowing number of times, C is integer, the sampling number at expression Short Time Fourier Transform windowing interval, general C<N
Fft, and K
c=N
Fft/ C is integer, and what that is to say employing is the Short Time Fourier Transform of overlapping windowing.
In embodiments of the present invention, in step S103, to the time-frequency domain data
Carry out preliminary treatment, when obtaining vectorial a (p, q), specifically comprise following two steps:
The first step is right
Go low-yield preliminary treatment, namely at each sampling instant p, will
Amplitude sets to 0 less than the value of thresholding ε, obtains
The setting of thresholding ε can be determined according to the average energy that receives signal;
Second step is found out the constantly time-frequency domain data of (p=0,1,2, L P-1) non-zero of p, uses
Expression, wherein
Constantly time-frequency response of expression p
Corresponding frequency indices when non-zero to these non-zero normalization preliminary treatment, obtains pretreated vectorial b (p, q)=[b
1(p, q), b
2(p, q), L, b
M(p, q)]
T, wherein
In embodiments of the present invention, in step S104, when utilizing clustering algorithm to estimate that the jumping moment of each jumping and each are jumped corresponding normalized hybrid matrix column vector, frequency hopping frequency, may further comprise the steps:
The first step, p (p=0,1,2 ... P-1) constantly, right
The frequency values of expression carries out cluster, the cluster centre number that obtains
The carrier frequency number that expression p exists constantly,
Individual cluster centre then represents the size of carrier frequency, uses respectively
Expression;
Second step to each sampling instant p, utilizes clustering algorithm pair
Carry out cluster, can obtain equally
Individual cluster centre is used
Expression;
In the 4th step, find out
The moment, use p
hExpression is to the p of the continuous value of each section
hAsk intermediate value, use
Represent the l section p that links to each other
hIntermediate value, then
Represent the estimation constantly of the 1st frequency hopping;
The 5th step is according to what estimate in the second step to obtain
P ≠ p
hAnd the 4th estimate that the frequency hopping obtain constantly estimates each and jumps corresponding in the step
Individual hybrid matrix column vector
Concrete formula is:
Here
It is corresponding to represent that l jumps
Individual hybrid matrix column vector estimated value;
The 6th step, estimate the carrier frequency frequency that each jumping is corresponding, use
It is corresponding to represent that l jumps
Individual frequency estimation, computing formula is as follows:
In embodiments of the present invention, in step S105, estimate the time-frequency domain source signal according to the normalization hybrid matrix column vector of estimating among the step S104 to obtain, concrete steps are as follows:
The first step judges that to all sampling instant index p which this moment index belongs to and jump, and concrete grammar is: if
Represent that then p belongs to the l jumping constantly; If
Represent that then p belongs to the 1st jumping constantly;
Second step is to l (l=1,2, all moment p that L) jump
l, estimate that this jumps the time-frequency domain data of each source signal, computing formula is as follows:
In embodiments of the present invention, in step S106, the frequency hopping synthesizer signal between the different frequency hopping points is spliced, concrete steps are:
The first step, it is corresponding to estimate that l jumps
Individual incident angle is used
Represent that l jumps n the incident angle that source signal is corresponding,
Computing formula as follows:
Second step, judge l (l=2,3 ...) and source signal and first that jump to estimate jumps the corresponding relation between the source signal of estimation, judgment formula is as follows:
M wherein
n (l)Represent that l jumps the m that estimates
n (l)Individual signal and first is jumped n the signal of estimating and is belonged to same source signal;
The 3rd the step, with different frequency hopping point estimation to the signal that belongs to same source signal be stitched together, as final time-frequency domain frequency hopping synthesizer Signal estimation, use Y
nN source signal of (p, q) expression the time time-frequency domain estimated value on frequency (p, q), p=0,1,2 ...., P, q=0,1,2 ..., N
Fft-1, namely
In embodiments of the present invention, in step S107, when recovering time domain frequency hopping synthesizer signal according to frequency hopping synthesizer signal time-frequency domain estimated value, concrete steps are as follows:
The first step, to each sampling instant p (p=0,1,2 ...) and frequency domain data Y
n(p, q), q=0,1,2, L, N
Fft-1 is N
FftThe IFFT conversion of point obtains time domain source signal corresponding to p sampling instant, uses y
n(p, q
t) (q
t=0,1,2, L, N
Fft-1) expression;
Second step is to above-mentioned all time domain source signal y that constantly obtain
n(p, q
t) merge processing, obtain final time domain frequency hopping synthesizer Signal estimation, concrete formula is as follows:
Here K
c=N
Fft/ C, the sampling number at C windowing interval.
Below in conjunction with drawings and the specific embodiments application principle of the present invention is further described.
The object of the invention is to owe surely blind separation problem for Frequency Hopping Signal, propose a kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster, the method comprises:
(1) utilization contains the aerial array reception of M array element from the Frequency Hopping Signal of a plurality of synchronized orthogonal frequency hopping radio sets, each road frequency-hopping mixing signal is sampled the discrete time-domain frequency-hopping mixing signal after obtaining sampling
M=1,2, L, M.
(2) to M road discrete time-domain Frequency Hopping Signal
Carry out overlapping windowing Short Time Fourier Transform (employing rectangular window function), obtain the time-frequency domain matrix of M mixed signal
P=0,1, L P-1, q=0,1, L N
Fft-1.Wherein (p, q) represents the time-frequency index, and concrete time-frequency value is
Here N
FftThe length of expression FFT conversion, P represents the windowing number of times, C is integer, the sampling number at expression Short Time Fourier Transform windowing interval, general C<N
Fft, and K
c=N
Fft/ C is integer, and what that is to say employing is the Short Time Fourier Transform of overlapping windowing.
(3) to mixed signal time-frequency domain matrix
Carry out preliminary treatment, obtain vectorial a (p, q), specifically comprise following two steps:
3.1) right
Go low-yield preliminary treatment, namely at each sampling instant p, will
Amplitude sets to 0 less than the value of thresholding ε, obtains
The setting of thresholding can be determined according to the average energy that receives signal.
3.2) find out the constantly time-frequency domain data of (p=0,1,2, L P-1) non-zero of p, use
Expression, wherein
Constantly time-frequency response of expression p
Corresponding frequency indices when non-zero to these non-zero normalization preliminary treatment, obtains pretreated vectorial b (p, q)=[b
1(p, q), b
2(p, q), L, b
M(p, q)]
T, wherein
(4) utilize clustering algorithm to estimate jumping moment and corresponding normalized hybrid matrix column vector, the frequency hopping frequency of each jumping of each jumping.May further comprise the steps:
4-1) p (p=0,1,2 ... P-1) constantly, right
The frequency values of expression carries out cluster, the cluster centre number that obtains
The carrier frequency number that expression p exists constantly,
Individual cluster centre then represents the size of carrier frequency, uses respectively
Expression;
4-2) to each sampling instant p, utilize clustering algorithm pair
Carry out cluster, can obtain equally
Individual cluster centre is used
Expression.
4-4) find out
The moment, use p
hExpression is to the p of the continuous value of each section
hAsk intermediate value, use
Represent the l section p that links to each other
hIntermediate value, then
Represent the estimation constantly of l frequency hopping.
4-5) obtain according to estimation in the step (4-2)
P ≠ p
hAnd estimate in the step (4-3) that the frequency hopping obtain constantly estimates each and jumps corresponding
Individual hybrid matrix column vector
Concrete formula is as follows:
Here
It is corresponding to represent that l jumps
Individual hybrid matrix column vector estimated value.
4-6) estimate the carrier frequency frequency that each jumping is corresponding, use
It is corresponding to represent that l jumps
Individual frequency estimation, computing formula is as follows:
(5) estimate the time-frequency domain source signal according to the normalization hybrid matrix column vector of estimating in the step (4) to obtain, concrete steps are as follows:
5-1) all sampling instant index p are judged which this moment index belongs to and jump, concrete grammar is as follows:
If
Represent that then p belongs to the l jumping constantly;
5-2) to l (l=1,2, all of L) jumping are pl constantly, utilize formula (5) to estimate that this jumps the time-frequency domain data of each source signal, computing formula is as follows:
(6) splicing of source signal between the different frequency hopping points, concrete steps are as follows:
6-1) estimation l jumps correspondence
Individual incident angle is used
Represent that l jumps n the incident angle that source signal is corresponding,
Computing formula as follows:
6-2) utilize formula (7) judge l (l=2,3 ...) and jump the source signal and first of estimating and jump corresponding relation between the source signal of estimation, as follows
M wherein
n (l)Represent that l jumps the m that estimates
n (l)Individual signal and first is jumped n the signal of estimating and is belonged to same source signal.
6-3) with different frequency hopping point estimation to the signal that belongs to same source signal be stitched together, estimate as final time-frequency domain source signal, use Y
nN source signal of (p, q) expression the time time-frequency domain estimated value on frequency (p, q), p=0,1,2 ..., P, q=0,1,2 ..., N
Fft-1, namely
(7) root signal time-frequency domain estimated value is recovered the time domain source signal, and concrete steps are as follows:
7-1) to each sampling instant p (p=0,1,2 ...) and frequency domain data Y
n(p, q), q=0,1,2, L, N
Fft-1 is N
FftThe IFFT conversion of point obtains time domain source signal corresponding to p sampling instant, uses y
n(p, q
t) (q
t=0,1,2, L, N
Fft-1) expression.
7-2) to above-mentioned all time domain source signal y that constantly obtain
n(p, q
t) merge processing, obtain final time domain source signal and estimate.Concrete formula is as follows:
Here K
c=N
Fft/ C, the sampling number at C windowing interval.
Before method of the present invention was elaborated, at first when adopting even linear antenna arrays to receive Frequency Hopping Signal, the characteristic of hybrid matrix described.
Suppose to have N source signal by evenly linear antenna arrays reception, the reception element number of array in this array is M.The incidence angle of supposing N source signal arrival aerial array is respectively θ
1, θ
2, L, θ
N, then evenly the response matrix of line style aerial array is:
Wherein, f
n(t) be that n source signal is in t carrier frequency constantly, c=3 * 10
8M/s is the light velocity, and d represents array element distance.Can see a certain moment, n source signal to the transmission coefficient between m the reception antenna is
This coefficient is only relevant with this moment signal carrier frequency and incident angle.
M road received signal vector x (t)=[x then
1(t) x
2(t) Lx
M(t)]
TWith original signal N frequency hopping synthesizer signal vector s (t)=[
s1 (t) s
2(t) L s
N(t)]
TBetween the pass be:
x(t)=H(t)s(t)+n(t) (11)
N (t)=[n wherein
1(t) n
2(t) L n
M(t)]
TBe the noise that M road antenna reception arrives, be generally white Gaussian noise.Since the Frequency Hopping Signal frequency at set intervals can saltus step once, so the hybrid matrix in this system is at each frequency hopping point in the duration, hybrid matrix remains unchanged, and when the Frequency generated saltus step, saltus step also occurs hybrid matrix thereupon.
By (10) as can be known
H wherein
M, n(t) (m, n) individual element of expression hybrid matrix H (t) can be seen the arbitrary column vector to hybrid matrix from following formula (12), and the ratio of its adjacent element equates that namely the hybrid matrix of even linear array has the Fan Demeng characteristic.
Embodiment one,
Saltus step characteristic and the Fan Demeng characteristic of Frequency Hopping Signal hybrid matrix when utilizing even linear array to receive, with reference to figure 1, the step of the inventive method is as follows:
Step 1), utilize and to contain the linear array array antenna received of M array element from the signal of a plurality of synchronized orthogonal frequency hopping radio sets, M the signal that receives sampled, sample frequency is f
s, obtain M discrete data, namely
Wherein k represents k sampling instant index.T
s=1/f
sThe expression sampling interval.
Step 2), the rear resulting multichannel data of sampling is carried out overlapping windowing Short Time Fourier Transform, wherein windowed function adopts rectangular window, obtains the time-frequency domain matrix of M observation signal, namely
N wherein
FftFourier transform length, N
FftLarger, frequency resolution is lower, and vice versa.C is integer, the sampling number at expression Short Time Fourier Transform windowing interval, and C is less, and temporal resolution is higher, and vice versa.General C≤N
Fft, and K
c=N
Fft/ C is integer.
Represent m road observation signal
At p moment q amplitude-frequency response that Frequency point is corresponding.Frequency index during (p, q) expression, corresponding time-frequency value is
We know that for the frequency hopping synthesizer signal multiple source signals has the probability of value very little simultaneously on synchronization, same Frequency point, especially when adopting the synchronized orthogonal networking mode, at a time on a certain Frequency point, only have at most a non-zero signal.Suppose at p
nConstantly, q
nOn the frequency, only have n signal that value is arranged, when not considering noise, by formula (11), (13)-(14) as can be known:
According to the analysis of front as can be known, n source signal is at moment p
nTransmission coefficient relevant with this signal carrier frequency constantly, and the carrier frequency of Frequency Hopping Signal is every T
h(T
hBe the frequency hopping time interval) the time saltus step once, namely at frequency hopping point b in the duration
m(p
n, q
n) be constant, namely
b
n(p
n, q
n)=[b
1(p
n, q
n), b
2(p
n, q
n), b
, M(p
n, q
n)]
TIt is constant within a certain jumping duration.Then in a certain jumping duration, the normalized vector on all Frequency points
Mainly concentrate on N the value.
Step 3), to the time-frequency domain data matrix
(m=1,2, L, M) does preliminary treatment and obtains vectorial a (p, q), specifically comprises following two steps:
3.1) right
Go low-yield preliminary treatment, namely at each sampling instant p, will
Amplitude sets to 0 less than the value of thresholding ε, obtains
Here the setting of thresholding can be determined according to the average energy that receives signal.
3.2) find out the constantly time-frequency domain data of (p=0,1,2, L, P-1) non-zero of p, use
Expression, wherein
Constantly time-frequency response of expression p
Corresponding frequency indices when non-zero to these non-zero normalization preliminary treatment, obtains pretreated vectorial b (p, q)=[b
1(p, q), b
2(p, q), L, b
M(p, q)]
T, namely
Step 4), utilize clustering algorithm to estimate jumping moment and corresponding normalized hybrid matrix column vector, the frequency hopping frequency of each jumping of each jumping.May further comprise the steps:
4-1) p (p=0,1,2 ..., P-1) constantly, right
The frequency values of expression carries out cluster, the cluster centre number that obtains
The carrier frequency number that expression p exists constantly,
Individual cluster centre then represents the size of carrier frequency, uses respectively
Expression;
4-2) to each sampling instant p, utilize clustering algorithm pair
Carry out cluster, can obtain equally
Individual cluster centre is used
Expression.
4-4) find out
The moment, use p
hExpression is to the p of the continuous value of each section
hAsk intermediate value, use
Represent the l section p that links to each other
hIntermediate value, then
Represent the estimation constantly of l frequency hopping.
4-5) obtain according to estimation in the step (4-2)
P ≠ p
hAnd estimate in the step (4-4) that the frequency hopping obtain constantly estimates each and jumps corresponding
Individual hybrid matrix column vector
Concrete formula is as follows:
Here
It is corresponding to represent that l jumps
Individual hybrid matrix column vector estimated value.
4-6) estimate the carrier frequency frequency that each jumping is corresponding, use
It is corresponding to represent that l jumps
Individual frequency estimation, computing formula is as follows:
Step 5), estimate the time-frequency domain source signal according to the normalization hybrid matrix column vector of estimating in the step (4) to obtain, concrete steps are as follows:
5-1) all sampling instant index p are judged which this moment index belongs to and jump, method is as follows:
If
Represent that then p belongs to the l jumping constantly;
5-2) to l (l=1,2, all moment p that L) jump
l, utilize formula (20) to estimate that this jumps the time-frequency domain data of each source signal, computing formula is as follows:
Step 6) splicing of source signal between the different frequency hopping point, concrete steps are as follows:
6-1) estimation l jumps correspondence
Individual incident angle is used
Represent that l jumps n the incident angle that source signal is corresponding,
Computing formula as follows:
Here
M the element that represents n the hybrid matrix column vector that l jumping estimation obtains,
6-2) utilize formula (22) judge l (l=2,3 ...) and jump the source signal and first of estimating and jump corresponding relation between the source signal of estimation
Here m
n (l)Represent that l jumps the m that estimates
n (l)Individual signal and first is jumped n the signal of estimating and is belonged to same source signal.
6-3) with different frequency hopping point estimation to the signal that belongs to same source signal be stitched together, estimate as final time-frequency domain source signal, use Y
nN source signal of (p, q) expression the time time-frequency domain estimated value on frequency (p, q), p=1,2 ...., q=1,2 ..., N
Fft, namely
Step 7) root signal time-frequency domain estimated value is recovered the time domain source signal.Concrete steps are as follows
7-1) to each sampling instant p (p=0,1,2 ..., frequency domain data Y P-1)
n(p, q), q=0,1,2, L, N
Fft-1 is N
FftThe IFFT conversion of point obtains time domain source signal corresponding to p sampling instant, uses y
n(p, q
t) (q
t=0,1,2, L, N
Fft-1) expression.
7-2) to above-mentioned all time domain source signal y that constantly obtain
n(p, q
t) merge processing, obtain final time domain source signal and estimate.Concrete formula is as follows:
Here K
c=N
Fft/ C, the sampling number at C windowing interval.
Embodiment two:
In order to verify that the present invention estimates the accuracy of the Frequency Hopping Signal that obtains, the Frequency Hopping Signal that the present invention is proposed is owed surely blind separation algorithm and is carried out emulation.Suppose to have 4 frequency hopping synthesizer signals, adopt the even linear array of 2 array elements to receive frequency hopping frequency sets [1:0.5:8] MHz, sample frequency fs=20MHz, the BPSK modulation system, character rate is respectively 20kbps, 40kbps, 25kbps, 50kbps, incident angle is respectively-20 °, and 60 °, 80 °,-40 °, since the restriction of calculator memory, each 3 frequency hopping points of simulation process, FFT counts 4096, at every turn 256 sampled points in sliding window interval.
Fig. 2-5 has provided respectively Frequency Hopping Signal that the present invention estimates and the time domain waveform figure of source signal.The Frequency Hopping Signal time domain waveform that as can be seen from the figure estimates and original Frequency Hopping Signal time domain waveform are basically identical.Owe surely blind separation algorithm performance in order quantitatively to weigh the Frequency Hopping Signal that the present invention provides, Fig. 6 has provided estimation obtains under the different signal to noise ratio conditions Frequency Hopping Signal and the coefficient correlation between the original Frequency Hopping Signal.The signal that estimates of the present invention and source signal are when signal to noise ratio 0dB as can see from Figure 6, and the time domain coefficient correlation reaches more than 0.9.
The invention has the advantages that:
1. the present invention can not know under the condition of any channel information, only according to the mixed signal of a plurality of Frequency Hopping Signals that receive, estimates the frequency hopping synthesizer signal;
2. the present invention can under the condition of reception antenna number less than the source signal number, carry out blind estimation to a plurality of Frequency Hopping Signals;
3. the present invention has only utilized Short Time Fourier Transform, and amount of calculation is little, realizes easily;
4. when the present invention just carries out blind separation to Frequency Hopping Signal, can also estimate partial parameters.
The synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster that the embodiment of the invention provides, do not knowing under the condition of any channel information, only according to the mixed signal of a plurality of Frequency Hopping Signals of receiving, estimate the frequency hopping synthesizer signal, can be under the condition of reception antenna number less than the source signal number, a plurality of Frequency Hopping Signals are carried out blind estimation, only utilized Short Time Fourier Transform, amount of calculation is little, realize that easily the method can also be estimated partial parameters when Frequency Hopping Signal is carried out blind separation, practical, have stronger propagation and employment and be worth.
The above only is preferred embodiment of the present invention, not in order to limiting the present invention, all any modifications of doing within the spirit and principles in the present invention, is equal to and replaces and improvement etc., all should be included within protection scope of the present invention.
Claims (7)
1. synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster is characterized in that the method may further comprise the steps:
Step 1 is utilized and is contained the array antenna received of M array element from the Frequency Hopping Signal of a plurality of synchronized orthogonal frequency hopping radio sets, each road is received signal sample the M road discrete time-domain mixed signal after obtaining sampling
M=1,2, L, M;
Step 2 is carried out overlapping windowing Short Time Fourier Transform to M road discrete time-domain mixed signal, obtains the time-frequency domain matrix of M mixed signal
p=0,1,L P-1,q=0,1,L N
fft-1;
Step 3 is to the frequency-hopping mixing signal time-frequency domain matrix that obtains in the step 2
Carry out preliminary treatment;
Step 4 utilizes clustering algorithm to estimate jumping moment and corresponding normalized hybrid matrix column vector, the frequency hopping frequency of each jumping of each jumping;
Step 5, the normalization hybrid matrix column vector that estimation obtains according to step 4 is estimated time-frequency domain frequency hopping synthesizer signal;
Step 6 is spliced the time-frequency domain frequency hopping synthesizer signal between the different frequency hopping points;
Step 7 according to source signal time-frequency domain estimated value, is recovered the time domain source signal.
2. the method for claim 1 is characterized in that, in step 2, and (p, q) expression time-frequency index, concrete time-frequency value is
Here Nfft represents the length of FFT conversion, and P represents the windowing number of times, and C is integer, the sampling number at expression Short Time Fourier Transform windowing interval, general C<N
Fft, and K
c=N
Fft/ C is integer, and what that is to say employing is the Short Time Fourier Transform of overlapping windowing.
3. the method for claim 1 is characterized in that, in step 3, to frequency-hopping mixing signal time-frequency domain matrix
Carry out preliminary treatment, specifically comprise following two steps:
The first step is right
Go low-yield preliminary treatment, namely at each sampling instant p, will
Amplitude sets to 0 less than the value of thresholding ε, obtains
The setting of thresholding ε can be determined according to the average energy that receives signal;
Second step is found out the constantly time-frequency domain data of (p=0,1,2, L P-1) non-zero of p, uses
Expression, wherein
Constantly time-frequency response of expression p
Corresponding frequency indices when non-zero to these non-zero normalization preliminary treatment, obtains pretreated vectorial b (p, q)=[b
1(p, q), b
2(p, q), L, b
M(p, q)]
T, wherein
4. the method for claim 1 is characterized in that, in step 4, when utilizing clustering algorithm to estimate that the jumping moment of each jumping and each are jumped corresponding normalized hybrid matrix column vector, frequency hopping frequency, may further comprise the steps:
The first step, p (p=0,1,2 ... P-1) constantly, right
The frequency values of expression carries out cluster, the cluster centre number that obtains
The carrier frequency number that expression p exists constantly,
Individual cluster centre then represents the size of carrier frequency, uses respectively
Expression;
Second step, to each sampling instant p (p=0,1,2 ... P-1), utilize clustering algorithm pair
Carry out cluster, can obtain equally
Individual cluster centre is used
Expression;
In the 4th step, find out
The moment, use p
hExpression is to the p of the continuous value of each section
hAsk intermediate value, use
Represent the l section p that links to each other
hIntermediate value, then
Represent the estimation constantly of l frequency hopping;
The 5th step is according to what estimate in the second step to obtain
P ≠ p
hAnd the 4th estimate that the frequency hopping obtain constantly estimates each and jumps corresponding in the step
Individual hybrid matrix column vector
Concrete formula is:
Here
It is corresponding to represent that l jumps
Individual hybrid matrix column vector estimated value;
The 6th step, estimate the carrier frequency frequency that each jumping is corresponding, use
It is corresponding to represent that l jumps
Individual frequency estimation, computing formula is as follows:
5. the method for claim 1 is characterized in that, in step 5, estimates time-frequency domain frequency hopping synthesizer signal according to the normalization hybrid matrix column vector of estimating in the step 4 to obtain, and concrete steps are as follows:
The first step judges that to all sampling instant index p which this moment index belongs to and jump, and concrete grammar is: if
Represent that then p belongs to the l jumping constantly; If
Represent that then p belongs to the 1st jumping constantly;
Second step, to l (l=1,2, all of L) jumping are pl constantly, estimate that this jumps the time-frequency domain data of each frequency hopping synthesizer signal, computing formula is as follows:
6. the method for claim 1 is characterized in that, in step 6, the time-frequency domain frequency hopping synthesizer signal between the different frequency hopping points is spliced, and concrete steps are as follows:
The first step, it is corresponding to estimate that l jumps
Individual incident angle is used
Represent that l jumps n the incident angle that source signal is corresponding,
Computing formula as follows:
Second step, judge l (l=2,3 ...) and source signal and first that jump to estimate jumps the corresponding relation between the source signal of estimation, judgment formula is as follows:
M wherein
n (l)Represent that l jumps the m that estimates
n (l)Individual signal and first is jumped n the signal of estimating and is belonged to same source signal;
The 3rd step, with different frequency hopping point estimation to the signal that belongs to same source signal be stitched together, estimate as final time-frequency domain source signal, use Y
nN source signal of (p, q) expression the time time-frequency domain estimated value on frequency (p, q), p=0,1,2 ...., P, q=0,1,2 ..., N
Fft-1, namely
7. the method for claim 1 is characterized in that, in step 7, when recovering the time domain source signal according to source signal time-frequency domain estimated value, concrete steps are as follows:
The first step, to each sampling instant p (p=0,1,2 ...) and frequency domain data Y
n(p, q), q=0,1,2, L, N
Fft-1 is N
FftThe IFFT conversion of point obtains time domain source signal corresponding to p sampling instant, uses y
n(p, q
t) (qt=0,1,2, L, N
Fft-1) expression;
Second step is to above-mentioned all time domain source signal y that constantly obtain
n(p, q
t) merge processing, obtain final time domain source signal and estimate, concrete formula is as follows:
Here K
c=N
Fft/ C, the sampling number at C windowing interval.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210525064.4A CN103051367B (en) | 2012-11-27 | 2012-11-27 | A kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210525064.4A CN103051367B (en) | 2012-11-27 | 2012-11-27 | A kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103051367A true CN103051367A (en) | 2013-04-17 |
CN103051367B CN103051367B (en) | 2015-08-26 |
Family
ID=48063881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210525064.4A Expired - Fee Related CN103051367B (en) | 2012-11-27 | 2012-11-27 | A kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103051367B (en) |
Cited By (71)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103346989A (en) * | 2013-05-06 | 2013-10-09 | 中国人民解放军重庆通信学院 | Multi-path frequency hopping-based single-channel blind source separation anti-interference communication system and method thereof |
CN103746722A (en) * | 2014-01-02 | 2014-04-23 | 东南大学 | Method for estimating jump cycle and take-off time of frequency hopping signal |
CN103871422A (en) * | 2014-01-08 | 2014-06-18 | 西安电子科技大学 | Time-domain aliasing blind signal separation method based on three-factor iteration joint block diagonalization |
CN104485979A (en) * | 2014-12-09 | 2015-04-01 | 西安电子科技大学 | Blind estimation method for underdetermined hybrid frequency hopping parameters based on time frequency diagram correction |
CN105337636A (en) * | 2015-10-08 | 2016-02-17 | 西安电子科技大学 | Asynchronous frequency hopping signal parameter blind evaluation method based on frequency splicing |
CN105844094A (en) * | 2016-03-22 | 2016-08-10 | 西安电子科技大学 | Gradient descent method and Newton method based underdetermined blind source separation source signal recovery method |
CN106202756A (en) * | 2016-07-15 | 2016-12-07 | 西安电子科技大学 | Based on monolayer perceptron owing determines blind source separating source signal restoration methods |
CN106203549A (en) * | 2016-07-05 | 2016-12-07 | 广东省心血管病研究所 | A kind of cord blood stem cell management system |
CN106226829A (en) * | 2016-09-21 | 2016-12-14 | 刘俊玉 | A kind of geologic structure face occurrence determinator |
CN106291751A (en) * | 2016-09-18 | 2017-01-04 | 王光栋 | A kind of detecting system of geologic(al) formation |
CN106339570A (en) * | 2016-06-20 | 2017-01-18 | 丽水学院 | Traditional Chinese medicine intelligent assistant prescription system |
CN106355968A (en) * | 2016-11-16 | 2017-01-25 | 北华大学 | Japanese teaching system based on intelligent terminal |
CN106355538A (en) * | 2016-08-26 | 2017-01-25 | 中国地质科学院探矿工艺研究所 | Collapse prevention and control engineering system |
CN106344354A (en) * | 2016-10-28 | 2017-01-25 | 张立沼 | Orthopedic walking aid |
CN106373018A (en) * | 2016-08-29 | 2017-02-01 | 巴中市方圆环保科技发展有限责任公司 | Internet management system used for konjac planting |
CN106383349A (en) * | 2016-08-31 | 2017-02-08 | 贵州省江口县气象局 | Rainfall estimating system and method based on X-waveband Doppler radar |
CN106446250A (en) * | 2016-10-14 | 2017-02-22 | 张成良 | Forecasting system for individual online news comment emotion based on cloud computing |
CN106483902A (en) * | 2016-11-23 | 2017-03-08 | 北华大学 | A kind of automatically controlled control system of the barotor based on internet |
CN106510988A (en) * | 2016-12-21 | 2017-03-22 | 王秀峰 | Intelligent wheelchair supporting intelligent terminal mechanical structure |
CN106510991A (en) * | 2016-12-21 | 2017-03-22 | 王秀峰 | Multifunctional wheelchair |
CN106571073A (en) * | 2016-07-27 | 2017-04-19 | 张法春 | English word root shape tone memorizing method |
CN106563262A (en) * | 2016-07-12 | 2017-04-19 | 山东医学高等专科学校 | Cloud computing-based sports detection system |
CN106658921A (en) * | 2017-01-11 | 2017-05-10 | 河南工程学院 | Lamp control switch system |
CN106647514A (en) * | 2016-12-28 | 2017-05-10 | 安徽工程大学 | Cement enterprise carbon emission real-time on-line monitoring management system |
CN106646661A (en) * | 2016-09-22 | 2017-05-10 | 张春志 | Comprehensive survey system for hydrogeology of mineral deposit |
CN106647861A (en) * | 2017-01-23 | 2017-05-10 | 阜阳师范学院 | Heating device based on mobile terminal control for fur-bearing animal nest box |
CN106686006A (en) * | 2017-03-01 | 2017-05-17 | 湖南财政经济学院 | Big data implementation control system based on Internet network security |
CN106691485A (en) * | 2016-07-18 | 2017-05-24 | 山东省肿瘤防治研究院 | Lung neoplasm diagnosis and treatment device under radiography guidance |
CN106781868A (en) * | 2017-01-04 | 2017-05-31 | 汤应 | A kind of tutoring system for College Computer Education |
CN106730213A (en) * | 2016-12-15 | 2017-05-31 | 湖北民族学院附属民大医院 | One kind digitlization anesthesia control system |
CN106803834A (en) * | 2016-11-16 | 2017-06-06 | 北华大学 | A kind of direct sequence spread spectrum train-ground communication control system based on internet |
CN106821425A (en) * | 2017-01-13 | 2017-06-13 | 贵定县人民医院 | A kind of binary channels bladder lithotriptoscepe system |
CN106887038A (en) * | 2017-02-20 | 2017-06-23 | 张光辉 | A kind of orthopaedics inside-fixture formation system and method |
CN106897558A (en) * | 2017-02-24 | 2017-06-27 | 王雯钰 | A kind of Chinese medicine nursing percussion device for rehabilitation control system |
CN106907761A (en) * | 2017-01-23 | 2017-06-30 | 房嵩然 | A kind of new type solar energy heating floor |
CN106901941A (en) * | 2017-02-21 | 2017-06-30 | 哈尔滨医科大学 | A kind of joint of vertebral column depressurized system |
CN106911767A (en) * | 2017-01-23 | 2017-06-30 | 常梅 | A kind of obstetrical and gynecological disease Medical treatment control system based on mobile terminal control |
CN106933984A (en) * | 2017-02-20 | 2017-07-07 | 周长英 | The dispatching method and system of a kind of distributed file system |
CN106952527A (en) * | 2017-05-09 | 2017-07-14 | 王学斌 | A kind of Modern digital electronic technological synthesis real training control system based on internet |
CN106992588A (en) * | 2017-05-02 | 2017-07-28 | 湖南工程学院 | A kind of integrated high intelligence control system of solar wind-energy based on Internet of Things |
CN107049464A (en) * | 2017-06-07 | 2017-08-18 | 任小宝 | A kind of cable formula inner fixing device for orthopaedics |
CN107050550A (en) * | 2017-03-30 | 2017-08-18 | 汪可 | A kind of department of general surgery's wound-cleaning device |
CN107103180A (en) * | 2017-03-13 | 2017-08-29 | 王靖 | A kind of multi-functional EICU nursing systems |
CN107145608A (en) * | 2017-06-12 | 2017-09-08 | 武汉梦之蓝科技有限公司 | A kind of intelligent environmental sanitation platform application is in the processing management control system of garbage apparatus |
CN107155161A (en) * | 2017-05-02 | 2017-09-12 | 贝尔顺(中国)科技股份有限公司 | A kind of bone conduction earphone hearing aid control system and control method |
CN107149470A (en) * | 2017-01-23 | 2017-09-12 | 常梅 | A kind of gynemetrics's electronic pulse meter control system based on internet |
CN107168057A (en) * | 2017-05-24 | 2017-09-15 | 颜振刚 | A kind of intelligence infection ward thimerosal configuration system |
CN107448075A (en) * | 2017-08-09 | 2017-12-08 | 王怀玲 | A kind of electric back door lock external lock signal instruction system |
CN107451522A (en) * | 2017-04-28 | 2017-12-08 | 山东省农业可持续发展研究所 | A kind of agricultural arid monitoring and early alarming and forecasting method |
CN107491037A (en) * | 2017-09-29 | 2017-12-19 | 华北科技学院 | A kind of Digit Control Machine Tool WeChat ID monitoring system |
CN107550471A (en) * | 2017-04-05 | 2018-01-09 | 李霄艳 | A kind of GI Medicine UGB nursing system |
CN107595288A (en) * | 2017-09-19 | 2018-01-19 | 华克胜 | A kind of enhancing respiratory function tester |
CN107622294A (en) * | 2017-09-18 | 2018-01-23 | 重庆科技学院 | A kind of Statistics Application principle and planar bar code technology control system |
CN107650713A (en) * | 2017-09-25 | 2018-02-02 | 盐城工学院 | A kind of method and system of automobile wireless charging transmitter module |
CN107707445A (en) * | 2017-11-24 | 2018-02-16 | 安徽省雨龙家具有限公司 | A kind of information-pushing method based on intelligent domestic system |
CN107787816A (en) * | 2017-10-30 | 2018-03-13 | 汤立志 | A kind of automatic garden irrigation system for greenbelt |
CN107860781A (en) * | 2017-10-12 | 2018-03-30 | 兰州职业技术学院 | A kind of detection method of food microorganisms qualitative and quantitative |
CN107886324A (en) * | 2017-11-29 | 2018-04-06 | 重庆科技学院 | A kind of small amount payment terminal of city smart cards control method based on 5G |
CN107951491A (en) * | 2017-09-26 | 2018-04-24 | 齐鲁师范学院 | A kind of novel visual motion tracking training system |
CN107961087A (en) * | 2016-10-19 | 2018-04-27 | 宋雨泽 | A kind of Intelligent orthopaedic operation auxiliary frame |
CN108272746A (en) * | 2018-03-06 | 2018-07-13 | 佳木斯大学附属第二医院(佳木斯大学附属口腔医院) | A kind of doctor trained in Western medicine surgery skin analgesia cream and preparation method thereof |
CN108297437A (en) * | 2018-02-05 | 2018-07-20 | 广东水利电力职业技术学院(广东省水利电力技工学校) | A kind of control system of 3D printer |
CN108389444A (en) * | 2018-03-29 | 2018-08-10 | 湖南城市学院 | A kind of English language tutoring system and teaching application method |
CN108755897A (en) * | 2018-05-28 | 2018-11-06 | 湖南城市学院 | A kind of energy saving building rain water collecting system |
CN108784708A (en) * | 2018-04-03 | 2018-11-13 | 长江师范学院 | A kind of human motion overall analysis system based on big data |
CN108877831A (en) * | 2018-08-28 | 2018-11-23 | 山东大学 | Blind source separating fast method and system based on multi-standard fusion frequency point screening |
CN109214469A (en) * | 2018-10-24 | 2019-01-15 | 西安交通大学 | A kind of source signal separation method based on non-negative tensor resolution |
CN109542146A (en) * | 2018-05-29 | 2019-03-29 | 菏泽学院 | The long-range plant physiological ecology monitor control system in greenhouse based on big data analysis |
CN109700457A (en) * | 2019-03-05 | 2019-05-03 | 邢台医学高等专科学校 | A kind of ECG acquisition platform |
CN109861721A (en) * | 2019-03-06 | 2019-06-07 | 重庆邮电大学 | A kind of compressed sensing based Frequency Hopping Signal sparse decomposition method |
CN113452409A (en) * | 2021-06-24 | 2021-09-28 | 湖南基石通信技术有限公司 | Frequency hopping implementation method of asynchronous ad hoc network |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105496702A (en) * | 2016-02-25 | 2016-04-20 | 滕艳玲 | Surgical nursing dressing-changing device |
CN107875591A (en) * | 2017-12-08 | 2018-04-06 | 湖南城市学院 | A kind of interactive physical training system |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102546116A (en) * | 2012-02-14 | 2012-07-04 | 西安电子科技大学 | Compressive-sensing-based frequency hopping signal blind-detection method |
-
2012
- 2012-11-27 CN CN201210525064.4A patent/CN103051367B/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102546116A (en) * | 2012-02-14 | 2012-07-04 | 西安电子科技大学 | Compressive-sensing-based frequency hopping signal blind-detection method |
Non-Patent Citations (2)
Title |
---|
WEIHONG FU 等: "Step-size Optimization EASI Algorithm for Blind Source Separation", 《IEEE COMPUTER SOCIETY》 * |
付卫红 等: "基于盲源分离的CDMA多用户检测与伪码估计", 《电子学报》 * |
Cited By (78)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103346989B (en) * | 2013-05-06 | 2016-05-25 | 中国人民解放军重庆通信学院 | The blind source of single channel based on multichannel frequency hopping separates jam-resistant communication system and method |
CN103346989A (en) * | 2013-05-06 | 2013-10-09 | 中国人民解放军重庆通信学院 | Multi-path frequency hopping-based single-channel blind source separation anti-interference communication system and method thereof |
CN103746722A (en) * | 2014-01-02 | 2014-04-23 | 东南大学 | Method for estimating jump cycle and take-off time of frequency hopping signal |
CN103746722B (en) * | 2014-01-02 | 2015-07-01 | 东南大学 | Method for estimating jump cycle and take-off time of frequency hopping signal |
CN103871422A (en) * | 2014-01-08 | 2014-06-18 | 西安电子科技大学 | Time-domain aliasing blind signal separation method based on three-factor iteration joint block diagonalization |
CN104485979A (en) * | 2014-12-09 | 2015-04-01 | 西安电子科技大学 | Blind estimation method for underdetermined hybrid frequency hopping parameters based on time frequency diagram correction |
CN105337636A (en) * | 2015-10-08 | 2016-02-17 | 西安电子科技大学 | Asynchronous frequency hopping signal parameter blind evaluation method based on frequency splicing |
CN105337636B (en) * | 2015-10-08 | 2017-11-21 | 西安电子科技大学 | Asynchronous Frequency Hopping Signal Blind Parameter Estimation based on frequency splicing |
CN105844094A (en) * | 2016-03-22 | 2016-08-10 | 西安电子科技大学 | Gradient descent method and Newton method based underdetermined blind source separation source signal recovery method |
CN105844094B (en) * | 2016-03-22 | 2018-03-06 | 西安电子科技大学 | Blind source separating source signal restoration methods are determined based on gradient descent method and the deficient of Newton method |
CN106339570A (en) * | 2016-06-20 | 2017-01-18 | 丽水学院 | Traditional Chinese medicine intelligent assistant prescription system |
CN106203549B (en) * | 2016-07-05 | 2018-09-14 | 广东省心血管病研究所 | A kind of cord blood stem cell management system |
CN106203549A (en) * | 2016-07-05 | 2016-12-07 | 广东省心血管病研究所 | A kind of cord blood stem cell management system |
CN106563262A (en) * | 2016-07-12 | 2017-04-19 | 山东医学高等专科学校 | Cloud computing-based sports detection system |
CN106202756A (en) * | 2016-07-15 | 2016-12-07 | 西安电子科技大学 | Based on monolayer perceptron owing determines blind source separating source signal restoration methods |
CN106202756B (en) * | 2016-07-15 | 2019-06-18 | 西安电子科技大学 | Deficient based on single layer perceptron determines blind source separating source signal restoration methods |
CN106691485A (en) * | 2016-07-18 | 2017-05-24 | 山东省肿瘤防治研究院 | Lung neoplasm diagnosis and treatment device under radiography guidance |
CN106571073A (en) * | 2016-07-27 | 2017-04-19 | 张法春 | English word root shape tone memorizing method |
CN106355538A (en) * | 2016-08-26 | 2017-01-25 | 中国地质科学院探矿工艺研究所 | Collapse prevention and control engineering system |
CN106373018A (en) * | 2016-08-29 | 2017-02-01 | 巴中市方圆环保科技发展有限责任公司 | Internet management system used for konjac planting |
CN106383349A (en) * | 2016-08-31 | 2017-02-08 | 贵州省江口县气象局 | Rainfall estimating system and method based on X-waveband Doppler radar |
CN106291751A (en) * | 2016-09-18 | 2017-01-04 | 王光栋 | A kind of detecting system of geologic(al) formation |
CN106226829A (en) * | 2016-09-21 | 2016-12-14 | 刘俊玉 | A kind of geologic structure face occurrence determinator |
CN106646661A (en) * | 2016-09-22 | 2017-05-10 | 张春志 | Comprehensive survey system for hydrogeology of mineral deposit |
CN106446250A (en) * | 2016-10-14 | 2017-02-22 | 张成良 | Forecasting system for individual online news comment emotion based on cloud computing |
CN107961087A (en) * | 2016-10-19 | 2018-04-27 | 宋雨泽 | A kind of Intelligent orthopaedic operation auxiliary frame |
CN106344354A (en) * | 2016-10-28 | 2017-01-25 | 张立沼 | Orthopedic walking aid |
CN106803834A (en) * | 2016-11-16 | 2017-06-06 | 北华大学 | A kind of direct sequence spread spectrum train-ground communication control system based on internet |
CN106355968A (en) * | 2016-11-16 | 2017-01-25 | 北华大学 | Japanese teaching system based on intelligent terminal |
CN106483902B (en) * | 2016-11-23 | 2019-07-05 | 北华大学 | A kind of automatically controlled control system of barotor Internet-based |
CN106483902A (en) * | 2016-11-23 | 2017-03-08 | 北华大学 | A kind of automatically controlled control system of the barotor based on internet |
CN106730213A (en) * | 2016-12-15 | 2017-05-31 | 湖北民族学院附属民大医院 | One kind digitlization anesthesia control system |
CN106510988A (en) * | 2016-12-21 | 2017-03-22 | 王秀峰 | Intelligent wheelchair supporting intelligent terminal mechanical structure |
CN106510991A (en) * | 2016-12-21 | 2017-03-22 | 王秀峰 | Multifunctional wheelchair |
CN106647514A (en) * | 2016-12-28 | 2017-05-10 | 安徽工程大学 | Cement enterprise carbon emission real-time on-line monitoring management system |
CN106781868A (en) * | 2017-01-04 | 2017-05-31 | 汤应 | A kind of tutoring system for College Computer Education |
CN106658921A (en) * | 2017-01-11 | 2017-05-10 | 河南工程学院 | Lamp control switch system |
CN106821425A (en) * | 2017-01-13 | 2017-06-13 | 贵定县人民医院 | A kind of binary channels bladder lithotriptoscepe system |
CN106647861A (en) * | 2017-01-23 | 2017-05-10 | 阜阳师范学院 | Heating device based on mobile terminal control for fur-bearing animal nest box |
CN106907761A (en) * | 2017-01-23 | 2017-06-30 | 房嵩然 | A kind of new type solar energy heating floor |
CN106911767A (en) * | 2017-01-23 | 2017-06-30 | 常梅 | A kind of obstetrical and gynecological disease Medical treatment control system based on mobile terminal control |
CN107149470A (en) * | 2017-01-23 | 2017-09-12 | 常梅 | A kind of gynemetrics's electronic pulse meter control system based on internet |
CN106933984A (en) * | 2017-02-20 | 2017-07-07 | 周长英 | The dispatching method and system of a kind of distributed file system |
CN106887038A (en) * | 2017-02-20 | 2017-06-23 | 张光辉 | A kind of orthopaedics inside-fixture formation system and method |
CN106901941A (en) * | 2017-02-21 | 2017-06-30 | 哈尔滨医科大学 | A kind of joint of vertebral column depressurized system |
CN106897558A (en) * | 2017-02-24 | 2017-06-27 | 王雯钰 | A kind of Chinese medicine nursing percussion device for rehabilitation control system |
CN106686006A (en) * | 2017-03-01 | 2017-05-17 | 湖南财政经济学院 | Big data implementation control system based on Internet network security |
CN107103180A (en) * | 2017-03-13 | 2017-08-29 | 王靖 | A kind of multi-functional EICU nursing systems |
CN107050550A (en) * | 2017-03-30 | 2017-08-18 | 汪可 | A kind of department of general surgery's wound-cleaning device |
CN107550471A (en) * | 2017-04-05 | 2018-01-09 | 李霄艳 | A kind of GI Medicine UGB nursing system |
CN107451522A (en) * | 2017-04-28 | 2017-12-08 | 山东省农业可持续发展研究所 | A kind of agricultural arid monitoring and early alarming and forecasting method |
CN107155161A (en) * | 2017-05-02 | 2017-09-12 | 贝尔顺(中国)科技股份有限公司 | A kind of bone conduction earphone hearing aid control system and control method |
CN106992588A (en) * | 2017-05-02 | 2017-07-28 | 湖南工程学院 | A kind of integrated high intelligence control system of solar wind-energy based on Internet of Things |
CN106952527A (en) * | 2017-05-09 | 2017-07-14 | 王学斌 | A kind of Modern digital electronic technological synthesis real training control system based on internet |
CN107168057A (en) * | 2017-05-24 | 2017-09-15 | 颜振刚 | A kind of intelligence infection ward thimerosal configuration system |
CN107049464A (en) * | 2017-06-07 | 2017-08-18 | 任小宝 | A kind of cable formula inner fixing device for orthopaedics |
CN107145608A (en) * | 2017-06-12 | 2017-09-08 | 武汉梦之蓝科技有限公司 | A kind of intelligent environmental sanitation platform application is in the processing management control system of garbage apparatus |
CN107448075A (en) * | 2017-08-09 | 2017-12-08 | 王怀玲 | A kind of electric back door lock external lock signal instruction system |
CN107622294A (en) * | 2017-09-18 | 2018-01-23 | 重庆科技学院 | A kind of Statistics Application principle and planar bar code technology control system |
CN107595288A (en) * | 2017-09-19 | 2018-01-19 | 华克胜 | A kind of enhancing respiratory function tester |
CN107650713A (en) * | 2017-09-25 | 2018-02-02 | 盐城工学院 | A kind of method and system of automobile wireless charging transmitter module |
CN107951491A (en) * | 2017-09-26 | 2018-04-24 | 齐鲁师范学院 | A kind of novel visual motion tracking training system |
CN107491037A (en) * | 2017-09-29 | 2017-12-19 | 华北科技学院 | A kind of Digit Control Machine Tool WeChat ID monitoring system |
CN107860781A (en) * | 2017-10-12 | 2018-03-30 | 兰州职业技术学院 | A kind of detection method of food microorganisms qualitative and quantitative |
CN107787816A (en) * | 2017-10-30 | 2018-03-13 | 汤立志 | A kind of automatic garden irrigation system for greenbelt |
CN107707445A (en) * | 2017-11-24 | 2018-02-16 | 安徽省雨龙家具有限公司 | A kind of information-pushing method based on intelligent domestic system |
CN107886324A (en) * | 2017-11-29 | 2018-04-06 | 重庆科技学院 | A kind of small amount payment terminal of city smart cards control method based on 5G |
CN108297437A (en) * | 2018-02-05 | 2018-07-20 | 广东水利电力职业技术学院(广东省水利电力技工学校) | A kind of control system of 3D printer |
CN108272746A (en) * | 2018-03-06 | 2018-07-13 | 佳木斯大学附属第二医院(佳木斯大学附属口腔医院) | A kind of doctor trained in Western medicine surgery skin analgesia cream and preparation method thereof |
CN108389444A (en) * | 2018-03-29 | 2018-08-10 | 湖南城市学院 | A kind of English language tutoring system and teaching application method |
CN108784708A (en) * | 2018-04-03 | 2018-11-13 | 长江师范学院 | A kind of human motion overall analysis system based on big data |
CN108755897A (en) * | 2018-05-28 | 2018-11-06 | 湖南城市学院 | A kind of energy saving building rain water collecting system |
CN109542146A (en) * | 2018-05-29 | 2019-03-29 | 菏泽学院 | The long-range plant physiological ecology monitor control system in greenhouse based on big data analysis |
CN108877831A (en) * | 2018-08-28 | 2018-11-23 | 山东大学 | Blind source separating fast method and system based on multi-standard fusion frequency point screening |
CN109214469A (en) * | 2018-10-24 | 2019-01-15 | 西安交通大学 | A kind of source signal separation method based on non-negative tensor resolution |
CN109700457A (en) * | 2019-03-05 | 2019-05-03 | 邢台医学高等专科学校 | A kind of ECG acquisition platform |
CN109861721A (en) * | 2019-03-06 | 2019-06-07 | 重庆邮电大学 | A kind of compressed sensing based Frequency Hopping Signal sparse decomposition method |
CN113452409A (en) * | 2021-06-24 | 2021-09-28 | 湖南基石通信技术有限公司 | Frequency hopping implementation method of asynchronous ad hoc network |
Also Published As
Publication number | Publication date |
---|---|
CN103051367B (en) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103051367B (en) | A kind of synchronized orthogonal Frequency Hopping Signal blind source separation method based on cluster | |
CN108462509B (en) | Asynchronous frequency hopping net platform method for separating based on time-frequency figure information | |
CN103607361B (en) | The parameter prediction method of time-frequency overlapped signal under Alpha Stable distritation noise | |
CN106772219A (en) | Indoor orientation method based on CSI signals | |
CN104410590A (en) | Short-wave OFDM (Orthogonal Frequency Division Multiplexing) interference suppression joint channel estimation method based on compressed sensing | |
CN103475602B (en) | MIMO-OFDM channel estimation method considering same frequency interference | |
CN103701730B (en) | Channel estimation method and device based on channel time-domain correlation and low-complexity compressed sensing | |
CN103235294A (en) | Method for estimating weak signal separation on basis of positioning for external radiation sources | |
CN101808053A (en) | Channel estimation method and device based on OFDM (Orthogonal Frequency Division Multiplexing) | |
CN102857996A (en) | Cell search timing synchronization method | |
WO2014063275A1 (en) | Method for determining remote same-frequency interference source and locating method therefor | |
CN104735014B (en) | A kind of time synchronization method related based on leading symbol difference | |
CN103107969A (en) | Incremental iterative time-varying channel evaluation and inter carrier interference (ICI) elimination method of fast orthogonal frequency division multiplexing (OFDM) system | |
CN103210603A (en) | Estimation of a multipath signal in a wireless communication system | |
CN106597381A (en) | Full coherent full polarization MIMO radar four-channel integrated target detecting method | |
CN103095638A (en) | Orthogonal Frequency Division Multiplexing (OFDM) system sampling frequency shift blind estimation method under multipath fading channel | |
CN103763222A (en) | Channel ambiguity removing method in MIMO signal blind detection process | |
WO2009142407A3 (en) | Mimo detection method and system based on channel correlation | |
CN105337906A (en) | Channel estimation method and device | |
CN105187352A (en) | Integer frequency offset estimation method based on OFDM preamble | |
CN104301281A (en) | Transmitting antenna number estimation method for MIMO-OFDM system under frequency selective fading channel | |
CN102215184B (en) | Method and system for estimating uplink timing error | |
CN106510988A (en) | Intelligent wheelchair supporting intelligent terminal mechanical structure | |
CN101378372B (en) | Estimation method for channel frequency response | |
CN101459455B (en) | Channel estimation method, signal detection method and apparatus |
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 | ||
DD01 | Delivery of document by public notice |
Addressee: Patent of Xi'an University of Electronic Science and technology The person in charge Document name: payment instructions |
|
DD01 | Delivery of document by public notice | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20211127 |
|
CF01 | Termination of patent right due to non-payment of annual fee |