# FFT-based selective harmonic current control of AC filter capacitors for audible noise analysis

J. Han, S. Sun, S. Zheng, Y. Gao, Yide Wang, T. Tang

► To cite this version:

J. Han, S. Sun, S. Zheng, Y. Gao, Yide Wang, et al., FFT-based selective harmonic current control of AC filter capacitors for audible noise analysis. Measurement and Control, 2021, 54 (3-4), pp.408-416. 10.1177/00202940211000078. hal-03195867

HAL Id: hal-03195867

https://hal.science/hal-03195867v1

Submitted on 12 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# FFT-based selective harmonic current control of AC filter capacitors for audible noise analysis

Jingang Han<fnref n="1" /> ID, Shuoshuo Sun<fnref n="1" /> ID, Shouzhi Zheng<fnref n="1" />, Yang Gao<fnref n="1" />, Yide Wang<fnref n="1" />,<fnref n="2" /> and TianHao Tang<fnref n="1" />

Measurement and Control
2021, Vol. 54(3-4) 408-416

© The Author(s) 2021

Article reuse guidelines:
sagepub.com/journals-permissions
DOI: 10.1177/00202940211000078
journals.sagepub.com/home/mac

<figure><img src="image_2.png" /></figure>

SAGE

## Abstract

With the development of High Voltage Direct Current (HVDC) and ultra HVDC, a growing number of high voltage film filter capacitors have been installed for harmonic suppression. But the capacitors generate severe audible noise in the converters, which has become more and more serious. Then it is important to imitate the practical working situation of capacitors for the research of noise characteristics to improve the capacitive performance in lab. In this paper, a FFT-based harmonic current control method is proposed for audible noise analysis of AC filter capacitor. Firstly, mathematical model of the Harmonic Current Injection System (HCIS) is built based on MATLAB. Then a Fast Fourier Transform (FFT) harmonics detection method is analyzed and used to detect the amplitude and phase of each harmonic considering a number of harmonic currents flowing through AC filter capacitors. Since the amplitude and phase spectrum of capacitor current are obtained, the Proportional Integral (PI) regulator is applied to adjust the specified harmonic amplitude to achieve zero-error control. Then accurate harmonic current can be injected to the tested capacitor and it will help to acquire the noise characteristics of the tested plant. At last, the simulation and experimental results validate the proposed strategy.

## Keywords

Harmonic detection, AC filter capacitors, current control, FFT algorithm

Date received: 30 January 2021; accepted: 13 February 2021

## Introduction

With the increasing of energy demand and the development of new energy generation technology, HVDC transmission technology becomes one of the development directions of power grid in the future in order to supplement the existing AC transmission system. The technology has significant advantages in long-distance, high-volume power delivery, and interconnection of the modern power grids.<fnref n="1" /> Among them, converter station is the key part of HVDC transmission system, and a large number of low-order harmonics are generated. In order to compensate reactive power and filter harmonics, the converter station needs a large number of AC filter capacitors which will generate audible noise. With the application of HVDC technology, the problem of noise pollution in converter station has becomes more and more serious. The audible noise created by AC filter capacitors in converter station may be over 100 dB (A) when capacitor currents contain multiple harmonics.<fnref n="2" /> The audible noise of AC filter capacitors is mainly caused by harmonic or superposition of harmonics, in which 11th and 13th harmonic

currents produce the most serious noise. When harmonic currents flow through the internal of capacitor, it will cause the vibration of capacitor inner core, and finally, the vibration can transfer from the inside to capacitor cases to produce the noise.

In the last decade, many scholars have been focusing on the audible noise characteristics of capacitors and noise suppression. Li et al.³ and Cox and Guan⁴ have investigated the relationship between the audible noise and the vibration characteristics. The acoustic noise of AC filter capacitor is affected by the attributes of harmonic currents, such as the amplitude, relative phase, and the harmonic order.⁵

## Corresponding author:

Jingang Han, Institute of Electric Drives and Control Systems, Shanghai Maritime University, 1550 Haigang Avenue, Shanghai 201306, China.

Email: jghan@shmtu.edu.cn

<footnote n="1">¹Institute of Electric Drives and Control Systems, Shanghai Maritime University, Shanghai, China</footnote>

<footnote n="2">²IETR-UMR CNRS 6164, l'Université de Nantes/Polytech Nantes, Nantes, Pays de la Loire, France</footnote>
Unfortunately, due to the limitation of the working condition, it is impossible to analyze exactly the acous- tic noise characteristics on site. So, a test bench, which could simulate the practical work circumstance of the capacitors, is necessary in laboratory. Then the rela- tively accurate noise characteristic of capacitor can be acquired and the results can be used to improve the capacitor.

For the sake of the adequate implementation of this work, two main parts are vital important: (1) harmonic detection algorithm to get the specific attributes of har- monic currents (frequency, phase, and amplitude) and (2) harmonic current controller to regulate the ampli- tude to meet reference.

For now, numerous harmonic detection algorithms have been proposed. In general, those algorithms fall into two categories: time-domain methods and frequency-domain methods.

Time-domain harmonic detection methods are based on the instantaneous reactive power theory, such as *p-q* method<sup>6,7</sup> and *d-q* method.<sup>8,9</sup> In addition to the instan- taneous reactive power, the filter algorithms have to be used to separate the fundamental and harmonic com- ponents, such as Adaptive notch filter (ANF)<sup>10</sup> or Kalman filter.<sup>11</sup> There has a common feature about the time-domain harmonic detection algorithms: after a number of mathematical manipulations, the harmonic components are obtained by the use of the correspond- ing filter. To a great extent, the filter determines the performance of the harmonic detection algorithm. As a result, the time-domain harmonic detection method can compute one harmonic component rapidly. But the computation burden of this detection algorithm will be greatly increased with the increase of injected harmonic number. Moreover, multiple filters used in time-domain can also increase the computational burden of the sys- tem considering that each harmonic needs a unique filter.

The frequency-domain harmonic detection algo- rithm is more appropriate for this work compared to time-domain methods because the phase and amplitude information of all harmonics can acquired during one calculation procedure. The these methods originate from the Fourier transform, such as DFT,<sup>12</sup> FFT,<sup>13</sup> and SDFT.<sup>14</sup> FFT is an efficient way to implement DFT. And it can decompose an N-point DFT into two N/2-point DFT to accelerate the operational speed and reduce the operation burden.

In this paper, in addition to detect the harmonic cur- rents exactly, regulation of the harmonic current is also an important task. In time domain, with the conversion of specified harmonic into DC component, PI control- ler for each harmonic can be designed to realize zero steady-state error tracking control.<sup>15,16</sup> According to the Internal Model Principle,<sup>17</sup> the Proportional Resonant controllers can also be used to control the AC reference components. Besides the above linear

controllers, some nonlinear controllers have been applied to harmonic control, such as Model predictive control (MPC), Repetitive control (RC) and Sliding Mode control (SMC), etc.<sup>15,18-20</sup>

Since the amplitude spectrum extracted by FFT algorithm is DC component, the discrete PI regulator can also be used to track the reference amplitude with zero steady-state error. On the other hand, the system inherent delay would result in the phase deviation, caused by the digital sample delay and system impe- dance. And the phase deviations will be more serious with the increase of harmonic order.²¹ Thus, the phase deviations should be compensated separately for each harmonic to ensure the performance of the overall sys- tem. Once the amplitude of harmonics is detected and controlled, and the phase deviations are compensated, the reference signal can be reconstructed on the time field with Inverse Fast Fourier Transform (IFFT). Then the output reference is converted into analog sig- nal by digital to analog converter.

In this paper, FFT-based selective harmonic current control is proposed for research on audible noise char- acteristics of AC capacitors. Firstly, in section 2, the mathematical model of HCIS is presented. In section 3, a detailed presentation about the detection and control strategy for harmonic currents will be given. Finally, the simulation and experimental results of the HCIS based on FFT are shown in section 4 and 5, based on the MATLAB/Simulink and RT-LAB, respectively, to validate the significance of the proposed strategy.

## System configuration

### Harmonic currents injection system

The HCIS shown in Figure 1 is composed of the refer- ence signal generator, harmonic detection module and current controller, transformer, tested capacitor, and switch mode power amplifier. The fundamental and selected harmonic signals are all generated by the refer- ence signal generator with the specific parameters such

<figure><img src="image_3.png" /><figcaption>Figure 1. Harmonic currents injection system.</figcaption></figure>
<table><caption>**Table 1.** Harmonic currents distribution.</caption><thead><tr><th>Harmonic current order</th><th>Current/A</th></tr></thead><tbody><tr><td>1</td><td>74.99</td></tr><tr><td>5</td><td>8.625</td></tr><tr><td>7</td><td>28.55</td></tr><tr><td>9</td><td>5.55</td></tr><tr><td>11</td><td>9.2</td></tr><tr><td>13</td><td>91.3</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** The simplified capacitor model.</figcaption></figure>

as the frequency, amplitude, and phase angle. The switch mode power amplifier can amplify the input signal linearly. And the amplified harmonic voltage can be injected into the power capacitors by the transformer.

The capacitor harmonic currents in HVDC converter station vary greatly under different working conditions. One of the actual measured harmonic currents of power capacitor (6.2 kV/38.5 uF) in one converter station is shown in Table 1. It can be found that the 13th harmonic is 91.3 A which will cause serious audible noise.

### Model of the HCIS

The power amplifier with inner voltage control can be approximately considered as a first-order inertial link, and its transfer function is expressed as:

$$G_{AMP} = \frac{k_{AMP}}{1 + \tau_s} \qquad (1)$$

where $k_{AMP}$ is the magnification factor, $\tau$ is the time constant considering the switching frequency.

The capacitor is simplified as an ideal capacitor C with an equivalent resistor $r$ as illustrated in Figure 2.

The transfer function of the output current $I_C(s)$ to voltage $U_C(s)$ is given as:

$$\frac{I_C(s)}{U_C(s)} = \frac{Cs}{Cr_s + 1} \qquad (2)$$

## Harmonic currents detection and control strategy

### Harmonics detection algorithm

In this section, the principles and processes of FFT detection algorithm will be presented in detail. FFT algorithm is an efficient method for DFT calculation, and DFT is the discrete Fourier analysis with finite length.

If a periodic and continuous signal $x(t)$ satisfies the Dirichlet conditions, it can be represented by its FS:

$$x(t) = \sum_{k=-\infty}^{\infty} X(k)e^{ik\omega_0 t}, \omega_0 = 2\pi/T \qquad (3)$$

where $T$ is the period of the signal. In real control system, the signal both in time and frequency is discrete and has limited sample length, thus N-point DFT is realized with a sampling rate $f_s$, where $f_s = N/T$, and the corresponding DFT will be expressed as:

$$X[k] = \sum_{n=0}^{N-1} x[n] W_N^{nk} \quad k = 0, 1, \dots, N-1 \qquad (4)$$

where $W_N^{nk} = e^{-j\frac{2\pi nk}{N}}$ is called the twiddle factor, and $X[k]$ is the spectrum of $x(n)$. Also there exists an Inverse Discrete Fourier Transform (IDFT), which is used to convert $X[k]$, $k = 0, 1, \dots, N-1$ in frequency-domain back into $x[n]$, $n = 0, 1, \dots, N-1$ in time-domain. The IDFT is expressed as:

$$x[n] = \frac{1}{N} \sum_{k=0}^{N-1} X[k] e^{j(2\pi/N)kn} = 0, 1, \dots, N-1 \qquad (5)$$

Usually, an N-point DFT operation requires $N^2$ times multiplication and $N^2$ times addition.

By the aid of the mathematical characteristics of the twiddle factor, the FFT algorithm decomposes the N-point DFT into smaller DFTs to reduce the computational burden of DFT.

$$\left\{ \begin{array}{l} W_N^0 = 1, W_N^{N/2} = -1 \\ W_N^{km + (N/2)} = -W_N^{km}, (W_N^{nk})^* = W_N^{-nk} \\ W_N^{(k+N)m} = W_N^{k(m+N)} = W_N^{km}, W_N^{km} = W_{nN}^{nkm} \end{array} \right. \qquad (6)$$

As shown in Figure 3, every point means a complex addition and arrows mean complex multiplication. The N-point FFT will be transformed into two N/2-point DFTs according to recursive decomposition, then the N/2-point DFT will be transformed into the N/4 DFTs. This operation will be repeated, until the final 2-point DFT is extracted.

Generally, the decomposition number of an N-point FFT is $\log_2 N - 1$, where N is a power of two. And there are two ways to decompose the sequence, one is the Decimation-In-Time FFT (DIT-FFT), the other is the Decimation-In-Frequency FFT (DIF-FFT).

Take the DIT-FFT for example, $x[n]$ is the sum of even sequence $x_1[n]$ and odd sequence $x_2[n]$, then
<figure><img src="image_5.png" /><figcaption>Figure 3. Principle of 8-point FFT algorithm.</figcaption></figure>

$$
\begin{align*}
x[n] &= x_1[n] + x_2[n] \\
X[k] &= \sum_{n=0}^{\frac{N}{2}-1} x_1[n] W_N^{2kn} + \sum_{n=0}^{N/2} x_2[n] W_N^{k(2n+1)} \tag{7} \\
&= \sum_{n=0}^{\frac{N}{2}-1} x_1[n] W_{N/2}^{kn} + W_N^k \sum_{n=0}^{\frac{N}{2}-1} x_2[n] W_{N/2}^{kn} \\
&= X_1[k] + W_N^k X_2[k]
\end{align*}
$$

where $X_1[k]$ and $X_2[k]$ are respectively $N/2$-point DFT of $x_1[k]$ and $x_2[k]$, and then, (7) can be transformed as:

$$
\begin{align}
X[k] &= X_1[k] + W_N^k X_2[k] \nonumber \\
X[k + \frac{N}{2}] &= X_1[k] - W_N^k X_2[k] \tag{8} \\
k &= 0, 1, \dots, N-1 \nonumber
\end{align}
$$

As shown in Figure 4, this basic operation unit of
FFT is called butterfly operation.

All this leads up to the basic principle of FFT: a long
sequence of DFT can be transformed into two short
sequences of DFT by properties of DFT. The computa-
tional complexity of N-point FFT is O[Nlog₂N] com-
pared to O[N²] for N-point DFT. It will reduce the
mathematics calculation and accelerate the operational
speed.

In the rest of this paper, some features about the
spectrum will be shown. It is clear that the output ele-
ments of FFT would be limited to the number of N.
And the frequency resolution is Δf = f_s/N, where f_s is
sampling frequency. Then the Nyquist frequency of the
system equals to f_s/2, which is the maximum positive
frequency of the spectrum. A frame of N points must
embody at least a complete period of the sampled
signal.

<figure><img src="image_6.png" /><figcaption>Figure 4. Butterfly operation.</figcaption></figure>

For the N complex numbers obtained by FFT, the
frequency of the first item is f = 0, and the frequencies
of the remaining items will multiply by Δf, which can be
written as: Δf, 2Δf, ..., (N/2 - 1)Δf. Through the mathe-
matical computation, the information of each harmonic
will be extracted from the spectrum, and the Fourier
coefficients are scaled by N/2.

If an exact integer period of signals is included dur-
ing the sampling window at the signal acquisition stage,
the final results of FFT algorithm will be precise. The
well-known result of “harmonic leakage” is caused by
incompleteness within the data acquisition window,
and the harmonics “leakage” will result in the inaccu-
rate detection. Under these circumstances, some rele-
vant measures have been proposed to compensate the
impact of non-cycle sampling, such as: adding the win-
dow function or zero padding.

When the amplitude and phase angle of a specific
harmonic are extracted from the data frames, the corre-
sponding harmonic current regulator will be adopted
to ensure the accurateness of the amplitude, and the
related phase compensation angle is applied to ensure
the synchronization. When the specific harmonic cur-
rent has been controlled, the instruction signal can be
obtained again in time domain by IFFT or other recon-
figuration algorithms.
<figure><img src="image_7.png" /><figcaption>Figure 5. Configuration of the proposed HCIS.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Figure 6. Current closed-loop control model.</figcaption></figure>

In general, the number of input and output points of FFT algorithm is equal. To reduce the computational burden, some specific points such as X(3), X(5), and X(7) are selected for calculating the specified butterfly of FFT.

## Harmonic control method

The configuration of the proposed HCIS is shown in Figure 5. The tested current measured by hall sensor is converted to digital data with the help of AD converter. Then the magnitude and phase of a specific harmonic current are detected based FFT method. And a PI regulator for specific harmonic is used to track the reference with phase compensation. At last, the reference signal for power amplifier is reconstructed in time domain by IFFT.

A PI regulator is applied to control the harmonic current amplitude to achieve zero-error control.²² The transfer function of PI controller for the fundamental current is expressed as:

$$H_1(s) = K_{p1} + \frac{K_{i1}}{s} \qquad (9)$$

where $K_{p1}$ is the proportional gain and $K_{i1}$ is integral coefficient for the fundamental current.

Let $U_1^*$ be the output voltage of PI regulator for the fundamental current, the transfer function is:

$$U_1^* \frac{k_{AMP}}{s\tau + 1} = I_1(s) \frac{sCr + 1}{sC} \qquad (10)$$

where $I_1(s)$ is the fundamental current. The current control model for a specific harmonic (here for the fundamental current) is shown in Figure 6. Then the open-
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Operation flow chart of HCIS.</figcaption></figure>

loop transfer function from $I_1(s)$ to $I_1^*(s)$ for the fundamental current is obtained as:

$$G(s) = \frac{K_{p1}(s + K_{i1}/K_{p1})}{s} \frac{k_{AMP}C_s}{rC\tau s^2 + (rC + \tau)s + 1} \quad (11)$$

Based on the open-loop transfer function (11), suitable cross frequency and phase margins are selected to meet the requirements of robustness and dynamic response. Then the values of the $K_{p1}$ and $K_{i1}$ can be determined with the help of Bode diagram.²³,²⁴ Assuming $k_{AMP} = 180/3$, $C = 600\mu F$ and $r = 0.01\Omega$, the obtained results are $K_{p1} = 0.001$, $K_{i1} = 169$. Similar to the fundamental current controller, the nth harmonic controller can be designed.

Considering the influence of the system inherent delay, such as the delay of AD, DA, and power amplifier, the phase deviation of each harmonic should be tested and compensated for control. A phase compensation is needed in this system to eliminate the impact of the deviation. Assuming the system parameters are generally invariant during the testing process, the angle offset of a specified harmonic order is basically constant. Thus, a preprocessing procedure can be used to detect the phase deviation of the relevant harmonic order.

Figure 7 shows the operation flow chart of the system. In the pretreatment segment, the phase delay of a given harmonic is detected and the angles is stored for the future experiments. Finally, the phase compensation can be directly used in the reference spectrum.

## Simulation results

In order to verify the effectiveness of the harmonic current detection and control strategy of the proposed

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Simulation results: (a) current waveform of the tested capacitor and (b) harmonic spectrum of injected current.</figcaption></figure>

<table><caption>**Table 2.** Parameters of simulation.</caption><thead><tr><td>Harmonic order</td><td>Amplitude (A)</td><td>Phase angle (°)</td></tr></thead><tbody><tr><td>1</td><td>75</td><td>0</td></tr><tr><td>7</td><td>29</td><td>0</td></tr><tr><td>13</td><td>91</td><td>0</td></tr></tbody></table>

HCIS, simulation model has been built based on MATLAB/Simulink. Table 2 shows the harmonics reference. The current is sampled with 2560 Hz sampling rate, and 256 points are used for FFT calculation. The frequency resolution ratio of this detection is 10 Hz, so the frame has integer number of cycles, and the spectrum leakage can be avoided. Figure 8(a) shows the current waveform of tested capacitor. As shown in Figure 8(b), the fundamental current is 75.01 A, the 7th harmonic is 29.0 A and the 13th harmonic is 90.99 A. It can be found that the harmonic components of the injected current are consistent with the reference current values shown in Table 2.

## Experimental results

In this section, experiments about HCIS are carried out based on the RT-LAB platform. The RT-LAB developed by Opal-RT Technologies has been widely employed in various domains. RT-LAB software can immediately compile and import mathematical model built with MATLAB/Simulink into real-time simulation platform for control and test. It is equipped with 16 channels for 16-bit ADC and the input ranges up to ±20 V. It also has 16 channels for 16-bit DAC and the output ranges up to ±15 V. The gap between the simulation clock and natural clock is only nanoseconds.
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Experimental test bench.</figcaption></figure>

<table><caption>**Table 3.** Experimental parameters of HCIS.</caption><thead><tr><td>Parameters</td><td>Value</td></tr></thead><tbody><tr><td>Model step size</td><td>1/12,800</td></tr><tr><td>Fundamental frequency</td><td>50 Hz</td></tr><tr><td>Filter capacitance</td><td>33 μF</td></tr><tr><td>Filter inductance</td><td>64 μH</td></tr><tr><td>V<sub>DC</sub> of power amplifier</td><td>180 V</td></tr><tr><td>Tested capacitor</td><td>200 μF</td></tr></tbody></table>

<table><caption>**Table 4.** Parameters of typical harmonic currents.</caption><thead><tr><td>Harmonic order</td><td>Amplitude (A)</td><td>Phase (°)</td></tr></thead><tbody><tr><td>1</td><td>5</td><td>0</td></tr><tr><td>5</td><td>1.22</td><td>0</td></tr><tr><td>7</td><td>4</td><td>0</td></tr><tr><td>13</td><td>13</td><td>0</td></tr></tbody></table>

<table><caption>**Table 5.** Phase angle compensation of typical harmonic currents.</caption><thead><tr><td>Harmonic order</td><td>Compensation angle (°)</td></tr></thead><tbody><tr><td>1</td><td>72.5</td></tr><tr><td>5</td><td>11</td></tr><tr><td>7</td><td>9</td></tr><tr><td>13</td><td>52</td></tr></tbody></table>

The experimental test bench is shown in Figure 9. Due to the limits of power level in lab, the transformer is not used and the compensation reactor is applied to reduce the active output power.

In this system, considering the sampling delay and the calculation burden of FFT algorithm, the sampling points is 1280 points and the model step size is 1/12,800. Table 3 shows the critical experimental parameters of HCIS. Both the harmonic current detection and control algorithms are implemented on the RT-LAB. And the

<figure><img src="image_12.png" /><figcaption>**Figure 10.** Reference current waveform of the system.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 11.** Injected current waveform of the tested capacitors.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Spectrum of the injected harmonic current.</figcaption></figure>

variable quantities of the system can be set in real time in the host computer.

Table 4 shows the parameters of several typical harmonic currents, which can be regarded as the reference values of HCIS. And the compensation phases for different harmonic orders can be calculated in the pre-processing stage, which are shown in Table 5.

After the phase compensation, the system can work normally with the harmonic current control strategy. As a comparison, Figure 10 shows the reference harmonic current waveform of the system. The injected harmonic current flowing through the tested capacitor is shown in Figure 11. The actual current is the same as the reference current. The spectrum of the injected harmonic current is calculated in Figure 12. As shown in Figure 12, the fundamental current is 5.02 A, the 5th
harmonic is 1.21 A, the 7th harmonic is 4.01 A, and the 13th harmonic is 12.98 A. It can be found that the harmonic components of the injected current are consistent with the reference current values shown in Table 4. Therefore, the injected harmonic current generated by the proposed HCIS can meet the requirements of audible noise testing for AC filter capacitors.

## Conclusion

In this paper, a FFT-based selective harmonic current control is proposed and applied to a HCIS, which can simulate the actual working conditions of AC filter capacitors in a converter station and help to analyze accurately the audible noise characteristics of AC filter capacitor in laboratory.

The theoretical analysis about harmonic detection algorithm based on FFT is presented and analyzed. Although the frequency leakage problem of FFT algorithm will affect the accuracy of measurement due to aperiodic sampling and cutting off sampling. This problem is not affecting so much. Since the amplitude and phase of each harmonic can be determined during one calculation procedure, the PI controller is designed to regulate the harmonic amplitude. A preprocessing procedure is applied to detect and compensate the angle offset. Once the amplitude of harmonics is detected and controlled, and the phase deviations are compensated, the reference signal can be reconstructed on the time field with IFFT and be used to control the power amplifier. Then the desired harmonic currents can be injected into the tested capacitor. The simulation and experimental results validate the proposed strategy.

## Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

## Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is supported by National Natural Science Foundation of China (Project Number:61673260)

## ORCID iDs

Jingang Han https://orcid.org/0000-0002-6347-1972

Shuoshuo Sun https://orcid.org/0000-0003-2634-0305

## References

1. Barnes M, Van Hertem D, Teeuwsen SP, et al. HVDC system in smart grids. Proc IEEE 2017; 105(11): 2082-2098.

2. High Voltage Direct Current (HVDC) substation audible noise, document IEC 22F/83/NP: 30-33, 2009.

3. Li J, Zhu L, Ji S, et al. Vibration analysis of power capacitors and vibration reduction with damping structures.

In: Proc IEEE int power modulator high voltage conf. (IPMHVC), San Francisco, CA, USA, 6-9 Jul 2016, pp.107-112. New York: IEEE.

4. Cox MD and Guan HH. Vibration and audible noise of capacitors subjected to nonsinusoidal waveforms. IEEE T Power Deliver 1994; 9(2): 856-862.

5. Zhu L, Li J, Shi Y, et al. Audible noise characteristics of filter capacitors used in HVDC converter stations. IEEE T Power Deliver 2017; 32(5): 2263-2271.

6. Jiang W, Hu Y, Zhang Y, et al. Different control objectives for grid-connected converter under unbalanced grid voltage using forgotten iterative filter as phase lock loop. IEEE T Power Electr 2015; 8(9): 1798-1807.

7. Zheng Peng F, Ott GW and Adams DJ. Harmonic and reactive power compensation based on the generalized instantaneous reactive power theory for three-phase four-wire systems. IEEE T Power Electr 1998; 13(6): 1174-1181.

8. Bahrani B, Kenzelmann S and Rufer A. Multivariable-PI based d-q current control of voltage source converters with superior axis decoupling capability. IEEE T Ind Electron 2011; 58(7): 3016-3026.

9. Lee T and Liu J. Modeling and control of a three-phase four-switch PWM voltage-source rectifier in d-q synchronous frame. IEEE T Power Electr 2011; 26(9): 2476-2489.

10. Yazdani D, Mojiri M, Bakhshai A, et al. A fast and accurate synchronization technique for extraction of symmetrical components. IEEE T Power Electr 2009; 24(3): 674-684.

11. Yu KKC, Watson NR and Arrillaga J. An adaptive Kalman filter for dynamic harmonic state estimation and harmonic injection tracking. IEEE T Power Deliver 2005; 20(2): 1577-1584.

12. Vodyakho O and Mi CC. Three-level inverter-based shunt active power filter in three-phase three-wire and four-wire systems. IEEE T Power Electr 2009; 24(5): 1350-1363.

13. Zheng S, Han J, Yao G, et al. Frequency-domian current control of harmonic current injection system for power capacitor. In: Proc IEEE Int Power Electron Appl Conf Exposit (PEAC), Shenzhen, CHN, 1-6 Nov 2018, pp.2009-2014. New York: IEEE.

14. Maza-Ortega JM, Rosendo-Macias JA, Gómez-Expósito A, et al. Reference current computation for active power filters by running DFT techniques. IEEE T Power Deliver 2010; 25(3): 1986-1995.

15. Xu S and Xu J. A current decoupling parallel control strategy of single phase inverter with voltage and current dual closed-loop feedback. IEEE T Ind Electron 2013; 60(4): 1306-1313.

16. Chen S and Huang H. Design of fractional order proportional integral controller using stability and robustness criteria in time delay system. Meas Control 2019; 52(9-10): 1552-1566.

17. Zou C, Liu B, Duan S, et al. Stationary frame equivalent model of proportional-integral controller in d-q synchronous frame. IEEE T Power Electr 2014; 29(9): 4461-4465.

18. Mattavelli P and Marafao FP. Repetitive-based control for selective harmonic compensation in active power filters. IEEE T Ind Electron 2004; 51(5): 1018-1024.

19. Tarisciotti L, Formentini A, Gaeta A, et al. Model predictive control for shunt active filters with fixed switching frequency. IEEE T Ind Appl 2017; 53(1): 296-304.
20. Chinnappan R, Logamani P and Ramasubbu R. Fixed-and variable-frequency sliding mode controller-maximum power point tracking converter for two-stage grid-integrated photovoltaic system employing nonlinear loads with power quality improvement features. *Meas Control* 2019; 52(7–8): 896–912.

21. Zhang Y, Li YW, Zargari NR, et al. Improved selective harmonics elimination scheme with online harmonic compensation for high-power PWM converters. *IEEE T Power Electr* 2015; 30(7): 3508–3517.

22. Han J, Zheng S, Yao G, et al. Multiple harmonic current injection system for audible noise analysis of AC filter capacitors in converter stations. *IEEE Access* 2020; 8: 94024–94032.

23. HO WK, Hang CC and Zhou J. Self-tuning PID control of a plant with under-damped response with specifications in gain and phase margins. *IEEE T Control Syst Technol* 1997; 5(4): 446–452.

24. Deshmukh AS and Vaishnav SR. A comparative performance analysis of PID tuning techniques based on frequency response specification. In: *Proc 2nd int conf emerg trends eng technol (ICETET)*, Nagpur, India, 16–18 Dec 2009, pp.971–974. New York: IEEE.