Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1582  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

Application of Time Recursive Processing for the 
Development of a Time/Phase Shifter 

 

Djurdje M. Perisic  
Faculty of Information Technology  

Slobomir P University, SPU 
Slobomir, Bjeljina, Republic of Srpska, BIH 

djurdje@beotel.rs  

Miroslav Bojovic  
Faculty of Electrical Engineering 

University of Belgrade 
Belgrade, Serbia 

mbojovic@etf.bg.ac.rs 
 

 
Abstract—This paper describes a powerful digital time and/or 
phase shifter of pulse rates, based on Time Recursive Processing. 
It can function either as a time shifter or as a phase shifter. The 
circuit can generate precise shifting for a pulse rate with a 
constant input period, but it is also capable to shift a pulse rate 
whose period is a ramp function. Besides that, the shifter can be 
used in tracking and prediction applications. The shifter is 
described by recursive equations as a linear discrete system. All 
mathematical analyses are made using the Z transform. 
Computer simulations are employed to prove the correctness of 
the mathematical analysis. The realization of the shifter is 
described and, to demonstrate the shifter functioning, actual 
oscilloscope screenshots are presented.    

Keywords- time shifter; phase shifter; digital circuit; PLL   

I. INTRODUCTION  
The field of Time Recursive Processing (TRP) has been 

described in detail in previous works such as [1-9]. The TRP 
approach is based on the measurement and processing of the 
input and output periods and time differences between them. 
Although the first idea was simply to discover a new kind of 
Phase Locked Loop (PLL) and Frequency Locked Loop (FLL), 
it turned out that TRP is suitable for a much broader range of 
applications. Namely, most of the circuits described in [1-9], 
possess some of the properties of either classic Phase Locked 
Loop (PLL) or classic Frequency Locked Loop (FLL). 
However, due to their properties, Time Recursive PLL (TR 
PLL) and Time Recursive FLL (TR FLL) exceed the 
applicability and functionality of classic PLL and FLL, 
improving the efficiency and expanding the scope and diversity 
of their applications. Different applications include phase 
shifting in [1, 2], the development of digital frequency 
synthesizers [3], noise rejection applications [4-6] and a wide 
range of tracking and prediction [4, 5, 7]. All the algorithms 
described are very suitable for usage in a software form, like 
for example a software predictor [8]. Although a very complex 
system, consisting of a lot of subsystems, it can also function as 
a FLL, whose realization is based on the same technique [9]. 

TR PLL and TR FLL represent a completely new approach 
both in theory and in implementation and application. For 
example, the measurement of phase, frequency or amplitude, 
which are used in classic PLL and FLL, are changed by the 

measurement of time, which is much more convenient and 
more precise. Let us remember that the classic PLL tends to 
equalize both frequency and phase between the input and the 
output signals, reducing the phase difference to zero. Unlike 
PLL, classic FLL tends to equalize the frequency, without 
taking care of the phase difference. When FLL reaches the 
stable state, the phase difference between the input and output 
signals depends on the initial variable conditions and, possibly, 
it depends on the system parameters. In other words, the phase 
difference is the random variable. According to [1, 2] some of 
TR PLLs possess the ability to regulate the phase difference by 
either system parameters or by output control. This is an 
important feature of the TR PLL, which can be utilized to 
generate phase differences of any value. Note that the phase 
difference of any PLL cannot depend on the initial conditions 
of the variables. 

In this paper, a new model of a time/phase shifter, based on 
TRP, is described. It provides phase shifting in a band without 
limitations which is controlled by the clock frequency. Under 
some system conditions, it can also provide time shifting 
linearly dependent on the value of an outside control. In 
comparison with the previous phase shifters of this type, 
described in [1, 2], this phase shifter can generate the proper 
shifting even in case the input period is a ramp function.  

II. MATHEMATICAL DESCRIPTION OF THE SHIFTER 
One general case of the time relation between an input 

signal Sin and an output signal Sop of the shifter, is shown in 
Figure 1. The periods TI0, TI1,….TIk, TIk+1, and TO0, 
TO1,….TOk, TOk+1, as well as the time differences 0, 1, 
2…..k, k+1, occur at discrete times respectively t0, t1, t2,…..tk, 
tk+1. The discrete times t0, t1, t2,….tk, tk+1 are defined by the 
falling edges of the pulses of Sop in Figure 1. Note that, unlike 
the classic PLL, all variables are distributed in time in Figure 1. 
The natural recursive relation (1), between the variables, yields 
from Figure 1. It was supposed that the time difference  is 
positive if the output signal leads the input signal. 
Correspondingly, in case that the input signal leads, time 
difference  will take a negative value through the analysis. 
The main algorithm of the circuit is presented in (2), where 
"Tk", as an outside control word, represents a step function and 
"a" and "m" are the system parameters. According to (1) and 



Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1583  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

(2), the circuit has two output variables, which describe its 
behavior: 

                           1k k k kTO TI                              (1) 
                    1 1a kk k kTI T mTO                            (2) 
 

    
Fig. 1.  Time relations between all variables of the circuit. 

The output variables are (k+1) = f[TI(k)] and TO(k+1) = 
f[TI(k)]. Note that, because of simplicity, TO(k), TI(k) and (k) 
are denoted in the article as TOk, TIk and k. It is now necessary 
to find the Z transform of the output variables (k+1) and 
TO(k+1) in order to analyze the properties of the shifter 
described. The Z transform of (1) and (2) are given by 
respectively (3) and (4). Time constants TO0 and 0 are the 
initial values of the output variables TO(k) and (k), which 
appear at discrete time t0 = 0.    

               0( ) ( ) ( ) ( )z z z z TO z TI z                         (3) 
   0 0a( ) ( ) ( ) ( )zTO z zTO TI z T z z z z m                 (4) 

Calculating (z) from (3) and using it into (4), gives: 

0 0z mz zz z
( ) (z 1)z(a m) a

( ) ( )
z(z m 1) z(z m 1)

T TO
TO TI

    
 

   
   (5) 

In the same way, changing TO(z) from (5) to (3), gives: 

0 0(z) (z )a( ) ( )
(z m 1) (z m 1)

T z m TO zz
z TI z

z z



   

  
   

     (6) 

Two transfer functions, presented by (7) and (8), describing 
the output variables of the shifter in dependence on the input 
period, can be defined from respectively (5) and (6). 

                      
( ) z(a m) a

( )
( ) z(z m 1)TO

TO z
H z

TI z
 

 
 

                   (7) 

                       
( ) a

( )
( ) (z m 1)
z z

H z
TI z z
 

  
 

                   (8) 

A. Step analyzes of the circuit 
Step analysis will discover under which conditions the 

circuit can possess the properties of a time/phase shifter. 
Providing that the control time function T(k)=T=const., let us 
suppose that a step function, TI(k)=TI=const., is applied to the 
input. If we change the Z transforms T(z) = T·z/(z-1) and TI(z) 
= TI·z/(z-1) into (5) and (6), we will get the Z transform of the 
output variables for this input.  Changing T(z) and TI(z) into  
(5) and using the final value theorem, it is possible to find the 

final value of the output period in the time domain as TO = 
lim TO(k) if k, using TO(z), as shown in (9). Applying (9), 
TO is calculated as shown in (10). To complete the 
information about the system properties, it is necessary to 
determine ∞, i.e. the final value of (k). If T(z) and TI(z)  are 
changed into (6), using the final value theorem in the same way 
like for TO, ∞ can be determined using (11). Applying (11), 
∞ is calculated and shown in (12). 

                       1lim[( 1) ( )]zTO z TO z                         (9) 
                       lim[ ( )]kTO TO k TI                         (10) 
              1lim ( ) lim[( 1) ( )]k zk z z                      (11) 
                       (1 ) / /TI a m T m                              (12) 

The expressions (10) and (12) are valued only if the shifter 
is the stable system i.e. if |z1| < 1 and |z2| < 1, where z1 and z2 
are the poles of the transfer function HTO(z) or H(z), given by 
(7) and (8). Since z1 = 0, and z2 = 1+m, it yields that the shifter 
is the stable system if (13) is satisfied. According to (10) and 
(12), it follows that the described shifter possesses the 
properties of a PLL. This conclusion comes out from the facts 
that the output period TO, for the stable system, equalizes the 
input period and that the time difference  does not depend on 
the initial conditions 0 and TO0. 

                                         2 0m                                  (13) 

B. The time and/or phase shifting 
Looking at (12), we can notice that total value of  

consists of two parts. The first one is TI(1-a)/m. The term (1-
a)/m  indicates the portion of a period TI, or the number of 
periods TI, for which the input pulse rate is phase shifted. 
According to its nature, the first part represents the phase 
shifting. The system stability does not depend on the parameter 
"a", so that the value of either positive or negative phase 
shifting is unlimited. The second part is equal to "-T/m". It 
represents the pure time shifting, which is controlled by "T". It 
can also be ether positive or negative. Since control word "T" 
does not affect on the system stability, this kind of shifting can 
be also unlimited. So, two kinds of shifting can be generated 
simultaneously. But they can be also generated separately. If 
a=1, the first part in (12) does not take part in the shifting. At 
the same time, if control time function T=0, the second part in 
(12) has no effect on the shifting. 

To prove the correctness of the previous analysis, 
simulations of the shifter functioning are presented. At the 
same time, the simulations are employed to enable a further 
insight into the procedure and into the physical meaning of the 
variables described. The simulations are also used to discover 
additional properties and possible efficient applications of the 
shifter. All discrete values in simulations were merged to form 
continuous curves. Note that all variables in the following 
diagrams were presented in time units. The time unit can be, 
sec, msec or any other, but assuming the same time units for 
TI, TO,  and T. It was more suitable to use just “time unit” or 
abbreviated “t.u.” in the text. It was more convenient to omit 
the indication “t.u.” in diagrams. All simulations were made by 
Matlab, using (1) and (2). 



Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1584  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

Let us first consider the time shifting, taking a=1. The 
locking procedure of the shifter is considered for three cases in 
Figure 2. The input is a step function TI = 10 t.u. The parameter 
"m" and control word "T" are chosen carefully in order to 
demonstrate the different properties of the shifter. All values of 
"m" and "T", as well as the values of the initial conditions TO0 
and 0  are presented in Figure 2. It can be seen in Figure 2a, 
that for any value of "m" satisfying (13), the output period 
tends to the input period, but the transition time depends on the 
value of "m". For case Nr. 2, (m=-1) the shifter is extremely 
fast. It reaches the stable state for only two steps. For case Nr. 
1, (m=-0.85), and case Nr. 3 (m=-1.2), the shifter takes longer 
time to reach the stable state. The conclusion is that the greater 
deviation of parameter "m", with respect to m=-1, provides 
longer transition time of the shifter. Note that for all three 
cases, the calculation of ∞=-T/m agrees with ∞ which can be 
seen on the simulated curves. For instance for case Nr. 2, ∞=-
T/m=-(-3)/-1=-3. The same result can be seen in Figure 2a. 
This agreement confirms the correctness of the simulation 
results and the correctness of the whole approach.  

 

 
Fig. 2.  Time shifting (a=1 for all of three cases): a. The output period 
tends to the input period and the final value of time difference ∞ = -T/m,  b.  
The case Nr. 3, in Figure 2a, is presented in the real time form. All periods of 
Sin, Sop and time differences k are calculated, showing the real time relation 
between them. For this case,  T=0 (∞=0) and the shifter functions as a classic 
PLL.   

The merged discrete values for case Nr. 3, representing the 
output variables TO3 and 3 in Figure 2a, are shown in real 
time, in Figure 2b. This presentation, in forms of pulse rates, 
exhibits the real time relations between all input and output 
variables. It also helps to better understand the physical mining 
of the simulated curves. The calculated values of all variables 
are shown for every step. All values are taken from the 
simulation results, but they can be calculated manually, step by 
step, using (1) and (2). According to the mathematical analyzes 
and simulation results obtained, the shifter functions as classic 

PLL for T=0 and a=1. The negative values of the time 
differences k are shaded in Figure 2b. 

Let us now consider the properties of the circuit in case that 
it functions as the phase shifter. In this case, according to (12),  
T=0 and ∞= TI(1-a)/m. The locking procedure of the shifter is 
simulated for three cases in Figure 3. The input is the step 
function TI = 10 t.u.  All values of "m" and "a", as well as the 
values of the initial conditions TO0 and 0  are presented in 
Figure 3. It can be seen in Figure 3a, that the role of parameter 
"m" is not changed. For any value of parameter "m" satisfying 
(13), the output period tends to the input period, but the 
transition time depends on the value of the parameter "m". For 
case Nr. 2, (m=-1) the shifter reaches the stable state for only 
two steps. The conclusion is the same as in the case of the time 
shifter. The greater deviation of parameter "m", with respect to 
m=-1, provides longer transition time of the shifter. Although 
the parameter "m" affects ∞, it is much more convenient to use 
"a" for the phase shifting. For case Nr. 2, a=1 and 2∞=0, 
regardless of the values of "m". In this case phase shifter 
functions as a classic PLL. For case Nr. 1, a=1.16 (a> 1) and 
1∞=TI(1-a)/m = 2 t.u. The same result for 1∞ can be seen in 
Figure 3b. For case Nr. 3, a=0.75 (a<1) and 3∞=TI(1-a)/m = - 2 
t.u. The same result can be also seen in Figure 3b.  

 
 

 
Fig. 3.  Phase shifting (T=0 for all of three cases) in the function of 
parameters "m" and "a": a. The output period tends to the input period,  b. The 
final value of time difference ∞ = TI(1-a)/m,  c. The final value of the phase 
shifting is Ph∞ =2π·∞/TO∞ [rad] = 2π·(1-a)/m [rad]. 

This agreement confirms the correctness of the simulation 
results and the correctness of the mathematical approach. We 
can now reach the general conclusion about the influence of 
parameter "a" to the phase shifting. If a>1, the phase shifting is 
positive. If however a<1, the phase shifting is negative. Since 
the parameter "a" does not affect the stability of the system, 
both of them, the positive and the negative shifting can be 
unlimited. The only limitation is the capacity of the circuits 
which are built in the shifter. The general expression for the 
phase shifting, for both, the time and the phase shifter, can be 



Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1585  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

expressed as Ph∞ = 2π·∞/TO∞ [rad]. Changing ∞=TI(1-a)/m 
for the phase shifter, and taking in account that for the stable 
shifter TI=TO∞, the expression turns into Ph∞ = 2π·(1-a)/m 
[rad]. The transition state, as well the stable state, for all of 
three cases of the phase shifting are presented in Figure 3c.  
According to the previous expression,  Ph1∞=2π·(1-1.16)/(-
0.8)=1.25 [rad],  Ph2∞=2π·(1-1)/(-1)=0 [rad] and Ph3∞=2π·(1-
0.75)/(-1.25)=-1.25 [rad]. All of calculated results agree with 
the simulated ones, shown in Figure 3c. 

C. The shifting in case when the input is a ramp function 
Besides the described shifting for the case when the input 

period is a step function, the time shifter (a=1) can be used for 
the shifting of an input pulse rate, whose period is a ramp 
function. To illustrate this ability, let us determine the well 
known velocity error, providing that the input period is the 
ramp function TI(k)=TIV(k)=p∙k, where “p” is a time constant. 
The definition of velocity error is KV=lim[TOV(k)-TIV(k)] for 
k→∞. One more convenient expression for velocity error is KV 
= lim TIV(k)[HTO(k)-1] for k→∞, where HTO(k) is the transfer 
function given by (7). Taking a=1 and using the previous 
expression and Z transform of velocity function TI(z)=TIV(z)= 
Z(p·k)=pz/(z-1)2, it can be found out that KV = 0. This means 
that the output of the shifter is capable to track the ramp input, 
without an error.  

In order to consider the shifting abilities of the circuit, let us 
now determine the behaviour of τV(k) for the velocity (ramp) 
input, if k→∞. Using the final value theorem, τV∞=lim τV(k)k→∞ 
is calculated using τV(z) and shown in (14). The expression 
τV(z) is get out by the substitution of TI(z)=TIV(z) in (6) and 
taking a=1. According to (14), τV∞ is the time constant, which 
does not depend on the initial conditions. Besides the 
parameters "m" and "T", τV∞ depends on the time constant "p",  

               1lim[( 1) ( )] ( ) /V V zz z p T m                   (14) 

which is the slope of the ramp input function. It means the time 
shifting of the output pulses can be controlled by the control 
word "T". If T>p, τV∞>0 and the output pulse leads in time, just 
like in Figure 1. If however T<p, τV∞ is negative and the output 
pulse is delayed in time, in comparison with the corresponding 
input pulse. It comes out that if T=p, τV∞ = 0, i.e. the shifter 
functions as classic PLL, reducing the phase difference to zero, 
even for the ramp input function. To check the correctness of 
the obtained results, the simulations of TO(k) and (k) for the 
velocity input TIk=(10+4·k) t.u. are shown in Figure 4. Three 
cases for different parameters “m” and “T” are presented. All 
initial conditions and final values are shown in Figure 4. For all 
of cases in Figure 4a, the output periods track the input periods 
without an error. That means KV = 0, which agrees with the 
previous theoretical conclusion. Note that the shifter takes only 
two steps to reach the stable state for m=-1, regardless of the 
fact that the input is the ramp function. For all cases in Figure 
4b, time differences tends to the time constant "(p-T)/m", given 
by (14), proving the correctness of the mathematical analyses. 
For case Nr. 2, T=p=4 t.u. and according to (14), V2∞ = 0. That 
means the shifter predicts and generates the output period 
which is completely synchronized with the input pulses, even 
for  the case when the input is a ramp function. For case Nr. 1, 

V1∞=(p-T)/m=(4-7.75)/(-0.75)=5 t.u. For case Nr. 3, V3∞=(p-
T)/m=[4-(-2)]/(-1.25)=-4.8 t.u. The calculated results agree 
with the simulated ones, shown in Figure 4b. Cases Nr. 1 and 
Nr. 3 prove that the shifter is able to generate both positive and 
negative shifting of the pulses of a ramp input, depending on 
the parameter "T". For case Nr.1 T>p (V∞>0), the output pulses 
lead for V∞+ in comparison with the input pulses. For the case 
Nr. 3, T<p (V∞ <0), the output pulses are delayed for V∞ in 
comparison with the input pulses. 

 

 
Fig. 4.  The time shifting of the ramp input function (a=1): a. The output 
periods, for all of three cases, track the input periods without an error, Kv=0,  
b. The time shiftings tend to the time differences V∞, given by eq. (14).  

III. REALIZATION OF THE SHIFTER 
According to the results obtained, the optimum value of 

"m" for the shifting applications, is m=-1. Firstly, the shifter is 
the fastest for m=-1 and secondly the realization of the shifter 
is much simpler if m=-1. For some other applications, where it 
is necessary to adapt the locking speed, different values of "m" 
can be required. Taking m=-1, and changing a=fa/fc, where fa 
and fc are the clock frequencies,  (2) turns into (15). According 
to (15), Tk and k+1 are measured by the clock period tc=1/fc, TIk 
is measured by the clock period ta = 1/fa and TOk+1 is generated 
by the clock period tc=1/fc. The principal scheme of the 
realization, corresponding to (15) is shown in Figure 5.  

       1 1c a c k ck k kf TI f T ff TO                     (15) 
                                                       

    
Fig. 5.  The principal scheme of the phase shifter   

The shifter consists of a Recursive Calculation Module 
(RCM) and a Programmable Period Generator (PPG). RCM 
calculates TOk+1 and PPG generates the output period at the 
next step. PPG was described in detail in [1-9]. The realization 
of RCM depends on the type of algorithm, but the realization of 



Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1586  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

RCM according to (15), can be easily performed using the 
technique described in [1-9]. According to the previous 
conclusions, if fa=fc (a=1), the circuit functions as time shifter. 
If however fa is not equal to fc and Tk=0, the circuit functions as 
the phase shifter. One eight bits shifter was realized by 
standard integrated circuits for a=1 (fa=fc), T=0 and m=-1. The 
input period, shown in Figure 6, is the step function TI=1 [ms], 
corresponding to the frequency of 1[kHz]. The clock frequency 
fc=6 [kHz], so that the ratio TI/tc=6. The oscilloscope picture of 
the voltage waveforms, shown in Figure 6, was taken when the 
shifter was in the stable state. In accordance with the chosen 
parameters, the shifter functions as a classic PLL, striving to 
reduce the time difference to zero. However, that is not 
possible because the ratio TI/tc is very small. The ratio TI/tc is 
chosen intentionally small to enable the visibility of ∞+ and ∞- 
and to provide the evidence of the whole physical process. In 
accordance with the chosen parameters, the time difference ∞ 
is changing the sign at every step. If TI/tc increases, the widths 
of ∞+ and ∞- would decrease and tend to zero. It can also be 
seen in Figure 6 that whenever two pulses of Sop occur during 
a period TI, it means that TI > TO, and the coming “∞” is 
negative. If two pulses of Sin comes during a period TO, it 
means that TO > TI, and the coming “∞” is positive. The 
previous explanations, at the same time, contribute to the 
understanding of the principles of the shifter functioning. 

 

 
Fig. 6.  The shifter is in the stable state. It functions as PLL (a=1, m=1, 
T=0). The ratio TI/tc is only about six and the time differences ∞+ and ∞- are 
visible. If TI/tc increases, the time differences tend to zero. 

IV. CONCLUSION 
The description and illustration of the realized shifter 

represent one additional contribution to the shifters described in 
[1, 2]. The presented theoretical approach in whole, also 
represents the contribution to the field of TRP. It was 
demonstrated that besides the processing of the input and 
output periods and the time differences between them, an 
additional control function can take part in the calculation of 
the output period. In dependence on the type of the control 
function, the shifter takes new properties. For instance, it was 
shown that if the additional control function is a step function, 
the shifter is able to provide the proper shifting even in case the 
input period is a ramp function. This example indicates that, 
using this approach, it is possible to develop new additional 
applications of the circuits and systems based on TRP.  

In comparison with the shifters described in [1, 2], this 
shifter provide important advantages consisting of the better 
functionality and the wider applicability. Firstly, unlike the 
mentioned shifters, this time shifter provides both the time 
shifting and the phase shifting. Secondly, the time shifting does 

not depend on the input period, so that the desired value of 
either positive or negative shifting is directly proportional to 
the outside control word. This property significantly facilitates 
its practical realization and its practical usage. Thirdly, this 
shifter can generate either a positive or negative shifting of the 
input pulse signal, even in case the input period is a ramp 
function. This feature provides a much wider field of 
application.  

It was shown that if T = 0, the shifter can function exactly 
as classic PLL, which has already found a number of 
applications. The analysis and simulations showed that this 
phase shifter is extremely fast if m=-1. In this case, the shifter 
takes only two steps to reach the stable state. Changing m in 
range of -2 to 0, it is possible to adapt a transient time of the 
shifter to the required application. The precision of the phase 
shifting can be extremely high. It is limited however by the 
speed of the used integrated circuits, since it depends on 
relation TI/tc. A higher TI/tc relation provides a higher shifting 
resolution. Generally, the shifting band can be as large as 
necessary. However, the capacity of the built-in up-down 
counters must be large enough to satisfy the requirements for 
the shifting band and the shifter resolution. The identity of the 
analytical and simulation results, in every step, proves the 
correctness of the entire theoretical approach, as well as the 
validity of the obtained results. 

REFERENCES 
[1] D. Perisic, M. Perisic, S. Rankov, “Phase shifter based on a Recursive 

Phase Locked Loop of the second order”, Rev. Roum. Sci. Techn.–
Electrotechn. et Energ., Vol. 59, No. 4, pp. 391–400, 2014 

[2] D. Perisic, A. Zoric, D. Babic, D. Perisic, “Recursive PLL of the First 
Order”, Przegląd Elektrotechniczny, Vol. 89, No.7, pp. 50-53, 2013 

[3] D. Perisic, A. Zoric, S. Obradovic, D. Perisic, “FLL as digital period 
synthesizer based on Binary Rate Multiplier control”, Przegląd 
Elektrotechniczny, Vol. 89, No.  1a, pp. 145-148, 2013 

[4] D. Perisic, M. Bojovic, “Multipurpose Time Recursive PLL ”, Rev. 
Roum. Sci. Techn.–Electrotechn. et Energ., Vol. 61, No. 3, pp. 283–288,  
2016 

[5] D. Perisic, A. Zoric, M. Perisic, D. Mitic, “Analysis and Application of 
FLL based on the Processing of the Input and Output Periods”, 
Automatika, Vol. 57, pp. 230–238, 2016 

[6] D. Perisic, M. Perisic, D. Mitic, M. Vasic “Time Recursive Frequency 
Locked Loop for the tracking applications”, Rev. Roum. Sci. Techn.–
Electrotechn. et Energ., Vol. 60, No. 2, pp. 195–203, 2015 

[7] D. Perisic, A. Zoric, M. Perisic, V. Arsenovic, L. Lazic “Recursive PLL 
based on the Measurement and Processing of Time”, Elektronika ir 
Elektrotechnika, Vol. 20, No. 5, pp. 33-36, 2014 

[8] D. Perisic, A. Zoric, D. Babic, D. Perisic, “Decoding and Prediction of 
Energy state in Consumption control”, Rev. Roum. Sci. Techn.–
Electrotechn. et Energ., Vol. 58, No. 3, pp. 263–272, 2013 

[9] D. Perisic, A. Zoric, S. Obradovic, S. Spalevic, “Application of 
Frequency Locked Loop in Consumption Peak Load Control”, Przegląd 
Elektrotechniczny, Vol. 88, No. 1b, pp. 264-267, 2012 

[10] R. N. Dean, A. K. Rane, “Digital Frequency-Locked Loop System for 
Capacitance Measurement”, IEEE Transactions on Instrumentation and 
Measurement, Vol. 62, No. 4, pp. 777-784, 2013 

[11] M. Abdul-Niby, M. Alameen, H.Baitie, "A Simple Phase Shifting 
Technique for an Injection Locked Oscillator" Engineering, Technology 
& Applied Science Research, Vol. 6, No. 6, pp. 1303-1306, 2016 

[12] D.  Jovcic, “Phase locked loop system for FACTS”, IEEE Transactions 
on Power Systems, Vol. 18, No. 3, pp. 2116-2124, 2003 



Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1582-1587 1587  
  

www.etasr.com Perisic and Bojovic: Application of Time Recursive Processing for the Development of a … 
 

[13] D. Abramovich, “Phase-locked loops: a control centric tutorial”, 
American Control Conference, Vol 1, pp. 1-15, USA, November 07, 
2002 

[14]  A. S. N. Mokhtar, M. B. I. Raez, M. Marufuzzaman, M. A. M. Ali, 
“Hardware implementation of a high speed inverse park transformation 
using cordic and pll for foc brushless servo drive”, Elektronika ir 
Elektrotechnika,  Vol. 19, No. 3, pp. 23-26, 2013 

[15] M. Pastor, J. Dudric, “Grid-tied 15-level Cascade Inverter with 
Predictive Current Control”, Elektronika ir Elektrotechnika, Vol 18, No. 
9, pp. 19-22, 2012 

[16] G. Bianchi, Phase-Locked Loop Synthesizer Simulation, Nc-Hill, Inc. 
New York, USA, 2005 

[17] R. Vich, Z Transform Theory and Application (Mathematics and 
Applications), Springer, 1987 

[18] D. B. Tablot, Frequency Acquisition Techniques for Phase Locked 
Loops, Wiley-IEEE Press, 2012        

ACKNOWLEDGEMENT 
This article was supported by the Ministry of Science and 
Technology of the Republic of Serbia within the project TR 32047. 
 

AUTHORS PROFILE 

Djurdje Perisic was born in Atenica near Cacak, Serbia in 
1943. He graduated from the Faculty of Electrical Engineering 
in Belgrade in 1967. Through an English government 
scholarship, he received an M.Sc. from the Cranfield Institute 
of Technology, England in 1976. He received a Ph.D. at the 
Technical Academy in Zagreb in 1983. After that he was 
awarded twice the Fulbright and once the Humboldt 
scholarship. He was awarded the "Tesla’s award", the highest 

Yugoslav recognition for contributions to the field of natural 
sciences, from the Serbian Academy of Sciences in 1996. In 
2003, as the dean of the faculty, he founded the first study 
programs for the Faculty of Information Technology at SPU. 
In 2013 he was elected a professor emeritus. His current 
research interests include the development of new approaches 
to electronic circuits and systems combining electronics, 
system theory and other fields. During last five years he has 
published fifteen papers in SCI list journals and about 40 
additional papers of international importance. 
 
Miroslav Bojovic was born in 1957. All his formal education, 
including the Ph.D. (1989) are from the Faculty of Electrical 
Engineering, University of Belgrade, Serbia. He realized the 
specialization at the faculty of Computer Science Department, 
UCLA, Los Angeles, USA, during 1988 and 1989. As a Head 
of the Computer Science Department at the Faculty, he is 
currently engaged as a professor at the Faculty of Electrical 
Engineering in Belgrade in the field of Database menagement 
System and Software Engineering and as a project leader of 
many projects in Software engineering. One of his projects 
"MobilePDR" was awarded as the best software product in 
USA in the field of medical information systems in 2004. As a 
software product, it has been referenced in more than 100,000 
different publications, all over the world.