HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS

Similar documents
Measurement of RMS values of non-coherently sampled signals. Martin Novotny 1, Milos Sedlacek 2

ON THE VALIDITY OF THE NOISE MODEL OF QUANTIZATION FOR THE FREQUENCY-DOMAIN AMPLITUDE ESTIMATION OF LOW-LEVEL SINE WAVES

Measurement of Amplitude Ratio and Phase Shift between Sinusoidal Voltages with Superimposed Gaussian Noise. Pawel Rochninski and Marian Kampik

POWER QUALITY MEASUREMENT SYSTEM BASED ON A DIGITAL SIGNAL PROCESSOR. Aleksandar Prodic and Predrag Pejovic*

An improvement for dual channel sampling wattmeter

ADC Based Measurements: a Common Basis for the Uncertainty Estimation. Ciro Spataro

Programming PC-Based Multi-Function Oscilloscopes: A User-Friendly Approach to Rapid Prototyping of Spectral Analysis Code

Harmonic Signal Processing Method Based on the Windowing Interpolated DFT Algorithm *

FAST Fourier Transform (FFT) and Digital Filtering Using LabVIEW

Metrol. Meas. Syst., Vol. XIX (2012), No. 4, pp METROLOGY AND MEASUREMENT SYSTEMS. Index , ISSN

ON THE BIAS OF TERMINAL BASED GAIN AND OFFSET ESTIMATION USING THE ADC HISTOGRAM TEST METHOD

ME scope Application Note 01 The FFT, Leakage, and Windowing

Discrete Fourier Transform (DFT)

METHOD OF TESTING AND CORRECTING SIGNAL AMPLIFIERS TRANSFER FUNCTION USING PRONY ANALYSIS

THE DESIGN OF RURAL POWER NETWORK POWER QUALITY MONITORING AND ANALYSIS PLATFORM ON LABVIEW

A POWER QUALITY INSTRUMENT FOR HARMONICS INTERHARMONICS AND AMPLITUDE DISTURBANCES MEASUREMENTS

Application of Digital Sampling Method for Voltage Transformer Test Set Calibrations. Hüseyin Çaycı

Detection Probability of Harmonics in Power Systems Affected by Frequency Fluctuation

Metrological Characterization of Low Power Instrument Transformer Integrated in MV Recloser

A Faster Method for Accurate Spectral Testing without Requiring Coherent Sampling

A NEW MOTOR SPEED MEASUREMENT ALGORITHM BASED ON ACCURATE SLOT HARMONIC SPECTRAL ANALYSIS

User-friendly Matlab tool for easy ADC testing

Research on Flicker Measurement Algorithm Based on FFT

LAB #7: Digital Signal Processing

IMPLEMENTATION OF VLSI BASED ARCHITECTURE FOR KAISER-BESSEL WINDOW USING MANTISSA IN SPECTRAL ANALYSIS

Compensation of Analog-to-Digital Converter Nonlinearities using Dither

USB Dynamic Signal Acquisition

Grid Power Quality Analysis of 3-Phase System Using Low Cost Digital Signal Processor

HARMONIC and interharmonic measurements are one of

Dynamic DAC Testing by Registering the Input Code when the DAC output matches a Reference Signal

ni.com Sensor Measurement Fundamentals Series

SITE-TO-SITE REPRODUCIBILITY IN CONDUCTED IMMUNITY TESTS ON PC-BASED DATA ACQUISITION SYSTEMS

Lab VIEW Programming for Vibration Analysis

FOURIER analysis is a well-known method for nonparametric

VARIABLE-FREQUENCY PRONY METHOD IN THE ANALYSIS OF ELECTRICAL POWER QUALITY

CHARACTERIZATION and modeling of large-signal

IMPROVING KWH-METER PERFORMANCE AT PV ON GRID SYSTEM BY MULTIPLYING THE NUMBER OF SAMPLING SIGNAL

Data acquisition and instrumentation. Data acquisition

Chapter 5 Window Functions. periodic with a period of N (number of samples). This is observed in table (3.1).

Low distortion signal generator based on direct digital synthesis for ADC characterization

Butterworth Window for Power Spectral Density Estimation

MULTIFUNCTION POWER QUALITY MONITORING SYSTEM

Direct Digital Synthesis Primer

Design of the Chaotic Signal Generator Based on LABVIEW

Data Acquisition Systems. Signal DAQ System The Answer?

Power Quality Measurements the Importance of Traceable Calibration

Frequency Domain Representation of Signals

The Hardware Design of Power Quality Online Monitoring Device Based. on MCU

ENGINEERING FOR RURAL DEVELOPMENT Jelgava, EDUCATION METHODS OF ANALOGUE TO DIGITAL CONVERTERS TESTING AT FE CULS

EE 435. Lecture 34. Spectral Performance Windowing Quantization Noise

Real-Time Selective Harmonic Minimization in Cascaded Multilevel Inverters with Varying DC Sources

Sub/super-synchronous harmonics measurement method based on PMUs

Virtual Instrument for FPGA based Spectrum Analyzer

MASSACHUSETTS INSTITUTE OF TECHNOLOGY /6.071 Introduction to Electronics, Signals and Measurement Spring 2006

Open Access Sparse Representation Based Dielectric Loss Angle Measurement

A novel power harmonic analysis method based on Nuttall-Kaiser combination window double spectrum interpolated FFT algorithm

Advanced Lab LAB 6: Signal Acquisition & Spectrum Analysis Using VirtualBench DSA Equipment: Objectives:

Laboratory Experiment #1 Introduction to Spectral Analysis

TIMA Lab. Research Reports

Agilent Technologies 3000 Series Oscilloscopes

Lab 12 Laboratory 12 Data Acquisition Required Special Equipment: 12.1 Objectives 12.2 Introduction 12.3 A/D basics

2.1 Make Accurate Low-Level Measurements with High-Resolution Instrumentation. Incorporating Best Practices for Ensuring Product Quality

Data Sheet. Digital Storage Oscilloscope. Features & Benefits. Applications. Ease-of-Use Feature DSO5202BMT DSO5102BMT DSO5062BMT

EE 422G - Signals and Systems Laboratory

Introduction. sig. ref. sig

Current Rebuilding Concept Applied to Boost CCM for PF Correction

Carrier Frequency Offset Estimation in WCDMA Systems Using a Modified FFT-Based Algorithm

A COMPARATIVE ANALYSIS IN TERMS OF CONDUCTED SUSCEPTIBILITY OF PC-BASED DATA ACQUISITION SYSTEMS

SHORT TIME ALGORITHM OF POWER WAVEFORMS FUNDAMENTAL HARMONIC ESTIMATION WITH USE OF PRONY S METHODS

Virtual FFT Analyser for identification of harmonics and inter-harmonics metrological aspects

Windows and Leakage Brief Overview

Chapter 2 Analog-to-Digital Conversion...

Measurement, Sensors, and Data Acquisition in the Two-Can System

Monitoring of Power Quality in Industry

A COMPLEX CURRENT RATIO DEVICE FOR THE CALIBRATION OF CURRENT TRANSFORMER TEST SETS

Design of Frequency Characteristic Test Instrument Based on USB

U1604A Handheld Oscilloscopes, 40 MHz

Control and Signal Processing in a Structural Laboratory

Available online at ScienceDirect. Measurement of Power Frequency with Higher Accuracy using PIC Microcontroller

A new scheme based on correlation technique for generator stator fault detection-part π

The Fundamentals of Mixed Signal Testing

Realization and characterization of a smart meter for smart grid application

Mel Spectrum Analysis of Speech Recognition using Single Microphone

Fourier Theory & Practice, Part I: Theory (HP Product Note )

A Compressed Sampling-Based Method Compliant with IEC for Harmonic and Interharmonic Measurements

Analyzing A/D and D/A converters

6 Sampling. Sampling. The principles of sampling, especially the benefits of coherent sampling

Gentec-EO USA. T-RAD-USB Users Manual. T-Rad-USB Operating Instructions /15/2010 Page 1 of 24

Capacitive MEMS accelerometer for condition monitoring

Using dspace in the Shunt Static Compensators Control

Module 1: Introduction to Experimental Techniques Lecture 2: Sources of error. The Lecture Contains: Sources of Error in Measurement

SAMPLING THEORY. Representing continuous signals with discrete numbers

U1571A Ni-MH Battery Pack for U1600A Handheld Oscilloscopes

The Fundamentals of FFT-Based Signal Analysis and Measurement Michael Cerna and Audrey F. Harvey

Transfer Function (TRF)

Publication List. 1. List of Relevant Articles (all are in ISI foreign journals): 2. Doctoral Thesis:

THE APPLICATION WAVELET TRANSFORM ALGORITHM IN TESTING ADC EFFECTIVE NUMBER OF BITS

Fourier Signal Analysis

PXIe Contents CALIBRATION PROCEDURE. Reconfigurable 6 GHz RF Vector Signal Transceiver with 200 MHz Bandwidth

Wavelet Transform Based Islanding Characterization Method for Distributed Generation

Transcription:

Metrol. Meas. Syst., Vol. XVIII (2011), No. 2, pp. 315-322 METROLOGY AND MEASUREMENT SYSTEMS Index 330930, ISSN 0860-8229 www.metrology.pg.gda.pl HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS Yue Tong 1), Guoxiong Ye 1), Keqin Guo 1), Hongbin Li 1,2) 1) State Grid Electronic Power Research Institut 2) Huazhong University of Science& Technology, Institute of Electrical & Electronic Engineering.Chian ( cherry8386@126.com) Abstract A high accurate electronic instrument transformer calibration system is introduced in this paper. The system uses the fourth-order convolution window algorithm for the error calculation method. Compared with Fast Fourier Transform, which is recommended by standard IEC-60044-8 (Electronic current transformers), it has higher accuracy. The relative measuring errors caused by asynchronous sampling could be reduced effectively without any special hardware technique adopted. The results show that the ratio error caused by asynchronous sampling can be reduced to 10-4, and the phase error can be reduced to 10-3 degrees when the deviation of frequency is within ±0.5 Hz. The present method of measurement processing is achieved by a high-accuracy USB multifunction data acquisition (DAQ) card and virtual measurement devices, with low cost, short exploitation period and high stability. Keywords: Electronic instrument transformers, Digital output, Calibration system, Virtual instrument, Fourthorder Convolution Window 2011 Polish Academy of Sciences. All rights reserved 1. Introduction As is well known, electronic transformers have been used in power systems because of their advantages of wide dynamic range, small size and adaptability for the development of digitization. As one of the most important equipments in power system, the accuracy of transformers is the basic guarantee for the metering, protection, monitoring and analysis of the power system. Take electronic instrument current transformers for example, according to the existing standards [1, 2], the electronic instrument current transformers for metering are classified into six classes: 0.1, 0.2, 0.5, 1, 3, 5; special purpose used instrument current transformers are classified into 0.2 s and 0.5 s. The ratio error and phase error in 1 % of rated current are also demanded. The accuracy of a relevant calibration system must be two grades higher. It means that the error of the calibration system must be less than 0.02 % for ratio error and 2 for phase error. Therefore highly accurate calibration of electronic transformers is very necessary. The requirement of electronic instrument voltage transformers is quite the same [3]. There are three main error sources of calibration system: the standard CT/VT, the digital sampling system and the error calculation method. Using the conventional CT and VT as standard transformer is easy to achieve a high accuracy, where the uncertainty is less than 0.01% for the ratio error and 0.3 for the phase error. The error of the digital sampling system includes two parts. One is the quantization error of the A/D converter and the other is the synchronous sample error. This system uses a NI 6281 for the digital sampling system. NI (National Instrument) Multifunction Series high-accuracy USB data acquisition (DAQ) devices are optimized for high accuracy at fast sampling rates with an outside trigger sample. The last and most important part of the main error resource is the error calculation method. Article history: received on Aug. 25, 2010; accepted on Apr. 05, 2011; available online on May 20, 2011; DOI: 10.2478/v10178-011-0013-x.

Yue Tong, Hongbin Li et al.: HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS The electronic apparatus and high voltage equipment are widely used in the power system which causes harmonic and inter-harmonic distortion in the power line. According to IEC- 60044-8, the range of standard frequency is 99 %fr ~101 % fr (fr is the rated frequency for metering), that is from 49.5 Hz to 50.5 Hz, when the rated frequency is 50 Hz. Even if the sampling frequency satisfies the Nyquist sampling criterion, the classic FFT, which is recommended by standard IEC-60044-8, cannot avoid the long-range spectral leakage, caused by an asynchronous sample. There must be some special signal processing for error calculation to meet the high accuracy requirement. In this paper, a highly accurate calibration system using a multifunction USB data acquisition card for the digital sampling system and a virtual instrument achieving the fourthorder convolution window of DFT (Discrete Fourier Transform) algorithm is presented. Labview is used as the unified software platform for computing, storage, calling, and displaying. The form of virtual instrument that software is machines provides a convenient way for the present calibration system [4]. 2. PRINCIPLE OF OPERATION The block diagram of the calibration system is shown in Fig. 1. The calibration system consists of a standard transformer, a digital sampling system (DDS) and calibration software. The calibrated ECT/EVT consists of primary converter, secondary converter (SC) and MU (merge units) with the digital output from Ethernet. A sampling sync circuit (SSC) is used for two signal synchronous sampling the output from DDS and MU respectively. Calibrated ECT/EVT EVT Standard CT/VT busbar SC I/V V or V/VV converter Calibration system MU SSC DDS Ethernet PC USB Fig.1 Block diagram of the calibration system for ECTs/EVTs 2.1. Standard CT/VT The conventional CT/VT is used as the standard transformer in this system. It can be easy to achieve a high accuracy of uncertainty less than 0.01% for ratio error and 0.3 for phase error. With the I/V or V/V converter, the signal output from the standard CT/VT can convert into a voltage signal in order to appropriately access to the DDS. 2.2. Digital sampling system The presented DDS is a high accuracy USB multifunction data acquisition (DAQ) card, NI-USB6281. The analog input is at 18-bits, 625 ks/s. Combined with advanced external synchronization and isolation it can start sampling under the external trigger. This highaccuracy device for USB delivers high-performance data acquisition in an easy-to-use and

Metrol. Meas. Syst., Vol. XVIII (2011), No. 2, pp. 315-322 portable form factor through USB ports on a personal computer and other portable computing platforms. The 18-bit analog-to-digital converter (ADC) and available filtering provide an increase in resolution and higher measurement sensitivity. With fast sampling rates and a low noise floor, these devices can accurately acquire a dynamic signal. 2.3 Method of error calculation algorithm The asynchronous sample caused by grid frequency fluctuation, results in leakage of the spectrum. A fundamental method of reducing long-range spectral leakage is selecting an appropriate window for sampling sequences. Some famous windows have been widely used, for example, a rectangular window, Hanning window, Blackman-Harris window, Rife-Vincent window and so on. In periodic signal measurement, the yardstick of windows is: minimum main lobe width, rapid side lobe attenuation. In recent years a new type of window has been raised the convolution window. The convolution window is a kind of windows which is constructed by convolutions via a few rectangular windows with same time width. When the sampling synchronous deviation is slight, just as the power network fluctuation, the convolution window is better than the other famous windows in decreasing the error of the asynchronous sample [5-9]. Rectangular window function and its Fourier transform are as follows: Rectangular window with time width T 0 can be written as: 1 t 1/ T0 t / T0 1/ 2 w1 ( t) = rect( ) = T0 T0 0 t / T0 > 1/ 2. The Fourier transform of a rectangular window is given by: j ft sin( πt f ) W( f ) w( t) e dt sinc( T f ). (2) 2π 0 1 = 1 = 0 = πt0 f The k-th order convolution window Wk(t) is defined as the convolutions of k rectangular windows with width T 0 : w ( t) = w ( t) w ( t)... w ( t). (3) k 1 1 1 According to the law of convolution, the Fourier transform of a k-th order convolution window in the frequency domain is given by: sin( πt f ) W f = T f = (4) k 0 k k ( ) [sinc( 0 )] [ ]. πt0 f The higher order window can accelerate the attenuation of the side lobe and it improves the precision of DFT. On the other hand, the width of the main lobe is increasing, which reduces the spectrum resolution. According to the needs of the calibration system, a fourth-order convolution window is chosen in this paper. Omitting the process of derivation, the fourth-order convolution window expression in the time domain is presented to be: 0 x 2 1 3 w4 ( t) = (2 x) 1 x < 2, (5) 3! T0 3 3 (2 x) 4(1 x) x < 1 where x = t / T0. (1)

Yue Tong, Hongbin Li et al.: HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS The accuracy of the algorithm is affinitive to the zero point value of the window function. And the fourth-order convolution window s third-order derivative is zero, which means that the fourth-order convolution window has an excellent performance in harmonic analysis and frequency fluctuation. 2.4. Design of calibration software The calibration system based on virtual instrument replaces the traditional calibrator with personal computer. That is easy to use, graphical user interface to programmer remarkably shortens exploitation period and reduces the need for high level in language programming and debugging. The Labview is multithreading and multitasking program. The main flow chart is shown in Fig. 2. The whole program includes receiving data from DAQ and network adapter, digital signal processing and error calculation [10]. Data from standard transducer is sampled by DAQ provided in Labview through NI- USB6281. It can be called expediently to specify the parameters like sample rate, samples to read and sample start with trigger. Receiving data from calibrated transducer through MU is much more complicated. The output interface of electronic instrument transformers is based on IEEE 802.3, which is beyond the offer of the Labview existing module. It needs an external dynamic link library with WinPcap function to achieve this task. The digital signal processing program is the core of the whole system. It achieves the fourth order convolution window algorithm whose time domain expression is proposed above. Then the samples are multiplied by the windows in time domain. Through discrete Fourier transform the parameters of fundamental frequency signal such as altitude and phase can be abstracted in frequency domain. Finally, the ratio error and phase error are calculated. The calculation results and signal waves are presented in front panel, while the program is presented in back panel. And the software also can do storage, playback, and document management. START Initialize the network adapter and DAQ Establish the communication Read data Fourth order convolution window for DFT The number of read n=n+1 n Sampled values multiple with window in time domain no n>1024? yes Frequency domain analysis extracting foundation frequency Calculate the ratio error and phase error Results show and record no Ask for stop? yes Clear up the memory End Fig. 2. The flow chart of software

Metrol. Meas. Syst., Vol. XVIII (2011), No. 2, pp. 315-322 3. EXPERIMENT The first experiment is to verify that the fourth-order convolution window of the DFT algorithm can obtain a more accurate representation of error estimate of electrical parameters. This experiment is simulated in Labview, using a virtual signal generator to simulate the output of an electronic instrument transformer. Generally speaking, the fluctuation of frequency in the grid is caused by harmonics. The odd harmonic is more harmful than the even harmonic. In a balanced three-phase system, the even harmonic is cancelled and only the odd harmonic left. So in this experiment, a signal with the 3 rd, 5 th, 7 th harmonic is simulated as the resource. The simulation resource in Labview is shown in Table 1: Table 1 The simulation resource Harmonic order 1 3 5 7 Amplitude(V) 100 30 20 10 phase() 0 20 130 150 Under the same conditions that the sample rate is 6.4 khz, the sampling number is 1024, and the simulation result is shown in Table 2 Table 2 The measuring error of harmonic test Harmonic order Window function A1(%) A3(%) A5(%) A7(%) Blackman-Harris 2.45E-4 1.74E-4 7.12E-5 4.12E-4 fourth-order convolution window 3.14e-8 5.11e-6 1.45e-6 1.37e-6 φ1( ) φ3( ) φ5( ) φ7( ) Blackman-Harris 4.87E-4 2.16E-4 2.56E-5 3.42E-4 fourth-order convolution window 1.86E-8 2.12E-6 8.45E-7 5.68E-6 The regular harmonic simulation result shows that the fourth-order convolution window algorithm can reduce the harmonic calculation error effectively. The second experiment is confirming the error caused by asynchronous sampling with a frequency fluctuation from 49.5 Hz to 50.5 Hz. It is carried out by comparing two algorithms with the same resource whose frequency varies from 49.5 Hz to 50.5 Hz in the same amplitude and phase. The result is shown in Table 3: Table 3 The measuring error of frequency fluctuation test Window function Frequency A49.5(%) A49.8(%) A50(%) A50.2(%) A50.5(%) Blackman-Harris -0.18-0.08 0-0.06-0.26 fourth-order convolution window -3.85E-4-4.26E-5 0-2.46E-5-2.7E-4 φ49.5( ) φ49.8( ) φ50( ) φ50.2( ) φ50.5( ) Blackman-Harris -0.27-0.09 0-0.073-0.25 fourth-order convolution window -7.5E-3-2.5E-4 0-1.75E-4-8.3E-3 Furthermore, another experiment is taken to compare with the calibration system in the State Grid Electric Power research Institute (SGEPI), where the highest national standard is established. A sample electronic instrument transformer is tested under two calibration systems. The sample electronic instrument transformer is a Gas-Insulated Metal-Enclosed Switchgear (GIS) 220 kv combined voltage and current transformer. The current transformer

Yue Tong, Hongbin Li et al.: HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS is of class 0.2 s, and the voltage transformer is class 0.1 for metering and class 3P for protecting. The results are shown in Fig. 3 and Fig. 4. Two groups of test data show basic coincidence. Fig.3 The normalized ratio error and phase error results of the current transformer comparison tests 4. Performance evaluation Fig.4 The normalized ratio error and phase error results of the voltage transformer comparison tests Every segment of the calibration system is traceable. The standard CT/VT uncertainty can be measured by comparing with a much more accurate reference transducer. The DDS, the

Metrol. Meas. Syst., Vol. XVIII (2011), No. 2, pp. 315-322 NI-USB-6281 device also has a two-year calibration interval in NI company services to help maintain the ongoing accuracy of the measurement hardware. Reference CT/VT Standard CT/VT busbar High-accuracy I/V V or V/VV converter I/V V or V/VV converter HP3458 3458A SSC DDS Calibration system PC Fig.5 Block diagram of the test setup for verification of the calibration Besides that, the calibration system can be evaluated in a laboratory too. The block diagram of the test setup is shown in Fig. 5. The test set includes a reference CT/VT, highaccuracy I/V or V/V converter and HP3458A. The reference CT/VT and the high-accuracy converter must be of at least few times higher accuracy than the standard CT/VT and I/V or V/V converter tested. A precise digital voltmeter HP3458A is used as a reference A/D converter which is used in numerous high precision measurements. It provides an effective precision of 24 bits with the least significant bit weighted 2-23. The error of the test setup is small enough to be ignored when it is used for verification of the calibration. Under this test setup, the linearity of the presented calibration system limits in 0.02% for the ratio error and 1 in the phase error, when the dynamic range is from 1% Ie~120% I e (I e is the rated current value) and from 80%U e ~120%U e (U e is the rated voltage value). It is qualified to calibrate class 0.2 s electronic current instrument transformers and class 0.1 electronic voltage instrument transformers. 5. Conclusions A high-accuracy calibration system for electronic instrument transformers is presented in this paper. A virtual instrument as a platform to achieve a fourth-order rectangular convolution window of the DFT algorithm can effectively reduce the measuring error caused by an asynchronous sample. The digital sampling system which is based on an USB data acquisition card can avoid the noise caused by hardware, and decrease the source of the measuring error. The system is portable and excellent in both flexibility and computing speed of digital signal processing as well. The realization with a virtual instrument is low fabrication costs and a short development cycle. This system has achieved a high accuracy of electronic instrument transformers measurement, which meets the requirement of calibrating class 0.2s electronic current instrument transformers and class 0.1 electronic voltage instrument transformers. References [1] IEC60044-7. (2002). Instrument Transformers-part 7. Electrical Voltage Transducers. [2] IEC60044-8. (2002). Instrument Transformers-part 8. Electrical Current Transducers.

Yue Tong, Hongbin Li et al.: HIGHLY ACCURATE CALIBRATION SYSTEM FOR ELECTRONIC INSTRUMENT TRANSFORMERS [3] Naumovic-Vukovic, D., Skundric, S., Kovacevic, D., Milosavljevic, S. (September 2009) Calibration of high accuracy class standard current transformers. XIX IMEKO World Congress Fundamental and Applied Metrology. Lisbon. [4] Bilski, P., Winiecki, W. (2002). Virtual Spectrum Analyzer Based on Data Acquisition Card, IEEE Transaction on instrumentation and measurement, 51(1), 82-87. [5] Zhang Jieqiu, Liang Changhong, Chen Yanpu (2004). Power System Harmonic Theory Analysis and Algorithm Based on Convolution Windows. Proceedings of the CSEE, 21(11), 48-52. [6] Nuttall, A.H. (1981). Some window with very good sidelobe behavior. IEEE Trans. On ASSP, 29(1), 84-91. [7] Zhao Fusheng, Geng Zhongxing, Ge Yaozhong (2001). FFT algorithm with high accuracy for harmonic analysis in the electric machine. Proceedings of the CSEE, 21(12), 83-87. [8] Harris, F.J. (1978). On the use of windows for harmonic analysis with the discrete Fourier transform. Proc. IEEE, 66(1), 51-83. [9] Jain, V.K., Collins, W.L., Davis, D.C. (1983). High-accuracy analog measurement via interpolated FFT. IEEE Trans on Instrum. Meas., 32(6), 350-355. [10] Teng J.H., Chan S.Y. (2000). A Labview Based Virtual instrument for Power Analyzers. IEEE Power System Technology, (1), 179-184.