This paper presents a single-phase energy metering chip with built-in calibration function to measure electric power quantities. The entire chip consists of an analog front end, a filter block, a computation engine, a calibration engine, and an external interface block. The key design issues are how to reduce the implementation costs of the computation engine from repeatedly used arithmetic operations and how to simplify calibration procedure and reduce calibration time. The proposed energy metering chip simplifies the computation engine using time-division multiplexed arithmetic units. It also provides a simple and fast calibration scheme by using integrated digital calibration functionality. The chip is fabricated with 0.18-μm six-layer metal CMOS process and housed in a 32-pin quad-flat no-leads (QFN) package. It operates at a clock speed of 4096 kHz and consumes 9.84 mW in 3.3 V supply.
s the energy consumption in a household continues to increase, the need to precisely measure and control alternating current (AC) power has increased. Recently, energy metering technology has been applied to various smart home devices and appliances in order to monitor and control energy usage
. An energy metering chip measures 110 V or 220 V AC voltage and current applied to a load, and it calculates the corresponding power consumption and energy. In addition, it also monitors a variety of conditions, such voltage sag or swell, over-current and so on.
The energy metering chip is interface to both voltage and current sensing elements. A resistor divider is widely used as the voltage sensor. A low resistance shunt or a current transformer (CT) is commonly used as the current sensor. The current shunt offers good accuracy at low cost, and the measurement is simple. However, its parasitic inductance and self-heating problem makes it unsuitable for large current measurement
. The CT can measure up to very high current and consumes little power. However, it typically has a small phase shift of between 0.1˚ and 0.3˚ because of the magnetizing current
. These sensing elements and an associated RC interface circuits affect the magnitude and phase of the sensed signals, and results in the degradation of measurement accuracy. In addition, component tolerances, residual analog-to-digital converter (ADC) offset, and system noise also degrade the performance. Therefore, the various system calibrations, such as offset correction, phase calibration, and gain calibration, are required for precise measurements.
The energy metering chip can be implemented with a microcontroller or in hardwired logic. The microcontroller-based implementation of the energy metering chip is easier and faster compared to the fully hardwired implementation. The ease of addition and modification of functionality through embedded S/W upgrade is another advantage of the microcontroller-based implementation, making it suitable for various applications. However, the power consumption and implementation costs of the microcontroller-based solution are known to be higher than its hardwired counterpart. Depending on the type of application, the energy metering chip needs to operate non-stop for long hours, making low-energy consumption essential. In addition, systems using the energy metering chip require system calibration for accurate energy measurement. The chip’s system calibration is a time consuming and complicated process which affects mass production. Therefore, an effective system calibration method is necessary.
There are several issues in the design and implementation of the fully hardwired energy metering chip. First, the arithmetic operations such as multiplication and division are repeatedly used to calculate the electric power quantities, which complicate the chip’s computation block and increase hardware (H/W) implementation costs. Therefore, an efficient architecture is required to simplify the computation block using time-division multiplexed arithmetic units. Second, an accurate source is commonly required for the system calibrations to be performed accurately. The source generates any desired voltage, current, and phase shifts between the voltage and current. However, it is usually expensive and is not easily available. Calibration method using only a purely resistive load can provide simplicity and convenience, although it has less accuracy compared with that using an accurate source. Last, the system calibrations are normally performed with a personal computer (PC) based calibration controller. The calibration parameters such as the gain, offset, and phase are calculated in the calibration controller and are delivered to the energy metering chip. This method increases calibration time and complicates calibration procedure. Hence, a built-in calibration function is required for a simple and fast calibration.
In this paper, an efficient energy metering chip architecture to achieve low H/W complexity as well as simple and fast calibration is proposed. The remainder of this paper is organized as follows. The system overview and power definitions are described in Section 2, and the computation engine architecture and the calibration scheme are described in Section 3. Implementation results and conclusions are given in Sections 4 and 5, respectively.
2. System Description
In this section, a brief description of a single-phase energy meter is given and power definitions are described in order to understand the architecture and calibration scheme proposed. In
, the list of symbols and definitions used across the paper is presented for reference purposes.
Symbols and definitions
The single-phase energy meter is composed of a voltage sensor, a current sensor, an energy metering chip, and a micro-controller as shown in
. The AC voltage and current applied to the load are input to the energy metering chip through the sensing elements. The chip converts the analog signals to digital samples, and then performs the process of calculating the root mean square (RMS) and power using the sensed samples. It also stores the calculated RMS and power measurements into the registers and communicates with the external micro-controller via the external interface block.
Single-phase energy meter
The IEEE standard 1459 shows some definitions to measure electric power quantities in a single-phase system under sinusoidal condition
. The instantaneous voltage and current can be expressed as shown in (1) and (2), respectively.
are the RMS values of the voltage and current signals, respectively.
is the angular frequency 2
is the line frequency, e.g. 50 or 60 Hz, and
is the phase angle between the current and voltage signals and
is the time.
The active power is the average value of the instantaneous power which is the product of the instantaneous voltage and the instantaneous current during the observation time interval from
, and it is given by
is the moment when the time starts,
is a positive integer number, and
is the line cycle time.
The reactive power is the average value of the instantaneous reactive power which is obtained by multiplying the instantaneous current by the instantaneous quadrature voltage. The instantaneous quadrature voltage is obtained by phase shifting instantaneous voltage
/2 using an integrator with
gain. The reactive power is defined by
The apparent power is defined as the maximum power that can be delivered to a load, and it is the product of the RMS voltage and the RMS current as shown in (5).
The power factor is the ratio between the energy transmitted to the load over the maximum energy. It is defined as the active power divided by the apparent power as shown in (6).
The RMS value of the voltage is defined as shown in (7), and the RMS value of the current is obtained by (7).
3. Energy Metering Chip Architecture
In this section, the proposed energy metering chip architecture is described. A brief description of the architecture is given, and a computation engine architecture and a calibration scheme are proposed.
shows the block diagram of the proposed energy metering chip. It includes an analog front end (AFE), a filter block, a computation engine (CE), a calibration engine (CAL), internal registers, and an external interface block. The AFE incorporates two programmable gain amplifiers (PGAs) and two third-order ΣΔ modulators. The PGAs are interface to the voltage and current sensors and independently amplify the sensed voltage and current, respectively. Both ΣΔ modulators simultaneously convert the analog voltage and current signals to single-bit digital data streams, and sample at a rate of 1024 kHz. The ΣΔ modulator is implemented using a cascaded integrator feed-forward (CIFF) architecture.
Block diagram of energy metering chip
The filter block consists of two decimation filters, two high-pass filters (HPFs), and a low-pass filter (LPF). The single-bit ΣΔ modulator output data streams are widened to 24-bit words and down-sampled to 8 kHz with both decimation filters. The decimation filter includes a fourth-order cascaded integrator-comb (CIC) filter and a compensation filter. The CIC filter is widely used as the decimation filter because it requires no multiplication or coefficient storage but rather only additions and subtractions
. The CIC and compensation filters are used to achieve sampling rate decrease and compensate for the amplitude roll-off of the CIC filter, respectively. The compensation filter is implemented using a finite impulse response (FIR) filter. Both HPFs remove any DC component from the voltage and current signals. The LPF removes noise and harmonic components included in the voltage signal, and it results in a stable and accurate frequency measurement tolerant of noise and harmonics. The HPF and LPF are implemented using first-order infinite impulse response (IIR) filters with a cutoff frequency of 4.45 Hz and 140 Hz, respectively.
The CE performs power measurements, line-voltage period measurement, and RMS calculation on the voltage and current. The CAL performs the precise estimate of errors in the gain and phase of the sensed voltage and current signals. The measured values generated from the CE are stored in registers. The external interface block communicates with the external micro-controller via such as an Inter-Integrated Circuit (I2C), a universal asynchronous receiver/transmitter (UART), or a Serial Peripheral Interface (SPI).
- 3.1 Computation Engine
From (7), the RMS value of the voltage or current signal is calculated by squaring the signal, taking the average, and obtaining the square root. The integration in continuous time is equivalent to the discrete time accumulation or summation. Therefore, the RMS calculation requires four arithmetic operations such as multiplication, accumulation, division, and square root operation. Arithmetic operations for the active, reactive, and apparent powers and the power factor can be obtained from (3) to (6), respectively. A line frequency measurement based on the zero-crossing detection of the voltage signal is given by
is the measurement count,
is the sampling frequency,
) is the number of zero-crossings for the voltage signal, and
) is the sample count during the zero-crossing detection. The line frequency measurement requires a multiplication and a division since 2
) is obtained by shifting 1-bit to the left. The arithmetic operations for the CE require six multiplications, four accumulations, six divisions, and two square roots, as shown in
Arithmetic operation for computation engine
Arithmetic operation for computation engine
The proposed CE architecture using time-division multiplexed arithmetic units is shown in
. It requires two multipliers, four accumulators, a divider, a square root block, seven multiplexers (MUXs), and temporary registers. The two multipliers, the divider, and the square root block are time-division multiplexed to perform RMS and power calculations. The values generated from one of the arithmetic units are temporarily stored in registers, and then are used in the other arithmetic unit.
Computation engine architecture
The proposed CE operates as follows. An instantaneous voltage, an instantaneous current, and an instantaneous quadrature voltage are sequentially input to the multiplier1 through two MUXs every sampling clock. The multiplier1 sequentially calculates four values: squared voltage, squared current, instantaneous active power, and instantaneous reactive power. These four values are independently accumulated in four accumulators during the averaging duration which is the sample count
) divided by the sampling frequency
, and then four accumulated samples are sequentially input to the divider once every averaging period. All accumulators are reset when the new averaging operation begins. The divider generates an averaged squared voltage, an averaged squared current, active power, and reactive power using the output of all accumulators and the sample count
) generated from the sample and zero-crossing counter. The square root block generates the RMS values of the voltage and current using the output of the divider. The multiplier2 generates the apparent power using the output of the square root block. The divider generates the power factor using the active power stored in the register and the output of the multiplier2. It also generates a temporary value
) for line frequency. The multiplier2 generates the line frequency.
For low H/W implementation costs, the multiplier, the divider, and the square root block are implemented using a shift-and-add multiplication algorithm, a non-restoring division algorithm, and a non-restoring square root algorithm, respectively
. The conventional CE using individual arithmetic units requires six multipliers, four accumulators, six dividers, and two square roots, whereas the proposed CE requires two multipliers, four accumulators, a divider, a square root block, seven multiplexers (MUXs), and temporary registers. The implementation result shows 33.99 % reduction in equivalent gate count compared with the conventional CE. Details of the result are given in
Comparison of conventional and proposed computation engine
1Gate refers to a 2-input NAND gate
- 3.2 Calibration Engine
A power factor is commonly used to calibrate the phase shift caused by a current sensor since it is equal to the cosine of the phase angle between the voltage and current under sinusoidal condition. The power factor is insensitive to phase changes at high power factor. Therefore, if the phase calibration is performed at a power factor of 1, the accuracy may degrade more compared to that performed at a power factor of 0.5. For example, a power factor error of 0.1% results in phase error of 2.5626˚ at a power factor of almost 1 because cos
(0.999) is equal to 2.5626˚, whereas the phase error is equal to 0.0661˚, that is cos
(0.499), at a power factor of 0.5. However, there may be more power factor errors at low power compared to that at high power factor, since the value of the active power is small at low power factor. For this reason, the phase calibration is usually performed at a power factor of 0.5 with an accurate source that is able to generate any desired voltage, current, and phase shifts between the voltage and current.
The proposed calibration scheme can be performed with a purely resistive load at a power factor of 1.
shows the proposed calibration scheme. It requires an additional variable phase shifter to change phase shifts between the voltage and current. The CAL calculates the calibrating gains and delays for both voltage and current channels using the measured values generated from the CE. The calibrating gains and delays are used to compensate both the voltage and current signals.
Gain and phase calibration scheme
The proposed phase calibration procedure is as follows. First, the variable phase shifter is used to set the phase of the instantaneous voltage to
degree, e.g. 60˚. After CE operation, a measured power factor and a measured line frequency are stored in the register. The CAL estimates the phase difference between the voltage and current and calculates the corresponding delay using the measured power factor and the measured line frequency as shown in (9) and (10), respectively.
is an iteration count,
) is the
-th measured power factor,
is the phase shift generated by the variable phase shifter,
) is the
-th measured line frequency, and
is the sampling frequency used for the AFE, i.e. 1024 kHz. If the absolute value of the estimated phase is greater than the predefined tolerance value and the iteration count is less than the predefined maximum iteration value, the iteration count is increased by one. And if the estimated phase is a positive value, the calibrating delay for the voltage channel is updated. Otherwise, the calibrating delay for the current channel is updated. Both delays are updated as follows:
where initial delays,
(0), are equal to 0 or the last calibrating delays. Then, the calibration procedure is performed repeatedly until the absolute value of the estimated phase is less than the tolerance value or the iteration count is equal to the predefined maximum iteration value.
The gain calibration is similar to the phase calibration. After CE operation, the measured voltage RMS is stored in the register. The CAL estimates the voltage RMS error and calculates the calibrating gain as shown in (13) and (14), respectively.
is the expected RMS value, and
) is the
-th measured voltage RMS value. An initial gain g
(0) is equal to 1 or the last calibrating gain. The gain calibration procedure is performed repeatedly until the estimated error is less than the tolerance value or the iteration count is equal to the predefined maximum iteration value. The current gain calibration is simultaneously performed with the same method of the gain calibration.
The output of the variable phase shifter is obtained by
is the phase shift generated by the variable phase shifter,
is the line frequency, and
is the time. The variable phase shifter includes two multipliers, an adder, and a π/2 phase shifter as shown in
. The π/2 phase shifter is implemented using the Simpson integrator as shown in
and its transfer function is given by (16)
) is the
-th estimated line frequency and
is the sampling frequency which is a fixed value of 8 kHz.
) is obtained by the CE.
Variable phase shifter
Simpson integrator based π/2 phase shifter
The proposed CAL architecture consists of an arccosine block, a multiplier, a divider, adders, MUXs, registers, and a control block as shown in
. The arccosine block is implemented using a coordinate rotation digital computer (CORDIC) algorithm
. The multiplier and the divider are time-division multiplexed to calculate the calibration parameters such as the gains and delays. These parameters are stored in registers and then used for gain and phase compensations. The implementation results shows 27.15 % reduction in equivalent gate count compared with the coneventional CAL using individual arithmetic units. Details of the results are given in
. The proposed calibration method can be performed only with a purely resistive load, which provides simplicity and convenience, although it has less accuracy compared with that using an accurate source. In addition, the calibration parameters such as the gain and phase are calculated in the built-in calibration engine, which decreases calibration time and simplifies calibration procedure.
Calibration engine architecture
Comparison of conventional and proposed calibration engine
1Gate refers to a 2-input NAND gate
4. Implementation Results and Performance Analysis
- 4.1 Implementation Results
The chip micrograph is shown in
. The chip is fabricated with 0.18-μm six-layer metal CMOS process and housed in a 32-pin QFN package. It operates at a clock speed of 4096 kHz, and consumes 9.84 mW in 3.3 V supply. The implementation details and results are given in
. The proposed energy metering chips is composed of an AFE and a digital block. The AFE is implemented with two PGAs, two ΣΔ modulators, and a low-dropout (LDO) regulator. The digital block includes a filter block, a CE, a CAL, internal registers, a UART, an I2C, and a SPI. The total gate count is 135 k. The gate count ratios of the CE and the CAL are 15.1 % and 10.7 % of the whole gate count, respectively.
shows the implementation details for the digital block.
Implementation results for energy metering chip
Implementation results for energy metering chip
Implementation results for digital block
1Equivalent Gate refers to a 2-input NAND gate
A signal-phase energy meter is implemented using the proposed chip as shown in
. It consists of a CT, a voltage transformer (VT), the proposed energy metering chip, a switched mode power supply (SMPS), a regulator, and an 8-pin connector. The CT and the VT are used as the current and the voltage sensors, respectively. The SMPS and the regulator are used to generate 3.3 V supply. The chip communicates with an external Wi-Fi module through the 8-pin connector. To verify the functionality and reliability of the proposed energy metering chip, it is tested in the laboratory as shown in
. The energy meter measures 220 V AC voltage and current applied to the loads, which are two 100 W and a 30 W light bulbs, and calculates the corresponding electric qualities. The values calculated from the energy meter are transported to the monitoring computer via Wi-Fi. The proposed energy metering chip was showed successful operation.
Implementation of the energy meter using the proposed chip
Test environment for the proposed energy metering chip
- 4.2 Performance Analysis
Simulations are performed to demonstrate the effectiveness of the proposed phase calibration scheme using the variable phase shifter. A 60 Hz frequency voltage and current signal is generated. The phase difference between the voltage and current is included. The voltage and current are sampled at 8 kHz. Taking the signal values sampled at 8 kHz for 1 second, the phase difference between the voltage and current is estimated using (9).
shows an absolute value of estimated phase errors (EPEs) at 0.1, 0.5, 1.0, and 2.0 phase differences (PDs) between the voltage and current. In this figure, solid lines and dashed lines denote the EPE performances of the conventional and proposed schemes as a function of the SNRs, respectively. The proposed scheme uses the variable phase shifter to rotate the phase of the voltage to 60 ˚. The EPE of the proposed scheme slightly increases with low SNR regardless of PD while the EPE of the conventional scheme dramatically increases in less than 70 dB SNR at low PD. Thus, the proposed scheme is less sensitive to PD and noise. The conventional scheme becomes more sensitive to noise as PD to be estimated becomes smaller.
Simulated phase error at various phase differences
To verify the accuracy of the proposed energy metering chip, relative errors of the current RMS, active power, and reactive power are measured. The relative error is defined by
is the reference value and
is the measured value. The measurement accuracy is tested on the energy metering chip calibrated with the proposed gain and phase calibration scheme. A programmable signal generator is used to supply the energy metering chip with known voltages and currents at various power factors. The signal generator creates 2 signals: a 60 Hz sine wave for the voltage channel and a 60 Hz sine wave for the current channel. The relative error for the current RMS, active power, and reactive power measured by the energy metering chip is calculated by varying the size of the sine wave for the current channel.
shows the measured relative error of the current RMS. The measured error is less 0.1% over a 2000:1 dynamic range. The measured relative error of the active power is shown in
. The solid line represents the relative error at a power factor of 1 and the measured error is less 0.1% over a 2000:1 dynamic range. The dash-dot and dashed lines show errors at lagging and leading power factors of 0.5, respectively. The measured errors are less 0.1% over a 1000:1 dynamic range and within 0.2% over a 2000:1 dynamic range.
shows the measured error of the reactive power. The solid line indicates the relative error at a power factor of 0. The dash-dot and dashed lines show errors at lagging and leading power factors of 0.87, respectively. The measured errors are less 0.1% over a 2000:1 dynamic range.
show the measured power spectral density and the measured signal-to-noise and distortion ratio (SNDR) of the ΣΔ modulator, respectively. The ΣΔ modulator achieves signal-to-noise ratio (SNR) of 88 dB and peak SNDR of 78 dB. Measurement results show that the proposed energy metering chip guarantees 0.1% measurement error precision for the current RMS, active power, and reactive power in current dynamic range of 1000:1. The wide current dynamic range allows a large variety of choices for the current sensor that interfaces with the chip, allowing simplification of the interface circuit.
Measured current RMS accuracy
Measured active power accuracy
Measured reactive power accuracy
Measured power spectral density
In this paper, a single-phase energy metering chip with built-in calibration function was proposed. The proposed energy metering chip is fully implemented in hardwired logic and simplifies the computation engine using time-division multiplexed arithmetic units. In addition, it provides a simple and fast calibration scheme by using integrated digital calibration functionality. The chip was fabricated with 0.18-μm six-layer metal CMOS process and packaged in a 32-pin QFN package. It operates at a clock speed of 4096 kHz and consumes 9.84 mW in 3.3 V supply. Complexity analysis results show that the proposed computation engine and calibration engine achieved 33.99% and 27.15% reduction in equivalent gate count, respectively, compared to the computation engine and calibration engine that use individual arithmetic units. In addition, measurement results show that the relative errors of the current RMS, active power, and reactive power are less than 0.1% over a 1000:1 current dynamic range. As a future work, efforts will be made in hardware implementation cost reduction by merging the calibration engine into the computation engine using time-division multiplexed arithmetic units. In addition, the performance of sigma-delta modulator and PGA needs to be improved to secure the energy metering chip’s precision and wide current dynamic range.
Youn-Sung Lee received the B.S. and M.S. degrees in computer engineering from Kwangwoon University, Seoul, Korea, in 1996 and 1998, respectively. From 1998 to 2003, he worked for Setri Microsystems Corp., Seongnam, Korea. Since 2003, he has been with Advanced Mobile Technology Research Center in Korea Electronics Technology Institute (KETI), Seoul, Korea. He is also currently working toward to Ph.D. degree in the school of Electrical and Electronic Engineering, Yonsei University, Seoul, Korea. His research interests are in the areas of low-power, high-performance VLSI architectures and circuit design for digital communications and digital signal processing.
Jeongwook Seo received the B.S. and M.S. degrees from the Department of Telecommunication and Information Engineering, Korea Aerospace University, Gyeonggi, Korea, in 1999 and 2001, respectively, and the Ph.D. degree from the School of Electrical and Electronic Engineering, Yonsei University, Seoul, Korea, in 2010. He is currently an Assistant Professor at the Department of Information and Communication Engineering, Namseoul University, Chungcheongnam-do, Korea. From 2001 to 2013, he was with Network Convergence Research Center in Korea Electronics Technology Institute, Seoul, Korea. His research interests include statistical signal processing for communications and networking, protocol design for Machine-to-Machine/Internet of Things, and software development of next-generation broadcasting and communication systems.
Jungwook Wee received the B.S. and M.S. degrees in electrical engineering from Chung-Ang University, Seoul, Korea, in 1999 and 2001, respectively, and Ph.D. degrees in the school of electrical and electronics engineering from Chung-Ang University in 2011. Since 2001, he has been with Contents Convergence Research Center in Korea Electronics Technology Institute (KETI), Seoul, Korea. His research interests are in the areas of signal processing for digital communication and wireless/wired communication systems.
Mingoo Kang received the B.S., M.S., and Ph.D. degrees in electronic engineering from Yonsei University, Seoul, Korea, in 1986, 1989 and 1994, respectively. He was a research engineer at Samsung Electronics from 1985 to 1997. He is in Hanshin University, Osan, Korea, from 2000 as a Professor in the Division of Inform. & Telecomm. His research interests include wireless communication devices, and Smart UX for Mobile TV & DTV.
Dong Ku Kim received the B. Eng. Degree from Hankuk Aviation University, Korea, in 1983 and the M.Eng and Ph.D. degrees from the University of Southern California, Los Angeles, in 1985 and 1992, respectively. He was a Research Engineer with the Cellular Infrastructure Group, Motorola by 1994, and he has been a professor in the school of Electrical and Electronic Engineering, Yonsei University, Seoul, since 1994. He is a director of Radio Communication Research Center at Yonsei University. He has been also a director of Qualcomm Yonsei CDMA Joint Research Lab since 1999. His main research interests are multi-user MIMO, OFDM, scheduling schemes, mobile multihop relay, and location finding algorithms.
Chin W. H.
“Smart grid communications: Overview of research challenges, solutions, and standardization activities,”
IEEE Commun. Surveys Tuts.
DOI : 10.1109/SURV.2011.122211.00021
“Design and Implementation of a High-Fidelity AC Metering Network,”
Information Processing in Sensor Networks
“Current sensing for energy metering,”
in Proc. of IIC-China Conference
IEEE Standard Definitions for the Measurement of Electric Power Quantities Under Sinusoidal, sinusoidal, Balanced, or Unbalanced Conditions
IEEE Std 1459-2010
Hogenauer E. B.
“An Economical Class of Digital Filters for Decimation and Interpolation,”
IEEE Trans. Acoust., Speech, Signal Processing
DOI : 10.1109/TASSP.1981.1163535
“A Low-Power Low-Area Multiplier Based on Shift-and-Add Architecture,”
IEEE Tran. Very Large Scale Integr. (VLSI) Syst.
DOI : 10.1109/TVLSI.2008.2004544
Marimuthu C. N.
Thangaraj D. P.
“Low Power Shift and Add Multiplier Design,”
Int’l J. Computer Science and Information Technology
“Division and Square Root: Digit-Recurrence Algorithms and Implementations,”
Oberman S. F.
Flynn M. J.
"Division Algorithms and Implementations,"
IEEE Trans. Computers
DOI : 10.1109/12.609274
“A New Non-Restoring Square Root Algorithm and Its VLSI Implementations,”
IEEE International Conference on Computer Design
“Implementation of Single Precision Floating Point Square Root on FPGAs,”
Proc. of the IEEE Symposium on FCCM
Al-Alaoui M. A.
“Novel digital integrator and differentiator,”
DOI : 10.1049/el:19930253
Bihan J. L.
“Novel class of digital integrators and differentiators,”
DOI : 10.1049/el:19930647
Hu Y. H.
“CORDIC-based VLSI Architectures for Digital Signal Processing,”
IEEE Signal Processing Magazine
DOI : 10.1109/79.143467