An Innovative Calibration Scheme for Interharmonic Analyzers in Power Systems under Asynchronous Sampling

Power-quality analyzers are commonly used in power systems to estimate waveform distortion, including the parameters of harmonics/interharmonics. In our study, a calibration scheme was developed and verified. This scheme is capable of calibrating the interharmonics specification of power-quality analyzers under asynchronous sampling. In our scheme, the hardware structure is composed of an interharmonic signal source, a wide-frequency resistive voltage divider, a broadband current shunt, and a data acquisition system. A new algorithm, based on discrete Fourier transform and interpolation, is presented. The procedure is implemented by LabVIEW software to process the sampling data and obtain the final interharmonic parameters. The test results of the amplitudes of the interharmonic current and voltage indicate that the calibration accuracy is 3.0‰ (16 Hz–6 kHz) and 6.8‰ (6 kHz–9 kHz) for the voltage signal, and 3.5‰ (16 Hz–6 kHz) and 6.5‰ (6 kHz–9 kHz) for the current signal. This index is higher than that acquired by the recommended methods in the International Electrotechnical Commission (IEC) standard.


Introduction
The International Electrotechnical Commission (IEC)-61000-2-1 standard firstly defined the term "interharmonic", and it was renewed with the subsequent IEC-61000-2-2 standard [1,2].According to the definition, interharmonic frequency refers to any frequency which is not a multiple of the fundamental frequency of a power supply or a frequency converter.Interharmonics attracted much attention in the past few years [3,4].In particular, electrical engineers who study current and voltage waveform distortion in power systems are interested in interharmonics [5,6].With nonlinear loads used increasingly in power systems, problems with interharmonics become more severe [7].The amplitude of the interharmonic and its ratio to the amplitude of the fundamental component are important parameters for evaluating the power-quality [8].Power-quality analyzers are commonly equipped in the power system to monitor its operation, including frequency deviation, voltage fluctuation and flicker, three-phase unbalance, and waveform distortion [9].In terms of waveform distortion, measuring the amplitude of harmonics/interharmonics of the waveform is a significant issue [10].To guarantee the accuracy of power-quality analyzers in interharmonic measurement, it is necessary to send them to legal metrological verification organizations of each country, such as the National Institute of Metrology in China or the Laboratoire National de Métrologie et d'Essais (LNE) Steers French Metrology in France, during a specified time interval.There are two verification methods, standard source method (SSM) and master meter method (MMM) [11,12].SSM employs a standard source to directly calibrate the instruments.The standard source should reach a corresponding level of accuracy.In some countries, such as America and Canada, Fluke 6100A, with enough accuracy (±5%) but high cost, is usually adopted as an interharmonic signal source for SSM.However, most interharmonic signal sources in China do not meet the accuracy requirement, and some even have no index for interharmonic parameters; thus, SSM is less used in general metrology organizations despite its convenience.MMM uses a meter as a reference, and this meter should be verified to reach a certain level of accuracy.When in use, MMM may be a bit complicated, but it does not have a strict requirement for input signal source.In this paper, based on MMM, we establish a new interharmonic calibration system.
The calibration algorithm largely decides the final measurement accuracy for any harmonic/ interharmonic calibration system.The approach recommended by IEC-61000-4-7 standard (hereafter called the IEC-approach) [13] is based on discrete Fourier transform (DFT).The IEC-approach calculates the amplitude of the interharmonics using a root-mean-square value of all interharmonic components between two adjacent harmonics.When all frequency components of the signal are synchronously sampled, the IEC-approach can obtain the interharmonic's amplitude without any error theoretically.However, to realize synchronous sampling, a highly elaborated measuring system set-up is required.In particular, the regular power-quality analyzers do not consider the situation where the fundamental frequency of power grid has some deviations, i.e., not always exactly 50 Hz/60 Hz.In this case, it is difficult for the data acquisition (DAQ) system to synchronously sample the signal [14].When the signal is asynchronously sampled, it is hard to estimate the measurement accuracy using the IEC-approach due to spectral leakage and the picket-fence effect [15,16].To reduce these errors, various harmonic/interharmonic measurement algorithms were put forward.In Reference [17], the method recommended by IEC was improved for asynchronous sampling conditions, but it requires that the fundamental component of the original signal is synchronously sampled.In Reference [18], a two-stage DFT-based algorithm was proposed to optimize the IEC-approach for asynchronously sampled signals.It can be used to estimate the frequency of each component, but not the amplitude parameter.In Reference [19], an interpolation DFT algorithm adding a Hanning window [20] was presented to measure the multi-frequency signal, but it mainly applies to harmonic, not interharmonic measurements.In Reference [21], an algorithm using fractional delay was proposed to handle asynchronous situations, but it was realized through massive calculations and a complex iterative process.There are other methods not adopting DFT, such as the methods based on a support vector machine [22], the Prony model [23], and wavelet packet transform [24].The methods mentioned above are often complex, and not appropriate for real-time engineering measurements.In view that the signal to be tested is only composed of one fundamental and one interharmonic component during the interharmonic calibration, a simple algorithm is preferred.This paper proposes a new algorithm based on DFT and interpolation, and it has a fast processing speed on sampled data and is easier to implement for real-time measurements.Via this algorithm, self-developed new voltage dividers, current shunts, and other devices, a low-cost calibration system on interharmonics was developed.Our interharmonic calibration system can be reproduced by general metrology laboratories.The developed voltage divider and the new measurement algorithm in this system can also be widely applied to measure other high-voltage signals and harmonic/interharmonic signals.
The rest of this paper is organized as follows: Section 2 gives a detailed description of the interharmonic calibration scheme, Section 3 tests the calibration system and evaluates the measurement error and uncertainty, and Section 4 summarizes the study.

Interharmonic Calibration Scheme
As shown in Figure 1, in our scheme, the interharmonic calibration system mainly consists of an interharmonic signal source, a wide-frequency resistive voltage divider (a combination of a divider and a buffer amplifier), a broadband current shunt, and a high-speed DAQ system.We describe each part in more detail below.

Interharmonic Signal Source
As a commercial instrument, the Fluke 6105A (Fluke Corporation, Everett, WA, USA) was used to produce the standard interharmonic voltage and current signals for the calibration.The Fluke 6105A can produce signals with a fundamental frequency between 16 Hz and 850 Hz and interharmonic frequencies between 16 Hz and 9 kHz.The accuracy (error) of the interharmonic's amplitude for Fluke 6105A was 1% (16 Hz to 6 kHz) and 4% (>6 kHz to 9 kHz) for both voltage and current.Fluke 6105A was actually used to first verify the accuracy of our scheme; then, it could be replaced by interharmonic signal sources that do not necessarily reach a certain accuracy.

Wide-Frequency Resistive Voltage Divider
The wide-frequency resistive voltage divider consisted of a divider and a buffer amplifier.We designed resistive voltage dividers with three rated ratios (75:1, 150:1, and 300:1).When the alternating current (AC) voltage ranges from 30 V to 240 V, the full-scale output voltage can be limited to 0.8 V using a suitable divider.Figure 2 shows the structure of a divider.Resistances with high precision (0.01%) were used to divide the voltage, and the stray parameters were compensated for by capacitances and shield rings.
When measuring high-frequency signals, the input impedance of the sampling device decreases rapidly.Because directly connecting the resistive voltage divider with the DAQ system causes a large error for the ratio, a buffer amplifier was required.The three-level bootstrap unit gain buffer amplifier, which was proposed by Budovsky at the Conference on Precision Electromagnetic Measurements (CPEM; 2010) [25], was used in our design.It was placed between the resistive voltage divider and the DAQ system to enable impedance matching.Figure 3 shows the circuit for the three-level buffer amplifier with output-current amplification.

Interharmonic Signal Source
As a commercial instrument, the Fluke 6105A (Fluke Corporation, Everett, WA, USA) was used to produce the standard interharmonic voltage and current signals for the calibration.The Fluke 6105A can produce signals with a fundamental frequency between 16 Hz and 850 Hz and interharmonic frequencies between 16 Hz and 9 kHz.The accuracy (error) of the interharmonic's amplitude for Fluke 6105A was 1% (16 Hz to 6 kHz) and 4% (>6 kHz to 9 kHz) for both voltage and current.Fluke 6105A was actually used to first verify the accuracy of our scheme; then, it could be replaced by interharmonic signal sources that do not necessarily reach a certain accuracy.

Wide-Frequency Resistive Voltage Divider
The wide-frequency resistive voltage divider consisted of a divider and a buffer amplifier.We designed resistive voltage dividers with three rated ratios (75:1, 150:1, and 300:1).When the alternating current (AC) voltage ranges from 30 V to 240 V, the full-scale output voltage can be limited to 0.8 V using a suitable divider.Figure 2 shows the structure of a divider.Resistances with high precision (0.01%) were used to divide the voltage, and the stray parameters were compensated for by capacitances and shield rings.
When measuring high-frequency signals, the input impedance of the sampling device decreases rapidly.Because directly connecting the resistive voltage divider with the DAQ system causes a large error for the ratio, a buffer amplifier was required.The three-level bootstrap unit gain buffer amplifier, which was proposed by Budovsky at the Conference on Precision Electromagnetic Measurements (CPEM; 2010) [25], was used in our design.It was placed between the resistive voltage divider and the DAQ system to enable impedance matching.Figure 3 shows the circuit for the three-level buffer amplifier with output-current amplification.

Interharmonic Signal Source
As a commercial instrument, the Fluke 6105A (Fluke Corporation, Everett, WA, USA) was used to produce the standard interharmonic voltage and current signals for the calibration.The Fluke 6105A can produce signals with a fundamental frequency between 16 Hz and 850 Hz and interharmonic frequencies between 16 Hz and 9 kHz.The accuracy (error) of the interharmonic's amplitude for Fluke 6105A was 1% (16 Hz to 6 kHz) and 4% (>6 kHz to 9 kHz) for both voltage and current.Fluke 6105A was actually used to first verify the accuracy of our scheme; then, it could be replaced by interharmonic signal sources that do not necessarily reach a certain accuracy.

Wide-Frequency Resistive Voltage Divider
The wide-frequency resistive voltage divider consisted of a divider and a buffer amplifier.We designed resistive voltage dividers with three rated ratios (75:1, 150:1, and 300:1).When the alternating current (AC) voltage ranges from 30 V to 240 V, the full-scale output voltage can be limited to 0.8 V using a suitable divider.Figure 2 shows the structure of a divider.Resistances with high precision (0.01%) were used to divide the voltage, and the stray parameters were compensated for by capacitances and shield rings.
When measuring high-frequency signals, the input impedance of the sampling device decreases rapidly.Because directly connecting the resistive voltage divider with the DAQ system causes a large error for the ratio, a buffer amplifier was required.The three-level bootstrap unit gain buffer amplifier, which was proposed by Budovsky at the Conference on Precision Electromagnetic Measurements (CPEM; 2010) [25], was used in our design.It was placed between the resistive voltage divider and the DAQ system to enable impedance matching.Figure 3 shows the circuit for the three-level buffer amplifier with output-current amplification.To verify the performance of the whole resistive voltage divider from 1 Hz to 10 kHz, we performed three groups of experiments as follows: a. Input voltages 37.5 V, 60 V, and 75 V for 75:1; b.Input voltages 75 V, 120 V, and 150 V for 150:1; c.Input voltages 150 V, 240 V, and 300 V for 300:1.
The relative errors for the rated voltage ratios are shown in Figure 4.For relatively low frequencies (<1 kHz), the compensation effect of the capacitances was not clear and the error increased with the frequency.However, for medium frequencies (1 kHz to 6 kHz), the capacitances could compensate for the stray parameters effectively and the error decreased.In addition, for high frequencies (>6 kHz), the effect from stray parameters exceeded the compensation, and the error increased.The overall relative error was below 0.01% (1 Hz to 6 kHz) and 0.04% (>6 kHz to 10 kHz).

Broadband Current Shunt
The Fluke A40B series (Fluke Corporation, Everett, WA, USA) [26] could be used as the broadband current shunt to measure current with a rated value from 0.1 mA to 100 A and frequencies ranging from direct current (DC) to 100 kHz (see Figure 1).We used shunts with four ranges (200 mA, 2 A, 10 A, and 20 A).The resistance of the shunt is very accurate and stable, with a low self-heating power coefficient and a tiny additional phase-shift.Therefore, it could be directly used to measure both DC and AC, with accuracies of 20 ppm and 23 ppm, respectively.

DAQ System
Two sets of Agilent 3458A were employed as digital multimeters to enable accurate sampling of voltage and current channels.The multimeters have three sampling modes: DC voltage (DCV) digitizing, direct sampling, and sub-sampling.In our calibration scheme, we used the DCV digitizing The relative errors for the rated voltage ratios are shown in Figure 4.For relatively low frequencies (<1 kHz), the compensation effect of the capacitances was not clear and the error increased with the frequency.However, for medium frequencies (1 kHz to 6 kHz), the capacitances could compensate for the stray parameters effectively and the error decreased.In addition, for high frequencies (>6 kHz), the effect from stray parameters exceeded the compensation, and the error increased.The overall relative error was below 0.01% (1 Hz to 6 kHz) and 0.04% (>6 kHz to 10 kHz).To verify the performance of the whole resistive voltage divider from 1 Hz to 10 kHz, we performed three groups of experiments as follows: a. Input voltages 37.5 V, 60 V, and 75 V for 75:1; b.Input voltages 75 V, 120 V, and 150 V for 150:1; c.Input voltages 150 V, 240 V, and 300 V for 300:1.
The relative errors for the rated voltage ratios are shown in Figure 4.For relatively low frequencies (<1 kHz), the compensation effect of the capacitances was not clear and the error increased with the frequency.However, for medium frequencies (1 kHz to 6 kHz), the capacitances could compensate for the stray parameters effectively and the error decreased.In addition, for high frequencies (>6 kHz), the effect from stray parameters exceeded the compensation, and the error increased.The overall relative error was below 0.01% (1 Hz to 6 kHz) and 0.04% (>6 kHz to 10 kHz).

Broadband Current Shunt
The Fluke A40B series (Fluke Corporation, Everett, WA, USA) [26] could be used as the broadband current shunt to measure current with a rated value from 0.1 mA to 100 A and frequencies ranging from direct current (DC) to 100 kHz (see Figure 1).We used shunts with four ranges (200 mA, 2 A, 10 A, and 20 A).The resistance of the shunt is very accurate and stable, with a low self-heating power coefficient and a tiny additional phase-shift.Therefore, it could be directly used to measure both DC and AC, with accuracies of 20 ppm and 23 ppm, respectively.

DAQ System
Two sets of Agilent 3458A were employed as digital multimeters to enable accurate sampling of voltage and current channels.The multimeters have three sampling modes: DC voltage (DCV) digitizing, direct sampling, and sub-sampling.In our calibration scheme, we used the DCV digitizing

Broadband Current Shunt
The Fluke A40B series (Fluke Corporation, Everett, WA, USA) [26] could be used as the broadband current shunt to measure current with a rated value from 0.1 mA to 100 A and frequencies ranging from direct current (DC) to 100 kHz (see Figure 1).We used shunts with four ranges (200 mA, 2 A, 10 A, and 20 A).The resistance of the shunt is very accurate and stable, with a low self-heating power coefficient and a tiny additional phase-shift.Therefore, it could be directly used to measure both DC and AC, with accuracies of 20 ppm and 23 ppm, respectively.

DAQ System
Two sets of Agilent 3458A were employed as digital multimeters to enable accurate sampling of voltage and current channels.The multimeters have three sampling modes: DC voltage (DCV) digitizing, direct sampling, and sub-sampling.In our calibration scheme, we used the DCV digitizing mode.This mode uses the DC measurement channel with a sampling frequency between 6 kHz and 100 kHz, and a sampling time between 500 ns and 1 s.The accuracies for the AC voltage and AC current could reach 100 ppm and 500 ppm, respectively.In addition, the accuracy of the frequency was 0.01%, and the time jitter was below 100 ps.During testing, the sampling frequency was set to 40 kHz, and the sample number was 5120.The sampled data from multimeters were transmitted to the computer and analyzed using the LabVIEW software.

Measurement of the Interharmonic
When calibrating the interharmonics of the power-quality analyzer, the signal generally contains both a fundamental component and the other interharmonic component.If the DAQ system can sample the signal synchronously, i.e., the sampling frequency is an integer multiple of the frequency of each component in the signal, we can obtain the amplitude and phase of the signal with minimal error.However, in most cases, the fundamental frequency and its interharmonic frequencies are unknown.Then, it is difficult to control the DAQ system and implement synchronous sampling.Using asynchronous sampling, the direct application of DFT leads to spectral leakage and the picket-fence effect.As a result, a significant error occurs.The existing methods based on fast Fourier transform (FFT), recommended by IEC 61000-4-7 and Institute of Electrical and Electronics Engineers (IEEE) 1459-2010 standards, can acquire accurate interharmonic parameters only if the signal is sampled synchronously.Under the asynchronous sampling condition, the errors are uncontrolled.In this paper, we propose a new interpolation algorithm to eliminate the error.
The signal is supposed to be sampled by the DAQ system with a sampling frequency f s and a sample number N; and it has the form shown in Equation (1).
where A 1 , f 1 , and ϕ 1 are the amplitude, frequency, and phase of the fundamental component; A 2 , f 2 , and ϕ 2 are the parameters of the interharmonic component, supposing f 1 < f 2 (convenient for discussion); and f s satisfies Shannon Theory (f s > 2 max{f 1 , f 2 }).f 0 = f s /N is defined as the spectral resolution, and we define parameters r 1 = f 1 / f 0 and r 2 = f 2 / f 0 to distinguish the sampling modes.When the signal is synchronously sampled, r 1 and r 2 are both integers.If the signal is asynchronously sampled, at least one of r 1 and r 2 is not an integer.For synchronous-sample mode, the interharmonic's parameters can be directly calculated by DFT.For asynchronous-sample mode, it can be divided into two cases: only one of the two frequency components is asynchronously sampled, or both of the components are asynchronously sampled.To improve the calibration accuracy under the condition of asynchronous sampling, we propose a new algorithm based on DFT and interpolation.The sampled data are processed using DFT first, and then the interpolation algorithm is applied to eliminate the errors caused by spectral leakage and the picket-fence effects.
For the case that neither of the two components is synchronously sampled, two lines with peak amplitudes in the DFT spectrogram are firstly searched, which are assumed to be the l 1 -th and the l 2 -th spectral lines.Secondly, near the l 1 -th line, another spectral line with its amplitude next to the l 1 -th line is found, and it is set as the l 3 -th line.Similarly, the fourth line is found near the l 2 -th line, and set as the l 4 -th line.As a result, a total of four lines in the spectrogram are obtained, with their values described as Y l m = R l m + jI l m (m = 1, 2, 3, 4), where R l m and I l m represent the real part and the imaginary part of Y l m .Thirdly, we utilize I l m (m = 1, 2, 3, 4) to establish Equation (2).On its left-hand side, the imaginary parts I l m are known.On the right-hand side, G m = 2πl m /N is known, while q s (s = 1, 2) and H s = 2πr s /N (s = 1, 2) are unknown and can be obtained from Equation (2).
Fourthly, we use R l m (m = 1, 2) to establish Equation (3).On its left-hand side, R l m is known.On the right-hand side, only p s (s = 1, 2) is unknown and can be obtained by solving two equations in Equation (3).
Finally, we acquire the amplitude for each component as follows: For the case that only one of the two components is asynchronously sampled, it is due to the sampling condition in IEC-approach being able to warrant for the synchronous sampling of fundamental component 50 Hz/60 Hz, but being unable to do so for the interharmonic component.We search two lines with peak amplitudes in the DFT spectrogram, which are assumed to be the l 1 -th (synchronously sampled component) and the l 2 -th (synchronously sampled component) spectral line.The frequency of the synchronously sampled component can be directly obtained by f 1 = l 1 f s /N.Similar to the interpolation approach proposed above, near the l 2 -th line, another spectral line with its amplitude next to the l 2 -th line is found, and set as the l 3 -th line.A total of three lines in the spectrogram are expressed as X l m = R l m + jI l m (m = 1, 2, 3), and I l m (m = 1, 2, 3) is used to establish equations as shown in Equation ( 5).On its left-hand side, the imaginary parts I l m are known.On the right-hand side, G m = 2πl m /N is known, H 1 = 2πr 1 /N is known, while q s (s = 1, 2) and H 2 = 2πr 2 /N are unknown and can be solved by Equation (5).
R p m (m = 1, 2) is used to establish equations as shown in Equation (6).On its left-hand side, R p i is known.On the right-hand side, only p s (s = 1, 2) is unknown and can be solved by Equation (6).
Finally, we acquire the amplitude for each component as follows: The proposed algorithm above is implemented by software LabVIEW in the host computer of the calibrating system.

Error Analysis
Figure 5 shows the hardware implementation of our interharmonic calibration scheme when calibrating the interharmonic voltage.When the wide-frequency resistive voltage divider is replaced by a broadband current shunt, it can be used to calibrate the interharmonic current.The program in LabVIEW has two functions, one is to control the two Agilent 3458A through a general purpose interface bus (GPIB) and, thus, determines the sampling frequency (30 kHz to 40 kHz) and sample number (5120) of DAQ system; the other is to implement the data-processing algorithm.To compare the measurement accuracy of the new algorithm and the IEC-approach, we programed both of them into LabVIEW.
Fifty-eight trials were run, and Figure 6 shows the relative errors of the interharmonics' amplitude for voltages and currents at different interharmonic frequencies with the range from 16 Hz to 9 kHz.The tests include the measurement of interharmonic voltage and current signals.We set the fundamental frequency as 50 Hz (working frequency in China).The fundamental voltage and current amplitudes were set as 30 V and 1 A, and the interharmonics' amplitudes were 10% of the fundamental amplitude.The interharmonics' frequencies were selected according to the IEC standard [27].Using the algorithm proposed in this paper, we obtained the maximum error of interharmonic amplitudes: for voltage signal, 3.0‰ for 16 Hz to 6 kHz, and 6.8‰ for 6 kHz to 9 kHz; for current signal, 3.5‰ for 16 Hz to 6 kHz, and 6.5‰ for 6 kHz to 9 kHz.However, when the IEC-approach was applied in the system, the maximum errors increased by up to 10 times either for voltage or for current, which indicates that the measurement accuracy of our scheme is higher than that required by the IEC standard.
amplitudes were set as 30 V and 1 A, and the interharmonics' amplitudes were 10% of the fundamental amplitude.The interharmonics' frequencies were selected according to the IEC standard [27].Using the algorithm proposed in this paper, we obtained the maximum error of interharmonic amplitudes: for voltage signal, 3.0‰ for 16 Hz to 6 kHz, and 6.8‰ for 6 kHz to 9 kHz; for current signal, 3.5‰ for 16 Hz to 6 kHz, and 6.5‰ for 6 kHz to 9 kHz.However, when the IEC-approach was applied in the system, the maximum errors increased by up to 10 times either for voltage or for current, which indicates that the measurement accuracy of our scheme is higher than that required by the IEC standard.

Evaluation of the Uncertainty
The interharmonic voltage was measured using a resistive voltage divider and sampled by a DAQ system.The model of the input voltage for the divider is described by Equation (8).amplitudes were set as 30 V and 1 A, and the interharmonics' amplitudes were 10% of the fundamental amplitude.The interharmonics' frequencies were selected according to the IEC standard [27].Using the algorithm proposed in this paper, we obtained the maximum error of interharmonic amplitudes: for voltage signal, 3.0‰ for 16 Hz to 6 kHz, and 6.8‰ for 6 kHz to 9 kHz; for current signal, 3.5‰ for 16 Hz to 6 kHz, and 6.5‰ for 6 kHz to 9 kHz.However, when the IEC-approach was applied in the system, the maximum errors increased by up to 10 times either for voltage or for current, which indicates that the measurement accuracy of our scheme is higher than that required by the IEC standard.

Evaluation of the Uncertainty
The interharmonic voltage was measured using a resistive voltage divider and sampled by a DAQ system.The model of the input voltage for the divider is described by Equation ( 8).

Evaluation of the Uncertainty
The interharmonic voltage was measured using a resistive voltage divider and sampled by a DAQ system.The model of the input voltage for the divider is described by Equation (8).
where U in is the amplitude of the voltage that is delivered to the divider, U out is the output voltage of the divider and the input voltage of the DAQ, K V is the rated ratio of the divider, and ε V is the error of the ratio.Five sources of uncertainty were considered (see Table 1).We combined them to obtain the standard uncertainty.For the voltage range (30 V to 240 V), the coverage factor was assumed to be two, and the expanded uncertainty was 0.5% (from 16 Hz to 6 kHz) and 2.0% (>6 kHz to 9 kHz), when calibrating the interharmonic voltage.
The interharmonic current was sampled using a current shunt and DAQ system, and the model of the input current for the current shunt can be described by Equation (9).
where I in is the amplitude of the current delivered to the shunt, U out is the output voltage of the shunt and the input voltage of the DAQ, K I is the rated ratio of the shunt, and ε I is the error of the ratio.Similar to the voltage calibration, five sources of uncertainty are listed in Table 1.We let the coverage factor be two and obtained the expanded uncertainty 0.3% (16 Hz to 6 kHz) and 1.2% (>6 kHz to 9 kHz) for this current range (100 mA to 10 A).
Table 1.Uncertainty sources and their calculated values.DAQ-data acquisition.

Interharmonic Current Measurement
The instability of the interharmonic source 0.23% 0.1% The inaccuracy of the rated ratio of divider 0.03% -The inaccuracy of the rated ratio of the shunt -0.002%The inaccuracy of the DAQ system 0.00023% 0.00023% The inaccuracy of the measurement algorithm 0.01% 0.01% The repeatability of the calibration system 0.12% 0.12%

Conclusions
We introduced and investigated an interharmonic calibration scheme for power-quality analyzers in power systems.In the hardware structure of the scheme, we used a Fluke 6105A as the standard interharmonic signal source, a set of Agilent 3458A as the high-speed data acquisition (DAQ) system, and Fluke A40B as broadband current shunt.We also developed wide-frequency resistive dividers (combination of a divider and buffer amplifier).After the signal was fed into the computer, a new algorithm based on both DFT and interpolation was used to analyze the signal under asynchronous sampling.Our results indicate that it can effectively reduce the error for asynchronous sampling.
For the calibration of the interharmonic voltage amplitudes, the interharmonic frequencies ranged from 16 Hz to 9 kHz, the error of amplitude was below 3‰ (16 Hz to 6 kHz) and 7‰ (>6 kHz to 9 kHz), and the expanded uncertainty was 0.5% (from 16 Hz to 6 kHz) and 2.0% (>6 kHz to 9 kHz).For the calibration of the interharmonic current, the interharmonic frequencies ranged from 16 Hz to 9 kHz, the error for the amplitude was below 4‰ (16 Hz to 6 kHz) and 7‰ (>6 kHz to 9 kHz), and the expanded uncertainties were 0.3% (16 Hz to 6 kHz) and 1.2% (>6 kHz to 9 kHz).The accuracy was improved using a thoroughly designed hardware and a new measurement algorithm that is different from the one described in the IEC standard.

Figure 1 .
Figure 1.Diagram of the interharmonic calibration system.

Figure 2 .
Figure 2. Structure of the resistive voltage divider compensated for by capacitances and shielding rings.

Figure 1 .
Figure 1.Diagram of the interharmonic calibration system.

Figure 1 .
Figure 1.Diagram of the interharmonic calibration system.

Figure 2 .
Figure 2. Structure of the resistive voltage divider compensated for by capacitances and shielding rings.

Figure 2 .
Figure 2. Structure of the resistive voltage divider compensated for by capacitances and shielding rings.

Figure 5 .
Figure 5. Hardware implementation for the interharmonic calibration scheme.

Figure 6 .
Figure 6.Test results of voltage and current interharmonics' amplitude.

Figure 5 .
Figure 5. Hardware implementation for the interharmonic calibration scheme.

Figure 5 .
Figure 5. Hardware implementation for the interharmonic calibration scheme.

Figure 6 .
Figure 6.Test results of voltage and current interharmonics' amplitude.

Figure 6 .
Figure 6.Test results of voltage and current interharmonics' amplitude.