RU2611102C1  Method for spectral analysis of polyharmonic signals  Google Patents
Method for spectral analysis of polyharmonic signals Download PDFInfo
 Publication number
 RU2611102C1 RU2611102C1 RU2015153710A RU2015153710A RU2611102C1 RU 2611102 C1 RU2611102 C1 RU 2611102C1 RU 2015153710 A RU2015153710 A RU 2015153710A RU 2015153710 A RU2015153710 A RU 2015153710A RU 2611102 C1 RU2611102 C1 RU 2611102C1
 Authority
 RU
 Russia
 Prior art keywords
 signal
 spectrum
 spectral analysis
 samples
 zero
 Prior art date
Links
 238000010183 spectrum analysis Methods 0 abstract title 3
 238000005516 engineering processes Methods 0 abstract 2
 238000004891 communication Methods 0 abstract 1
 230000000694 effects Effects 0 abstract 1
 238000001228 spectrum Methods 0 abstract 1
 239000000126 substances Substances 0 abstract 1
 230000001502 supplementation Effects 0 abstract 1
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
 G01R23/16—Spectrum analysis; Fourier analysis
Abstract
Description
The invention relates to digital signal processing, in particular to spectral analysis of signals in the Fourier basis, and can be used in radar, radio communications and measuring equipment.
A known method of spectral analysis using weighted windows of Hanning, Kaiser, DolphChebyshev and others (L. Rabiner, B. Gold. Theory and application of digital signal processing. M., 1978)
where n = 0, ... N1, is the number of time reference,
m = 0, ... M1, is the number of the spectral component, M = N,
Signal weighting is used to reduce side lobe levels. The side lobes not only affect the spreading of the spectrum, in which the monoharmonic signal is reproduced by a whole group of spectral components; their level means the transmission coefficient in the suppression band, and the gain in the signaltonoise ratio obtained by spectral transformation of the signal depends on its magnitude.
The best of the "classical" window functions, such as Kaiser, Hanning, DolphChebyshev, atomic functions of Kravchenko, etc., provide suppression of the "tails" of the frequency response to levels (80 ÷ 100) dB. However, suppression of the side lobes is achieved at the cost of a twofold decrease in resolution due to the expansion of the main lobe of the frequency response of the spectrum analyzer. This is a common drawback for all window functions: whatever the shape of the weight window, it narrows the already short time interval of the signal, and narrowing the interval in the time domain inevitably leads to the expansion of the spectrum.
Twofold expansion of the main lobe not only degrades the resolution, but also doubles the noise level. The efficiency of isolating a narrowband signal from noise by a window function is usually estimated by an equivalent noise band (ESR) equal to the bandwidth of a conventional rectangular frequency filter, whose area under the amplitude frequency response curve is equal to the corresponding area of the window function, and the amplitude is equal to the effective value of the transmission coefficient of its main lobe :
ESR is measured in bins  the number of steps in frequency. For an ideal (rectangular) frequency response ESR = 1. For a rectangular time window ESR = 4.1, for a window of Hanning ESR = 4.01. Such a slight difference shows the inefficiency of using window functions to suppress noise. For comparison, the proposed method gives ESR within 1.1 ÷ 1.5.
In FIG. Figure 1 shows the frequency characteristics of two adjacent channels of the fast Fourier transform (FFT)  a spectrum analyzer. The bold line shows the desired (ideal) characteristics. In the ideal case, the squareness of the characteristic ensures the independence of the transmission coefficient from the frequency and guarantees the absence of a signal skip at the boundary between the channels in the absence of their overlap. The more gentle the borders, the more overlap is necessary for reproducing the “boundary” frequencies. For the FFT method, the real characteristic of the first zeros is twice as wide as the ideal one. When using the Hanning window, the characteristic becomes twice as wider. Not much smaller expansion of the main petal is provided by the windows of DolphChebyshev and KravchenkoGauss.
It is also known that in the case where resolution is the main parameter, you have to abandon the method of direct Fourier transform of the original signal, and use such "superresolution" methods as the Pisarenko, Proni, and Kapon methods (Modern methods of spectral analysis: Review. S.M. Kay, S.L. Marple, TIIER, vol. 69, No. 11, November 1981). However, these methods are characterized by the appearance of false frequencies and a change in amplitude from the number of analyzed signals.
There are also known methods for eliminating the spreading effect of the spectrum by adjusting the sampling frequency of the signal to the fundamental frequency when analyzing polyharmonic signals. So, in the method of spectral analysis of polyharmonic signals and a device for its implementation (RU 2363005), such adjustment is carried out by multiplying the basic Fourier functions by the calculated coefficient. This method is applicable only for large signaltonoise ratios, when the fundamental frequency can be measured without resorting to spectral processing.
There is also a known method of spectral analysis of polyharmonic signals (Middleton V. Introduction to the statistical theory of communication. T. 2.  M .: Sov. Radio, 1962), in which for the analysis of polyharmonic signals in the frequency domain, the nonamplitude spectrum is often used  the complex samples module obtained by the FFT procedure, and the instantaneous power spectral density (SPM):
Given the complex nature of S (m), the operations of calculating the PSD are represented as the sum of the products of vectors of the real and imaginary parts (we multiply the complex and imaginary parts of the signal spectra elementwise, respectively):
Multiplying the material and imaginary parts of the samples by their copies does not eliminate the side lobes of the frequency response, however, in terms of the combination of signal conversion operations, it is the closest analogue to the prototype of the proposed method.
The technical result of the proposed method is to reduce the level of the side lobes without compromising the resolution of the spectral analysis, which entails an increase in the signaltonoise ratio.
The indicated technical result is achieved by the fact that the Fourier transform is performed on the sample of the analyzed signal, the obtained real and imaginary parts of the complex spectrum samples are multiplied by their copies and the products are summed, while the original signal sample is supplemented with zero intervals in front and behind, and the real and imaginary parts of the spectrum samples multiplied by the same parts of the neighboring sample and after adding the product vectors, the total vector is multiplied by minus one and all the samples are reset to zero.
The essence of the claimed technical solution is illustrated in FIG. 28.
In FIG. 2, the solid line shows the frequency response of two adjacent channels obtained by the proposed method, and the dotted line shows the frequency response of the FFT channel. Black color shows the overlap zone, in which the spectrum of the monoharmonic signal will contain two nonzero samples instead of one. The width of the main lobe is almost twice as narrow as the FFT lobe.
In FIG. 3 shows the frequency characteristics of one FFT channel (dotted line) and obtained by the proposed method. The width of the main lobe can be adjusted by changing the number of zero samples.
In FIG. 4a shows the spectrum of a sinusoidal signal obtained by the proposed method (Fig. 4a) before the operation of zeroing negative samples of the spectrum, and in FIG. 4b shows the spectrum obtained by the FFT method.
In FIG. Figure 5 shows the phase change in the spectrum of a monoharmonic signal of frequency 11.5 when noise is applied (FFT with a rectangular window)  even in noise, the phase difference between 11 and 12 samples remains equal to π.
In FIG. 6a shows the amplitude spectrum of a noisy signal; FIG. 6b is a spectrum obtained by the proposed method from the same FFT samples.
In FIG. 7 shows the frequency characteristics of the Hanning weighing window and PS.
The sequence of signal conversion operations of the proposed method is presented in paragraphs 14:
1. The signal sample x (n), n = 1, 2, ... N, is supplemented on both sides by k zeros:
2. The Fourier transform (FFT) is performed:
3. The real and imaginary parts of each mth spectral sample are multiplied by the same parts of the next (m + 1) th sample, and the products are summed:
S (m) = a (m) + jb (m);
P (m) = a (m) ⋅ a (m + 1) + b (m) ⋅ b (m + 1).
4. The sign of the vector P is reversed, after which the samples with a negative sign are assigned a zero value or the value of the selected small value in the case of recalculation of P in decibels:
Comparing the result with relation (4), one can write the operations listed in the form:
We call the obtained characteristic "the product of neighboring samples of the spectrum (PS)."
The proposed method is based on the following regularity of the phasefrequency characteristic of the main lobe of any channel of the Fourier spectrum: on the entire frequency interval of the main lobe, the phase varies linearly from π to + π. A monoharmonic signal that has got into two adjacent channels of a spectrum analyzer due to overlap will form two spectral coefficients, which can be represented as vectors on the complex plane. The linearity of the phasefrequency characteristic provides the opposite direction of these vectors  a constant phase shift of 180 degrees. The most efficient way to isolate such vectors is their scalar product, which is formed by the proposed method.
Regardless of the phase and frequency of the monoharmonic signal, both the real and imaginary parts of the complex coefficients of adjacent spectrum channels will have opposite signs. Therefore, the product of their material parts, as well as the product of their imaginary parts, will be negative. Their sum will also be negative. To bring the result to its usual form, we multiply the sum by 1. Now all samples with a minus sign can be discarded as not belonging to the main lobe.
This method can be implemented using computing devices.
The main advantages of the proposed method are:
1. The width of the main lobe is the smallest of all known window functions and even 1.42 times narrower than a rectangular window. This minimizes spreading of the spectrum.
2. Low level of side lobes  most of them are zeros, and the total weight of all side lobes is 0.16% of the main lobe.
3. The rectangularity of the main lobe (the ratio of the width of the lobe at 0.5 to the width at the level of zeros) is significantly higher than that of existing windows (Table 1).
4. The location of the maximum of the main lobe corresponds to the middle of the interval between the samples of the spectrum, which is most convenient for the combined FFTPS processing, which guarantees a twofold increase in resolution.
5. The width of the main lobe can be adjusted by the number of zeros added, gradually increasing it to ΔF = 2 (Fig. 5). Since the addition of zeros leads to an increase in the time base of the analysis, at the same time, the resolution can be doubled.
6. This is the most significant property of the proposed method: it gives an increase in signaltonoise ratio by 1020 dB ;.
This result is due to the combined influence of the following factors:
 narrowing the main lobe almost twice and a sharp decrease in the side lobes gives a total gain, which can be estimated by the ratio of the ESR of the rectangular window and the PS:
 the basis of the method of PS formation is the scalar product of two vectors a * b =  a  b  cosϕ, which allows you to sharply increase the amplitude of the inphase (antiphase) factors and throw out (zero) about 50% of the noise due to the use of the difference in the shapes of the laws of distribution of difference phases of adjacent samples of the noise spectrum and the monoharmonic signal;
 like any other quadratic characteristic (dispersion, PSD, correlation), the PS contrasts even the subtle difference in amplitude values.
The efficiency of signal noise cleaning allows us to estimate the spectra in FIG. 6a and 6b. In FIG. 6a shows the amplitude spectrum of a noisy signal
in FIG. 6b is a spectrum obtained by the proposed method from the same FFT samples.Even more striking is the difference in the spectra on a logarithmic scale. In FIG. 7 shows the superposition of the amplitude spectrum of the monoharmonic signal weighted by the Hanning window and the PS spectrum.
Claims (1)
 The method of spectral analysis of polyharmonic signals, namely, that the Fourier transform is performed on the sample of the analyzed signal, the obtained real and imaginary parts of the complex spectrum samples are multiplied by their copies and the products are summarized, characterized in that the original signal sample is supplemented with zero intervals in front and behind, and the real and imaginary parts of the spectrum samples are multiplied by the same parts of the neighboring sample, and after adding the product vectors, the total vector is multiplied by minus one in and zero all samples less than zero.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

RU2015153710A RU2611102C1 (en)  20151214  20151214  Method for spectral analysis of polyharmonic signals 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

RU2015153710A RU2611102C1 (en)  20151214  20151214  Method for spectral analysis of polyharmonic signals 
Publications (1)
Publication Number  Publication Date 

RU2611102C1 true RU2611102C1 (en)  20170221 
Family
ID=58458876
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

RU2015153710A RU2611102C1 (en)  20151214  20151214  Method for spectral analysis of polyharmonic signals 
Country Status (1)
Country  Link 

RU (1)  RU2611102C1 (en) 
Cited By (2)
Publication number  Priority date  Publication date  Assignee  Title 

RU2685972C1 (en) *  20171201  20190423  Акционерное общество "Всероссийский научноисследовательский институт радиотехники"  Method and device for filtering frequencymodulated signals 
RU2702917C1 (en) *  20190311  20191014  Федеральное государственное бюджетное учреждение науки Институт проблем морских технологий Дальневосточного отделения Российской академии наук (ИПМТ ДВО РАН)  Method of detecting noisy objects in shallow and deep sea 
Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

SU1256044A1 (en) *  19850219  19860907  Рязанский Радиотехнический Институт  Digital spectrum analyzer 
SU1691852A1 (en) *  19861205  19911115  Предприятие П/Я В8251  The spectrum analyzer 
US6882947B2 (en) *  20011231  20050419  Teradyne, Inc.  Discrete fourier transform (DFT) leakage removal 
RU2363005C1 (en) *  20080125  20090727  Открытое акционерное общество ордена Трудового Красного Знамени Всесоюзный научноисследовательский проектноконструкторский институт "Тяжпромэлектропроект" имени Ф.Б. Якубовского  Method of spectral analysis of polyharmonic signals and device to this end 
US8321491B2 (en) *  20070626  20121127  The United States Of America As Represented By The Secretary Of The Army  System and method for detecting a weak signal in a noisy environment 

2015
 20151214 RU RU2015153710A patent/RU2611102C1/en active
Patent Citations (5)
Publication number  Priority date  Publication date  Assignee  Title 

SU1256044A1 (en) *  19850219  19860907  Рязанский Радиотехнический Институт  Digital spectrum analyzer 
SU1691852A1 (en) *  19861205  19911115  Предприятие П/Я В8251  The spectrum analyzer 
US6882947B2 (en) *  20011231  20050419  Teradyne, Inc.  Discrete fourier transform (DFT) leakage removal 
US8321491B2 (en) *  20070626  20121127  The United States Of America As Represented By The Secretary Of The Army  System and method for detecting a weak signal in a noisy environment 
RU2363005C1 (en) *  20080125  20090727  Открытое акционерное общество ордена Трудового Красного Знамени Всесоюзный научноисследовательский проектноконструкторский институт "Тяжпромэлектропроект" имени Ф.Б. Якубовского  Method of spectral analysis of polyharmonic signals and device to this end 
Cited By (2)
Publication number  Priority date  Publication date  Assignee  Title 

RU2685972C1 (en) *  20171201  20190423  Акционерное общество "Всероссийский научноисследовательский институт радиотехники"  Method and device for filtering frequencymodulated signals 
RU2702917C1 (en) *  20190311  20191014  Федеральное государственное бюджетное учреждение науки Институт проблем морских технологий Дальневосточного отделения Российской академии наук (ИПМТ ДВО РАН)  Method of detecting noisy objects in shallow and deep sea 
Similar Documents
Publication  Publication Date  Title 

Howard et al.  Autoregressive modeling of wideband indoor radio propagation  
Tepedelenlioglu et al.  The Ricean K factor: estimation and performance analysis  
Irahhauten et al.  An overview of ultra wide band indoor channel measurements and modeling  
Wilkinson et al.  Steppedfrequency processing by reconstruction of target reflectivity spectrum  
Grandke  Interpolation algorithms for discrete Fourier transforms of weighted signals  
Benesty et al.  A generalized MVDR spectrum  
US8223989B2 (en)  Signal processing system and method for calibrating channel signals supplied from an array of sensors having different operating characteristics  
Shi et al.  A novel fractional wavelet transform and its applications  
US20140126624A1 (en)  LowAmplitude Echo Estimation for a Received Communication Signal  
Napolitano  Cyclostationarity: New trends and applications  
US7088765B1 (en)  Vector calibration system  
Turin  On the estimation in the presence of noise of the impulse response of a random, linear filter  
Walden  A unified view of multitaper multivariate spectral estimation  
Shen et al.  Lowcomplexity directionofarrival estimation based on wideband coprime arrays  
EP0986933A1 (en)  Continuous frequency dynamic range audio compressor  
Hewitt et al.  An autoregressive approach to the identification of multipath ray parameters from field measurements  
Baransky et al.  SubNyquist radar prototype: Hardware and algorithm  
Jakobsson et al.  Computationally efficient twodimensional Capon spectrum analysis  
Li et al.  Convolutional compressed sensing using deterministic sequences  
Myers et al.  A fast gridded method for the estimation of the power spectrum of the cosmic microwave background from interferometer data with application to the cosmic background imager  
Glentis et al.  Efficient implementation of iterative adaptive approach spectral estimation techniques  
US9118401B1 (en)  Method of adaptive interference mitigation in wide band spectrum  
US5874916A (en)  Frequency selective TDOA/FDOA crosscorrelation  
US9160310B2 (en)  Linearity compensator for removing nonlinear distortion  
Thomson  Jackknifing multitaper spectrum estimates 