This PR adds a number of improvements for frequency plots (Bode, Nyquist, etc). The changes are motivated by my experience in using python-control for a course, where I found a few things that were harder to do than they should have been.

Changes in this PR:

• Allowlabel keyword in frequency response commands to override default label generation.

• Fix bug in processingindent_radius keyword whennyquist_plot is passed a system (keyword was not getting sent through tonyquist_response).

• Restore functionality that allows omega to be specified as a list of 2 elements (indicating a range) in all frequency response/plotting routines. This used to work fornyquist but got removed at some point. It now works for all frequency response commands.

  • Note: this could break some code if someone meant to evaluate the frequency response at exactly two frequencies. This case seems rare, but it can be fixed by passing the two elements as an array.

• Fix up theax keyword processing to allow arrays or lists + uniform processing in all freqplot routines.

• Fix uprcParams processing (uniformity + remove unnecessary code). (This needs more work in a future pass to eliminate the rcParams keywords that are part of frequency plotting functions and to pull separate time and frequency plot parameters into a single place. There are some TODOs in the code marking what needs to be done.)

• Added new suptitle() function to add titles that are better centered (on axes instead of figure). Using this for frequency domain plots for now; will update for time domain plots as well.

• Set upfrd as factory function with keywords, including setting the signal/system names. Calls toct.FRD should how be replaces withct.frd, though the former will still work.

• Allow Bode and Nyquist plots for FRD systems with different omega vectors (Fix nyquist plotting from FrequencyResponseData #996) as well as mixtures of FRD and other LTI systems.

• Added a notebookexamples/cds110-bode_nyquist.ipynb that shows a lot of the new functionality in the context of frequency domain analysis.

• Minor: changed the way thatfreq_label (frequency axis label) is implemented so that it can be more easily overriden (use fstring with '{units}' instead of '%' => can put in custom units).

For the most part these changes just make it easier to do things that were already possible. Example of some new code that makes use of these features (fromexamples/cds110-bode_nyquist.ipynb):

import control as ctP = ct.tf([1], [1, 0.1, -1])  # normalized inverted pendulumC = ct.tf([2, 10], [1])       # PD controllerL = P * C                     # Loop transfer function# Bode and Nyquist plotsplt.figure(figsize=[7, 4])ax1 = plt.subplot(2, 2, 1)plt.title("Bode plot for L", size='medium')ax2 = plt.subplot(2, 2, 3)ct.bode_plot(    [P, C, L], label=['P', 'C', 'L'], ax=[ax1, ax2],    initial_phase=0)ax3 = plt.subplot(1, 2, 2)ct.nyquist_plot(L, ax=ax3)plt.title("Nyquist plot for L", size='medium')ct.suptitle("Loop analysis for inverted pendulum")plt.tight_layout()

For the mixed FRD functionality (motivated by comments I made in#996), here's an example showing how Nyquist plots can be combined in ways that weren't possible before:

s = ct.tf('s')sys = ct.tf(    (0.02 * s**3 - 0.1 * s) / (s**4 + s**3 + s**2 + 0.25 * s + 0.04),    name='tf')sys1 = ct.frd(sys, np.logspace(-1, 1, 15), name='frd1')sys2 = ct.frd(sys, np.logspace(-2, 2, 20), name='frd2')ct.nyquist_plot([sys, sys1, sys2])ct.suptitle("Mixed FRD, tf data")