SimulationofaPhaseLockedLoo SimulationofaPhaseLockedLoopUsingLabVI pUsingLabVIEW EW.
T.J.Moir MasseyUniversityatAlbany InstituteofInformationandMathem InstituteofInformationandMathematicalSci aticalSciences ences Albany Auckland NewZealand Email:
[email protected]
2 Abstract
A meth ethod for sim simulat ulatiing an anal analog ogue ue phas phasee-llocke ocked d loop (PLL (PLL) ) is show shownn. The The programminglanguagegivenisLabVIEW(trademarkNationalInstruments)butthe techn techniiques ques used used are qui quite gen general eral and and appl applyy equa equallly wel well to other other progr program amm ming languag languages esand and packages.LabVIEW packages.LabVIEW seems seemsanunl anunlik ikely ely candid candidate ateto to simu simulatea lateaPLL PLL as it is mo more re often often associ associated ated with with process process control control system systems s and and virtual virtual instrum instrumen entati tation on ratherthancommunicationsystems.HoweveritwillbeshownthatinfactLabVIEW givesapowerfulsolutionwhichisbotheasytoimplementandwithausergraphical interface interface comparabl comparable e to any mod modern ern programmi programming ng languag language. e. The simul simulation ation is full fully y interact interactivean iveanddemodul ddemodulatesordinaryFM atesordinaryFMli likearadioreceiver. kearadioreceiver. Keywords:
Phase-LockedLoop(PLL),Feedbacksystem Phase-LockedLoop(PLL),Feedb acksystem,LabVIEW ,LabVIEW
1.Introduction
ThePhase-LockedLoop(PLL)isoneofthemostcommonlyusedintegratedcircuits (ICs) (ICs) in use in mo moder dern n comm communic unicati ations ons system systems[1 s[1]. ]. Althou Although gh perha perhaps ps surpri surprisi sing ngly ly firstinventedasearlyas1932byBellescise[2]itnevergainpopularityuntiltheearly 1970 1970s s when when chea cheap p ICs ICs were were read readiily avai availlabl able. It qui quickl ckly foun ound appl appliicati cation on as a precisi precision on FMdem FM demodul odulator ator as a replacem replacement ent for Foster-Seely Foster-Seely discri discrimi minators. nators.Dig Digital ital commu communi nication cation system systems s quickl quickly yfol follow lowed ed and the PLL has found found applica application tionin in such area reas as mode odems, mobi obile com omm munication tionss, satel tellite receivers and tel television electronics.ThePLLisusedextensivelyinmodernelectronicsystemsbutitsdesignis often met met with with trepidation trepidation . This This is perhaps perhapsund understan erstandab dable lesi since nceto to fully fully understand understand the the operat operatiion of a PLL PLL requi requires res some some knowl knowled edge ge of comm communi unicati cation on syste system ms and and controlsystems,twosubjectswhicharetreatedinisolation.Forexampleseldomare PLLs PLLs covered covered in a taught taughtcon control trol course courseat at under undergra gradua duate te leve levelwhi lwhils lstt feedb feedback ack and and stabilityisonlybrieflymentionedwhencoveringPLLsinacommunicationcourse. The The purp purpos ose e of thi this pape paper r is to show show how a real realiisti stic PLL PLL can can be sim simulate ulated. d. The The sim simulat ulatiion uses uses Lab LabVIEW VIEW as it is easy easy to buil uild a grap graphhical cal interf terfac ace e whi which is interactiv interactive. e. Although Although such packages packages as MATLA MATLAB B and MATRIXx MATRIXx [3] could could equally equally wel well be used used the the resul result t wil will not not be as intera nteracti ctivve. For For exam exampl ple e in [3] [3] alth althou ough gh the the simu simulati lation onis is realis realisticit ticitisoft isofthe‘one-sh he‘one-shot’ ot’vari varietywherethe etywheretheprogram program mu must stbere-run bere-run toshowtheeffectofanychangesandcannotbemadetoruncontinuouslyinpseudo ‘real-tim ‘real-time’ e’li like kethe theapproac approachusedhere. husedhere.Thefin Thefinalgoalisa algoalisa single single PLL program which which canthenbeappliedtoawholerangeofpossibleapplicationssuchasdemodulationof FM, FM, Frequ Frequen ency cy Shif Shiftt Key Keying ing (FSK) (FSK) or to furthe further r investi nvestigate gate such such proble problem ms such such as multi-pathandadditivenoise(thresholding).
3
2.TheBasicPLL
TheblockdiagramofagenericPLLisshowninFigure1below.
Figure1.BlockDiagramofGenericPLL
We consi consideran deran inputto inputto the PLL PLLwh whic ichis hisFM FM..Wh Whil ilstt stthe hereare reare other othertyp typesof esofin input put whichcanbeused,thisapproachgivestheopportunityofmodulatingvariousdifferent basebandsignalswhichcanbelaterusedtotestthePLL.Forexampleastepinputora sinu sinusoi soidal dal frequ frequen ency cy respons response e are comm commonly only used used for testing testing all all feedb feedback ack control control syste system ms. The The FM sign signal al whi which must ust be sim simulate ulated d is an FM modul odulated ated sin sine wave wave (al (althou though gh squa square re wav waves and and other other wav wavefor eform ms are are also also con conside sidere red) d) and and has the the analogueformf(t)where f ( t ) = cos[
c
( t ) +
sin(
m
t )]
Intheabove = ∆ / m isdefinedtobetheusualFMmodulationindexwith c and res respecti ectivvely ely the the carr carriier and base aseband freq requenc encies in rad rad/s. The The depth of m modulationisgivenby ∆ rad/s. InFigure1thePLLcomprisesaphasedetector(PD),avoltage-controlledoscillator (VCO)andafilter.ThePLLwilldemodulatetheFMandgiveandoutputwhichisthe ori origin ginal baseb aseban and d sign signal al.. The The PLL PLL woul would d normal ormallly oper operat ate e on the the interm termed ediiate ate frequency(IF)wavef frequency(IF)waveformofaradioreceiv ormofaradioreceiver. er. The The theory theory of the PLL is well well docum documented ented [1,4] [1,4] and and only only an outlin outline e will will be give given n here. ere. The The operat operatiion is easi easies est t to see see when when there there is no FM and and onl only a fixed xed carri carrier er frequ frequen ency cy is presen presentt (ie (ie ∆ω = 0 ). When in lock ock the the outp output ut of the the VCO VCO wil will be a
4
wave wavefform whi which is in phas phase-q e-qua uadra dratur ture e with with the the incom ncomiing wave wavefform. orm. The The phas phase e detector detector (PD) (PD) is a line linear ar multip ultipli lier er (for (for our exam exampl ple) e) and and give gives s an output output whic which h is approx approxim imatel atelyy zero with with an additi additive ve term at twice twicethe thecarr carrie ier r frequ frequen ency cy 2ω c . When When properly properly designe designed d (iethe (ie the bandwi bandwidth dthis ischose chosen napprop appropriatel riately) y)the the termat term at 2ω c willbe attenuatedbythefilter(apartfromsomeresidualleft-over)andthePLLoutputwillsit at zero. zero.Wh When en aba a baseb seban and dsi sign gnal al is presen presented(ie ted(ieFM FM) )the the VCO output outputwi will ll track the variationinphaseoftheincomingFMandthePLLoutput(theinputtotheVCO)will be the the rate rate of chan change ge of phas phase e (ie (ie instan stanta tanneous eous frequ requen ency cy). ). The The VCO VCO thus thus acts acts dyna dynam mical cally as an integ ntegrat rator or and and thi this is importan portant t when when exam examiining ning the the openopen-lloop frequencyresponse. The Thebe beauty auty ofthe PLL PLLiistha s thatt it can be anal analyysed in a line linear arfo form rm indep indepen enden dentt ofthe carrierfrequency.TheblockdiagramoftheclosedloopPLLisshowninFigure2.The VCOtransferfunctionH(s)isshownasapureintegratorandthephasedetectorasa summi summing ngju juncti nction on providi providing ngnega negativ tive e feedba feedback. ck. It remain remains s to find findthe thefi filter lterdy dynam namics ics F(s). in
+
Filter
PD
F(s)
DemodulatedSignal
o
K v s
VCO H(s) Figure2.LinearPLL
Thefilterdynamicsarefoundbydrawingtheopen-loopBodeplotwhichshouldhave asimil asimilarformtotheoneshownin arformtotheoneshowninFigure3.(although Figure3.(althoughotherformsarepossibl otherformsarepossible) e)
5 dBGain -40dB/decade
f 1 0
-20dB/decade
f
Frequency(Hz)
f 2 -40dB/decade
Phase(degrees)
PhaseMargin Frequency(Hz) o
-180
Figure3.OpenLoopPLLBodePlot
This This typeof typeofPL PLLissome Lissometim timesknow esknownas nasa athi thirdorder rdorder typeII typeII PLLastherearetwo PLLastherearetwo integratorswithintheloop.ThefirstintegratoristheVCOandthesecondisanadded elec electro tronnic integra tegrator tor.. Sin Since two two integr tegrat ator ors s with with negati egativve feedb eedbac ack k resu resullts in an oscill oscillator, ator, apha a phase selea lead d(adv (advanc ance) e)stabi stabili lisation sation is needed. needed.Hen Hencet cetheoveral heoverallBode lBodepl plot ot has has the form form shown shown in Figur Figure e 3. This This particul particular ar desig design nis is prefe preferred rred asit as itha has s better better trackingabilitiesthanatypeIPLL.Thehigherthegainatlowfrequenciesresultsin goodtrackingandhencelowerror.(forthecontrolsystemdetailsseereference[5]) The The uni unity gain ain bandw andwiidth of the the PLL PLL shoul ould be chos chosen en high enou enough gh to trac trackk adequa adequatel telyy but but nottoo high highso so astole asto letttoomuc toomuch h 2 c throug through. h. Byexper Byexperiienceithas enceithas beenfoundthataunitygainbandwidthof
givesgoodresults.
f
=
2 f c 10
6
3.LabVIEWSimulation
Thegraph The graphical ical userin user interface terface of LabVIEW LabVIEW together with with thebl the block ockdi diagram agram approach approach makesitagoodcandidateforsimulatingaPLL.Thesimulationisdividedintoseveral steps.Thegen steps.The generation eration ofFM, simu simulati latinga nga first firstorder orderli lineartime-i neartime-inv nvaria ariant ntsy system, stem, the VCO and and phase ase detec etecto tor r and and the the com composi osite des design. gn. The The spec speciificati cation on for the the simulationisasfollows: sampl samplin ing g frequ frequen ency cy 10kHz 10kHz,car ,carri rier er frequ frequen ency cy 2kHz, 2kHz, mo modul dulati ation on frequ frequen ency cy as a +/- percentageofthecarrierfrequencyuptosay+/-10%(+/-200Hz),unity-gaincrossover frequency400Hz,phase-m frequency400Hz,phase-marginaround55 arginaround55degrees. degrees. 3.1SimulationofFrequencyModulationandtheVCO
TheFMsignalgenerationiscoveredinoneofthelabVIEWexamplesandsoisoneof theeasiestofthestepstofoll theeasiestofthestepsto follow.Theblockdi ow.Theblockdiagrami agramisshowni sshowninFi nFigure4. gure4.
Figure4FMSimulationinLabVIEW
In common common with with thewhol the whole e PLL simul simulation ation theFMgenerator the FMgenerator (which (whichis is only onlya a slig slight ht modifi odifica cati tion on of a Lab LabVIEW VIEW exam exampl ple e progr program am) ) uses uses scal scalar ar quan quanti titi ties es rathe rather r than than arrays.Figure4showsasinewaveasthemodulatingsignalandasinewaveasthe carrier carrier signal signal.. Thecase statement statementenab enables les diff different erentki kinds nds ofbase of baseban band dsig signal nals snam namely ely,, square square,, trian triangl gle e and and sawtooth. sawtooth. The The vari various ous wavef waveform orm generato generators rs defaul defaultt to a vector vector outputandsotheymustbeconvertedtoscalarformbyindexingthearrayandtaking the the zero zeroth th value alue (the (the first rst sam sample) ple).. The The FM gen generat eratiion is sim similar to how FM is generatedusingtwosinusoidalgenerators.Theoutputofthebasebandsignalgenerator isfedintothe‘sweep’inputofthesecondwhichactsasthecarrierfrequency.Thesine generatorforthecarrieronlyacceptsscalarinputsforfrequencyandhencetheneed for for a scal scalar sim simulati ulation on.. Normal Normalis ised ed frequ frequen ency cy is used used through throughout out defin defined ed as actual actual frequenc frequency y(Hz)/Sa (Hz)/Samp mpli ling ngFrequ Frequency ency (Hz). Thesam The sampli plingfrequen ngfrequency cy of10kH of 10kHz z usedin thePLLsimulationwaschosenasitismorethantentimestheunitygainbandwidthof 400Hzandabouteighttimestheupperbreak-frequencyoftheopen-loopBodeplot. Thecontinuouslychangingamplitudeofthefirstsinewavegeneratorismultipliedby the the fixed xed carri carrier er freque requenc ncyy of 2kHz 2kHz/1 /10k 0kHz Hz and and thi this chan change ges s the the freque requenncy of the the secondoscillatorandproducesFM.Adcoffsetatthefirstoscillatoroutputisrequired sothatwhenthereisnobasebandsignal,unityismultipliedbythecarrierfrequencyto give give a contin continuou uous s wavef waveform orm.. The The ampl amplitu itude de of the first first oscil oscilla lator tor (baseb (baseban and d sign signal al) )
7
determin determines es the depth of mod modul ulation ation.. Ifthe If the baseban baseband d signal signal had ampl amplitud itude e unity unitythen then byoff-settingitsoutputbyunityandmultiplyingby2kHz/10kHzthesecondoscillator willsweepfrom0hzto4kHzwhichrepresents100%FMmodulation.Theamplitude ofthebasebandosci ofthebasebandoscill llator atorX100 X100 thereforerepresentspercen thereforerepresentspercentagedepthofmodul tagedepthofmodulation ation.. It is nom nomin inal ally ly setto set to 10% through throughout out whic which hrep represe resents nts a carrie carrier r frequ frequen ency cy whic which his is centredon2kHzandsweeps+/-200Hz. ItisworthmentioningtheVCOoperationinthissectionasitisnearlyidenticaltothe aboveFMgeneration.Figure5showstheLAbVIEWblockdiagramoftheVCOand thephasedetector.
Figure5VCOandPhaseDetector
The The VCO sim simulati ulation on is the FM generator generator with with noba no baseb seban and d input. input. As with with the FM generato generator r the VCO VCOin inpu putt has hasan an offs offset et ofunity ofunity sothatwh sothat when en the VCO inputiszero, inputiszero, unity multip tiplies the the VCO freeree-rrunning freq requency (2k (2kHz/1 z/10kHz) and gives a continuouswaveform.ThescalingoftheVCOis2kHz/voltandhencetheVCOgain is K v = 12566 rad / s / volt .TheinputunitistakenasvoltstoconformtoarealVCO. Thi This gai gain is inheren herent t in the the VCO VCO and and must ust be accou account nted ed for when when cal calcul culatin ating the the overallgainoftheloop.ThephasedetectorisalinearmultiplierasshowninFigure5. ItstwoinputscomefromtheFMgeneratorandtheVCOoutput.TheVCOoutputisa sinewaveherebutcanbeeasilychangedtoasquarewaveasisthecaseinmostICs. Thephasedetectoroutputfeedsin Thephasedetectoro utputfeedsintot tothefi hefilterwhi lterwhichisdi chisdiscussednex scussednext. t. 3.2FilterSimulation
ExcludingthedynamicsoftheVCOwhichisoftheformofanintegrator,itremains tosimulatethefiltertransferfunctionwhichisoftheform
F ( s ) =
(1 + sT 1 ) s (1 + sT 2 ) K
forT1>T2.
8
Theaboveequationisa secondlineartim secondlineartime-in e-invari variant(LTI) ant(LTI) transferfuncti transferfunctionwhich onwhichcan can be spli split into into two first-ord first-order ertran transf sfer erfu func nction tions sin in cascad cascade. e.Lab LabVIEW VIEW in its basic basic form hasmany.viblockswhichareavailableforfilteringetcbutnonearedirectlyrelevant to implem plemen enti ting ng the above above equat equatiion. It was theref therefore ore deci decided ded to bui build a .vi .vi block block whichwouldimplementanyfirstorderLTIsystemandcascadethemtoconstructF(s). Itisamatterofchoiceastowhethertoimplementonesecond-orderLTIortwofirstorder order LTI syste system ms but but it was deci decide ded d to go for the sim simples plestt soluti solution. on. Consi Conside der r the generalfirst-ordersystemC(s)where g (b0 + b1s ) C ( s) = (s + a) OnepossiblesignalflowgraphforC(s)isshowninFigure6.
Figure6.Signalflowgraphofgenericfir Figure6.Signalflowgraphofgenericfirst-ordersystem st-ordersystem
Where here y is the the syst system em outpu output, t, u is the the syst system em input put and and x is a state state vari variab ablle. For For LabV LabVIEW IEW to implem plemen entt this this an integra ntegrati tion on algori algorithm thm is requi required red.. There There are many any techniquesforintegrationbutperhapsthemostpopularmethodistouseTrapezoidal integr ntegrati ation. on. A Trapez Trapezoi oidal dal integra ntegrator tor is repres represen ented ted by the Bil Bilinear near transf transform orm.. In differenceequationformitbecomes
y k
= y k −1 +
T
2
[uk + uk 1 ] −
whereTisthesamplinginterval(0.1ms)andhastheLabVIEWblockdiagramshown inFigure7.
9
Figure7TrapezoidalIntegrationusingLabVIEW
Inthedia Inthe diagram gram delta deltatisthe tisthestepsizeT,and stepsizeT,andgai gainis0.5set nis0.5setextern externall ally.Whe y.Whencombi ncombined ned with with the flow flow graph graph a gene general ral first-ord rst-order er syste system m is constr construc ucted ted and and has has the form form of Figure8.
Figure8.Generalfirst-orderLTIsystem.
Intheabovediagramtheblockwiththeintegrationsymbolrepresentsthepreviously disc discus usse sed d Trap Trapez ezoi oida dall integra tegrator tor.. The The While loops oops in Figu Figure re 7 and and 8 hav have thei their r condi conditi tion onse sett to Fal False so that that they they only only iterate terate one tim time for for every every loopof oop of the main ain
10
extern external al Whi While loop. loop. When When imp mple lem menting enting feedb feedback ack in LabVIEW LabVIEW a sign signal al cannot cannot be connected connected directl directlyy back backas asin in abl a block ockdi diagra agram m.Ins . Insteadit teaditm mustbe stored in aregi a register ster andthepreviousvaluefedbackasshownwiththestatevariablexinFigure8.Thisis becauseadigital becauseadigitalsystem systemcannotrespondinstan cannotrespondinstantaneousl taneouslyastheremustbeat yastheremustbeatlea leastao staone ne steptimecompuationaldelayforinformationtopassfrominputtooutput.Theabove LabVI abVIEW EW prog progra ram m can can be con converte erted d into a sub .vi .vi and and used used as many any tim times as necessaryprovideditisdefinedasre-entrant.Itcanbeusedasanintegrator,phaseleadorasalow-passfilter. 3.3.CompositePLLSimulation
TheLTIblockusedintheprevioussectioncanbeusedtoconstructanintegratorand phase-le phase-leadcompensa adcompensator torforthePLL(iet forthePLL(ietheFil heFilter).It ter).Itisnecessa isnecessaryto rytocal calcul culatethe atethegain gain val values and and any any param paramete eters. rs. One One possi possibl ble e desi design gn for a band bandwi width dth of 400Hz 400Hz give gives s a phase-leadL(s)of
L( s ) = 10
( s + 794.2) ( s + 7942 )
sothatfortheLTIblocka=7942,g=10,b0=794.2andb1=1.Theintegratorhasagain whi which can can be found ound to be appr approx oxiimatel ately 2X10 2X106 but but thi this does does not not accou accounnt for any any existin existing g(hi (hiddenor ddenorim impli plicit)gain cit)gainsalready salready in theloop.These theloop. These hiddengain hiddengain terms consist consist oftheVCOgain(12566)andthePhasede oftheVCOgain(12566)andthePhasedetectorgain(0.5),atotalof tectorgain(0.5),atotalof6283.2.Div 6283.2.Dividi iding ng this this value value into theover the overal alllgai gain ngi gives ves a remai remaini ning nggai gain nof of318 318.141 .141and and itis it is this this gain gain whichmustbeusedontheintegratorie318.14/s.ThisisshowninFigure9below.
Figure9.PartofthePLLshowingthefilterdynamics .
11
Thegainadjustparameterisnominallysettounitybutcanbevariedtoseetheeffect ofvaryingtheoverallgainoftheloop.ThecompletePLLFMdemodulatorisshowin Figure10below. The The trans transiient ent respon response se of the loop is found ound by FM modul odulating ating a squa square re wave. wave. An extern external al first-ord first-order erfi filte lter rsetat setat a900H a 900Hz z cut-offfrequ cut-offfrequenc encywas ywascons constructe tructed dusi usingthe ngthe LTIblockandisanintegralpartofthePLLblock.Thecut-offfrequencycanbevaried fromthefrontpanel.TheresultingtransientresponseisshowninFigure11belowfor a10Hzbaseband a10Hzbaseband signa signal.l. Adepthofmodulati Adepthofmodulationof onof+/-200Hz +/-200Hz was usedo usedonthe2kHz nthe2kHz carrier.TheFMmodulationindexforthisexampleistherefore = ∆ / m =20. Iftheexternalfil Iftheexternalfilterisset terissetto toahigh ahighvalue(4kHz)sothatitisin value(4kHz)sothatitisineff effectiv ectiveandthegain eandthegainof of theloopisincreasedsignificantlytheeffectoftwicethefrequencyofthecarrierfeedthroughcanbeillustratedinFigure12.
Figure10MainsimulationloopforPLLandFMgenerator.
12
Figure11TransientresponseofPLLtoa10Hzsquarewave .
Figure12Illustratesfeed-throughtermforhighbandwidth andnoexternalfiltering.
13
Itisimportanttonotethatintheabovesimulationsthereisnohardlimiterpresentas wouldbeina wouldbeinareal realradioreceive radioreceiver.SincetheFMhasa r.SincetheFMhasaconstantamp constantampli litudeofunitythere tudeofunitythere is no need. eed. Howe Howeve ver, r, if addi additi tive ve noi noise or co-ch co-chan anne nell inter nterfferen erence ce is adde added d to the the simulationitisthennecessarytoaddeitherahardlimiterorsomeformofAutomatic Gai Gain Control Control (AGC) (AGC).. Thi This is beca becaus use e the phas phase e detec detector tor is a linear near multi ultipl pliier and and hence hence any anycha change nge in ampl amplitu itude deof ofthe theFM FM signa signallwi will ll resultina resultina change change ofloop-gain ofloop-gain.. Thiscanleadtoinstabilityasthegaincanbeeithertoolowortoohigh.(seefigure3) 4.Conclusionsandfurtherwork
A sim simulat ulatiion of a PLL PLL usi using Lab LabVIEW VIEW has been een desc descri ribbed in som some dept depthh. The The sim simulati ulation on diff differs ers from fromprev previou ious s studie studies s in thatit that itiis full fullyy interac interactiv tive e and not a ‘one‘oneshot’ shot’ type type sim simulati ulation. on.Unl Unlik ike e a real PLL all allof ofthe thepara parame meters ters can bevar be varie ied din in realrealtime time to see theef the effe fect. ct. For exam exampl ple e thegai the gain nis is easil easilyy increa increased sed or decreas decreased edan and dthe the externalfi externalfiltercut-offcanbe ltercut-offcanbevari varied.At ed.At present presentforsim forsimpli plicity citytthePLLisdesi hePLLisdesignedfora gnedfora givenfixedcarrierfrequencyandbandwidthbutthattoocanbemadeinteractivewith thedesign equationsenteredas equationsenteredas equationsinLabVIEW equationsinLabVIEWandcontinuousl andcontinuouslyupdated.The yupdated.The PLLbehavesinan PLLbehavesinanide identic nticalfash alfashionto ionto anICPLLandassuchisid anICPLL andassuchisideal ealfor forin investi vestigatin gating g suchpropertiesasco-channelinterferenceinFM,thresholding,multi-pathandsoon.
5.References.
[1]Specialissueonphase-lockedloops.,IEEETrans.onCommunications,vol.COM20,No.10,Oct.1982. [2]H.deBellescise,Laréceptionsynchrone,OndeÉlectrique,vol.11,1932. [3] [3] T.J T.J.Moi .Moir, r, Sim Simula ulation tion of a Pha Phase-Loc -Lockked Loop oop usi using Matri trixx. Ele Electron troniic EngineeringJune1995,pp45-49 [4]F.M.Gardner,Phase-lockLoopTechniques,NewYorkWiley1979. [5] [5] T.J.M T.J.Moi oir, r, Digi Digital tal contro controll syst system em compe compennsati sation on usi using the the span span rati ratio, o, Jour Journa nall A,vol31,1,1990,pp65-69