costas loop is a classical phase - locked loop ( pll ) based circuit for carrier recovery and signal demodulation @xcite .
the pll is an automatic control system which is designed to generate an electrical signal ( voltage ) , the frequency of which is automatically tuned to the frequency of the input ( reference ) signal .
various pll based circuits are widely used in modern telecommunications , computer architectures , electromechanical systems ( see , e.g. @xcite ) . nowadays among the applications of costas loop
there are global positioning systems ( see , e.g. , @xcite ) , wireless communication ( see , e. g. , @xcite ) and others ( @xcite ) .
dynamic behavior of the pll and the costas loop has been described extensively in the literature @xcite , and a number of key parameters has been defined that describe its lock - in and lock - out characteristics . when the pll is initially out of lock ,
two different types of acquisition processes can occur , either the so - called lock - in process or the so - called pull - in process .
the first of those is a fast process , i.e. the acquisition takes place within at most one beat note of the difference between reference frequency @xmath0 and initial vco ( voltage controlled oscillator ) frequency @xmath1 , cf .
figure  [ f1 - 1 ] for signal denotations .
the frequency difference for which such a fast acquisition process takes place corresponds to the lock - in range @xmath2 , and the duration of the locking process is called lock time @xmath3 . when the difference between reference and vco frequency is larger than the lock - in range but less than the pull - in range @xmath4 , a slow acquisition process occurs .
the time required to get acquisition is called pull - in time @xmath5 . in case of the pll
all these acquisition parameters can be approximated by characteristic parameters of the pll , i.e. from natural frequency @xmath6 and damping factor @xmath7 .        to the authors knowledge
such acquisition parameters have not been analytically derived for the different types of costas loops .
it seems that most authors only described the static properties of the costas loop such as the derivation of the phase error in the locked state and the like .
based on methods developed earlier for the pll , the authors could now derive similar expressions for all relevant acquisition parameters of the costas loop .
this enables the designer to determine the lock - in and pull - in ranges , and to estimate the duration of the corresponding processes .
because the systems considered are highly non linear , exact computation of such parameters is very difficult or even impossible .
therefore it is necessary to introduce a number of simplifications .
this implies that the obtained results are only approximations , in some cases rather crude approximations .
as will be shown in the following sections there are different types of costas loops .
the very first of these loops has been described by j. costas in 1956 @xcite and was primarily used to demodulate amplitude - modulated signals with suppressed carrier ( dsb - am ) .
the same circuit was used later for the demodulation of bpsk signals ( binary phase shift keying ) @xcite . with the advent of qpsk ( quadrature phase shift keying )
this costas loop was extended to demodulate qpsk signals as well .
these two types of costas loop operated with real signals . in case of bpsk ,
the input signal @xmath8 is a sine carrier that was phase modulated by a binary signal , i.e. @xmath9 where @xmath0 is the ( radian ) carrier frequency , and @xmath10 can have two different values , either @xmath11 or @xmath12 , or two arbitrary equal and opposite values @xmath13 and @xmath14 , where @xmath15 can be any value . in case of qpsk ,
two quadrature carriers are modulated by two modulating signals , i.e.    @xmath16    where @xmath17 and @xmath18 can both have two equal and opposite values @xmath13 and @xmath14 .
it is obvious that in both cases the input signal is a real quantity . in the following these two types of costas loop
will be referred to as `` conventional costas loops '' .
much later , costas loops have been developed that operate not on real signals , but on pre - envelope signals @xcite .
these types of costas loops will be referred to as `` modified costas loop '' in the following sections .
the block diagram shown in figure  [ f1 - 2 ] explains how the pre - envelope signal is obtained .
the real input signal @xmath8 is applied to the input of a hilbert transformer [ 2 ] , [ 5 ] .
the output of the hilbert transformer @xmath19 is considered to be the imaginary part of the pre - envelope signal , i.e the pre - envelope signal is obtained from @xmath20 the costas loops operating with pre - envelope signals will be referred to as `` modified costas loops '' , cf .
sections [ s4 ] and [ s5 ] .        because there are different types of costas loops
the acquisition parameters must be derived separately for each of these types .
this will be performed in the following sections . in order to see how good or bad the obtained approximations
, we will develop simulink models for different types of costas loops and compare the results of the simulation with those predicted by theory .
the operation of the costas loop is considered first in the locked state with zero phase difference ( see figure  [ costas_locked ] ) , hence the frequency of the carrier is identical with the frequency of the vco .        by the input signal @xmath8 is the product of a transferred binary data and the harmonic carrier @xmath21 with a high frequency @xmath22 . since the costas loop is considered to be locked , the vco orthogonal output signals are synchronized with the carrier ( i.e. there is no phase difference between these signals ) .
the input signal is multiplied ( multiplier block ( @xmath23 ) ) by the corresponding vco signal on the upper branch and by the vco signal , shifted by @xmath24 , on the lower branch .
therefore on the multipliers outputs one has @xmath25    consider the low - pass filters ( lpf ) operation .    [ as - twice - frequency ] signals components , whose frequency is about twice the carrier frequency , do not affect the synchronization of the loop ( since they are suppressed by the low - pass filters ) .
[ as - lpf - initstate ] initial states of the low - pass filters do not affect the synchronization of the loop ( since for the properly designed filters , the impact of filter s initial state on its output decays exponentially with time ) .
[ as - lpf - data ] the data signal @xmath10 does not affect the synchronization of the loop .    assumptions  [ as - twice - frequency],[as - lpf - initstate ] , and [ as - lpf - data ] together lead to the concept of so - called _ ideal low - pass filter _ , which completely eliminates all frequencies above the cutoff frequency ( assumption  [ as - twice - frequency ] ) while passing those below unchanged ( assumptions  [ as - lpf - initstate],[as - lpf - data ] ) . in the classic engineering theory of the costas loop
it is assumed that the low - pass filters lpf are ideal low - pass filters .    since in figure
[ costas_locked ] the loop is in lock , i.e. the transient process is over and the synchronization is achieved , by assumptions  [ as - twice - frequency],[as - lpf - initstate ] , and [ as - lpf - data ] for the outputs @xmath26 and @xmath27 of the low - pass filters lpf one has @xmath28 thus , the upper branch works as a demodulator and the lower branch works as a phase - locked loop .
since after a transient process there is no phase difference , a control signal at the input of vco , which is used for vco frequency adjustment to the frequency of input carrier signal , has to be zero : @xmath29 . in the general case when the carrier frequency @xmath22 and a free - running frequency @xmath30 of the vco are different
, after a transient processes the control signal at the input of vco has to be non - zero constant : @xmath31 , and a constant phase difference @xmath32 may remain .
consider the costas loop before synchronization ( see figure  [ costas_before_sync ] ) . here
the phase difference @xmath33 varies over time , because the loop has not yet acquired lock ( frequencies or phases of the carrier and vco are different ) .        in this case , using assumption  [ as - twice - frequency ] , the signals @xmath34 and @xmath35 can be approximated as @xmath36 approximations depend on the phase difference of signals , i.e. two multiplier blocks ( @xmath23 ) on the upper and lower branches operate as phase detectors .
the obtained expressions with @xmath37 coincide with well - known ( see , e.g. , @xcite ) phase detector characteristic of the classic pll with multiplier / mixer phase - detector for sinusoidal signals .    by assumptions
[ as - lpf - initstate ] and [ as - lpf - data ] the low - pass filters outputs can be approximated as @xmath38 since @xmath39 , the input of the loop filter ( lf ) is @xmath40 such an approximation is called a _ phase detector characteristic of the costas loop_.    since an ideal low - pass filter is hardly realized , its use in the mathematical analysis requires additional justification .
thus , the impact of the low - pass filters on the lock acquisition process must be studied rigorously .
the relation between the input @xmath41 and the output @xmath42 of the loop filter has the form @xmath43 where @xmath44 is a constant matrix , the vector @xmath45 is the loop filter state , @xmath46 are constant vectors , h is a number .
the filter transfer function has the form : @xmath47 the control signal @xmath42 is used to adjust the vco frequency to the frequency of the input carrier signal @xmath48 here @xmath30 is the free - running frequency of the vco and @xmath49 is the vco gain . the solution of with initial data @xmath50 ( the loop filter output for the initial state @xmath50 ) is as follows @xmath51 where @xmath52 is the impulse response of the loop filter and @xmath53 is the zero input response of the loop filter , i.e. when the input of the loop filter is zero .
[ as - lf - initstate ] zero input response of loop filter @xmath54 does not affect the synchronization of the loop ( one of the reasons is that @xmath54 is an exponentially damped function for a stable matrix @xmath44 ) .
consider a constant frequency of the input carrier : @xmath55 and introduce notation @xmath56 then assumption  [ as - lf - initstate ] allows one to obtain the classic mathematical model of pll - based circuit        in signal s phase space ( see figure  [ f2 - 1b ] ) : @xmath57    for the locked state a linear pll model can be derived , which is shown in figure  [ linear - pll ] .
this model is useful for approximation of hold - in range .        in the locked state
both reference and vco frequencies are approximately the same , hence the input of the lowpass filter is a very low frequency signal .
therefore the lowpass filter can be ignored when setting up the linear model of the costas loop .
the linear model is made up of three blocks , the phase detector pd , the loop filter lf and the vco . in digital costas
loops the vco is replaced by a dco ( digital controlled oscillator ) .
this will be discussed in later sections . for these building blocks
the transfer functions are now defined as follows
.    * phase detector ( pd ) .
* in the locked state , the phase error @xmath58 is very small so by we can write @xmath59 with @xmath60 called phase detector gain .
@xmath61    note that the uppercase symbols are laplace transforms of the corresponding lower case signals .
* loop filter ( lf ) .
* for the loop filter we choose a pi ( proportional + integral ) filter whose transfer function has the from @xmath62 this filter type is the preferred one because it offers superior performance compared with lead - lag or lag filters .
* vco*. the transfer function of the vco is given by @xmath63 where @xmath49 is called vco gain .
consider another non linear model of costas loop in figure  [ delay - costas ] ( delay model ) .
here we use assumtions 13 ( initial states of filters are omitted , double - frequency terms are completely filtered by lpfs , and @xmath10 does nt affect synchronization ) and filters lpfs are replaced by the corresponding phase - delay blocks @xmath64 .
outputs of low - pass filters are @xmath65 where @xmath66 then after multiplication of @xmath26 and @xmath27 we have @xmath67 and the output @xmath42 of the loop filter satisfies the following equations @xmath68 equations of costas loop in this case are @xmath69 for lpf transfer functions @xmath70 phase shift is equal to @xmath71 .
therefore is equal to the following system @xmath72 where @xmath73 .
equation is hard to analyze both numerically and analytically , however this model is still useful . in the following discussion
it is used to approximate pull - in range and pull - in time . for this purpose
we need to simplify delay model shown in figure  [ delay - costas ] .
consider block diagram in figure  [ f2 - 4a ] .
the lowpass filters ( lpf ) used in both @xmath74 an @xmath75 branches are assumed to be first order filters having transfer function . as will be demonstrated later the corner frequency of these filters must be chosen such that the data signal i is recovered with sufficient accuracy , i.e. the corner frequency @xmath76 must be larger than the symbol rate .
typically it is chosen twice the symbol rate , i.e. @xmath77 with @xmath78 symbol rate and @xmath79 .
the output signal @xmath80 of the multiplier in the @xmath74 branch consists of two terms , one having the sum frequency @xmath81 and one having the difference frequency @xmath82 .
because the sum frequency term will be suppressed by the lowpass filter , only the difference term is considered . the same holds true for signal @xmath83 in the @xmath75 branch .
it will show up that the range of difference frequencies is markedly below the corner frequency @xmath76 of the lowpass filter .
hence the filter gain will be nearly 1 for the frequencies of interest .
as will also be shown later the phase at frequency @xmath84 can not be neglected .
the lowpass filter is therefore represented as a delay block whose transfer function has the value @xmath85 , where @xmath86 is the phase at frequency @xmath87 . the delayed signals @xmath88 and @xmath89 are now multiplied by the product block at the right in the block diagram .
consequently the output signal @xmath41 of this block will have a frequency of @xmath90 .
this signal is now applied to the input of the loop filter lf .
its transfer function has been defined in ( [ loop - filter - tf ] ) .
the corner frequency of this filter is @xmath91 .
because the phase of the loop filter can not be neglected , it is represented as a delay block characterized by @xmath92 where @xmath93 is the phase of the loop filter at frequency @xmath90 .
the analysis of dynamic behavior becomes easier when the order of some blocks in figure  [ f2 - 4a ] is reversed ( see figure  [ f2 - 4b ] ) , i.e. when we put the multiplying block before the lowpass filter .
because the frequency of signal @xmath41 in figure  [ f2 - 4a ] is twice the frequency of the signals @xmath88 and @xmath89 , the phase shift created by the lowpass filter at frequency @xmath90 is now twice the phase shift at frequency @xmath87 .
the lpf is therefore represented here by a delay block having transfer function @xmath94 .
we can simplify the block diagram even more by concatenating the lowpass filter and loop filter blocks .
the resulting block delays the phase by @xmath95 .
this is shown in figure  [ f2 - 4c ] .
the output signal @xmath42 of this delay block now modulates the frequency generated by the vco .        to compute pull - in time we need to consider costas loop model in figure  [ f2 - 1b ] with averaged signals of phase detector output @xmath96 and filter output @xmath97
( see figure  [ f2 - 7 ] ) .
the model is built from three blocks .
the first of these is labeled `` phase - frequency detector '' .
we have seen that in the locked state the output of the phase detector depends on the phase error @xmath58 . in the unlocked state ,
however , the average phase detector output signal @xmath98 is a function of frequency difference as will be shown in next section ( eq .  ( [ 20 ] ) ) , hence it is justified to call that block `` phase - frequency detector '' . as we will recognize the pull - in process
is a slow one , i.e. its frequency spectrum contains low frequencies only that are below the corner frequency @xmath99 of the loop filter , cf .
( [ loop - filter - tf ] ) . the loop filter can therefore be modeled as a simple integrator with transfer function @xmath100 therefore @xmath101 the frequency @xmath1 of the vco output signal is defined as @xmath102 where @xmath30 is the free running frequency and @xmath49 is the vco gain .
now we define the instantaneous frequency difference @xmath87 as @xmath103 substituting and ( [ 24 ] ) into ( [ 22 ] ) finally yields @xmath104      consider qpsk costas loop operation ( see figure  [ costas_after_sync ] ) for the sinusoidal carrier and vco in lock state for the same initial frequencies @xmath105 .        by , the input qpsk signal has the form @xmath106 where @xmath107 is the transmitted data , @xmath21 and @xmath108 are sinusoidal carriers , @xmath109  phase of input signal .
the outputs of the vco are @xmath110 and @xmath111 .
after multiplication of vco signals and the input signal by multiplier blocks ( @xmath23 ) on the upper @xmath74 branch one has @xmath112 on the lower branch the output signal of vco is multiplied by the input signal : @xmath113    here from an engineering point of view , the high - frequency terms @xmath114 and @xmath115 are removed by ideal low - pass filters lpfs ( see assumption [ as - twice - frequency ] in previews section ) . in this case , the signals @xmath26 and @xmath27 on the upper and lower branches can be approximated as @xmath116    apart from considered case there are two possible cases : 1 ) the frequencies are different or 2 ) the frequencies are the same but there is a constant phase difference .
consider costas loop before synchronization ( see figure  [ qpsk_before_sync ] )        in the case when the phase of the input carrier @xmath117 and the phase of vco @xmath118 are different : @xmath119    in this case , using assumption 1 , the signals @xmath26 and @xmath27 on the upper and lower branches can be approximated as @xmath120    after the filtration , both signals , @xmath34 and @xmath35 , pass through the limiters ( sgn blocks ) .
then the outputs of the limiters @xmath121 and @xmath122 are multiplied with @xmath27 and @xmath26 , respectively . by assumption  2 and corresponding formula ,
the difference of these signals @xmath123 can be approximated as @xmath124 with @xmath125 . here
@xmath126 is a piecewise - smooth function depends on @xmath127 at the points @xmath128 . ] shown in figure  [ qpsk - char ] .    ]    the resulting signal @xmath129 , after the filtration by the loop filter , forms the control signal @xmath42 for the vco .    to derive _ mathematical model in the signal space _ describing _ physical model _ of qpsk costas loop one takes into account and : @xmath130 however equations are nonlinear and non autonomous with discontinuous right - hand side , which are extremely hard to investigate .
therefore , the study of is outside of the scope of this work .    to derive linear model , we consider and the corresponding figure  [ qpsk - char ] .
the curve looks like a `` chopped '' sine wave .
the costas loop can get locked at four different values of @xmath58 , i.e. with @xmath131 , @xmath132 , or @xmath133 . to simplify the following analysis
, we can define the phase error to be zero wherever the loop gets locked .
moreover , in the locked state the phase error is small , so we can write @xmath134 i.e. the output signal of the adder block at the right of figure  [ qpsk_before_sync ] is considered to be the phase detector output signal @xmath96 .
the phase detector gain is then @xmath135 it is easily seen that the linear model for the locked state is identical with that of the costas loop for bpsk , cf .
figure  [ linear - pll ] . because only small frequency differences are considered here
, the lowpass filters can be discarded .
the transfer functions of the loop filter and of the vco are assumed to be the same as in the case of the costas loop for bpsk , hence these are given by eqs .
( [ loop - filter - tf ] ) and ( [ 6 ] ) .
similar to bpsk costas loop , it is reasonable to consider delay model of qpsk costas loop ( see figure  [ delay - qpsk - costas ] ) .
filters lpfs are replaced by the corresponding phase - delay blocks @xmath136 .
the outputs of low - pass filters are @xmath137 then @xmath41 can be approximated as @xmath138 consider the loop filter transfer function .
equations of delay model of qpsk costas loop in this case are @xmath139    the non linear model of the costas loop for qpsk is developed on the basis of the non linear model we derived for the costas loop for bpsk , cf .
figure  [ f2 - 4c ] . here
again the order of lowpass filters and the blocks shown at the right of figure  [ costas_after_sync ] is reversed .
this results in the model shown in figure  [ f3 - 4]a .        in the block labeled `` b '' the function blocks at the right of figure  [ f3 - 4]a have been integrated , cf . figure  [ f3 - 4]b .
the output signal @xmath96 of block b is the `` chopped '' sine wave as shown in figure  [ f3 - 2 ] .
its fundamental frequency is 4 times the frequency difference @xmath140 .
the lowpass filters and the loop filter have been concatenated in the block labeled `` lpf + lf '' at the right of figure  [ f3 - 4]a .
referring to figure  [ costas_after_sync ] signals @xmath80 and @xmath83 are passed through lowpass filters . as in the case of the costas loop for bpsk
we assume here again that the difference frequency @xmath87 is well below the corner frequency @xmath76 of the lowpass filters , hence the gain of the lowpass filters is nearly 1 at @xmath141 . because the phase shift must not be neglected , we represent the lowpass filter by a delay , i.e. its frequency response at @xmath141 is @xmath142 where @xmath86 is the phase of the lowpass filter . due to the arithmetic operations in block `` b '' ( cf .
figure  [ f3 - 4 ] ) the frequency of the @xmath96 is quadrupled , which implies that the phase shift at frequency @xmath143 becomes @xmath144 .
the frequency response of the loop filter at @xmath145 is given by @xmath146 where @xmath93 is the phase of the loop filter at frequency @xmath147 .
hence the cascade of lowpass filter and loop filter can be modeled by the transfer function @xmath148)$ ] as shown in figure  [ f3 - 4]a .
the block diagram of the modified costas loop for bpsk is shown in figure  [ f4 - 1 ] .
the input signal is given by @xmath149 where @xmath150 is initial phase .
the input signal is first converted into a pre - envelope signal , as explained in section [ ss1.1 ] .
the output signal of the hilbert transformer is @xmath151 = m_1(t)\sin(\omega_1t+\theta_1).\ ] ]    note that because the largest frequency of the spectrum of the data signal @xmath10 is much lower than the carrier frequency @xmath0 , the hilbert transform of the product @xmath152 $ ] equals @xmath153 $ ] [ 5 ] .
the pre - envelope signal is obtained now from @xmath154).\ ] ]    the exponential in eqn .
( [ 56 ] ) is referred to as a `` complex carrier '' . in figure
[ f4 - 1 ] complex signals are shown as double lines .
the solid line represents the real part , the dotted line represents the imaginary part .
to demodulate the bpsk signal , the pre - envelope signal is now multiplied with the output signal of the vco , which is here a complex carrier as well .
the complex output signal of the vco is defined as @xmath155).\ ] ]    in the locked state of the costas loop both frequencies @xmath0 and @xmath1 are equal , and we also have @xmath156 . hence the output signal of the multiplier @xmath157 is @xmath158)\approx m_1(t),\ ] ] i.e. the output of the multiplier is the demodulated data signal @xmath10 . to derive the linear model of this costas loop
, it is assumed that @xmath159 and @xmath160 .
the output signal of multiplier @xmath157 then becomes @xmath161).\ ] ]     in the complex plane ]    this is a phasor having magnitude @xmath162 and phase @xmath163 , as shown in figure  [ f4 - 2 ] .
two quantities are determined from the phase of phasor @xmath164 , i.e. the demodulated data signal i and the phase error @xmath58 .
the data signal is defined as @xmath165),\ ] ] i.e. when the phasor lies in quadrants i or iv , the data signal is considered to be + 1 , and when the phasor is in quadrants ii or iii , the data signal is considered to be -1 .
this means that i can be either a phasor with phase 0 or a phasor with phase @xmath132 .
these two phasors are plotted as thick lines in figure  [ f4 - 2 ] .
the phase error @xmath58 is now given by the difference of the phases of phasor @xmath164 and phasor @xmath74 , as shown in figure [ f4 - 2 ] , i.e. @xmath58 is determined from @xmath166    the product @xmath167 is computed by multiplier @xmath168 in figure  [ f4 - 2 ] .
the block labeled `` complex @xmath169 mag , phase '' is used to convert the complex signal delivered by @xmath168 into magnitude and phase .
the magnitude is not used in this case , but only the phase .
it follows from eqn .
( [ 61 ] ) that the phase output of this block is the phase error @xmath58 , hence the blocks @xmath170 , sgn , and complex @xmath169 mag , phase represent a phase detector with gain @xmath171 .
the phase output of block complex @xmath169 mag , phase is therefore labeled @xmath96 . figure  [ linear - pll ] shows the complete linear model of the modified costas loop for bpsk .
the transfer functions of the loop filter and vco have been defined in eqs .
( [ loop - filter - tf ] ) and ( [ 6 ] ) .
note that with this type of costas loop there is no additional lowpass filter , because the multiplication of the two complex carriers ( cf .
( [ 58 ] ) ) does not create the unwanted double frequency component as found with the conventional costas loops .
figure  [ f5 - 1 ] shows the block diagram of the modified costas loop for qpsk .
the reference signal @xmath8 is defined by @xmath172 the hilbert transformed signal is then given by @xmath173 and the pre - envelope signal then becomes @xmath174    this can be rewritten as @xmath175+j\sin[\omega_1t+\theta_1]=\\ & = ( m_1(t)+jm_2(t))\exp(j[\omega_1t+\theta_1 ] ) .
\end{aligned}\ ] ] herein the term @xmath176 is complex envelope , and the term @xmath177 is complex carrier .
the vco generates another complex carrier given by .
the multiplier @xmath157 creates signal @xmath164 that is given by @xmath178).\ ] ]    when the loop has acquired lock , @xmath179 , and @xmath180 , so we have @xmath181 hence the output of @xmath157 is the complex envelope . in the locked state ,
the complex envelope can take four positions , as shown in figure  [ f5 - 2 ] .
when there is a phase error , @xmath164 deviates from the ideal position , as demonstrated in the figure .
the phase error @xmath58 then is the angle between @xmath164 and the closest of the four possible positions .
when @xmath164 is in quadrant i , e.g. , phasor @xmath182 is considered the estimate of the complex envelope .
when @xmath164 is in quadrant ii , the estimate of the complex envelope is @xmath183 etc .
the estimates @xmath74 and @xmath75 are taken from the output of sgn blocks , cf .
figure  [ f5 - 1 ] .
the phase error is obtained from @xmath184\ ] ] where @xmath185 is the conjugate of the complex envelope .
multiplier @xmath168 delivers the product @xmath186 , and the block `` complex @xmath169 mag , phase '' is used to compute the phase of that complex quantity . note that the magnitude is not required .
the blocks @xmath157 , sgn , inverter , @xmath168 , and complex @xmath169 mag , phase form a phase detector having gain @xmath171 .
the phase output of block complex @xmath169 mag , phase is therefore labeled @xmath96 .     in the complex plane , scaledwidth=50.0% ]
figure  [ linear - pll ] shows the completed linear model of the modified costas loop for qpsk , which is the same as for bpsk .
the transfer functions of the loop filter and vco have been defined in eqs .
( [ loop - filter - tf ] ) and ( [ 6 ] ) .      in the classic books on phase - locked loops @xcite such concepts as hold - in pull - in lock - in and other frequency ranges for which pll can achieve lock were introduced .
usually in engineering literature non - rigorous definitions are given for these concepts . in the following we introduce definitions , based on rigorous discussion in @xcite .
* definition of hold - in range . * the largest interval @xmath187 of frequency deviations @xmath188 , such that the loop re - achieves locked state after small perturbations of the filters state , the phases and frequencies of vco , and the input signals ,
is called a hold - in range .
this effect is also called steady - state stability . in addition , for a frequency deviation within the hold - in range , the loop in a locked state tracks small changes in input frequency , i.e. achieves a new locked state ( _ tracking process _ ) @xcite .
assume that the loop power supply is initially switched off and then at @xmath189 the power is switched on , and assume that the initial frequency difference is sufficiently large .
the loop may not lock within one beat note , but the vco frequency will be slowly tuned toward the reference frequency ( acquisition process ) .
this effect is also called a transient stability . the pull - in range
is used to name such frequency deviations that make the acquisition process possible .
* definition of pull - in range . * the largest interval @xmath190 of frequency deviations @xmath188 , such that the loop achieves locked state for any initial states ( filters and initial phase of vco ) ,
is called a pull - in range @xcite . the largest frequency deviation @xmath191
is called a pull - in frequency@xcite .
* definition of lock - in range .
* lock - in range is a largest interval of frequency deviations @xmath192 inside pull - in range , such that after an abrupt change of @xmath0 within a lock - in range the pll re - acquires lock without cycle slipping , if it is not interrupted . here
@xmath2 is called a lock - in frequency@xcite .
finally , our definitions give @xmath193 @xmath194 which is in agreement with the classical consideration ( see , e.g. @xcite,@xcite,@xcite,@xcite,@xcite ) .
recall linear model of costas loop in phase space ( see figure  [ linear - pll ] ) . by , , and we can derive the open loop transfer function of the costas loop , which is defined by the ratio @xmath195 : @xmath196     ]    figure  [ f2 - 2 ] shows a bode plot of the magnitude of @xmath197 .
the plot is characterized by the corner frequency @xmath99 , which is defined by @xmath198 , and gain parameters @xmath60 and @xmath49 . at lower frequencies the magnitude rolls off with a slope of  40 db / decade . at frequency
@xmath99 the zero of the loop filter causes the magnitude to change its slope to  20 db / decade . to get a stable system , the magnitude curve should cut the 0 db line with a slope that is markedly less than  40 db / decade .
setting the parameters such that the gain is just 0 db at frequency @xmath99 provides a phase margin of 45 degrees , which assures stability [ 2 ] . from the open loop transfer function
we now can calculate the closed loop transfer function defined by @xmath199    after some mathematical manipulations we get @xmath200 it is customary to represent this transfer function in normalized form , i.e. @xmath201 with the substitutions @xmath202 where @xmath6 is called natural frequency and @xmath7 is called damping factor .
the linear model enables us to derive simple approximations for lock - in range @xmath2 and lock time @xmath3 .    for the following analysis
we assume that the loop is initially out of lock .
the frequency of the input signal ( figure  [ costas_before_sync ] ) is @xmath0 , and the frequency of the vco is @xmath1 .
the multiplier in the i branch therefore generates an output signal consisting of a sum frequency term @xmath203 and a difference frequency term @xmath204 .
the sum frequency term is removed by the lowpass filter , and the frequency of the difference term is assumed to be much below the corner frequency @xmath76 of the lowpass filter , hence the action of this filter can be neglected for this case . under this condition
the phase detector output signal @xmath41 will have the form ( cf .
( [ ud ] ) and ( [ 3 ] ) ) @xmath205    with @xmath206 .
@xmath41 is plotted in figure  [ f2 - 3 ] , left trace .
this signal passes through the loop filter . in most cases
the corner frequency @xmath207 is much lower than the lock - in range , hence we can approximate its transfer function by @xmath208        let us define the gain of this filter at higher frequencies by constant @xmath209 . now
the output signal of the loop filter is a sine wave having amplitude @xmath60 @xmath210 as shown by the middle trace in figure  [ f2 - 3 ] .
consequently the frequency of the vco will be modulated as shown in the right trace .
the modulation amplitude is given by @xmath211 .
in this figure the reference frequency and the initial frequency @xmath30 of the vco are plotted as horizontal lines . when @xmath0 and @xmath30 are such that the top of the sine wave just touches the @xmath0 line , the loop acquires lock suddenly , i.e. the lock - in range @xmath2 is nothing more than the modulation amplitude @xmath60 @xmath49 @xmath210 . making use of the substitutions ( [ omega zeta ] ) we finally get @xmath212    now the lock process is a damped oscillation whose frequency is the natural frequency .
because the loop is assumed to lock within at most one cycle of that frequency , the lock time can be approximated by the period of the natural frequency , i.e. we have @xmath213      we have seen that all signals found in this block diagram are sine functions , i.e. all of them seem to have zero average , hence do not show any @xmath214 component .
this would lead to the ( erroneous ) conclusion that a pull - in process would not be possible . in reality it will be recognized that some of the signals become asymmetrical
, i.e. the duration of the positive half wave is different from the duration of the negative one .
this creates a non zero @xmath214 component , and under suitable conditions acquisition can be obtained .
we are therefore going to analyze the characteristics of the signals in figure  [ f2 - 4c ] .
, scaledwidth=50.0% ]    all considered signals are plotted in figure  [ f2 - 5 ] . for signals @xmath80 and @xmath83 we obtain @xmath215 @xmath216 the sum frequency terms are discarded because they are removed by the lowpass filter .
the signal @xmath41 is the product of @xmath80 and @xmath83 and is given by . for small arguments @xmath217
this can be written as @xmath218 where @xmath219 .
because the phase detector gain is defined by @xmath220 we have @xmath221 .
next the loop filter output signal @xmath42 is plotted .
its amplitude is @xmath209 @xmath222 , and its phase is delayed by @xmath95 .
this signal modulates the frequency of the vco as shown in the bottom trace of figure  [ f2 - 5 ] .
the modulation amplitude is given by @xmath223 . in order to get an estimate for the non zero @xmath214 component of @xmath41 we will have to analyze the asymmetry of the signal waveforms .
it will be shown that @xmath98 ( the average of @xmath41 ) is a function of frequency difference @xmath87 and phase @xmath224 .
the analysis becomes easier when we first calculate @xmath98 for some special values for @xmath224 , i.e. for @xmath225 ; and @xmath226 .
let us start with @xmath227 , cf .
figure  [ f2 - 6a ] .     for @xmath227 , scaledwidth=50.0% ]    in figure
[ f2 - 6a ] the waveforms for @xmath41 and @xmath228 are shown .
the asymmetry of the signals is exaggerated in this plot . during the positive half cycle ( duration @xmath229 ) the average value of vco output frequency @xmath228 is increased , which means that the average difference frequency @xmath230 is lowered .
consequently the duration of the positive half wave becomes larger than half of a full cycle . during the negative half cycle ( duration @xmath231 ) ,
however , the average value of vco output frequency @xmath228 is decreased , which means that the average difference frequency @xmath232 is increased .
consequently the duration of the negative half wave becomes less than half of a full cycle .
next we are going to calculate the average frequency difference in both half cycles .
the average frequency difference during half cycle @xmath229 is denoted @xmath233 , the average frequency difference during half cycle @xmath231 is denoted @xmath234 .
we get @xmath235 @xmath236 for the durations @xmath229 and @xmath231 we obtain after some manipulations @xmath237 @xmath238 now the average value @xmath98 can be calculated from @xmath239 the average signal @xmath98 is seen to be inversely proportional to the frequency difference @xmath87 . because @xmath98 is positive , the instantaneous frequency @xmath228 is pulled in positive direction , i.e. versus @xmath0 , which means that a pull - in process will take place .     for @xmath240 , scaledwidth=50.0% ]    next we are going to analyze the dependence of @xmath98 on phase @xmath224 .
let us consider now the case for @xmath240 , cf .
figure  [ f2 - 6b ] .
we observe that in interval @xmath229 the instantaneous frequency @xmath228 is pulled in negative direction , hence the average difference frequency @xmath233 becomes larger . consequently interval @xmath229 becomes shorter . in interval @xmath231 , however , the reverse is true . here
the instantaneous frequency @xmath229 the pulled in positive direction , hence the average @xmath234 is reduced , and interval @xmath231 becomes longer .
the average @xmath98 is now equal and opposite to the value of @xmath98 for @xmath227 .
because it is negative under this condition , a pull - in process can not take place , because the frequency of the vco is `` pulled away '' in the wrong direction .
last we consider the case @xmath241 , cf . figure  [ f2 - 6c ] . in the first half of interval
@xmath229 the instantaneous frequency @xmath228 is lowered , but in the second half it is increased . consequently the average difference frequency @xmath233 does not change its value during @xmath229 .
the same happens in interval @xmath231 .
@xmath234 does not change either , and @xmath98 remains 0 .     for @xmath241 , scaledwidth=50.0% ]
it is easy to demonstrate that @xmath98 varies with @xmath242 , hence we have @xmath243 eq .
( [ 20 ] ) tells us that the pull - in range is finite . the pull - in range
can be found as the frequency difference for which phase @xmath241 .
an equation for the pull - in range will be derived in section [ ss2.4 ] .
we also will have to find an equation for the pull - in time .
the model shown in figure  [ f2 - 7 ] will enable us to obtain a differential equation for the average frequency difference @xmath244 as a function of time .
recall equations of filter output @xmath245 and frequency deviation @xmath246 eqs .
( [ 20 ] ) , ( [ 21b ] ) , and ( [ 25 ] ) enable us to compute the three variables @xmath98 , @xmath247 , and @xmath87 as a function of time .
this will be demonstrated in section [ ss2.4 ] .
the pull - in range can be computed using eq .
( [ 20 ] ) .
lock can only be obtained when the total phase shift @xmath224 is not more negative than @xmath248 .
this leads to an equation of the form @xmath249 according to eqs .
( [ loop - filter - tf ] ) and ( [ 14 ] ) @xmath86 and @xmath93 are given by @xmath250 @xmath251 with @xmath252 .
hence the pull - in range @xmath4 can be computed from the transcendental equation @xmath253    to solve this equation for @xmath4 we use the addition formula for the tangent function @xmath254 and can replace @xmath255 by @xmath256 .
( [ 27 ] ) can now be rewritten as @xmath257 .    when the arctg expressions on both sides of the equation are equal , their arguments must also be identical , which leads to @xmath258 hence we get for the pull - in range @xmath259 last an equation for the pull - in time @xmath5 will be derived .
( [ 20 ] ) , ( [ 21b ] ) , and ( [ 25 ] ) describe the behavior of the three building blocks in figure  [ f2 - 7 ] and enable us to compute the three variables @xmath98 , @xmath247 , and @xmath87 .
we only need to know the instantaneous @xmath87 vs. time , hence we eliminate @xmath98 and @xmath247 from eqs .
( [ 21b ] ) and ( [ 25 ] ) and obtain the differential equation @xmath260 this differential equation is non linear , but the variables @xmath87 and @xmath261 can be separated , which leads to an explicit solution . putting all terms containing @xmath87 to the left side and performing an integration
we get    @xmath262    the limits of integration are @xmath263 and @xmath2 on the left side , because the pull - in process starts with an initial frequency offset @xmath264 and ends when @xmath87 reaches the value @xmath2 , which is the lock - in range . following that instant a lock - in process will start
. the integration limits on the right side are @xmath265 and @xmath5 , respectively , which means that the pull - in process has duration @xmath5 , and after that interval ( fast ) lock - in process starts .
performing the integration on the left imposes some considerable problems , when we remember that @xmath266 is given by @xmath267 finding an explicit solution for the integral seems difficult if not impossible , but the @xmath268 term can be drastically simplified .
when we plot @xmath266 vs. @xmath87 we observe that within the range @xmath269 the term @xmath266 is an almost perfect straight line .
hence we can replace @xmath266 by @xmath270 inserting that substitution into eq .
yields a rational function of @xmath87 on the left side , which is easily integrated . after some mathematical procedures we obtain for the pull - in time @xmath5 @xmath271.\ ] ] making use of eqs .
and we have @xmath272 using these substitutions eq .
can be rewritten as @xmath273.\ ] ] this equation is valid for initial frequency offsets in the range @xmath274 . for lower frequency offsets ,
a fast pull - in process will occur , and eq . should be used .
an analog costas loop for bpsk shall be designed in this section .
it is assumed that a binary signal is modulated onto a carrier .
the carrier frequency is set to 400 khz , i.e. the costas loop will operate at a center frequency @xmath275 400000 = 2512000 rad @xmath276 . the symbol rate is assumed to be @xmath277 symbols / s .
now the parameters of the loop ( such as time constants @xmath278 and @xmath279 , corner frequencies @xmath99 and @xmath76 , and gain parameters such as @xmath280 ) must be determined .
( note that these parameters have been defined in eqs .
( [ 4 ] ) , ( [ loop - filter - tf ] ) , ( [ 6 ] ) and ( [ 13 ] ) ) .
the modulation amplitude is set @xmath281 . according to eq .
( [ 3 ] ) the phase detector gain is then @xmath171 .
it has proven advantageous to determine the remaining parameters by using the open loop transfer function @xmath282 of the loop [ 2 ] .
this is given by @xmath283        the magnitude @xmath284 ( bode diagram ) is plotted in figure  [ f2 - 9 ] .
the magnitude curve crosses the 0 db line at the so called transit frequency @xmath285 .
it is common practice to choose @xmath285 to be about @xmath286 . here
we set @xmath287 , i.e. @xmath288 rad @xmath276 .
furthermore we set corner frequency @xmath289 .
when doing so , the slope of the asymptotic magnitude curve changes from  40 db / decade to  20 db / decade at @xmath290 . under this condition
the phase of @xmath291 is -135@xmath292 at @xmath99 .
consequently the phase margin of the loop becomes 45@xmath292 , which provides sufficient stability . according to eq .
( [ loop - filter - tf ] ) @xmath279 becomes @xmath293 .
next corner frequency @xmath76 will be determined .
the corner frequency of the lowpass filter must be chosen such that the demodulated data signal ( i. e. the output of the lowpass filter in the i branch ) is recovered with high fidelity . to fulfill this requirement ,
@xmath76 should be chosen as large as possible . on the other hand
, the lowpass filter should suppress the double frequency component ( here at about 800 khz ) sufficiently , which means that @xmath76 should be markedly less than @xmath294 .
it s a good compromise to set corner frequency to twice the symbol rate , i.e. @xmath295 . last the remaining parameters @xmath278 and @xmath49 must be chosen .
they have to be specified such that the open loop gain becomes 1 at frequency @xmath290 . according to eq .
( [ 33 ] ) we can set @xmath296 because 2 parameters are still undetermined , one of those can be chosen arbitrarily , hence we set @xmath297 . finally from ( [ 34 ] )
we get @xmath298 .
the design of the costas loop is completed now , and we can compute the most important loop parameters . for the natural frequency and damping factor we get from @xmath299 from the lock - in range becomes @xmath300 and from ( [ 13 ] ) the lock time becomes @xmath301 next we want to compute the pull - in range .
( [ 28 ] ) yields @xmath302      to convert the analog loop into a digital one , we first must define a suitable sampling frequency @xmath303 ( or sampling interval @xmath304 ) . to satisfy the nyquist theorem , the sampling frequency must be higher than twice the highest frequency that exists in the loop . in our case
the highest frequency is found at the output of the multipliers in the @xmath74 and @xmath75 branches ( cf .
figure  [ costas_before_sync ] ) .
the sum frequency term is about twice the center frequency , hence @xmath303 must be greater than 4 times the center frequency . a suitable choice would be @xmath305 @xmath306 mhz .
next the transfer functions of the building block have to be converted into discrete transfer functions , i.e. @xmath307 . for best results it is preferable to use the bilinear @xmath308 transform .
given an analog transfer function @xmath309 , this can be converted into a discrete transfer function @xmath310 by replacing @xmath311 by @xmath312 now the bilinear z transform has the property that the analog frequency range from @xmath313 is compressed to the digital frequency range from @xmath314 .
to avoid undesired `` shrinking '' of the corner frequencies ( @xmath99 and @xmath76 ) , these must be `` prewarped '' accordingly , i.e. we must set @xmath315 @xmath316 where @xmath317 and @xmath318 are the prewarped corner frequencies .
now we can apply the bilinear z transform to the transfer functions of the lowpass filters ( cf .
( [ 14 ] ) ) and of the loop filter ( cf .
( [ loop - filter - tf ] ) ) and get @xmath319+\left[1-\frac{2}{\omega_{3,p}t}\right]z^{-1}}{1+z^{-1}},\ ] ] @xmath320+\left[1-\frac{2}{\omega_{c , p}t}\right]z^{-1 } } { \frac{2\tau_1}{t}-\frac{2\tau_1}{t}z^{-1}}.\ ] ] because the vco is a simple integrator , we can apply the discrete @xmath308 transform of an integrator , i.e. @xmath321 the digital costas loop is ready now for implementation .
a simulink model will be presented in section [ ss2.7 ] .      a simulink model of a costas loop for bpsk
is shown in figure  [ f2 - 10 ] .
a data signal is created by a random number generator at the left in the block diagram .
the other blocks are self explanatory .
the model is used now to check the validity of the approximations found for pull - in range and pull - in time .
( [ 28 ] ) predicts a pull - in range @xmath322 khz .
the simulation revealed a pull - in range of @xmath323 khz . , which shows that the theoretical result is a rather crude approximation .
a series of other simulation delivered results for the pull - in time @xmath324 .
the results are listed in table 2 - 1 .    [ cols="<,<,<,<",options="header " , ]     table 5 - 1 .
comparison of predicted and simulated results for the pull - in range    the predictions come very close to the results obtained from the simulation .
as demonstrated in figure  [ f5 - 1 ] the phase error signal @xmath96 was obtained from the phase output of block `` complex @xmath169 mag , phase '' .
the phase of the complex input signal to this block can be obtained from the @xmath325 function .
this imposes no problem when a processor is available .
this is the case in most digital implementations of the costas loop . as an alternative a phase error signal
can also be obtained directly from the imaginary part of multiplier @xmath168 ; this is shown in figure  [ f5 - 9 ] .
, scaledwidth=40.0% ]    it is easily seen that here @xmath96 is given by @xmath326    the blocks shown in figure  [ f5 - 9 ] therefore represent a phase detector having gain @xmath327 .
this must be taken in account when specifying the open loop transfer function , cf .
section [ ss5.5 ] . for this design
the pull - in time of the loop is given by @xmath328
this work was supported by russian science foundation ( project 14 - 21 - 00041 ) and saint - petersburg state university
one needs to be cautious using model in figure  [ f2 - 1b ] even for calculating hold - in range for bpsk costas .
consider an example : costas loop with lead - lag loop filter @xmath329 and low - pass filters lpfs @xmath330    in locked state phase error @xmath32 satisfies @xmath331 therefore we get a bound for the hold - in range @xmath332 in order to find hold - in range we need to find poles of the closed - loop transfer function ( roots of the characteristic polynomial ) for the linearized model ( small - signal model ) of the system on figure  [ costas_before_sync ] .
open - loop transfer function is @xmath333 @xmath334 phase error @xmath335 corresponds to hold - in range ( see ) if all roots of the polynomial have negative real parts ( i.e. polynomial is stable ) . applying routh - hurwitz criterion to study stability of the polynomial
, we get that for the following parameters @xmath336 polynomial is stable for all @xmath337
. however , if @xmath338 the following condition is necessary for stability @xmath339 then , taking into account static phase error equation , we get different hold - in ranges for different values of @xmath76 @xmath340              r.  lazzari , c.  parma , a.  de  marco , s.  bittanti , enabling a flexible exchange of energy of a photovoltaic plant with the grid by means of a controlled storage system , international journal of control 88  ( 7 ) ( 2015 ) 13531365 .
y.  a. sidorkina , v.  sizykh , b.  shakhtarin , v.  shevtsev , costas circuit under the action of additive harmonic interferences and wideband noise , journal of communications technology and electronics 61  ( 7 ) ( 2016 ) 807816 .
f.  bizzarri , a.  brambilla , g.  s. gajani , periodic small signal analysis of a wide class of type - ii phase locked loops through an exhaustive variational model , circuits and systems i : regular papers , ieee transactions on 59  ( 10 ) ( 2012 ) 22212231 .
g.  leonov , n.  kuznetsov , m.  yuldashev , r.  yuldashev , nonlinear dynamical model of costas loop and an approach to the analysis of its stability in the large , signal processing 108 ( 2015 ) 124135 . http://dx.doi.org/10.1016/j.sigpro.2014.08.033 [ ] .
r.  best , n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , simulation of analog costas loop circuits , international journal of automation and computing 11  ( 6 ) ( 2014 ) 571579 , 10.1007/s11633 - 014 - 0846-x .
n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , simulation of nonlinear models of modified bpsk costas loop for non sinusoidal waveforms in matlab simulink , in : 2014 6th international congress on ultra modern telecommunications and control systems and workshops ( icumt ) , vol . 2015-january , ieee , 2014 , pp .
http://dx.doi.org/10.1109/icumt.2014.7002084 [ ] .
n.  kuznetsov , o.  kuznetsova , g.  leonov , p.  neittaanmaki , m.  yuldashev , r.  yuldashev , simulation of nonlinear models of qpsk costas loop in matlab simulink , in : 2014 6th international congress on ultra modern telecommunications and control systems and workshops ( icumt ) , vol . 2015-january , ieee , 2014 , pp .
n.  kuznetsov , g.  leonov , p.  neittaanmki , s.  seledzhi , m.  yuldashev , r.  yuldashev , nonlinear mathematical models of costas loop for general waveform of input signal , in : ieee 4th international conference on nonlinear science and complexity , nsc 2012 - proceedings , 2012 , pp .
http://dx.doi.org/10.1109/nsc.2012.6304729 [ ] .
g.  leonov , n.  kuznetsov , m.  yuldahsev , r.  yuldashev , analytical method for computation of phase - detector characteristic , ieee transactions on circuits and systems - ii : express briefs 59  ( 10 ) ( 2012 ) 633647 . http://dx.doi.org/10.1109/tcsii.2012.2213362
n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , analytical methods for computation of phase - detector characteristics and pll design , in : isscs 2011 - international symposium on signals , circuits and systems , proceedings , ieee , 2011 , pp . 710 . http://dx.doi.org/10.1109/isscs.2011.5978639 [ ] .
n.  kuznetsov , o.  kuznetsova , g.  leonov , p.  neittaanmaki , m.  yuldashev , r.  yuldashev , limitations of the classical phase - locked loop analysis , proceedings - ieee international symposium on circuits and systems 2015-july ( 2015 ) 533536 . http://dx.doi.org/10.1109/iscas.2015.7168688 [ ] .
r.  best , n.  kuznetsov , o.  kuznetsova , g.  leonov , m.  yuldashev , r.  yuldashev , a short survey on nonlinear models of the classic costas loop : rigorous derivation and limitations of the classic analysis , in : proceedings of the american control conference , ieee , 2015 , pp . 12961302 , art .
7170912 , http://arxiv.org/pdf/1505.04288v1.pdf .
http://dx.doi.org/10.1109/acc.2015.7170912 [ ] .
n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , rigorous mathematical definitions of the hold - in and pull - in ranges for phase - locked loops , ifac - papersonline 48  ( 11 ) ( 2015 ) 710713 . http://dx.doi.org/10.1016/j.ifacol.2015.09.272 [ ] .
g.  leonov , n.  kuznetsov , m.  yuldashev , r.  yuldashev , hold - in , pull - in , and lock - in ranges of pll circuits : rigorous mathematical definitions and limitations of classical theory , ieee transactions on circuits and systems
i : regular papers 62  ( 10 ) ( 2015 ) 24542464 .
http://dx.doi.org/10.1109/tcsi.2015.2476295 [ ] .
g.  bianchi , n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , limitations of pll simulation : hidden oscillations in matlab and spice , international congress on ultra modern telecommunications and control systems and workshops ( icumt 2015 ) 2016-january ( 2016 ) 7984 .
http://dx.doi.org/10.1109/icumt.2015.7382409 [ ] .
n.  kuznetsov , o.  kuznetsova , g.  leonov , s.  seledzhi , m.  yuldashev , r.  yuldashev , bpsk costas loop : simulation of nonlinear models in matlab simulink , in : 2014 6th international congress on ultra modern telecommunications and control systems and workshops ( icumt ) , vol .
2015-january , ieee , 2014 , pp .
http://dx.doi.org/10.1109/icumt.2014.7002083 [ ] .
g.  leonov , n.  kuznetsov , hidden attractors in dynamical systems . from hidden oscillations in hilbert - kolmogorov ,
aizerman , and kalman problems to hidden chaotic attractors in chua circuits , international journal of bifurcation and chaos 23  ( 1 ) , art .
http://dx.doi.org/10.1142/s0218127413300024 [ ] .
n.  kuznetsov , g.  leonov , m.  yuldashev , r.  yuldashev , nonlinear analysis of classical phase - locked loops in signal s phase space , ifac - papersonline 47  ( 3 ) ( 2014 ) 82538258 . http://dx.doi.org/10.3182/20140824-6-za-1003.02772 [ ] .
e.  kudryashova , o.  kuznetsova , n.  kuznetsov , g.  leonov , s.  seledzhi , m.  yuldashev , r.  yuldashev , nonlinear models of bpsk costas loop , icinco 2014 - proceedings of the 11th international conference on informatics in control , automation and robotics 1 ( 2014 ) 704710 . http://dx.doi.org/10.5220/0005050707040710 [ ] .
g.  bianchi , n.  kuznetsov , g.  leonov , s.  seledzhi , m.  yuldashev , r.  yuldashev , http://www.sciencedirect.com/science/article/pii/s2405896316312617[hidden oscillations in spice simulation of two - phase costas loop with non - linear vco ] , ifac - papersonline 49  ( 14 ) ( 2016 ) 45  50 . http://dx.doi.org/10.1016/j.ifacol.2016.07.973 [ ] . http://www.sciencedirect.com/science/article/pii/s2405896316312617