This looks ugly, and I want to avoid the jumps in the phase (angle) when the phase passes through -180^\circ. I.e., instead of jumping to +180^\circ, I want it to continue to a negative value past -180^\circ.
Question: Is there a simple/efficient/smart way to do this?
Note: This is a somewhat complex transfer function from a system of PDEs, and does not fit into the ControlSystems framework. [Unless ControlSystems allows me to pass on frequency and magnitude and phase instead of the system, and then does it form me.]
Not quite clear to me what you mean. I compute some complex numbers as a function of omega. Then I compute the magnitude and phase angle of the complex numbers, and plot these as a function of omega (Bode plot).
Another common plot is the Nichols plot, where one plots magnitude as a function of phase angle with omega as parameter.
A third common plot is a Nyquist plot, where the imaginary value is plotted as a function of the real part, again with omega as parameter.
Do you prefer to a Nichols plot, a Nyquist plot, or something else?