Digital two way automatic communication system (TWACS) outbound receiver and method

BACKGROUND OF THE INVENTION

p-0002This invention relates to communications sent through a power distribution network and, more particularly, to a point-to-point communications system by which information is readily transmitted from any one location within the power distribution network to any other location within the network. Specifically, the present invention relates to digitally demodulating TWACS® (Two Way Automatic Communication System) outbound signals including outbound preamble detection.

p-0003Power line communications systems are known in the art. A typical system enables a utility to send messages over its power line to or from a central location such as a substation to most, if not all, of its customers connected to that site. The messages involve such things as current electrical usage at the customer's site, polling requests to determine whether or not an outage has occurred within a service area, or commands to reduce or shut-off the amount of power provided to a load at the customer's site during periods of peak electrical usage. Replies received from the various locations to which messages are sent enable the utility to determine its current operational status, as well as changes that may need to be made to reconfigure the power distribution system for changes (or prospective changes) in its operating circumstances.

p-0004Electrical usage has grown significantly in recent years so that, over time, the demands placed on utilities have greatly increased and many utilities are now hard pressed to maintain adequate levels of service to their customers. Similarly, the demands placed on current communications systems employed by these utilities to support their operations have also greatly increased to the point where it has become difficult for these systems to timely provide the information necessary for the utility to operate at the level at which it needs to operate. For example, the amount of information required by the utility, on an almost continuous basis, has expanded to the point where the information throughput (data transmission rates) required of communications systems is at, or near the limits of the communications system's capabilities.

p-0005TWACS outbound signals have been decoded using comparators and timers. In this way the digital processor is supplied a sequence of times at which the incident signal has taken on known values. This has the advantage of not requiring expensive hardware.

p-0006The present analog outbound receiver requires a large signal-to-noise ratio (SNR) in order to maintain a reasonable bit error ratio. It is designed to work with a specific outbound signal constellation. This constellation contains only two signals, each with a duration of two complete mains cycles. Thus, the data rate is only 1 bit every 1/30 seconds, or 30 bps. This constellation size can be increased to as many as 16 symbols for this signal duration, so that up to 4 bits may be transmitted every 2 cycles giving a maximum data rate of 120 bps. However, adapting the present receiver to the new constellation would be difficult if not impossible.

p-0007There is a need for an outbound receiver using analog-to-digital conversion that could be adapted to the new constellation and which could be implemented on general purpose microprocessors.

SUMMARY OF THE INVENTION

p-0008In one form, the invention comprises a receiver for a two-way automatic communications system (TWACS) used by an electrical utility in which outbound messages are sent from the utility to a consumer and inbound messages are sent from the consumer to the utility. The respective outbound and inbound messages are sent and received over the utility's power distribution system. The receiver detects the outbound messages and comprises an analog filter, an A/D converter and a digital processor. The analog filter component receives the analog outbound signal (AOS) r(t) including an outbound message ({circumflex over (m)}) to provide a filtered signal. The A/D converter converts the filtered signal to a corresponding digitized signal and provides the corresponding digitized signal. The digital processor executes a digital preprocessor component for receiving the corresponding digitized signal and for preprocessing the received corresponding digitized signal to provide a preprocessed digital signal z[k]. The processor also executes a digital symbol determination component for receiving the preprocessed digital signal z[k] and for digitally demodulating the received preprocessed digital signal z[k] to provide the outbound message ({circumflex over (m)}).

p-0009In another form, the invention comprises a two-way automatic communications system (TWACS) used by an electrical utility including a substation sending analog outbound messages from the utility to a consumer, a transceiver sending inbound messages from the consumer to the utility, and a receiver included in the transceiver for enabling the detection of the outbound messages. The receiver comprises an analog filter component, an A/D converter and a digital processor executing a digital preprocessor component and a digital symbol determination component.

p-0010In another form, the invention comprises a method in a two-way automatic communications system (TWACS) used by an electrical utility in which analog outbound messages are sent from the utility to a consumer and inbound messages are sent from the consumer to the utility, the respective outbound and inbound messages being sent and received over the utility's power distribution system. The method detects the outbound messages and comprises: <ul><li id="ul0001-0001" num="0000"><ul><li id="ul0002-0001" num="0010">receiving the analog outbound signal including an outbound message;</li><li id="ul0002-0002" num="0011">converting the received analog outbound signal to a corresponding digitized signal;</li><li id="ul0002-0003" num="0012">preprocessing the corresponding digitized signal to provide a preprocessed digital signal; and</li><li id="ul0002-0004" num="0013">digitally demodulating the preprocessed digital signal to provide the outbound message.</li></ul></li></ul>

p-0011Other objects and features will be in part apparent and in part pointed out hereinafter.

BRIEF DESCRIPTION OF THE DRAWINGS

p-0012<figref idrefs="DRAWINGS">FIG. 1</figref> is a simplified schematic of a point-to-point communications system incorporated into a utility's power distribution network. The box labeled X_Nrepresents one of a plurality of outbound signal transmitters. The box labeled Y_Nrepresents one of a plurality of outbound receivers. <figref idrefs="DRAWINGS">FIG. 1</figref> illustrates the outbound receiver on the LV (low voltage) side of the distribution transformer. It is also contemplated that the receiver may be placed instead or additionally on the MV (medium voltage) side (not illustrated in <figref idrefs="DRAWINGS">FIG. 1</figref>).

p-0013<figref idrefs="DRAWINGS">FIG. 2</figref> is a waveform diagram illustrating one embodiment of an outbound TWACS (Two Way Automatic Communication System) signal, which is digitally demodulated by a receiver and method of the invention.

p-0014<figref idrefs="DRAWINGS">FIG. 3</figref> is a flow chart detailing the interaction between the preamble detector and the receiver.

p-0015<figref idrefs="DRAWINGS">FIG. 4</figref> is a waveform diagram illustrating one embodiment of a 30 bps TWACS outbound signal set.

p-0016<figref idrefs="DRAWINGS">FIG. 5</figref> is a block diagram of one embodiment of a receiver according to the invention.

p-0017<figref idrefs="DRAWINGS">FIG. 6</figref> illustrates unfiltered preamble coefficients P_kas seen at the output of the FFT of <figref idrefs="DRAWINGS">FIG. 5</figref> and filtered preamble coefficients q[k] as seen at the output of FIR filter −1+z⁻²of <figref idrefs="DRAWINGS">FIG. 5</figref>.

p-0018<figref idrefs="DRAWINGS">FIG. 7</figref> illustrates a block diagram of one embodiment of the digital phase locked loop (PLL) of <figref idrefs="DRAWINGS">FIG. 5</figref>.

p-0019<figref idrefs="DRAWINGS">FIG. 8</figref> illustrates one embodiment of the frequency and impulse response of the loop filter H(z) of the digital PLL of <figref idrefs="DRAWINGS">FIG. 7</figref>.

p-0020Corresponding reference characters indicate corresponding parts throughout the drawings.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

p-0021As shown in <figref idrefs="DRAWINGS">FIG. 1</figref>, a point-to-point communications system <b>20</b> enables messages to be transmitted from any one location A within the power distribution system to any other location B within the system. Typically, messages are sent from a substation to one or more of the facilities and a reply message is separately sent from each facility back to the substation. It is a feature of the present invention, however, that a message can be sent from locations other than a substation to any other location (which could be, but is not necessarily a substation) within communications system <b>20</b>.

p-0022As shown in <figref idrefs="DRAWINGS">FIG. 1</figref>, a voltage waveform W_Ggenerated or propagated by the utility (i.e., the mains signal) is impressed across the primary windings of a high voltage transformer T_HV. Waveform W_Gis typically a 3 phase, usually between 35 KVAC and 160 KVAC, 60 Hz waveform; although it will be understood by those skilled in the art that communications system <b>20</b> works equally as well with other utility generated waveforms, for example, 50 Hz waveforms generated by utilities in many countries.

p-0023The secondary windings of transformer T_HVare, in turn, connected across the primary windings of transformers T_I-T_N. A transmitter X_Nof a transceiver indicated generally 12 is connected across the secondary or low voltage windings LV of transformer T_lat a location A. At a location B, a transceiver <b>13</b> includes a receiver Y_Nconnected to the power line for receiving and processing messages sent over the communications system <b>20</b>.

p-0024A signal waveform W_Mshown in <figref idrefs="DRAWINGS">FIG. 1</figref> represents the waveform W_Gincluding a modulated message added by transceiver <b>12</b> for receiving and decoding by transceiver <b>13</b>. TWACS transmissions are phase-locked to the mains signal MS as shown in <figref idrefs="DRAWINGS">FIG. 2</figref>. In any given half-cycle of this signal the transmitter has the option to fire its signal. The resulting current may be represented by the function g(t). The inbound and outbound signal constellation for any TWACS transmission is then summarized by the equation

p-0025<maths id="MATH-US-00001" num="00001"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>s</mi><mi>m</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>⁢</mo><mrow><msup><mrow><msub><mi>c</mi><mi>mk</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mrow><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mi>k</mi></msup><mo>⁢</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>1.1</mn><mo>⁢</mo><mi>.1</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where c_mkε{0,1}, g(t) is the TWACS pulse with support on [0, T_p], T_pis the duration of one half-cycle, and K is the signal duration in half-cycles. c_mk=1 indicates that the firing element is active during the kth half-cycle of the mth symbol, where c_mk=0 indicates that the element has not fired. g(t) is therefore the signal generated by the element firing.

p-0026Outbound TWACS uses binary signaling with four half-cycle symbol lengths. The constellation is depicted in <figref idrefs="DRAWINGS">FIG. 4</figref> where s₁(t) and s₂(t) correspond to binary 1 and 0, respectively. Thus, each symbol may be represented using (1.1.1) where <br /><i>c</i>₁=[1 0 0 0]^T<br /><i>c</i>₂=[0 0 1 0]^T  (2.1.1)

p-0027Code division multiple access (CDMA), although not frequently used explicitly for outbound TWACS signaling, is frequently used for inbound TWACS signaling.

p-0028According to one embodiment, the invention comprises a digital implementation of an additive white Gaussian noise (AWGN) receiver including a processor executing an algorithm for preamble detection and an algorithm for phase locking to the mains signal.

p-0029A receiver for dealing with AWGN is given by the following:

p-0030<maths id="MATH-US-00002" num="00002"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><mo>-</mo><msub><mi>r</mi><mn>0</mn></msub></mrow><mo>+</mo><msub><mi>r</mi><mn>2</mn></msub></mrow><mo>⁢</mo><munderover><munder><mo><</mo><mo>></mo></munder><mrow><mover><mi>m</mi><mo>^</mo></mover><mo>=</mo><mn>2</mn></mrow><mrow><mover><mi>m</mi><mo>^</mo></mover><mo>=</mo><mn>1</mn></mrow></munderover><mo>⁢</mo><mn>0</mn></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.2</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths>

p-0031The receiver statistic r_kin (2.1.2) is the integral

p-0032<maths id="MATH-US-00003" num="00003"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><msubsup><mo>∫</mo><mn>0</mn><msub><mi>T</mi><mi>s</mi></msub></msubsup><mo>⁢</mo><mrow><mrow><mi>r</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>⁢</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mrow><mi>t</mi><mo>-</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>⁢</mo><mstyle><mspace width="0.2em" height="0.2ex" /></mstyle><mo>⁢</mo><mrow><mo>ⅆ</mo><mi>t</mi></mrow></mrow></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><msubsup><mo>∫</mo><msub><mi>kT</mi><mi>p</mi></msub><mrow><msub><mi>T</mi><mi>s</mi></msub><mo>-</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow></msubsup><mo>⁢</mo><mrow><mrow><mi>r</mi><mo>⁡</mo><mrow><mo>(</mo><mrow><mi>τ</mi><mo>+</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>⁢</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow></mrow><mo>⁢</mo><mrow><mrow><mo>ⅆ</mo><mi>τ</mi></mrow><mo>.</mo></mrow></mrow></mrow></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.3</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> But since g(t) has support only on the interval [0,T_p]

p-0033<maths id="MATH-US-00004" num="00004"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><msubsup><mo>∫</mo><mn>0</mn><msub><mi>T</mi><mi>p</mi></msub></msubsup><mo>⁢</mo><mrow><mrow><mi>r</mi><mo>⁡</mo><mrow><mo>(</mo><mrow><mi>τ</mi><mo>-</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>⁢</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>τ</mi><mo>)</mo></mrow></mrow><mo>⁢</mo><mstyle><mspace width="0.2em" height="0.2ex" /></mstyle><mo>⁢</mo><mrow><mo>ⅆ</mo><mi>τ</mi></mrow></mrow></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.4</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Where we have defined a new signal

p-0034<maths id="MATH-US-00005" num="00005"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>-</mo><mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><mrow><mi>r</mi><mo>⁡</mo><mrow><mo>(</mo><mrow><mi>t</mi><mo>+</mo><msub><mi>kT</mi><mi>p</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>t</mi><mo>∈</mo><mrow><mo>[</mo><mrow><mn>0</mn><mo>,</mo><msub><mi>T</mi><mi>p</mi></msub></mrow><mo>]</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mi>else</mi></mtd></mtr></mtable><mo>.</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.5</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where r(t) is the received signal.

p-0035<figref idrefs="DRAWINGS">FIG. 3</figref> is a flow chart detailing the interaction between the preamble detector and the receiver. The preamble detection algorithm according to the embodiment illustrated in <figref idrefs="DRAWINGS">FIG. 3</figref> lends itself to noise covariance matrix estimation, which can, in turn, be used to increase the fidelity of the receiver. One embodiment of a system for performing this computation is demonstrated in <figref idrefs="DRAWINGS">FIG. 5</figref>. At <b>252</b>, digital processor <b>200</b> checks the detect flag (set in the detect preamble routine to reflect the detection of a preamble sequence). If this is false, the preamble detection algorithm is executed at <b>254</b>. If true, a half cycle counter hcl is checked at <b>256</b>. If not equal to 4, a count is added. If equal to 4, the signal is processed (received) at <b>258</b> by the digital processor <b>200</b> and the hcl count is reset to zero at <b>260</b>. At the end of the message, the detect counter is set to zero at <b>262</b>, thus ending the receiving process and enabling preamble detection once more.

p-0036Computing r_kis explained above; r_kis the inner product of the received signal with the signal waveform g(t−kT_p). For frequency-selective channels, this must be the TWACS signal as seen by the receiver. That is, it is the transmitted pulse operated on by the channel. The system in <figref idrefs="DRAWINGS">FIG. 5</figref> estimates the received pulse, as well as computing the necessary inner product.

p-0037We will call r_k(t) the kth half-cycle snapshot since it is the signal received during one half-cycle of the mains. If we let r_k+(t) be the analytic function of r_k(t) (2.1.4) becomes

p-0038<maths id="MATH-US-00006" num="00006"><math overflow="scroll"><mtable><mtr><mtd><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><mo>〈</mo><mrow><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>+</mo><msubsup><mi>r</mi><mi>k</mi><mo>*</mo></msubsup><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi /><mo>⁢</mo><mrow><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow></mrow><mo>+</mo><mrow><mo>〈</mo><mrow><mrow><msubsup><mi>r</mi><mi>k</mi><mo>*</mo></msubsup><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi /><mo>⁢</mo><mrow><mrow><mfrac><mn>1</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow></mrow><mo>+</mo><msup><mrow><mo>〈</mo><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>*</mo></msup></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi /><mo>⁢</mo><mrow><mfrac><mn>2</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mi>Re</mi><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mo>{</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>g</mi><mo>⁡</mo><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>}</mo></mrow></mrow></mrow></mtd></mtr></mtable></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.6</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Where we have used the fact that g(t) is a real signal. By Parseval's theorem (2.1.6) is

p-0039<maths id="MATH-US-00007" num="00007"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>=</mo><mrow><mfrac><mn>2</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mi>Re</mi><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mo>{</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>R</mi><mi>k</mi></msub><mo>+</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>G</mi><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>}</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.7</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Where R_k+(f) and G(f) are the Fourier transforms of r_k+(t) and g(t) respectively. The revised expression for the AWGN receiver is found by substitution of (2.1.7) for k=0 and k=2 into (2.1.2):

p-0040<maths id="MATH-US-00008" num="00008"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mrow><mrow><mrow><mo>-</mo><mfrac><mn>2</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac></mrow><mo>⁢</mo><mi>Re</mi><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mo>{</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>R</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>G</mi><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>}</mo></mrow></mrow><mo>+</mo><mrow><mfrac><mn>2</mn><msqrt><msub><mi>E</mi><mi>g</mi></msub></msqrt></mfrac><mo>⁢</mo><mi>Re</mi><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mo>{</mo><mrow><mo>〈</mo><mrow><mrow><msub><mi>R</mi><mrow><mn>2</mn><mo>+</mo></mrow></msub><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>,</mo><mrow><mi>G</mi><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>}</mo></mrow></mrow></mrow><mo>⁢</mo><munder><mover><mo>≶</mo><mrow><mover><mrow><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mi>m</mi></mrow><mo>^</mo></mover><mo>=</mo><mn>2</mn></mrow></mover><mrow><mover><mrow><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mi>m</mi></mrow><mo>^</mo></mover><mo>=</mo><mn>1</mn></mrow></munder><mo>⁢</mo><mn>0</mn></mrow><mo>⁢</mo><mstyle><mtext /></mstyle><mo>⁢</mo><mstyle><mspace width="4.4em" height="4.4ex" /></mstyle><mo>⁢</mo><mrow><mrow><mi>Re</mi><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mo>{</mo><mrow><mo>〈</mo><mrow><mrow><mrow><mo>-</mo><msub><mi>R</mi><mn>0</mn></msub></mrow><mo>+</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>+</mo><mrow><msub><mi>R</mi><mrow><mn>2</mn><mo>+</mo></mrow></msub><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mrow><mo>,</mo><mrow><mi>G</mi><mo>⁡</mo><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mrow><mo>〉</mo></mrow><mo>}</mo></mrow></mrow><mo>⁢</mo><munder><mover><mo>≶</mo><mrow><mover><mrow><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mi>m</mi></mrow><mo>^</mo></mover><mo>=</mo><mn>2</mn></mrow></mover><mrow><mover><mrow><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mi>m</mi></mrow><mo>^</mo></mover><mo>=</mo><mn>1</mn></mrow></munder><mo>⁢</mo><mn>0</mn></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.8</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Where the summation over k has been strategically placed inside the inner product. (2.1.8) is what must be computed in the receiver. We shall shift our focus now to how that is done.

p-0041Referring again to the block diagram in <figref idrefs="DRAWINGS">FIG. 5</figref>, the received signal, r(t), is operated on by a filter component <b>202</b> including a mains notch filter <b>204</b> to cancel out the 60 Hz (or 50 Hz) mains signal to provide an analog message signal (AMS). This is done to reduce the signal dynamic range and reduce quantization noise in the subsequent A/D conversion. The AMS is digitally sampled by A/D converter <b>206</b>B to provide a digital message signal (DMS) stored in a FIFO buffer <b>206</b> large enough to accommodate at least one mains half-cycle worth of data. The received signal r(t) is also low-pass filtered (LPF) prior to digitizing by A/D converter <b>206</b>A so that the mains signal is tracked without aliasing. The filtered analog outbound signal (FAOS) is converted by the A/D converter <b>206</b>A to a digitized filtered signal (DFS). The DFS is provided to a digital phase-locked loop (PLL) <b>210</b> which tracks the mains signal and provides half-cycle synchronization. The converters <b>206</b> may be part of the digital processor <b>200</b> or separate components. When a half-cycle boundary occurs, as illustrated by dotted lines in <figref idrefs="DRAWINGS">FIG. 4</figref>, the contents of the FIFO buffer are passed to the next stage of the algorithm as a column vector:

p-0042<maths id="MATH-US-00009" num="00009"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>r</mi><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>=</mo><mrow><mo>[</mo><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><msub><mi>T</mi><mi>Δ</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mo>*</mo></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mi>k</mi></msub><mo>⁡</mo><mrow><mo>(</mo><mrow><mrow><mo>(</mo><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow><mo>⁢</mo><msub><mi>T</mi><mi>Δ</mi></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>]</mo></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.9</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where T_Δ is the sample time and N is the half-cycle snapshot length.

p-0043Equation (2.1.8) requires half-cycle snapshots to be in the frequency domain. Thus, this column vector is operated on by an FFT <b>212</b>. Since only the Fourier transform of the analytic function, R_k+(f) is required, half of the output of the FFT may be discarded. In practice, many more than that may be discarded, provided the signal is sufficiently narrowband.

p-0044The FFT is used with a dimension reducing component DRC <b>214</b> for reducing the dimensionality of the vector prior to computing its inner product. The DRC eliminates unwanted signal components. Most of the elements of the vector output of the FFT can be discarded since the transmitted signal g(t) is bandlimited. Those elements not discarded should collectively represent that portion of the spectrum in which the SNR is significant. The resulting complex column vector, R₊[k], is a discretization of R_k+(f). With this discretization algorithm in place the receiver equation (2.1.8) is approximated by the matrix equation

p-0045<maths id="MATH-US-00010" num="00010"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mi>Re</mi><mo>⁢</mo><mrow><mo>{</mo><mrow><msup><mrow><mo>(</mo><mrow><mrow><mo>-</mo><mrow><msub><mi>R</mi><mo>+</mo></msub><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></mrow><mo>+</mo><mrow><msub><mi>R</mi><mo>+</mo></msub><mo>⁡</mo><mrow><mo>[</mo><mrow><mi>k</mi><mo>-</mo><mn>2</mn></mrow><mo>]</mo></mrow></mrow></mrow><mo>)</mo></mrow><mi>H</mi></msup><mo>⁢</mo><mi>G</mi></mrow><mo>}</mo></mrow></mrow><mo>⁢</mo><munder><mover><mo>≶</mo><mrow><mover><mi>m</mi><mo>.</mo></mover><mo>=</mo><mn>2</mn></mrow></mover><mrow><mover><mi>m</mi><mo>.</mo></mover><mo>=</mo><mn>1</mn></mrow></munder><mo>⁢</mo><mn>0.</mn></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.1</mn><mo>⁢</mo><mi>.10</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths>

p-0046G is a vectorization of the signal G(f) in the same way that R₊[k] is a vectorization of the signal R_+k(f). An estimate of G is provided by the preamble detection algorithm (see <figref idrefs="DRAWINGS">FIG. 3</figref>). Note that in <figref idrefs="DRAWINGS">FIG. 5</figref> the preamble detector output is denoted w rather than Ĝ. As noted below, the weight vector returned by the preamble detector may implement adaptive equalization. Also, as shown in <figref idrefs="DRAWINGS">FIG. 5</figref> the sum −R+[k]+R+[k−2] is computed using the FIR filter [−1+z−2]. The output of this is fed to the inner product calculator. The preamble detector also provides symbol synchronization information as described by <figref idrefs="DRAWINGS">FIG. 3</figref>. The synchronized stream is compared to 0 to determine the transmitted symbol.

p-0047The receiver <b>258</b> implements the following equation using the weight vector provided by the preamble detector: <br />γ=<i>Re{w</i>^H<i>z[k]}</i>

p-0048The bit estimate is obtained by thresholding this scalar:

p-0049<maths id="MATH-US-00011" num="00011"><math overflow="scroll"><mrow><mrow><mover><mrow><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mi>m</mi></mrow><mo>^</mo></mover><mo>⁢</mo><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mrow><mi>if</mi><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>γ</mi></mrow><mo>></mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>else</mi></mtd></mtr></mtable></mrow></mrow></math></maths>

p-0050The length field encodes the length of the variable fields of the outbound packet. The total length in bits is then an affine function of the value encoded in the length field. Prior to reading this field the length is set to a default value which should represent a maximum length. After the number of symbols retrieved corresponds to length the receiver is reset to search for preambles once more.

h-0005Preamble Detector

p-0051In summary, the digital processor <b>200</b> executes a digital preprocessor component <b>220</b> receiving the corresponding digitized signal r[k] and preprocessing the received corresponding digitized signal to provide a preprocessed digital signal z[k]. The processor <b>200</b> also executes a digital symbol determination component <b>222</b> receiving the preprocessed digital signal z[k] and digitally demodulating the received preprocessed digital signal z[k] to provide the outbound message ({circumflex over (m)}). The symbol determination component <b>222</b> includes the preamble detector <b>254</b> which provides a weight vector w to a digital symbol estimation component <b>224</b> generating a noise product signal γ which is selectively sampled by the synchronization signal to provide the outbound message ({circumflex over (m)}).

p-0052A filtered preamble sequence (see detect preamble <b>254</b> in <figref idrefs="DRAWINGS">FIG. 3</figref>) is formed: <br /><i>{circumflex over (q)}</i>⁽⁰⁾<i>[k</i>]←(−<i>p</i>_k<i>+p</i>_k−2)(−1)^k  (2.1.11)<br /><figref idrefs="DRAWINGS">FIG. 6</figref> illustrates unfiltered preamble coefficients P_kas seen at the output of the FFT of <figref idrefs="DRAWINGS">FIG. 5</figref> and filtered preamble coefficients q[k] as seen at the output of FIR filter <b>216</b> (−1+z−²) of <figref idrefs="DRAWINGS">FIG. 5</figref>.

p-0053The received symbols Ĝ⁽ⁿ⁾and steering vector {circumflex over (q)}⁽ⁿ⁾are estimated iteratively, using the following algorithm:

p-0054<maths id="MATH-US-00012" num="00012"><math overflow="scroll"><mrow><mrow><mi>for</mi><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>n</mi></mrow><mo>=</mo><mrow><mn>1</mn><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>to</mi><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>N</mi><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>do</mi></mrow></mrow></math></maths><maths id="MATH-US-00012-2" num="00012.2"><math overflow="scroll"><mrow><mstyle><mspace width="1.4em" height="1.4ex" /></mstyle><mo>⁢</mo><mrow><msup><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>←</mo><mfrac><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>⁢</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo>⁢</mo><mrow><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>⁢</mo><mrow><mi>z</mi><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></mrow></mrow><mrow><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo>⁢</mo><msup><mrow><mo></mo><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo></mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></mrow></math></maths><maths id="MATH-US-00012-3" num="00012.3"><math overflow="scroll"><mrow><mstyle><mspace width="1.4em" height="1.4ex" /></mstyle><mo>⁢</mo><mrow><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>←</mo><mfrac><mrow><msup><mrow><mo>(</mo><msup><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mi>H</mi></msup><mo>⁢</mo><mrow><mi>z</mi><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></mrow><mrow><msup><mrow><mo>(</mo><msup><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow><mi>H</mi></msup><mo>⁢</mo><msup><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup></mrow></mfrac></mrow></mrow></math></maths><maths id="MATH-US-00012-4" num="00012.4"><math overflow="scroll"><mrow><mi>end</mi><mo>⁢</mo><mstyle><mspace width="0.8em" height="0.8ex" /></mstyle><mo>⁢</mo><mi>for</mi></mrow></math></maths>

p-0055Next, the symbols are thresholded using:

p-0056<maths id="MATH-US-00013" num="00013"><math overflow="scroll"><mrow><mrow><msup><mi>q</mi><mi>′</mi></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>=</mo><mrow><mo>{</mo><mtable><mtr><mtd><mrow><mrow><mo>-</mo><mn>1</mn></mrow><mo>,</mo><mrow><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo><</mo><mrow><mo>-</mo><mn>0.5</mn></mrow></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><mrow><mrow><mo>-</mo><mn>0.5</mn></mrow><mo><</mo><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo><</mo><mn>0.5</mn></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mn>1</mn><mo>,</mo><mrow><mrow><msup><mover><mi>q</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>></mo><mn>0.5</mn></mrow></mrow></mtd></mtr></mtable></mrow></mrow></math></maths>

p-0057The resulting sequence is compared to the sequence in (2.1.11). If there are fewer than M errors, then detect is set to 1.

p-0058If optional adaptive equalization is not used, the weight vector is set to: <br /><i>w=Ĝ</i>^(N)

p-0059To implement adaptive equalization, the noise covariance matrix, S, is estimated using equation 2.3.1:

p-0060<maths id="MATH-US-00014" num="00014"><math overflow="scroll"><mtable><mtr><mtd><mrow><mover><mi>S</mi><mo>^</mo></mover><mo>=</mo><mrow><mfrac><mn>1</mn><mrow><mo></mo><msub><mi>K</mi><mn>0</mn></msub><mo></mo></mrow></mfrac><mo>⁢</mo><mrow><munder><mo>∑</mo><mrow><mi>k</mi><mo>∈</mo><msub><mi>K</mi><mn>0</mn></msub></mrow></munder><mo>⁢</mo><mrow><mrow><mi>z</mi><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow><mo>⁢</mo><mrow><msup><mi>z</mi><mi>H</mi></msup><mo>⁡</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></mrow></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.3</mn><mo>⁢</mo><mi>.1</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where k ε K_oif q[k]=0.

p-0061The weight vector returned is then: <br /><i>w=S</i>⁻¹<i>Ĝ</i>^(N)<br /> Phased Lock Loop (PLL)

p-0062<figref idrefs="DRAWINGS">FIG. 7</figref> illustrates a block diagram of one embodiment of the digital phase locked loop (PLL) <b>210</b> of <figref idrefs="DRAWINGS">FIG. 5</figref>. <figref idrefs="DRAWINGS">FIG. 8</figref> illustrates one embodiment of the frequency and impulse response of the loop filter H(z) of the digital PLL of <figref idrefs="DRAWINGS">FIG. 7</figref>.

p-0063A precision phase-locked loop (PLL) <b>210</b> locking the receiver to the mains signal is a part of the receiver illustrated in <figref idrefs="DRAWINGS">FIG. 5</figref>. Zero-crossing detection is one way of doing this, but is sensitive to noise. Doing so would require some kind of jitter mitigation algorithm. Since the ratio of the mains signal power to the noise power is so high, a digital version of an analog phase locked loop may be used. This algorithm is illustrated in the block diagram of <figref idrefs="DRAWINGS">FIG. 7</figref>. The loop consists of a phase detector <b>302</b>, a loop filter <b>304</b>, a digitally controlled oscillator <b>306</b>, and an anti-aliasing filter <b>308</b>.

p-0064The phase-locked loop <b>210</b> is designed for digital signals in which the lock-on component is at or very near a digital frequency of π/2. This requires a sampling rate four times the mains frequency, but in practice signals are sampled at much higher frequencies. <figref idrefs="DRAWINGS">FIG. 7</figref> illustrates several decimation blocks for downsampling to this lower frequency. If, for example, for inputs v[n] sampled at 12 KHz, operating on 60 Hz power systems, M equals 50. For sampling rates of 12.5 KHz, M equals 52. For 12 KHz sampling rates operating on 50 Hz power systems, M equals 60.

p-0065In general, the mains signal is not simply a sinusoid. G(z) is an anti-aliasing filter <b>308</b> whose purpose is to eliminate power system harmonics which would otherwise be aliased into the low-rate digital signal fed to the PLL. If the principal component after downsampling is at digital frequency π/4, then it is only necessary to attenuate harmonics 3 and above. A low pass filter (LPF) with a transition region between 70 and 170 Hz (on 60 Hz power systems) is one option. The signal amplitude must be between −1 and 1. The filter gain should be set accordingly. Optionally, an adaptive gain control stage that adapts to the mains signal strength may be inserted. The amplitude of the signal fed to the mixer is critical when considering lock-in and hold-in range.

p-0066The digitally controlled oscillator (DCO) <b>306</b> generates a sinusoid with a frequency that is an affine function of the input from the loop filter. However, the estimate of the mains phase, {circumflex over (φ)}[n], is provided at the same sampling frequency as the mains signal v[n]. If the input to the DCO is a constant x₀, then the phase estimate is

p-0067<maths id="MATH-US-00015" num="00015"><math overflow="scroll"><mtable><mtr><mtd><mrow><mrow><mover><mi>Φ</mi><mo>^</mo></mover><mo>⁡</mo><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mo>=</mo><mrow><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><msup><mi>z</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow></mfrac><mo>⁢</mo><mrow><mo>(</mo><mrow><msub><mi>x</mi><mn>0</mn></msub><mo>|</mo><mi>k</mi></mrow><mo>)</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.4</mn><mo>⁢</mo><mi>.1</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> Where k is the bias constant. In the time domain <br />{circumflex over (φ)}[<i>n]={circumflex over (φ)}[n−</i>1]+<i>x</i>₀<i>+k</i>  (2.4.2)<br />{circumflex over (φ)}[<i>n]=n</i>(<i>x</i>₀<i>+k</i>).  (2.4.3)<br /> So the DCO sends this signal to the phase detector: <br /><i>y[n</i>]=cos(2π<i>n</i>(<i>x</i>₀<i>+k</i>)).  (2.4.4)<br /> The bias constant k must be set so that when the input x0=0 the DCO sends a sinusoid with the mains frequency to the phase detector. Thus

p-0068<maths id="MATH-US-00016" num="00016"><math overflow="scroll"><mtable><mtr><mtd><mrow><mi>k</mi><mo>=</mo><mfrac><msub><mi>f</mi><mn>0</mn></msub><msub><mi>F</mi><mi>s</mi></msub></mfrac></mrow></mtd><mtd><mrow><mo>(</mo><mrow><mn>2.4</mn><mo>⁢</mo><mi>.5</mi></mrow><mo>)</mo></mrow></mtd></mtr></mtable></math></maths><br /> where f₀is the mains frequency and F_sis the sampling rate of v[n].

p-0069The phase detector <b>302</b> output will consist of the actual phase estimate with a frequency near 0 and nuisance components with frequencies near ±π. There may also be components with frequency near ±π/2 resulting from the mains harmonics. The loop filter exists to remove all but the phase estimate. Thus, zeros are placed at z=±j and z=−1 to mitigate harmonic output from the phase detector arising from the 0th and 2nd mains harmonics, which pass through the anti-aliasing filter. Poles are placed at z=0.5 to increase gain at low frequencies. However, care must be taken not to make the phase slope too steep. Analyses have suggested that hold-in range degrades as the filter delay increases.

h-0006The filter is an IIR filter of the form

p-0070<maths id="MATH-US-00017" num="00017"><math overflow="scroll"><mrow><mrow><mi>H</mi><mo>⁡</mo><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow></mrow><mo>=</mo><mfrac><mrow><munder><mo>∑</mo><mi>n</mi></munder><mo>⁢</mo><mrow><msub><mi>b</mi><mi>n</mi></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mi>n</mi></mrow></msup></mrow></mrow><mrow><munder><mo>∑</mo><mi>n</mi></munder><mo>⁢</mo><mrow><msub><mi>a</mi><mi>n</mi></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mi>n</mi></mrow></msup></mrow></mrow></mfrac></mrow></math></maths><br /> having the coefficients as shown in Table 1, below.

p-0071<tables id="TABLE-US-00001" num="00001"><table frame="none" colsep="0" rowsep="0"><tgroup align="left" colsep="0" rowsep="0" cols="1"><colspec colname="1" colwidth="217pt" align="center" /><thead><row><entry namest="1" nameend="1" rowsep="1">TABLE 1</entry></row></thead><tbody valign="top"><row><entry namest="1" nameend="1" align="center" rowsep="1" /></row><row><entry>Filter Coefficients</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="2"><colspec colname="offset" colwidth="28pt" align="left" /><colspec colname="1" colwidth="189pt" align="center" /><tbody valign="top"><row><entry /><entry>n</entry></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="6"><colspec colname="offset" colwidth="28pt" align="left" /><colspec colname="1" colwidth="42pt" align="center" /><colspec colname="2" colwidth="14pt" align="center" /><colspec colname="3" colwidth="63pt" align="center" /><colspec colname="4" colwidth="14pt" align="center" /><colspec colname="5" colwidth="56pt" align="center" /><tbody valign="top"><row><entry /><entry>0</entry><entry>1</entry><entry>2</entry><entry>3</entry><entry>4</entry></row><row><entry /><entry namest="offset" nameend="5" align="center" rowsep="1" /></row></tbody></tgroup><tgroup align="left" colsep="0" rowsep="0" cols="7"><colspec colname="offset" colwidth="14pt" align="left" /><colspec colname="1" colwidth="14pt" align="left" /><colspec colname="2" colwidth="42pt" align="center" /><colspec colname="3" colwidth="14pt" align="char" char="." /><colspec colname="4" colwidth="63pt" align="char" char="." /><colspec colname="5" colwidth="14pt" align="center" /><colspec colname="6" colwidth="56pt" align="center" /><tbody valign="top"><row><entry /><entry>a_n</entry><entry>1</entry><entry>−1</entry><entry>0.25</entry><entry>0</entry><entry>0</entry></row><row><entry /><entry>b_n</entry><entry>1</entry><entry>2</entry><entry>2</entry><entry>2</entry><entry>1</entry></row><row><entry /><entry namest="offset" nameend="6" align="center" rowsep="1" /></row></tbody></tgroup></table></tables><br /> Notice that all the filter coefficients are powers of 2. This has the potential of simplifying implementation since each multiplication is a simple bit shift. The resulting filter characteristics are illustrated in <figref idrefs="DRAWINGS">FIG. 8</figref>.

p-0072The loop gain, K, affects the lock-in range and the hold-in range. The choice of filter gain affects the hold-in range of the PLL and the phase accuracy. In general, increasing K increases the hold-in range but also increases phase estimate accuracy. Increasing K beyond about 2⁻¹³will cause the loop to abruptly stop functioning altogether.

p-0073For purposes of illustration, programs and other executable program components, such as the operating system, are illustrated herein as discrete blocks. It is recognized, however, that such programs and components reside at various times in different storage components of the computer, and are executed by the data processor(s) of the computer.

p-0074Although described in connection with an exemplary computing system environment, embodiments of the invention are operational with numerous other general purpose or special purpose computing system environments or configurations. The computing system environment is not intended to suggest any limitation as to the scope of use or functionality of any aspect of the invention. Moreover, the computing system environment should not be interpreted as having any dependency or requirement relating to any one or combination of components illustrated in the exemplary operating environment. Examples of well known computing systems, environments, and/or configurations that may be suitable for use with aspects of the invention include, but are not limited to, personal computers, server computers, hand-held or laptop devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable consumer electronics, mobile telephones, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.

p-0075Embodiments of the invention may be described in the general context of components, data and/or computer-executable instructions, such as program modules, stored one or more tangible non-transitory computer storage media and executed by one or more computers or other devices. Generally, program modules include, but are not limited to, routines, programs, objects, components, and data structures that perform particular tasks or implement particular abstract data types. Aspects of the invention may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.

p-0076In operation, computers and/or servers may execute the computer-executable instructions such as those illustrated herein to implement aspects of the invention.

p-0077The order of execution or performance of the operations in embodiments of the invention illustrated and described herein is not essential, unless otherwise specified. That is, the operations may be performed in any order, unless otherwise specified, and embodiments of the invention may include additional or fewer operations than those disclosed herein. For example, it is contemplated that executing or performing a particular operation before, contemporaneously with, or after another operation is within the scope of aspects of the invention.

p-0078Embodiments of the invention may be implemented with computer-executable instructions. The computer-executable instructions may be organized into one or more computer-executable components or modules on a tangible computer readable storage medium. Aspects of the invention may be implemented with any number and organization of such components or modules. For example, aspects of the invention are not limited to the specific computer-executable instructions or the specific components or modules illustrated in the figures and described herein. Other embodiments of the invention may include different computer-executable instructions or components having more or less functionality than illustrated and described herein.

p-0079When introducing elements of aspects of the invention or the embodiments thereof, the articles “a,” “an,” “the,” and “said” are intended to mean that there are one or more of the elements. The terms “comprising,” “including,” and “having” are intended to be inclusive and mean that there may be additional elements other than the listed elements.

p-0080In view of the above, it will be seen that several advantages of the invention are achieved and other advantageous results attained.

p-0081Not all of the depicted components illustrated or described may be required. In addition, some implementations and embodiments may include additional components. Variations in the arrangement and type of the components may be made without departing from the spirit or scope of the claims as set forth herein. Additional, different or fewer components may be provided and components may be combined. Alternatively or in addition, a component may be implemented by several components.

p-0082Having described aspects of the invention in detail, it will be apparent that modifications and variations are possible without departing from the scope of aspects of the invention as defined in the appended claims. As various changes could be made in the above constructions, products, and methods without departing from the scope of aspects of the invention, it is intended that all matter contained in the above description and shown in the accompanying drawings shall be interpreted as illustrative and not in a limiting sense.