How to unwrap the Phase plot of BodePlot?

Questions
, ,
1

Hello,

a=s^2+1/L1+1/(L2*C1)
b=s^3+1/(L1*L2)*s^2+1/(L1*L2*C1)*s
a/b//Together
a/b/.{L1->1,L2->1,C1->1}
test=TransferFunctionModel[a/b/.{L1->1,L2->1,C1->1},s]
BodePlot[test,GridLines->Automatic,PlotRange->{Automatic,{-270,0}}]
GainPhaseMargins[test]

1906Capture.PNG
1906Capture.PNG531Ă—665 35 KB

So, I want to unwrap the phase diagram of the Bode Plot. The Phase diagram of Bode Plot should look something like that after unwrapping:

9342Capture.PNG
9342Capture.PNG558Ă—348 14.7 KB

I tried this (from https://reference.wolfram.com/applications/insydes/ReferenceManual/BodePlot.html) but does not function:

test=TransferFunctionModel[a/b/.{L1->1,L2->1,C1->1},s]
BodePlot[test,GridLines->Automatic,PlotRange->{Automatic,{-270,0}}, UnwrapPhase -> True]
GainPhaseMargins[test]

9478Capture.PNG
9478Capture.PNG678Ă—252 12.9 KB

Thank you.

2

Which Mathematica are you using. The documentation you’re showing is a very old package.

Check this doc:

Specifically the PhaseRange option.

3

The documentation you’re showing is a very old package.

Yes, it looks like older documentation, but that’s all I found when I searched for UnwrapPhase for Bode Plot…

I’ll also look at PhaseRange (I tried now quickly by putting PhaseRange → {-Pi, Pi} (but also other combinations), but it didn’t work for me so far…I didn’t manage to find the right combination until now.
If anyone has found a solution then it is welcome.

4

I might be wrong, but your jump is 180 degrees, so it is a fundamentally different jump than a phase-wrap which is 360 degrees.

5

Other sugestions?
So, i want to be able to change the Phase Plot of Bode Plot in following manner from (0, 360) degrees, (90, 450) degrees, (180, 540) degrees, and so on, i think that it is the idea (or at least something like that)

6

Cornel,

I have no experience with BodePlot[] whatsoever and so I do not know whether this unwrapping makes sense in this case. But anyway - here is a simple approach:

bpl = BodePlot[test, GridLines -> Automatic, PlotRange -> {Automatic, {-270, 0}}, PlotLayout -> "Phase"];
(* extracting point values from the plot: *)
pts0 = First@Cases[bpl, Line[pts__] :> pts, All];
(* "true" x-value: *)
pts = {10^#1, #2} & @@@ pts0;
{{xmin, xmax}, {ymin, ymax}} = MinMax /@ ({xVals, yVals} = Transpose[pts]);
xGrid = Flatten[Range[10] # & /@ PowerRange[xmin, xmax]];

For the unwrapping there is an elegant algorithm from @Peter Karpov in this discussion :

unWrap[\[Tau]_][a_] := a - \[Tau] Prepend[Accumulate[Round[Differences[a]/\[Tau]]], 0]

So - we get:

(* unwrapped result: *)
ptsUnWrapped = Transpose[{xVals, unWrap[ymax - ymin][yVals]}];
ListLinePlot[ptsUnWrapped, ScalingFunctions -> {"Log", None}, GridLines -> {xGrid, Automatic}, Frame -> True]

wrappedbode.png
wrappedbode.png561Ă—351 10.1 KB

Does that help? Regards – Henrik

7

Hi Henrik,

Thank you for your response.

8

I am looking also for other solutions If exists, more simpler more better.

9

It is wired because for what I saw until now it is not possible to do unwrapping of phase of bode plot inside the BodePlot[…] command in a more simple manner in Mathematica, and only way is to extract the points from the PhasePlot of BodePlot and then doing something like Henrik did, which is something not so straight forward, mainly if someone is a beginner. And we are talking here that this unwrapping of phase plot of bode plot is something that someone who deals with control systems and signal processing will encounter very often, therefore unwrapping of phase of bode plot I think it should have already taken into account and have been implemented by the developers of wolfram mathematica…