Calculating gain of a phase detector

  • Thread starter ollie456
  • Start date
  • #1
9
0
I have a linear model of a phase detector (phase lock loop) represented by the block diagram below. It is designed to maintain zero difference in phase between the input carrier signal and a local voltage controlled oscillator.

Phase demand→Ka→F(s)→K/s→Phase output

F(s) = 10(s+10)/(s+1)(s+100)

We want to minimise the steady state error for a ramp change in the phase information signal.

Can anyone help determine the limiting value of the gain KaK = Kv in order to maintain a stable system.
 

Answers and Replies

  • #2
818
67
We usually work with the open-loop function of the system when determining its margins for closed-loop stability. I assume you're showing the closed-loop function since there's no feedback loop included. Maybe you could show a block diagram of your system before any reductions are made and we can go from there?

The system you have shown has infinite gain margin, be it in open or closed loop.
 
  • #3
9
0
I've attached a more accurate block diagram which has the feedback look included.
 

Attachments

  • Phase Detector Block Diagram.png
    Phase Detector Block Diagram.png
    4.6 KB · Views: 542
  • #4
818
67
Right, are you familiar with root-locus and/or frequency-response analysis?
 
  • #5
9
0
I have been trying to use the root locus approach, i'm not sure about frequency response. I worked the closed loop transfer function to 100K2/(s+1)(s+100)+K and then used the quadratic formula to try and obtain roots getting -50.5≠.5√9801-4k. The steady state error acceptable is 1° for a ramp signal of 100rad/s. I am confused at how to go about calculating gain for a ramp change compared with a step change.
 
  • #6
818
67
First thing, I think you took a wrong turn somewhere in calculating the closed-loop transfer function. You could show the details of your work if you want some pointers, but you're going to end up having to find the roots of a cubic polynomial. That's not an intractable problem per se, but root-locus analysis is, in large part, about finding the locus of the roots of the characteristic polynomial without having to actually solve it analytically (which can be hard to impossible).

It's difficult to see on this plot since I have no idea how to set the gain and axis limits manually, but hopefully it will give you the basic idea. The locus never enters the right-half plane, i.e. your system is stable for any gain.

I am confused at how to go about calculating gain for a ramp change compared with a step change.
It's not much different. Your system is type 1 since you have a free integrator so it will have a bounded [strike]output[/strike] error for a ramp input. Do you know what the error constant is for a ramp input?
 
Last edited:
  • #7
9
0
I know the error constant is the velocity error constant Kv, i now have the transfer function as k/(100s)(s+1)(s+100)
 
  • #8
818
67
I know the error constant is the velocity error constant Kv ...
Great :approve: You need to use the velocity error constant instead of the position error constant when your input is a ramp signal but otherwise the procedure is the same. Show your work and where you get stuck.

... i now have the transfer function as k/(100s)(s+1)(s+100)
That's still not correct. Reminder:

The open-loop function L(s) is given by:
L(s) = Ka*K*F(s)/s

The closed-loop function T(s) is given by:
T(s) = L(s)/(1 + L(s))

It would help a great deal if you could show each step of what you did in order to arrive at your result.
 
  • #9
9
0
Ive attached some workings of the open loop transfer function calculation, is this evaluated at s=0 to obtain Kv
 

Attachments

  • 2013-08-13 19.01.07 HDR.jpg
    2013-08-13 19.01.07 HDR.jpg
    21.4 KB · Views: 486
  • #10
818
67
Ok so, third line, there's a problem with what you're doing in the denominator:
[tex]
\frac{1}{10}(s + 1)(s + 100) \neq \left(\frac{s + 1}{10}\right) \left(\frac{s + 100}{10}\right)
[/tex]

You should be able to tell from:
[tex]
\left(\frac{s + 1}{10}\right) \left(\frac{s + 100}{10}\right) = \frac{1}{10}(s + 1)\left(\frac{s + 100}{10}\right) = \frac{1}{10}\frac{1}{10}(s + 1)(s + 100) = \frac{1}{100}(s + 1)(s + 100)
[/tex]

Also:
[tex]K_v = \lim_{s \rightarrow 0}sL(s)
[/tex]
 
  • #11
9
0
Thanks, i've attached my working to get a value of 0.1 for Kv is this correct?
 

Attachments

  • 2013-08-14 13.28.01.jpg
    2013-08-14 13.28.01.jpg
    21.7 KB · Views: 488
  • #12
818
67
The L(s) you started out with in your first attachment was correct. You've changed it around some in your second one.

So you have:
[tex]
L(s) = K_a K \frac{10(s + 10)}{s(s+1)(s+100)} \\
K_v = \lim_{s \rightarrow 0}s L(s) = \lim_{s \rightarrow 0}s K_a K \frac{10 (s + 10)}{s(s+1)(s+100)} = \lim_{s \rightarrow 0}K_a K \frac{10 (s + 10)}{(s+1)(s+100)} = K_a K \frac{10 \cdot 10}{1 \cdot 100} = K_a K
[/tex]
The steady-state error for your system is then given by:
[tex]
e_{ss} = \frac{\alpha}{K_v} = \frac{\alpha}{K_a K}
[/tex]
where ##\alpha## is the slope of your ramp signal.

Now you just need to come up with a value for your gain ##K_a K## such that your constraint on the steady-state error is satisfied.
 
  • #13
9
0
So this is the value of the gain Kv that satisfies a steady state error of 1° with a ramp signal (slope) of 100 rads? Is there a limiting value of gain Kv to keep the system stable, or is this it?
 

Attachments

  • 2013-08-14 17.12.23.jpg
    2013-08-14 17.12.23.jpg
    24.2 KB · Views: 487
  • #14
818
67
So this is the value of the gain Kv that satisfies a steady state error of 1° with a ramp signal (slope) of 100 rads?
Yes :smile: Usually you write it as an upper bound on the steady-state error, which gives a lower bound on the gain:
[tex]
e_{ss} = \frac{\alpha}{K_a K} \leq 1\frac{\pi}{180} \Leftrightarrow K_a K \geq \frac{180\alpha}{\pi}
[/tex]
Is there a limiting value of gain Kv to keep the system stable, or is this it?
As we talked about before, your system is stable for any gain.
 
  • #15
9
0
How would you go about determining the location of the poles for the system?
 
  • #16
818
67
How would you go about determining the location of the poles for the system?
I would use a numerical algorithm to find the roots of the characteristic polynomial of the closed-loop system:
http://en.wikipedia.org/wiki/Root-finding_algorithm

Which is basically what your computer does with varying loop gain when you ask it to plot the root locus for the system.

The open-loop function usually tells you all you need to know about the performance of the closed-loop system. If its characteristic polynomial admits an analytical solution, it rarely contributes significantly to your body of knowledge of the system.

It just occured to me that this might just be a homework problem. Just FYI, if that's the case, there's a whole other forum for it on PF, which you should have used instead.
 

Related Threads on Calculating gain of a phase detector

  • Last Post
Replies
22
Views
1K
  • Last Post
Replies
2
Views
18K
  • Last Post
Replies
9
Views
2K
  • Last Post
Replies
6
Views
5K
Replies
2
Views
4K
Replies
3
Views
6K
  • Last Post
Replies
5
Views
8K
  • Last Post
Replies
5
Views
6K
Replies
10
Views
6K
Replies
5
Views
3K
Top