Next Article in Journal
Shaping Perpendicular Magnetic Anisotropy of Co2MnGa Heusler Alloy Using Ion Irradiation for Magnetic Sensor Applications
Previous Article in Journal
Grad-MobileNet: A Gradient-Based Unsupervised Learning Method for Laser Welding Surface Defect Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Wide-Band Directional MUSIC Algorithm Using the Strength Proportion

1
Department of Automation, Shanghai Jiao Tong University, No. 800 Road Dongchuan, Shanghai 200240, China
2
Department of Electrical Engineering, Shanghai Jiao Tong University, No. 800 Road Dongchuan, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(9), 4562; https://doi.org/10.3390/s23094562
Submission received: 20 February 2023 / Revised: 29 April 2023 / Accepted: 3 May 2023 / Published: 8 May 2023
(This article belongs to the Section Internet of Things)

Abstract

:
The directional multiple signal classification (Dir-MUSIC) algorithm based on the antenna gain array manifold has been proposed to find the direction of the partial discharge (PD) source in substations. However, PD signals are wideband signals and the antenna gain pattern functions are always different at different frequencies; therefore, the accuracy can be improved using a wideband Dir-MUSIC algorithm. In this paper, wideband Dir-MUSIC algorithms are discussed and a novel wideband Dir-MUSIC algorithm using the strength proportion (DirSP) is proposed. This algorithm estimates a focusing PD signal at a certain frequency using the strength proportion among different directions, and then the Dir-MUSIC algorithm can process the focusing PD signal at this frequency. In simulations, when the antenna gain functions among different frequency bins are quite different, the Dir-MUSIC algorithm loses accuracy; meanwhile, DirDP performs very well. In the experiments, we deal with six sets of samples, and the mean error and standard deviation are both smaller than 4 ° better than other methods.

1. Introduction

Partial discharge (PD) detection and location systems are indispensable for performing quality assurance and fault identification in high-voltage apparatus. PD can be detected and located through PD signals (electromagnetic waves), especially the ultra-high-frequency (UHF) component of PD signals [1,2,3]. In recent years, PD detection and location methods based on UHF signals have been popular, with the advantages of strong anti-interference, high sensitivity and a stable transmitting speed. Thus, UHF detection methods, especially those using an omnidirectional UHF sensor array to detect and locate the PD signals in substations, have obtained fruitful results [4,5,6,7,8,9].
In some situations, the omnidirectional UHF sensor array is not a good choice because it need a small and light sensor array. For example, the intelligent inspection robot used for the substation inspection system can only carry a small sensor array [10,11]. However, PD detection methods using the time difference of arrival (TDOA) algorithm based on an omnidirectional UHF sensor array require a relatively large size sensor array to obtain an accurate TDOA. Therefore, a directional antenna array is a new choice in this situation. As a tentative research direction, the directional multiple signal classification (Dir-MUSIC) algorithm for direction of arrival (DOA) estimation based on signal strength, as represented by antenna gain array manifold, is proposed [12]. A miniaturized spiral directional antenna array has been designed and developed to receive the PD signal. Using the signal strength information, the Dir-MUSIC algorithm can successfully and precisely find the direction of the PD source. However, in the Dir-MUSIC algorithm, the antenna gain is assumed to be approximatively equal at different frequencies of interest. In fact, when the antenna gain at different frequencies is too different, the direction error will be unacceptable. Therefore, it is necessary to improve a wideband Dir-MUSIC algorithm.
It is difficult to deal with the wideband signals directly using a subspace DOA estimation algorithm. Most wideband subspace methods decompose the sensor outputs into several narrowband signals using a filter bank or the discrete Fourier transform (DFT). Normally, they nclude the incoherent signal subspace method (ISSM) [13] and the coherent signal subspace method (CSSM) [14]. ISSM is one of the simplest wideband methods. It process the decomposed narrowband signals independently, applying narrowband techniques, and the results of all the frequency bins are averaged to obtain the final DOA estimation. CSSM calculates a focusing matrix between different frequency bins and estimates a coherently averaged sample covariance matrix, which can be obtained by applying narrowband techniques. CSSM requires initial values to find the focusing matrices, and the estimation performance of CSSM is sensitive to these initial values.
Wideband subspace techniques have many applications. An efficient implementation of CSSM on the field-programmable gate array (FPGA) has been proposed to achieve the desired performance [15]. A generic subspace model has been proposed for characterizing a variety of radio frequency interference (RFI) types, which reveals a low-dimensional structure of the RFI subspace [16]. For PD location in substations in oil, a PD location method based on an ultrasonic phased array and wideband array signal processing is proposed [17].
In this paper, based on the Dir-MUSIC algorithm, we propose a novel wideband Dir-MUSIC algorithm using the strength proportion (DirSP), which can be regarded as another form of CSSM. In the proposed algorithm, PD signals will be decomposed into several narrowband signals using the DFT. We will choose some frequency bins to estimate an averaged sample covariance matrix at a certain frequency, and we will first estimate a focusing PD signal at this frequency. Then, we could apply the Dir-MUSIC algorithm on the averaged sample covariance matrix and the focusing PD signal.

2. Wide-Band Dir-MUSIC Algorithms

A PD signal is always a microsecond pulse signal, and the time interval between two pulses occurs randomly at the interval 50 us, 1 ms. Therefore, the probability of multi-source signal pulse superposition is very small, and this paper can only consider a single PD source in one time. We consider a uniform circular array of N elements (directional antenna), and the opening direction of each element is toward outside. The bandwidths of PD signals need not be identical, but there should be some frequency band [ ω L , ω H ] where PD signals’ bandwidths overlap.
Most wideband subspace methods decompose the sensor outputs into several narrowband signals using a filter bank or the DFT. Generally, PD signals are always measured at some frequency band of interest, such as the designed resonance points of antennas. If [ ω L , ω H ] contains all the frequency bands of interest, then the output of the filter bank or DFT module can be written in vector form at K frequencies, as follows:
X ( ω i ) = G ( ω i , θ ) S ( ω i ) + N ( ω i ) , i = 1 , , K ,
where
X ( ω i ) = [ X 1 ( ω i ) , , X N ( ω i ) ] T ,
where X j ( ω i ) is the output of the jth antenna at frequency ω i for j = 1 , , N , i = 1 , , K and ω L < ω i < ω H . S ( ω i ) is the component of the PD source and G ( ω i , θ ) is the antenna gain matrix at frequency ω i ( G ( ω i , θ ) is a vector in this paper):
G ( ω i , θ ) = [ g 1 ( ω i , θ 1 ) , , g N ( ω i , θ N ) ] T .
where g j ( ω i , θ j ) is the jth antenna gain when the direction angle of the PD source is θ j at frequency ω i , and θ j + 1 = θ j + 360 N holds.

2.1. The Normal Incoherent and Coherent Method

Next, a normal incoherent method is introduced. First, find the direction using the Dir-MUSIC algorithm at frequency ω i , as follows:
θ i = Dir-MUSIC ( X ( ω i ) , G ( ω i , θ ) ) .
Then take an arithmetic mean of θ i as the estimated direction:
θ e s t = i = 1 K θ i .
In the coherent method, an estimation of the coherently averaged sample covariance matrix is indispensable. Focusing matrices are effective in this operation, and can transform the constituent narrow-band components of the received signal into appropriate narrow-band representations at a certain given frequency. Next, we will introduce the rotational signal subspace (RSS) focusing matrix T ( ω i ) , which satisfies the following constrained minimization problem:
min T ( ω i , θ ) G ( ω 1 , θ ) T ( ω i , θ ) G ( ω i , θ ) F , i = 1 , , K ,
subject to
T T ( ω i , θ ) T ( ω i , θ ) = I ,
where F is the Frobenius matrix norm, and ω 1 is the centre frequency. One solution to (6) is
T ( ω i , θ ) = V ( ω i , θ ) U ( ω 1 , θ ) T ,
where the columns of U ( ω 1 , θ ) and V ( ω i , θ ) are the left and right singular vectors of G ( ω i , θ ) G ( ω 1 , θ ) T .
Then, the coherently averaged sample covariance matrix at ω 1 can be estimated by
R 1 = ( i = 1 K T ( ω i ) X ( ω i ) ) ( i = 1 K T ( ω i ) X ( ω i ) ) T .
It must be noted that directly applying the RSS focussing matrix is not a good choice. G ( ω 1 , θ ) and G ( ω i , θ ) in (6) always have different modules, so T can not have a theoretical solution that satisfies G ( ω 1 , θ ) T ( ω i , θ ) G ( ω i , θ ) F = 0 . In the traditional CSSM, the steering vectors have the same modules, and this characteristic ensures the effectiveness of the traditional CSSM. Therefore, when the antenna gain vectors have quite different modules, CSSM in this paper will lose accuracy.
As an alternative, a wideband Dir-MUSIC algorithm using the strength proportion is proposed in the next subsection.

2.2. The Wide-Band Dir-MUSIC Algorithm Using the Strength Proportion

In this paper, g j ( ω i , θ ) is only a real number when ω i , θ is given. Then, from (1) and (3), we can obtain
X j ( ω 1 ) = g j ( ω 1 , θ ) S ( ω 1 ) + N j ( ω 1 ) X j ( ω 2 ) = g j ( ω 2 , θ ) S ( ω 2 ) + N j ( ω 2 ) X j ( ω K ) = g j ( ω K , θ ) S ( ω K ) + N j ( ω K ) j = 1 , , N .
Let ω 1 be the reference frequency; (9) could be written as
X j ( ω 1 ) = g j ( ω 1 , θ ) S ( ω 1 ) + N j ( ω 1 ) g j ( ω 1 , θ ) g j ( ω 2 , θ ) X j ( ω 2 ) = g j ( ω 1 , θ ) S ( ω 2 ) + g j ( ω 1 , θ ) g j ( ω 2 , θ ) N j ( ω 2 ) g j ( ω 1 , θ ) g j ( ω K , θ ) X j ( ω K ) = g j ( ω 1 , θ ) S ( ω K ) + g j ( ω 1 , θ ) g j ( ω K , θ ) N j ( ω K ) j = 1 , , N .
Let X ˜ j be the sum of the left side in (10):
X ˜ j = i = 1 K g j ( ω 1 , θ ) g j ( ω i , θ ) X j ( ω i ) = g j ( ω 1 , θ ) S + i = 1 K g j ( ω 1 , θ ) g j ( ω i , θ ) N j ( ω i ) .
Then, the coherently averaged sample covariance matrix at ω 1 can be estimated by
R 1 = X ˜ X ˜ T ,
where X ˜ = [ X ˜ 1 , , X ˜ N ] T .
Next, we will explain that this method is equivalent to CSSM when it is applied in the normal array signal processing.
We consider a single source and a linear array, and the incident angle is θ . Then (9) and (10) convert to
X j ( ω 1 ) = e i ω 1 d j sin θ S ( ω 1 ) + N j ( ω 1 ) X j ( ω 2 ) = e i ω 2 d j sin θ S ( ω 2 ) + N j ( ω 2 ) X j ( ω K ) = e i ω K d j sin θ S ( ω K ) + N j ( ω K ) j = 1 , , N .
X j ( ω 1 ) = e i ω 1 d j sin θ S ( ω 1 ) + N j ( ω 1 ) e i ( ω 1 ω 2 ) d j sin θ X j ( ω 2 ) = e i ω 1 d j sin θ S ( ω 2 ) + e i ( ω 1 ω 2 ) d j sin θ N j ( ω 2 ) e i ( ω 1 ω N ) d j sin θ X j ( ω 2 ) = e i ω 1 d j sin θ S ( ω 2 ) + e i ( ω 1 ω N ) d j sin θ N j ( ω N ) j = 1 , , N .
Let
T ( ω i ) = d i a g { 1 , e i ( ω 1 ω 2 ) d j sin θ , , e i ( ω 1 ω N ) d j sin θ } .
The coherently averaged sample covariance matrix of this method and the CSSM are identical. Therefore, the proposed method is an application of CSSM in this paper.
Now, we can estimate the direction of the PD source using the Dir-MUSIC algorithm, and the steps are as follows.
We calculate the eigenvalues and eigenvectors of R 1 via eigenvalue decomposition. It is clear that R 1 is a real symmetric matrix, and the eigenvalues are all real numbers. S p a n { R 1 } also consists of signal subspace S p a n { R S } and noise subspace S p a n { R N } :
R 1 = G 1 R S G 1 T + R N ,
where G 1 is a simple expression of G ( ω 1 , θ ) . Therefore, R 1 is a positive definite. Since λ 1 λ 2 λ N are the eigenvalues, and ν 1 , ν 2 , , ν N are the related corresponding eigenvectors. Additionally, the eigenvectors are orthogonal to each other; that is,
ν i 1 T ν i 2 = 0 , 1 i 1 i 2 N .
This paper considers only one signal at one time; therefore, there is only one eigenvalue related to the signal. Certainly, λ 1 is the biggest eigenvalue related to the signal, and suppose λ i is one of the smaller eigenvalues. Then, we have
R 1 ν i = λ i ν i ( G 1 R S G 1 T + R N ) ν i = λ i ν i ( G 1 R S G 1 T + σ 2 I ) ν i = λ i ν i G 1 R S G 1 T ν i = 1 .
Because R S is positive definite, we have
G 1 T ν i = 0 , i = 2 , , N .
From (3), we can see that G 1 is also a vector function of θ ; therefore, we can search the direction of the PD source by traversing the possible value of θ . Let N 1 eigenvectors be related to the smaller N 1 eigenvalue to construct a noise matrix:
E n = [ ν 2 , ν 3 , , ν N ] .
Let the space spectrum be:
P m u ( θ ) = 1 g 1 T ( θ ) E n E n T g 1 ( θ ) = 1 E n T g 1 ( θ ) ,
where g 1 ( θ ) represents G 1 .
When we take the maximum value of P m u ( θ ) , it indicates that E n T g 1 ( θ ) is closest to 0 in all θ ; then, we decide this θ is the estimated direction angle.

3. Simulations

In this section, simulations will be carried out to verify the advantages of the wideband Dir-MUSIC algorithm. A relatively simple but effective wideband PD signal is designed for simulations.
A double exponential oscillation attenuation function is used to simulate the PD signal [18], and the specific expression is
f ( t ) = ( exp ( 1.3 t k 1 ) exp ( 2.2 t k 2 ) ) s i n ( 2 π t f ) ,
where the parameters can be k 1 = 1.5 × 10 8 , k 2 = 2.0 × 10 8 , f = 1.0 × 10 8 . This is a narrow-band signal in the approximate frequency band [ 90 M , 110 M ] . Consider other sets of parameters k 1 = 1.5 × 10 8 , k 2 = 2.0 × 10 8 , f = 2.0 × 10 8 and k 1 = 1.5 × 10 8 , k 2 = 2.0 × 10 8 , f = 3.0 × 10 8 ; the related simulated PD signals are written as f 2 ( t ) and f 3 ( t ) , and their approximate frequency bands are [ 190 M , 210 M ] and [ 290 M , 310 M ] . Meanwhile, the simulated PD signal related to the first set of parameters is written as f 1 ( t ) . Then, a wide-band PD signal is simulated as
f ( t ) = 1 3 ( f 1 ( t ) + f 2 ( t ) + f 3 ( t ) ) .
The sampled pulse signal of f ( t ) without noise at the sample rate 1 × e 9 is shown in Figure 1, and the frequency spectrum obtained from DFT is shown in Figure 2.
From Figure 1 and Figure 2, this simulated PD signal is a wideband signal that contains three bands. ω 1 = 100 M , ω 2 = 200 M , ω 3 = 300 M are the centre frequencies of the three bands. When the antenna gain is different among ω 1 , ω 2 and ω 3 , the Dir-MUSIC algorithm will lose its effectiveness.
Based on the actual antenna pattern, a linear combination of three Gaussian functions is enough to describe the two-dimensional antenna pattern function, as shown in
g ( θ ) = a 1 e ( θ b 1 ) 2 c 1 2 + a 2 e ( θ b 2 ) 2 c 2 2 + a 3 e ( θ b 3 ) 2 c 3 2 .
Suppose the set of parameters at ω 1 is a 1 = 0.5255 , b 1 = 218.1 , c 1 = 51.73 ,   a 2 = 0.3405 , b 2 = 304.5 , c 2 = 41 , a 3 = 0.6251 , b 3 = 156.1 , c 3 = 109.1 , and the other two set sof parameters are different at c 1 = 101.73 and c 1 = 251.73 . The normalized antenna gains of the three bands are shown in Figure 3.
We provide simulation steps in Figure 4.
Let the direction of the PD source θ change from 91 ° to 270 ° , and SNR is set at 5, 0, −5, and the direction results of DirSP compared to Dir-MUSIC are shown in Figure 5, Figure 6 and Figure 7.
From Figure 5, Figure 6 and Figure 7, we can see when the antenna pattern functions among different frequency bins are quite different; the Dir-MUSIC algorithm only using one function will lose accuracy. Meanwhile, the DirSP algorithm using three functions performs very well.
Let the direction of the PD source θ change from 91 ° to 270 ° ; SNR is set at 5, 0, 5 , and the direction results of the wide-band methods are shown in Figure 8, Figure 9 and Figure 10. Additionally, each result in the following figures is a mean of 50 random estimations.
From Figure 8, Figure 9 and Figure 10, we can see that the errors of ISSM and DirSP are similar, but much less than CSSM. Specifically, the mean errors of ISSM, CSSM and DirSP are 1.8130 ° , 5.1066 ° and 2.1403 ° at SNR = 5 , 0.9914 ° , 2.7778 ° and 1.1081 ° at SNR = 0 , 0.5532 ° , 1.4948 ° and 0.5677 ° at SNR = 5 . This is ineffective for CSSM when directly applied in the wide-band Dir-MUSIC algorithm. On the other hand, errors in ISSM and DirSP are acceptable even when SNR is 5 .
It seems that ISSM is a little better than the proposed DirSP. The reason for this is that measurements are generated from the ideal model perfectly, and every Dir-MUSIC estimation is a perfect estimation of ISSM. Additionally, in DirSP, the strength proportion in (10) is always not equal to 1, and the noise distribution changes in X ˜ in (12). Therefore, ISSM performs best in the ideal model. However, the relationship between the actual measurements and the measured antenna pattern are not as ideal. Moreover, the proportion of f 1 , f 2 and f 3 is always not equal in (23); having the wrong proportion will cause mistakes in direction estimation. Next, we simulate a new PD signal as
f ( t ) = 1 3 ( f 1 ( t ) + 2 f 2 ( t ) + 3 f 3 ( t ) ) .
Then, we repeat the last simulation with the new signal, and the direction results are shown in Figure 11, Figure 12 and Figure 13.
From Figure 11, Figure 12 and Figure 13, we can see that ISSM errors become higher when the proportion of frequency bins is wrong. Specifically, the mean errors of ISSM, CSSM and DirSP are 3.4939 ° , 3.7559 ° and 2.6799 ° at SNR = 5 , 1.6623 ° , 2.1102 ° and 1.3461 ° at SNR = 0 , 0.8794 ° , 1.1434 ° and 0.7001 ° at SNR = 5 . Another method could possibly be used to estimate the proportion; however, this will increase computational complexity and the accuracy will rely on the proportion estimation. Fortunately, errors in the new method are not influenced with changes in proportion and the errors of DirSP are quite acceptable.
At last, we should analyze why CSSM performs badly in the simulations. ω 0 and ω i in (6) always have different modules, so T could not have a theoretical solution. In the traditional CSSM, the steering vectors have the same modules, and this characteristic ensures the effectiveness of the traditional CSSM. Therefore, when the antenna gain vectors have quite different modules, CSSM in this paper will lose accuracy.

4. Experimental Data Processing

Experiments were carried out to verify the effectiveness of DirSP. The experimental platform consisted of a digital oscilloscope with storage function, four equal-length radio frequency coaxial cables, the developed uniform fan-shaped directional Vivaldi antenna array, and a lighter as a signal generator. Figure 14 and Figure 15 are two photos of the experiment environment.
To achieve a good direction, high gain and high angular resolution requirements, we used the developed uniform fan-shaped directional Vivaldi antenna array to receive signals. These four antennas were denoted from Antenna 1 to Antenna 4 in counter-clockwise order. The angle between the symmetry axes of Antenna 1 and Antenna 4 was 60 ° , and the angle between the symmetry axes of adjacent antennas was 20 ° . The location of the experimental site was calibrated by infrared range finders, and PD signals were generated by the lighter. We provide an example of the PD signal received from the lighter.
In Figure 16, the signal contains 50 pulses. We used a high sample rate of 10 GHz to verify our algorithm. We calculated the frequency spectrum of a single pulse using DFT, as shown in Figure 17. Normally, interference is stronger at low frequencies than at high frequencies. In this paper, we chose the frequency band 1–2 GHz for performance verification. Whether the frequency band 0.5–1 GHz can be used for the algorithms could be studied in subsequent research.
From Figure 18, we chose three frequency bands, 1.3–1.5 GHz, 1.5–1.7 GHz and 1.7–1.9 GHz, to run the wideband Dir-MUSIC algorithm. We measured the antenna patterns for all four antennas every 50 M in the band 1–2 GHz. Because it is a fan-shaped array, we only measured the antenna pattern every 0.5 ° , from 90 ° to 90 ° , and the 0 ° direction was the direction of the symmetry axis. The measured normalized patterns at 1.4 GHz, 1.6 GHz and 1.8 GHz of all four antennas are shown in Figure 18 and Figure 19.
Next, we will process the experimental data with the wideband Dir-MUSIC algorithm, and the mean error and standard deviation of the direction results will be shown in Table 1. Additionally, it should be noted that the location area is the fan-shaped area, which means that the true azimuth is in [ 0 ° , 60 ° ] .
In Table 1, the first column is the PD coordinate (distance and azimuth). At each position, angles in the first row are mean errors (between the average direction result and the true azimuth) and angles in the second row are standard deviation. The mean errors of DirSP are smaller than other methods in most situations. Only at (7.8 m, 30 ° ) are the mean errors of ISSM and CSSM smaller than DirSP. The standard deviation of DirSP is also smaller than other methods in most situations. More importantly, the mean error and standard deviation are both smaller than 4 ° . Bad results can be obtained using other methods. In general, the experimental data-processing results are consistent with simulations and DirSP performs better than other methods.

5. Conclusions

In this paper, wideband Dir-MUSIC algorithms to find PD directions are discussed, and a novel wideband Dir-MUSIC algorithm using the strength proportion is proposed. Conclusions are drawn as follows.
When the antenna gain pattern functions among different frequency bins are quite different, the Dir-MUSIC algorithm only using one function will lose accuracy. Meanwhile, DirSP using three functions performs very well.
It is non-effective for CSSM to be directly applied to the Dir-MUSIC algorithm because the steer vectors have different modules, and DirSP is another application of CSSM in this paper.
ISSM may perform better when the signals are generated in an ideal model. However, DirSP performs better when the strength proportion among the frequency bins is unknown. DirSP is widely applicable in most situations.
Wideband Dir-MUSIC algorithms can avoid the influence of interference signals.
Finally, experimental data-processing has verified the effectiveness of DirSP. Compared to Dir-MUSIC, ISSM and CSSM, DirSP performs better in terms of both mean errors and standard deviation. More importantly, DirSP is more stable than other methods, while some bad results are obtained using other methods. Therefore, DirSP is effective for PD location.

Author Contributions

Conceptualization, W.X. and Y.H.; methodology, W.X.; software, W.X.; validation, B.C., Y.H. and J.L.; formal analysis, W.X.; investigation, B.C.; resources, Y.H.; data curation, W.X. and B.C.; writing—original draft preparation, W.X.; writing—review and editing, W.X.; visualization, W.X.; supervision, Y.H. and J.L.; project administration, Y.H. and J.L.; funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

NSFC, Grant/Award Number: 51777122.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All data, models, or code generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Robles, G.; Fresno, J.M.; Sánchez-Fernández, M.; Martínez-Tarifa, J.M. Antenna Deployment for the Localization of Partial Discharges in Open-Air Substations. Sensors 2016, 16, 541. [Google Scholar] [CrossRef] [PubMed]
  2. Li, J.; Zhang, X.; Han, X.; Yao, X. A Partial Discharge Detection Method for SF6 Insulated Inverted Current Transformers Adopting Inner Shield Case as UHF Sensor. IEEE Trans. Power Deliv. 2018, 33, 3237–3239. [Google Scholar] [CrossRef]
  3. Li, X.; Tang, H.; Mu, S.; Song, K.; Liu, R.; Xu, G.; Li, Q. Partial discharge monitoring system for PD characteristics of typical defects in GIS using UHF method. In Proceedings of the 2012 International Conference on High Voltage Engineering and Application, Shanghai, China, 17–20 September 2012; pp. 625–628. [Google Scholar]
  4. Hou, H.; Sheng, G.; Miao, P. Spatial location of partial discharge in Substation Based on UHF electromagnetic wave. High Volt. Technol. 2012, 38, 1334–1340. [Google Scholar]
  5. Hu, Y.; Wang, H.; Tang, L.; Wang, J.; Xu, W.; Jiang, X. Research and practice of partial discharge monitoring in Substation Based on UHF antenna array. Trans. China Electromech. 2014, 29, 502–509. [Google Scholar]
  6. Zhang, G.; Zhu, M.; Wang, Y. Detection and Localization of Discharge Sources in Substation Area Using Movable UHF Antenna Array. Proc. CSEE 2017, 37, 2761–2773. [Google Scholar]
  7. Tang, L.; Hu, Y.; Wang, H. Time delay estimation of partial discharge UHF pulse signal based on interpolation correlation method. High Volt. Technol. 2015, 41, 3320–3325. [Google Scholar]
  8. Li, Z.; Luo, L.; Sheng, G.; Liu, Y.; Jiang, X. UHF partial discharge location method in substation based on dimension-reduced RSSI fingerprint. IET Gener. Transm. Distrib. 2018, 12, 398–405. [Google Scholar] [CrossRef]
  9. Liu, Q.; Zhu, M.X.; Wang, Y.B.; Deng, J.B.; Li, Y.; Zhang, G.J.; Shao, X.J.; Zhao, X.F. UHF antenna array arrangement optimization for partial discharge direction finding in air-insulted substation based on phased array theory. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 3657–3668. [Google Scholar] [CrossRef]
  10. Li, B.; Quan, Q.; Xie, B. Application of infrared thermal image monitoring robot in substation inspection. Electron. Des. Eng. 2020, 28, 175–179. [Google Scholar]
  11. Zhao, H.; Qin, J.; Xie, Y. Application of substation intelligent inspection robot in operation and maintenance. Electr. Energy 2020, 41, 579–581. [Google Scholar]
  12. Xu, W.; Chen, B.; Li, Y.; Hu, Y.; Li, J.; Zeng, Z. Dir-MUSIC Algorithm for DOA Estimation of Partial Discharge Based on Signal Strength Represented by Antenna Gain Array Manifold. Sensors 2022, 22, 5406. [Google Scholar] [CrossRef] [PubMed]
  13. Wax, M.; Shan, T.; Kailath, T. Spatio-temporal spectral analysis by eigenstructure methods. IEEE Trans. Acoust. Speech Signal Process. 1984, 32, 817–827. [Google Scholar] [CrossRef]
  14. Hung, H.; Kaveh, M. Focussing matrices for coherent signal-subspace processing. IEEE Trans. Acoust. Speech Signal Process. 1988, 36, 1272–1281. [Google Scholar] [CrossRef]
  15. Almomany, A.M.; Jarrah, A.; Alsobeh, A.; Alqudah, E. High-Performance Implementation of Wideband Coherent Signal-Subspace (CSS) based DOA Algorithm on FPGA. J. Circuits Syst. Comput. 2011, 30, 1–24. [Google Scholar]
  16. Yang, H.; Li, K.; Li, J.; Du, Y.; Yang, J. BSF: Block Subspace Filter for Removing Narrowband and Wideband Radio Interference Artifacts in Single-Look Complex SAR Images. IEEE Trans. Geosci. Remote. Sens. 2022, 30, 1–16. [Google Scholar] [CrossRef]
  17. Xie, Q.; Li, Y.Q.; Lü, F.C.; Li, C.R.; Wang, N.; Ding, Y.J. Method for PD location in oil combining ultrasonic phased array with wideband array signal processing. Proc. Chin. Soc. Electr. Eng. 2009, 29, 13–17. [Google Scholar]
  18. Li, A.Q.; Wang, N.; Yi, D. Numerical Analysis, 4th ed.; Tsinghua University Press: Beijing, China, 2001; pp. 52–58. [Google Scholar]
Figure 1. The sampled pulse signal.
Figure 1. The sampled pulse signal.
Sensors 23 04562 g001
Figure 2. The frequency spectrum of the sampled pulse signal.
Figure 2. The frequency spectrum of the sampled pulse signal.
Sensors 23 04562 g002
Figure 3. The normalized antenna gain of three bands.
Figure 3. The normalized antenna gain of three bands.
Sensors 23 04562 g003
Figure 4. Simulation steps of DirSP.
Figure 4. Simulation steps of DirSP.
Sensors 23 04562 g004
Figure 5. Direction errors of DirSP compared to Dir-MUSIC at SNR = 5.
Figure 5. Direction errors of DirSP compared to Dir-MUSIC at SNR = 5.
Sensors 23 04562 g005
Figure 6. Direction errors of DirSP compared to Dir-MUSIC at SNR = 0.
Figure 6. Direction errors of DirSP compared to Dir-MUSIC at SNR = 0.
Sensors 23 04562 g006
Figure 7. Direction errors of DirSP compared to Dir-MUSIC at SNR = −5.
Figure 7. Direction errors of DirSP compared to Dir-MUSIC at SNR = −5.
Sensors 23 04562 g007
Figure 8. Direction errors of wide-band algorithms at SNR = 5 .
Figure 8. Direction errors of wide-band algorithms at SNR = 5 .
Sensors 23 04562 g008
Figure 9. Direction errors of wide-band algorithms at SNR = 0 .
Figure 9. Direction errors of wide-band algorithms at SNR = 0 .
Sensors 23 04562 g009
Figure 10. Direction errors of wide-band algorithms at SNR = 5 .
Figure 10. Direction errors of wide-band algorithms at SNR = 5 .
Sensors 23 04562 g010
Figure 11. Direction errors of wide-band algorithms at SNR = 5 with a new signal.
Figure 11. Direction errors of wide-band algorithms at SNR = 5 with a new signal.
Sensors 23 04562 g011
Figure 12. Direction errors of wide-band algorithms at SNR = 0 with a new signal.
Figure 12. Direction errors of wide-band algorithms at SNR = 0 with a new signal.
Sensors 23 04562 g012
Figure 13. Direction errors of wide-band algorithms at SNR = 5 with a new signal.
Figure 13. Direction errors of wide-band algorithms at SNR = 5 with a new signal.
Sensors 23 04562 g013
Figure 14. A photo of experiment environment.
Figure 14. A photo of experiment environment.
Sensors 23 04562 g014
Figure 15. A photo of experiment environment.
Figure 15. A photo of experiment environment.
Sensors 23 04562 g015
Figure 16. The received signal generated by the lighter.
Figure 16. The received signal generated by the lighter.
Sensors 23 04562 g016
Figure 17. The frequency spectrum of The received signal.
Figure 17. The frequency spectrum of The received signal.
Sensors 23 04562 g017
Figure 18. The measured antenna gain of Antenna 1 and Antenna 2.
Figure 18. The measured antenna gain of Antenna 1 and Antenna 2.
Sensors 23 04562 g018
Figure 19. The measured antenna gain of Antenna 3 and Antenna 4.
Figure 19. The measured antenna gain of Antenna 3 and Antenna 4.
Sensors 23 04562 g019
Table 1. The mean errors and standard deviation of the experimental data processing.
Table 1. The mean errors and standard deviation of the experimental data processing.
PD CoordinatesDir-MUSICISSMCSSMDirSP
(7.8 m, 43 ° ) 4.94 ° 6.89 ° 1.88 ° 0.39 °
9.10 ° 3.71 ° 9.94 ° 2.96 °
(7.8 m, 22 ° ) 9.11 ° 3.37 ° 9.65 ° 3.18 °
1.68 ° 5.28 ° 11.56 ° 2.85 °
(7.8 m, 30 ° ) 3.01 ° 0.11 ° 0.88 ° 2.01 °
1.52 ° 4.10 ° 4.80 ° 2.98 °
(14.4 m, 27.6 ° ) 2.98 ° 0.38 ° 1.20 ° 0.12 °
2.95 ° 2.64 ° 10.02 ° 2.45 °
(14.4 m, 39.5 ° ) 3.75 ° 3.00 ° 4.65 ° 1.29 °
1.63 ° 2.63 ° 8.33 ° 3.46 °
(6.6 m, 45.3 ° ) 5.37 ° 10.37 ° 1.68 ° 1.01 °
5.13 ° 5.12 ° 8.19 ° 2.69 °
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, W.; Chen, B.; Hu, Y.; Li, J. A Novel Wide-Band Directional MUSIC Algorithm Using the Strength Proportion. Sensors 2023, 23, 4562. https://doi.org/10.3390/s23094562

AMA Style

Xu W, Chen B, Hu Y, Li J. A Novel Wide-Band Directional MUSIC Algorithm Using the Strength Proportion. Sensors. 2023; 23(9):4562. https://doi.org/10.3390/s23094562

Chicago/Turabian Style

Xu, Wencong, Bingshu Chen, Yue Hu, and Jianxun Li. 2023. "A Novel Wide-Band Directional MUSIC Algorithm Using the Strength Proportion" Sensors 23, no. 9: 4562. https://doi.org/10.3390/s23094562

APA Style

Xu, W., Chen, B., Hu, Y., & Li, J. (2023). A Novel Wide-Band Directional MUSIC Algorithm Using the Strength Proportion. Sensors, 23(9), 4562. https://doi.org/10.3390/s23094562

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop