Numerical Differential Protection of Power Transformer using Algorithm based on Fast Haar Wavelet Transform
|
|
- Sheila Griffith
- 5 years ago
- Views:
Transcription
1 IDIA ISTITUTE O TECHOLOGY, KHARAGPUR 7232, DECEMBER 27-29, umerical Differential Protection of Power Transformer using Algorithm based on ast Haar Wavelet Transform K. K. Gupta and D.. Vishwakarma Abstract-- A fast and simple numerical filtering algorithm is an important requirement for the efficient power system relaying. This paper presents a wavelet-based algorithm for numerical differential protection of power transformer, using harmonic restraint approach. This algorithm provides an accurate and computationally efficient tool for distinguishing internal faults from magnetizing inrush and over excitation inrush. This algorithm is essentially a numerical filter, which extract the fundamental frequency and second and fifth harmonic s of differential current, to provide operating and restraining signal respectively. This algorithm generates it s coefficient by additions and subtractions routines only and does not involve time consuming multiplication and division calculation. The computation complexity of the algorithm presented here is O() additions and subtractions as compared to DT based algorithm, whose computation complexity is O(log2) operations. Index Terms Differential Protection, Haar Transform, Wavelet Transform, ourier Transform. I. ITRODUCTIO A power transformer belongs to a class of very expensive and vital of electric power system and its protection is one of the most challenging problems in the area of power system relying. The frequency of occurrence faults in power transformer is less than on lines. But if a power transformer experiences a fault, it is necessary to take the transformer out of service as soon as possible so that damage is minimized. The cost associated with repairing a damaged transformer may be very high. An unplanned outage of a power transformer can cost electric utilities crores of rupees. Consequently, it is of great importance to minimize the frequency and duration of unwanted outages. Accordingly, high demands are imposed on power transformer protective relays. Requirements include Dependability (no missing operation) Security (no false tripping) Speed of operation (short fault clearing time) Differential protection scheme based on circulating current principle is widely used to protect the power transformer against internal faults [4], [9]. The authors are with the Institute of Technology, Banaras Hindu University, Varanasi - 225, India (telephone: , gupta_k_k@rediffmail.com, dnv@banaras.ernet.in). The differential protection scheme shows certain limitations. Detection of a differential current does not provide a clear distinction between internal faults and other conditions, e.g. Magnetizing inrush, over-excitation of the transformer core, external faults. The conventional approach to mitigate these problems is to apply percentage (biased) differential characteristic along with second and fifth harmonic restraints for inrush and over-excitation conditions, respectively. The conventional percentage differential relays used for the protection of power transformer against internal faults are either of electromagnetic or static type. Electromagnetic and static relays have several drawbacks. The concept of numerical protection, which evolved during the last two decades, shows much promise in providing improved performance. The main features of the umerical relays are their economy, reliability, compactness, flexibility and the possibility of integrating a umerical relay into the hierarchical computer system within the substation. II. AST HAAR WAVELET TRAORM (HWT) TECHIQUE The proposed algorithm is derived from the well-known Haar Transform based algorithm. In Haar Transform the Haar coefficients are calculated first and from Haar coefficients the sine and cosine ourier coefficients are calculated, using established relationship between Haar and ourier coefficients. The Haar transform can also be termed as Discrete Wavelet Transform (DWT), using Haar wavelet. The Discrete Wavelet (Haar) Transform requires log 2 operations to transform a sample vector. The DWT matrix is not sparse in general, so it has the same complexity issues as the discrete ourier transform. So it is solved in same way as for the T, by factoring the DWT into a product of a few sparse matrices using selfsimilarity properties. This result in an algorithm that requires only order operations to transform an -sample vector. This fast version of the DWT can be called as ast Wavelet Transform (WT). The WT using Haar Wavelet calculates Haar coefficients by 4(-) multiplications. Hence its computation complexity is O() operations. In the algorithm presented here these 4(-) multiplications are converted in 2(-) additions. Hence the computation complexity is reduced to O() additions. But the coefficients calculated by this modified algorithm will not be
2 5 ATIOAL POWER SYSTEMS COERECE, PSC 22 same as Haar coefficients. We can call them as Modified Haar coefficients. A relationship between the Modified Haar coefficient and ourier coefficients can be established and from that the sine and cosine ourier coefficients of the input signal can be calculated. A signal vector [ X ], whose length is an integer power of two, can expressed as [ X ] [ H 6 ] [ B ] () Where H 6 is the Haar matrix and B is the Haar coefficient vector, which can be calculated as [ B ] [ H 6 ]. [ X ] (2) The Haar matrix H 6 (3) (4. ).5.5 (4.4) Here blank entries signify zeros. Like the ast ourier Transform (T), the ast wavelet Transform (WT) is a fast, linear operation that operates on a data vector whose length is an integer power of two, transforming it into a numerically different vector of the same length. The WT consists of applying above matrix [ H ] hierarchically, first to full data vector [ X ] of length (i.e. multiplying [H ] to data vector [ X ]), next all the high frequency information is shifted to bottom of the vector by using a permutation matrix. This process is applied again to the vector of low frequency information of length /2 (i.e. multiplying [H 2 ] to data vector of low frequency information), then to the vector of length /4 and so on, until only two s remains. The procedure is sometimes called a pyramidal algorithm [5], for obvious reason. The equation below makes the procedure clear: Hence by using above matrix and equation (2), the Haar coefficients vector [ B ] can be calculated and by relating these coefficients with ourier coefficients, the sine and cosine ourier coefficients can be determined. This algorithm requires order log 2 operations. The same Haar coefficients can be computed with less number of mathematical operations by factorizing above matrix as follows [] (4.) (4.2) In the above procedure the last vector is the vector of Haar coefficients [B] we got before. This procedure requires only 4(-) multiplications and 2(-) additions. ow if above procedure is followed in reverse direction without Haar coefficient vector [B], we get a matrix [H] -, which is same as the matrix [H 6 ] - mentioned before. ow if all the non-zero coefficients (.5 & -.5) of matrix [H] are replaced by ones ( & - respectively) and above procedure is followed, we will get a vector of coefficients, like [ B ]. But the coefficients of this vector will be different then coefficients of vector [ B ]. Let call this vector as modified Haar coefficient vector [B m ]. ow it will take only 2(-) additions to get these coefficients.
3 IDIA ISTITUTE O TECHOLOGY, KHARAGPUR 7232, DECEMBER 27-29, 22 5 [ ] (6.) Η a [ ] (6.2) Η a 2 [ ] (6.3 ) Η a 2 [ ] (6.4 ) Η a 4 By following the reverse procedure without involving modified Haar coefficient vector [ B m ], we will get a x matrix [ H a ] -. Hence the input vector X(m) in terms of this matrix can be expressed as [ X ] [ H a ] [ B m ] (7) Using matrix [ Ha ], a relationship between the modified Haar coefficients and ourier coefficients can be established as follows [8] ourier expansion of the periodic signal x(t) in the interval (, T) is defined as X(t) + 2 sin wt cos wt sin 2wt cos 2wt sin 5wt + 2 cos 5wt +.. (8) Where w is the fundamental angular frequency which is given by w 2πf 2π/T, and,, 2, 9,,. are the ourier coefficients. A finite series expansion is obtained by truncating the above series of equation (8) and reducing the number of terms from infinite to, where is of the form 2 n. The ourier coefficients are defined in the discrete form as follows. X ( m) m 2 m 2πm X ( m) sin 2 2πm 2 X ( m) cos m 2 X ( m) sin πm m The above equations can be written in matrix form as follows. [ ] [ S ] [ X ] (9) Where [ ] [,, 2,.. - ] T is a vector whose s are the ourier coefficients, [ X ] is the data vector, and [ S ] is a x matrix forms by values of sine and cosine functions at different sampling points, along with the factor, 2/. Considering to be even and of the form 2 n, the elements of the matrix S can be obtained as follows. S km 2 ( k + ) mπ sin 2 k mπ cos for k foroddk, ie.. k,3, KK( ) forevenk, ie.. k 2,4, KK ( 2) () Where, m,,2 (-) and S km is the element of matrix S in row k and column m. Thus all the elements of [ S ] for 6 can be obtained easily. The ourier coefficients are the s of vector [ ] representing x(t) over the interval (,T). The square of the length of this vector, [ ] T [ ] 2 k, is equal to the k mean square value of x(t) over this interval. Similarly the modified Haar coefficients are the s of another vector [ B m ] that has the same length and represents x(t) in terms of a different basis. By substituting equation (7) into equation (9) [ ] [ S ] [ H a ] [ B m ] () Hence, by calculating the matrix [ S ] from Equ. () and matrix [ H a ], the ourier coefficient vector [ ] can be determined directly from the modified Haar coefficient vector [ B m ].
4 52 ATIOAL POWER SYSTEMS COERECE, PSC 22 III. ALGORITHM OR DIERETIAL PROTECTIO O POWER TRASORER The algorithm for extracting the fundamental frequency and the harmonic s from differential current signal is based on the ast Haar Wavelet Transform (HWT) technique. Based on the relative magnitude of these fundamental and harmonic s the operating condition of the relay is decided. This algorithm for the umerical Differential protection of power transformer can be applied in the following way: The differential transformer current, whether due to fault or inrush, is represented as [6] I ατ βτ m ( τ) IDCe + IkSin( kwτ + Φk ) e k (2) Where the time constant /α and /β are associated with decaying DC and harmonic s, respectively. The current signal is sampled at regular intervals. In order to prevent aliasing effect due to sampling, the sampling rate (samples per cycle) should be greater than (2n+). Where n is the order of harmonic to be detected. Therefore, to satisfy the computational requirement the sampling rate is chosen as 6 samples per cycle, i.e. a sampling frequency of 8Hz for the 5Hz power frequency. The selection of data window is an important factor in the computation of the threshold quantities. The full cycle data window is used for an accurate estimation of the fundamental and harmonic s and the rejection of the DC s. The modified Haar coefficients of sampled current vector are obtained by using equation (5) and equation (6). These coefficients are obtained mere addition and subtraction of the current samples. The sine and cosine ourier coefficients corresponding to modified Haar coefficients are obtained by the relationship between them, which is represented by the equation (). rom these ourier coefficients the fundamental and different harmonic s of input differential current are obtained. IV. RESULTS The algorithm based on ast Haar Wavelet Transform for the numerical differential protection of power transformer, to extract fundamental and harmonic s from the differential current signal, has been programmed in the MATLAB 5.3 and tested successfully. The flow chart to find the modified Haar coefficients and relating them with ourier coefficients, to find different ourier coefficients is given in the ig.. Software on MATLAB 5.3, based on this flow chart is given in Appendix I. This software is tested for the following type of post fault test signal: i(t) I m Sin (wt - φ ) + I m2 Sin (2wt - φ 2 ) + I m5 Sin (5wt - φ 5 ) + I d e (- t / τ) (3) Where the terms used are as follows: I m maximum value of fundamental current I m2 maximum value of second harmonic current I m5 maximum value of fifth harmonic current φ angle lag from reference for fundamental φ 2 angle lag from reference for second harmonic φ 5 angle lag from reference for fifth harmonic I d maximum value of dc decaying τ time constant associated with decaying dc The post fault current signal based on which results are calculated is i(t) 2 Sin (wt - 8 )+.5 Sin (2wt - 27 )+.3 Sin (5wt - 8 ) +.7 e (- t /.) (4) The value of fundamental current obtained using the algorithm based on ast Haar Wavelet transform. I.9938 The value of second harmonic I The value of fifth harmonic I 5.32 Calculated errors and frequency response The error in the fundamental frequency, second harmonic and fifth harmonic can be calculated from the formula given below. % Error ((Actual value Calculated value) / Actual value) * (5) or the fundamental frequency : % error ((2.9938) / 2)*.3 or second harmonic : % error ((.5.556) /.5)* -.2 or second harmonic : % error ((.3.32) /.3)* -.4 The frequency response is a measure of the filtering characteristics of algorithm. The frequency response for the umerical filter using proposed algorithm, for the test signal is shown in ig. 2. requency response is curve between normalized frequency and gain, where normalized frequency is assumed frequency divided by power frequency and gain is the calculated value of current at assumed frequency divided by actual current. The response to all undesired frequencies is
5 IDIA ISTITUTE O TECHOLOGY, KHARAGPUR 7232, DECEMBER 27-29, ig.. Program flow chart for computation of I m, I m2 and I m5 nonzero, but small. The efficacy of the filter is clear from the response shown in ig. 2., that for particular frequency, the corresponding frequency is extracted and other are almost rejected. V. COCLUSIO The algorithm based on fast Haar wavelet transform, generates its coefficients using additions and subtractions of the data samples more accurately than earlier schemes. The equations for extraction of fundamental frequency and harmonic s using ast Haar Wavelet transform are derived. Software has been developed to calculate modified Haar coefficients and ourier coefficients and hence magnitude of fundamental and second and fifth harmonic of differential current. The calculated errors for different frequency of differential current are very small and do not affect the accuracy of algorithm. The frequency responses of the fundamental frequency filter, second harmonic filter and fifth filer shows good filtering characteristics of the algorithm. During magnetizing inrush and over excitation, the algorithm is flexible enough to generate the necessary restrain signal. ig. 2. requency Response of HWT Based ilter VI. REERECES [] Z. Moravej, D.. Vishwakarma, and S. P. Singh, Digital filtering algorithms for the differential relaying of power transformer: An Overview, Electrical Machines and Power Systems, Vol. 28, pp , 2. [2] M. Gomez-Morgante and D. W. icollete, A wavelet based differential transformer protection, IEEE Transactions on Power Delivery, vol. 4, o. 4, pp , 999. [3] B. Kasztenny, and M. Kezunovic, Digital relays improve protection of large transformers, IEEE Computer Application in Power, pp , Oct [4] A. D. Guzman, S. Zocholl, G. Benmouyal, and H. J. Altuve, A current based solution for transformer differential protection Part-I; Problem statement, IEEE Transactions on Power Delivery, vol. 6, o. 4, pp , Oct. 2. [5] O. Rioul and M. Vetterli, Wavelets and Signal Processing, IEEE Signal Processing Magazine, pp. 4-38, Oct. 99. [6] H. K. Verma, and A. M. Basha, Transformer differential relay based on waveform approach, The institute of Engineers (India) Journal (EL), Vol. 67, pp , Dec [7] M. A. Rahman, B. so, M. R. Zaman, and M. A. Hoque, Testing of algorithm for a stand-alone digital relaying for power transformers, IEEE Transactions on Power Delivery, vol. 3, o. 2, pp , April 997. [8] D.. Vishwakarma and B. Ram, Power System Protection and Switch Gear, Tata McGraw-Hill Publishing Company Limited, ew Delhi, 995. [9] L. P. Singh, Digital Protection Protective Relaying rom Electro - Mechanical to Microprocessor, Willey Eastern Limited, 994. [] A. G. Phadke and J. S. Thorpe, Computer Relaying for Power Systems, Research Studies Press Ltd., England 988. []
New Algorithms for Removal of DC Offset and Subsynchronous. Resonance terms in the Current and Voltage Signals under Fault.
ew Algorithms for Removal of DC Offset and Subsynchronous Resonance terms in the Current and Voltage Signals under Fault Conditions KAALESH KUAR SHARA Department of Electronics & Communication Engineering
More informationSecuring Sequence-Current Differential Elements
Securing Sequence-Current Differential Elements Gabriel Benmouyal and Tony Lee Schweitzer Engineering Laboratories, Inc. Presented at the 31st Annual Western Protective Relay Conference Spokane, Washington
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationIdentification and Classification of High Impedance Faults using Wavelet Multiresolution Analysis
92 NATIONAL POWER SYSTEMS CONFERENCE, NPSC 2002 Identification Classification of High Impedance Faults using Wavelet Multiresolution Analysis D. Cha N. K. Kishore A. K. Sinha Abstract: This paper presents
More informationOptimal PMU Placement
Optimal PMU Placement S. A. Soman Department of Electrical Engineering Indian Institute of Technology Bombay Dec 2, 2011 PMU Numerical relays as PMU System Observability Control Center Architecture WAMS
More informationWAVELET TRANSFORMS IN TIME SERIES ANALYSIS
WAVELET TRANSFORMS IN TIME SERIES ANALYSIS R.C. SINGH 1 Abstract The existing methods based on statistical techniques for long range forecasts of Indian summer monsoon rainfall have shown reasonably accurate
More informationAdaptive under frequency load shedding using synchrophasor measurement
Adaptive under load shedding using synchrophasor measurement Abstract - The imbalance between the generation and the demand is the major factor that causes instability in a power system. Conventional Under
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationQ 1 Find the square root of 729. 6. Squares and Square Roots Q 2 Fill in the blank using the given pattern. 7 2 = 49 67 2 = 4489 667 2 = 444889 6667 2 = Q 3 Without adding find the sum of 1 + 3 + 5 + 7
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationFinite Element Modeling and Performance Analysis of a Permanent Magnet Synchronous Machine
Finite Element Modeling and Performance Analysis of a Permanent Magnet Synchronous Machine Saadoun 1,. Guersi,A. Ouari 3 1- Université Badji Mohtar,, Faculté des Sciences de l Ingénieur, saadoun_a@yahoo.fr
More informationAdaptive Distance Relaying Scheme for Power Swing Tripping Prevention
Adaptive Distance Relaying Scheme for Power Swing Tripping Prevention 1 NOR ZULAILY MOHAMAD, AHMAD FARID ABIDIN, 3 ISMAIL MUSIRIN Centre of Electrical Power Engineering Studies Universiti Teknologi MARA
More informationA Machine Intelligence Approach for Classification of Power Quality Disturbances
A Machine Intelligence Approach for Classification of Power Quality Disturbances B K Panigrahi 1, V. Ravi Kumar Pandi 1, Aith Abraham and Swagatam Das 1 Department of Electrical Engineering, IIT, Delhi,
More informationReactive Power Solutions
GE Digital Energy Reactive Power Solutions Effects of Series Capacitors on Line Protection Relaying Design and Settings Presented by: Paul Datka, GE Energy Consulting Luis Polanco, GE Energy Consulting
More informationAssignment #09 - Solution Manual
Assignment #09 - Solution Manual 1. Choose the correct statements about representation of a continuous signal using Haar wavelets. 1.5 points The signal is approximated using sin and cos functions. The
More informationA Novel Method of Voltage Sag and Swell Estimation for Power System Applications
American Journal of Engineering and Applied Sciences 6 (2): 233-240, 2013 ISSN: 1941-7020 2014 A. Vetrivel et al., This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0
More informationUNIT 1. SIGNALS AND SYSTEM
Page no: 1 UNIT 1. SIGNALS AND SYSTEM INTRODUCTION A SIGNAL is defined as any physical quantity that changes with time, distance, speed, position, pressure, temperature or some other quantity. A SIGNAL
More informationMathematical Modelling of Permanent Magnet Synchronous Motor with Rotor Frame of Reference
Mathematical Modelling of Permanent Magnet Synchronous Motor with Rotor Frame of Reference Mukesh C Chauhan 1, Hitesh R Khunt 2 1 P.G Student (Electrical),2 Electrical Department, AITS, rajkot 1 mcchauhan1@aits.edu.in
More informationANALOG AND DIGITAL SIGNAL PROCESSING ADSP - Chapter 8
ANALOG AND DIGITAL SIGNAL PROCESSING ADSP - Chapter 8 Fm n N fnt ( ) e j2mn N X() X() 2 X() X() 3 W Chap. 8 Discrete Fourier Transform (DFT), FFT Prof. J.-P. Sandoz, 2-2 W W 3 W W x () x () x () 2 x ()
More informationWavelet Transform in Speech Segmentation
Wavelet Transform in Speech Segmentation M. Ziółko, 1 J. Gałka 1 and T. Drwięga 2 1 Department of Electronics, AGH University of Science and Technology, Kraków, Poland, ziolko@agh.edu.pl, jgalka@agh.edu.pl
More informationChap 4. Software Reliability
Chap 4. Software Reliability 4.2 Reliability Growth 1. Introduction 2. Reliability Growth Models 3. The Basic Execution Model 4. Calendar Time Computation 5. Reliability Demonstration Testing 1. Introduction
More informationPART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review
More informationThe Solution of Linear Systems AX = B
Chapter 2 The Solution of Linear Systems AX = B 21 Upper-triangular Linear Systems We will now develop the back-substitution algorithm, which is useful for solving a linear system of equations that has
More informationWavelets. Lecture 28
Wavelets. Lecture 28 Just like the FFT, the wavelet transform is an operation that can be performed in a fast way. Operating on an input vector representing a sampled signal, it can be viewed, just like
More informationMitigating Subsynchronous resonance torques using dynamic braking resistor S. Helmy and Amged S. El-Wakeel M. Abdel Rahman and M. A. L.
Proceedings of the 14 th International Middle East Power Systems Conference (MEPCON 1), Cairo University, Egypt, December 19-21, 21, Paper ID 192. Mitigating Subsynchronous resonance torques using dynamic
More informationFair Exponential Smoothing with Small Alpha
Fair Exponential Smoothing with Small Alpha (fessa) Juha Reunanen juha.reunanen@tomaattinen.com October 12, 2015 Abstract This short paper shows that standard exponential smoothing with a small α may not
More informationWavelets Based Identification and Classification of Faults in Transmission Lines
Wavelets Based Identification and Classification of Faults in Transmission Lines 1 B Narsimha Reddy, 2 P Chandrasekar 1,2EEE Department, MGIT Abstract: - Optimal operation of a power system depends on
More informationHigh Impedance Faults Detection Technique Based on Wavelet Transform
Proceedings of the 11th WSEAS International Conference on SYSTEMS, Agios Nikolaos, Crete Island, Greece, July 23-25, 2007 142 High Impedance Faults Detection Technique Based on Wavelet Transform MING-TA
More informationACCURACY OF APPROXIMATE FORMULAS FOR INTERNAL IMPEDANCE OF TUBULAR CYLINDRICAL CONDUCTORS FOR LARGE PARAMETERS
Progress In Electromagnetics Research M, Vol. 16, 171 184, 2011 ACCURACY OF APPROXIMATE FORMULAS FOR INTERNAL IMPEDANCE OF TUBULAR CYLINDRICAL CONDUCTORS FOR LARGE PARAMETERS D. Lovrić Faculty of Electrical
More informationChapter 2: The Fourier Transform
EEE, EEE Part A : Digital Signal Processing Chapter Chapter : he Fourier ransform he Fourier ransform. Introduction he sampled Fourier transform of a periodic, discrete-time signal is nown as the discrete
More informationMethods for the Simulation of Dynamic Phasor Models in Power Systems
Using Frequency-Matched Numerical Integration Methods for the Simulation of Dynamic Phasor Models in Power Systems Turhan Demiray, Student, IEEE and Göran Andersson, Fellow, IEEE Abstract In this paper
More informationThe Discrete Fourier Transform
In [ ]: cd matlab pwd The Discrete Fourier Transform Scope and Background Reading This session introduces the z-transform which is used in the analysis of discrete time systems. As for the Fourier and
More informationLoss of excitation protection function (40Z)
Loss of excitation protection function (40Z) Budapest, April 2016 Loss of excitation protection function The loss of excitation protection function can be applied mainly for synchronous generators. On
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More informationTransient Analysis of Three Phase Squirrel Cage Induction Machine using Matlab
Transient Analysis of Three Phase Squirrel Cage Induction Machine using Matlab Mukesh Kumar Arya*, Dr.Sulochana Wadhwani** *( Department of Electrical Engineering, Madhav Institute of Technology & Science,
More informationarxiv: v1 [math.ca] 6 Feb 2015
The Fourier-Like and Hartley-Like Wavelet Analysis Based on Hilbert Transforms L. R. Soares H. M. de Oliveira R. J. Cintra Abstract arxiv:150.0049v1 [math.ca] 6 Feb 015 In continuous-time wavelet analysis,
More informationWavelet based feature extraction for classification of Power Quality Disturbances
European Association for the Development of Renewable Energies, Environment and Power Quality (EA4EPQ) International Conference on Renewable Energies and Power Quality (ICREPQ 11) Las Palmas de Gran Canaria
More informationPerformance Evaluation and Review of System Protection Scheme Design with the help of Synchrophasor Measurements
Performance Evaluation and Review of System Protection Scheme Design with the help of Synchrophasor Measurements P. Mukhopadhyay V. Pandey Rajkumar Anumasula Chandan Kumar Sunil Patil Srinivas Chitturi
More informationCAPACITOR PLACEMENT IN UNBALANCED POWER SYSTEMS
CAPACITOR PLACEMET I UBALACED POWER SSTEMS P. Varilone and G. Carpinelli A. Abur Dipartimento di Ingegneria Industriale Department of Electrical Engineering Universita degli Studi di Cassino Texas A&M
More informationFault Location in Distribution Feeders with Distributed Generation using Positive Sequence Apparent Impedance
Fault Location in Distribution Feeders with Distributed Generation using Positive Sequence Apparent Impedance ARTURO SUMAN BRETAS Federal University of Rio Grande do Sul Department of Electrical Engineering
More informationAP PHYSICS 2012 SCORING GUIDELINES
AP PHYSICS 2012 SCORING GUIDELINES General Notes About 2012 AP Physics Scoring Guidelines 1. The solutions contain the most common method of solving the free-response questions and the allocation of points
More informationSecured Busbar Differential Protection Using A Computationally Efficient Dot Product Technique
Secured Busbar Differential Protection Using A Computationally Efficient Dot Product Technique Krish Narendra, ERLPhase Power Technologies Dave Fedirchuk, ERLPhase Power Technologies Abstract The use of
More informationEDEXCEL NATIONAL CERTIFICATE UNIT 28 FURTHER MATHEMATICS FOR TECHNICIANS OUTCOME 3 TUTORIAL 1 - TRIGONOMETRICAL GRAPHS
EDEXCEL NATIONAL CERTIFICATE UNIT 28 FURTHER MATHEMATICS FOR TECHNICIANS OUTCOME 3 TUTORIAL 1 - TRIGONOMETRICAL GRAPHS CONTENTS 3 Be able to understand how to manipulate trigonometric expressions and apply
More informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationFourier Transform Fast Fourier Transform discovered by Carl Gauss ~1805 re- re invented by Cooley & Tukey in 1965
ourier Transform ast ourier Transform discovered by Carl Gauss ~85 re-invented by Cooley & Tukey in 965 wikipedia Next concepts Shannon s Theorem ourier analysis Complex notation Rotating vectors Angular
More information1.Chapter Objectives
LU Factorization INDEX 1.Chapter objectives 2.Overview of LU factorization 2.1GAUSS ELIMINATION AS LU FACTORIZATION 2.2LU Factorization with Pivoting 2.3 MATLAB Function: lu 3. CHOLESKY FACTORIZATION 3.1
More informationSparse Solutions of Linear Systems of Equations and Sparse Modeling of Signals and Images: Final Presentation
Sparse Solutions of Linear Systems of Equations and Sparse Modeling of Signals and Images: Final Presentation Alfredo Nava-Tudela John J. Benedetto, advisor 5/10/11 AMSC 663/664 1 Problem Let A be an n
More informationStator Thermal Time Constant
Stator Thermal Time Constant Jon Steinmetz, Subhash C. Patel, and Stanley E. Zocholl Schweitzer Engineering Laboratories, Inc. 13 IEEE. Personal use of this material is permitted. Permission from IEEE
More informationEddy Current Heating in Large Salient Pole Generators
Eddy Current Heating in Large Salient Pole Generators C.P.Riley and A.M. Michaelides Vector Fields Ltd., 24 Bankside, Kidlington, Oxford OX5 1JE, UK phone: (+44) 1865 370151, fax: (+44) 1865 370277 e-mail:
More informationParameter Prediction and Modelling Methods for Traction Motor of Hybrid Electric Vehicle
Page 359 World Electric Vehicle Journal Vol. 3 - ISSN 232-6653 - 29 AVERE Parameter Prediction and Modelling Methods for Traction Motor of Hybrid Electric Vehicle Tao Sun, Soon-O Kwon, Geun-Ho Lee, Jung-Pyo
More informationA New Novel of transverse differential protection Scheme
A New Novel of transverse differential protection Scheme Li Xiaohua, Yin Xianggen, Zhang Zhe, Chen Deshu Dept of Electrical Engineering, Huazhong University of science and technology, Wuhan Hubei, 430074,
More informationOrder Tracking Analysis
1. Introduction Order Tracking Analysis Jaafar Alsalaet College of Engineering-University of Basrah Mostly, dynamic forces excited in a machine are related to the rotation speed; hence, it is often preferred
More informationEEL3135: Homework #3 Solutions
EEL335: Homework #3 Solutions Problem : (a) Compute the CTFT for the following signal: xt () cos( πt) cos( 3t) + cos( 4πt). First, we use the trigonometric identity (easy to show by using the inverse Euler
More informationCritical clearing time evaluation of Nigerian 330kV transmission system
American Journal of Electrical Power and Energy Systems 2013; 2(6): 123-128 Published online October 20, 2013 (http://www.sciencepublishinggroup.com/j/epes) doi: 10.11648/j.epes.20130206.11 Critical clearing
More informationNonlinear system identification with the use of describing functions a case study
Nonlinear system identification with the use of describing functions a case study Zhongge Zhao 1, Chuanri Li 2, Kjell Ahlin 3, Huan Du 4 1, 2, 4 School of Reliability and System Engineering, Beihang University,
More informationTransient Stability Analysis of Single Machine Infinite Bus System by Numerical Methods
International Journal of Electrical and Electronics Research ISSN 348-6988 (online) Vol., Issue 3, pp: (58-66), Month: July - September 04, Available at: www.researchpublish.com Transient Stability Analysis
More informationPowerApps Optimal Power Flow Formulation
PowerApps Optimal Power Flow Formulation Page1 Table of Contents 1 OPF Problem Statement... 3 1.1 Vector u... 3 1.1.1 Costs Associated with Vector [u] for Economic Dispatch... 4 1.1.2 Costs Associated
More informationCoupled Electrical Oscillators Physics Advanced Physics Lab - Summer 2018 Don Heiman, Northeastern University, 5/24/2018
Coupled Electrical Oscillators Physics 3600 - Advanced Physics Lab - Summer 08 Don Heiman, Northeastern University, 5/4/08 I. INTRODUCTION The objectives of this experiment are: () explore the properties
More informationTrajectory Sensitivity Analysis as a Means of Performing Dynamic Load Sensitivity Studies in Power System Planning
21, rue d Artois, F-75008 PARIS CIGRE US National Committee http : //www.cigre.org 2014 Grid of the Future Symposium Trajectory Sensitivity Analysis as a Means of Performing Dynamic Load Sensitivity Studies
More information5.1 The Discrete Time Fourier Transform
32 33 5 The Discrete Time ourier Transform ourier (or frequency domain) analysis the last Complete the introduction and the development of the methods of ourier analysis Learn frequency-domain methods
More informationSolving Linear Systems Using Gaussian Elimination
Solving Linear Systems Using Gaussian Elimination DEFINITION: A linear equation in the variables x 1,..., x n is an equation that can be written in the form a 1 x 1 +...+a n x n = b, where a 1,...,a n
More informationFourier Transform. Find the Fourier series for a periodic waveform Determine the output of a filter when the input is a periodic function
Objectives: Be able to Fourier Transform Find the Fourier series for a periodic waveform Determine the output of a filter when the input is a periodic function Filters with a Single Sinusoidal Input: Suppose
More informationThen r (t) can be expanded into a linear combination of the complex exponential signals ( e j2π(kf 0)t ) k= as. c k e j2π(kf0)t + c k e j2π(kf 0)t
.3 ourier Series Definition.37. Exponential ourier series: Let the real or complex signal r t be a periodic signal with period. Suppose the following Dirichlet conditions are satisfied: a r t is absolutely
More informationI. Signals & Sinusoids
I. Signals & Sinusoids [p. 3] Signal definition Sinusoidal signal Plotting a sinusoid [p. 12] Signal operations Time shifting Time scaling Time reversal Combining time shifting & scaling [p. 17] Trigonometric
More informationStockbridge-Type Damper Effectiveness Evaluation: Part II The Influence of the Impedance Matrix Terms on the Energy Dissipated
1470 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 18, NO. 4, OCTOBER 2003 Stockbridge-Type Damper Effectiveness Evaluation: Part II The Influence of the Impedance Matrix Terms on the Energy Dissipated Giorgio
More informationHarmonic Modeling of Networks
Harmonic Modeling of Networks Thomas H. Ortmeyer ECE Dept. Clarkson University Potsdam, NY 13699-5720 M. Fayyaz Akram Dept. of Elec. Eng. Univ. of Engineering and Technology Lahore, Pakistan Takashi Hiyama
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More information1 Unified Power Flow Controller (UPFC)
Power flow control with UPFC Rusejla Sadikovic Internal report 1 Unified Power Flow Controller (UPFC) The UPFC can provide simultaneous control of all basic power system parameters ( transmission voltage,
More informationHARMONIC WAVELET TRANSFORM SIGNAL DECOMPOSITION AND MODIFIED GROUP DELAY FOR IMPROVED WIGNER- VILLE DISTRIBUTION
HARMONIC WAVELET TRANSFORM SIGNAL DECOMPOSITION AND MODIFIED GROUP DELAY FOR IMPROVED WIGNER- VILLE DISTRIBUTION IEEE 004. All rights reserved. This paper was published in Proceedings of International
More informationMATLAB SIMULINK Based DQ Modeling and Dynamic Characteristics of Three Phase Self Excited Induction Generator
628 Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 MATLAB SIMULINK Based DQ Modeling and Dynamic Characteristics of Three Phase Self Excited Induction Generator A. Kishore,
More informationStudy of Transient Behaviour of the Capacitor Voltage Transformer
Study of Transient Behaviour of the Capacitor Voltage Transformer Amit Kumar 1, Dr. A. G. Thosar 2, Vivek Moroney 3 PG Student [EPS], Dept. of EE, Government College of Engineering, Aurangabad, Maharashtra,
More informationEffects of Capacitor Bank Installation in a Medium Voltage (MV) Substation
Effects of Capacitor Bank Installation in a Medium Voltage (MV) Substation Adesina, Lambe Mutalub Department of Engineering & Standardization, Eko Electricity Distribution Plc, 24/25, Marina, Lagos Island,
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationSIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 6a. Dr David Corrigan 1. Electronic and Electrical Engineering Dept.
SIGNALS AND SYSTEMS: PAPER 3C HANDOUT 6a. Dr David Corrigan. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.mee.tcd.ie/ corrigad FOURIER SERIES Have seen how the behaviour of systems can
More informationCyclic fatigue demands on structures subjected to the Canterbury Earthquake Sequence
Cyclic fatigue demands on structures subjected to the -11Canterbury Earthquake Sequence J.B. Mander Texas A&M University, College Station, Texas, USA. G.W. Rodgers Dept. of Mechanical Engineering, University
More informationEvolutionary Power Spectrum Estimation Using Harmonic Wavelets
6 Evolutionary Power Spectrum Estimation Using Harmonic Wavelets Jale Tezcan Graduate Student, Civil and Environmental Engineering Department, Rice University Research Supervisor: Pol. D. Spanos, L.B.
More informationCahier Technique N 13 Principe de réduction des courants d enclenchement des transformateurs
Cahier Technique N 13 Principe de réduction des courants d enclenchement des transformateurs Numerical transformer inrush current minimizer Principle of the operation Rev 1.0 Document version information
More informationTHE MEASUREMENT OF SOLAR ULTRAVIOLET SPECTRAL IRRADIANCE PROBLEMS & SOLUTIONS
THE MEASUREMENT OF SOLAR ULTRAVIOLET SPECTRAL IRRADIANCE PROBLEMS & SOLUTIONS INTRODUCTION In recent years, researchers involved in many unrelated scientific disciplines have acquired an interest in accurately
More informationChapter 12 Variable Phase Interpolation
Chapter 12 Variable Phase Interpolation Contents Slide 1 Reason for Variable Phase Interpolation Slide 2 Another Need for Interpolation Slide 3 Ideal Impulse Sampling Slide 4 The Sampling Theorem Slide
More informationCIGRE US National Committee 2015 Grid of the Future Symposium. Setting-less Protection: Laboratory Experiments and Field Trials
http : //www.cigre.org CIGRE US National Committee 2015 Grid of the Future Symposium Setting-less Protection: Laboratory Experiments and Field Trials A. P. MELIOPOULOS, G.J. COKKINIDES Georgia Institute
More informationNASA Contractor Report. Application of FEM to Estimate Complex Permittivity of Dielectric Material at Microwave Frequency Using Waveguide Measurements
NASA Contractor Report Application of FEM to Estimate Complex Permittivity of Dielectric Material at Microwave Frequency Using Waveguide Measurements M. D.Deshpande VIGYAN Inc., Hampton, VA C. J. Reddy
More informationModule 7:Data Representation Lecture 35: Wavelets. The Lecture Contains: Wavelets. Discrete Wavelet Transform (DWT) Haar wavelets: Example
The Lecture Contains: Wavelets Discrete Wavelet Transform (DWT) Haar wavelets: Example Haar wavelets: Theory Matrix form Haar wavelet matrices Dimensionality reduction using Haar wavelets file:///c /Documents%20and%20Settings/iitkrana1/My%20Documents/Google%20Talk%20Received%20Files/ist_data/lecture35/35_1.htm[6/14/2012
More informationReal-Life Observations of Power System Dynamic Phenomena Some Interesting Aspects
Real-Life Observations of Power System Dynamic Phenomena Some Interesting Aspects A.M.Kulkarni Department of Electrical Engineering, IIT Bombay and Visiting Professor, University of Manitoba The Indian
More informationTransient analysis of the behaviour of grounding systems consisted by driven rods
Transient analysis of the behaviour of grounding systems consisted by driven rods I.F. GONOS M.K. ANTONIOU I.A. STATHOPULOS F.V. TOPALIS Department of Electrical and Computer Engineering, High Voltage
More informationDriven Harmonic Oscillator
Driven Harmonic Oscillator Physics 6B Lab Experiment 1 APPARATUS Computer and interface Mechanical vibrator and spring holder Stands, etc. to hold vibrator Motion sensor C-209 spring Weight holder and
More informationA New Predictive Control Strategy Dedicated to Salient Pole Synchronous Machines
A New Predictive Control Strategy Dedicated to Salient Pole Synchronous Machines Nicolas Patin Member IEEE University of Technology of Compiègne Laboratoire d Electromécanique de Compiègne Rue Personne
More informationPower System Stability and Control. Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India
Power System Stability and Control Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India Contents Chapter 1 Introduction to Power System Stability
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationModeling Buck Converter by Using Fourier Analysis
PIERS ONLINE, VOL. 6, NO. 8, 2010 705 Modeling Buck Converter by Using Fourier Analysis Mao Zhang 1, Weiping Zhang 2, and Zheng Zhang 2 1 School of Computing, Engineering and Physical Sciences, University
More informationThe Effect of the Harmonics, the Fault Location and the Fault Resistance on the Performance of the Impedance-Type Distance Relay
American Journal of Applied Sciences 6 (4): 788-796, 2009 SSN 1546-9239 2009 Science Publications The Effect of the Harmonics, the Fault Location and the Fault esistance on the Performance of the mpedance-type
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 1 Chapter 4 Roundoff and Truncation Errors PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationDigital Power System Protection
Digital Power System Protection S.R. Bhide Associate Professor of Electrical Engineering Visvesvaraya National Institute of Technology Nagpur Delhi 110092 2014 Digital power system protection S.R. Bhide
More informationThe Structure of Digital Imaging: The Haar Wavelet Transform. Wavelet Transformation and Image Compression
The Structure of Digital Imaging: The Haar Wavelet Transformation and Image Compression December 7, 2009 Abstract The Haar wavelet transform allows images to be compressed and sent over a network of computers
More informationModeling and Simulation of Air Insulated and Gas Insulated Substations
International Journal of Electrical Engineering. ISSN 0974-2158 Volume 11, Number 2 (2018), pp. 177-187 International Research Publication House http://www.irphouse.com Modeling and Simulation of Air Insulated
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More information