René's URL Explorer Experiment

Title: Better I/O systems support for phase plots by murrayrm · Pull Request #1001 · python-control/python-control · GitHub

Open Graph Title: Better I/O systems support for phase plots by murrayrm · Pull Request #1001 · python-control/python-control

X Title: Better I/O systems support for phase plots by murrayrm · Pull Request #1001 · python-control/python-control

Description: This PR provides better I/O systems support for making phase plots: Allow I/O systems with no inputs and/or no outputs to be specified. This is useful when generating phase plots, where you might just have a simple dynamical systems but still want to use the I/O system objects. Improve the error messages in input_output_response when the number of states, inputs, or outputs are incompatible with the system size by telling you which one didn't match. Generate warnings when simulations fail for individual initial conditions and drop individual trace (rather than terminating). This yields much nicer phase plane plots by allowing the system to generate lines when it can and ignoring the rest. Added Jupyter notebook analyzing open and closed loop dynamics of an inverted pendulum, showing off some of the capabilities of the package. Changed the way plot titles are created, using axes.set_title (centers title over axes) instead of fig.suptitle (centers over figure, which is good for multi-axes plots but otherwise looks funny). I'll propagate this same change to time plots and frequency plots in separate PRs. Did some tweaking on arrow placement so that very short lines have zero or one arrows. Could use more work, but this is better than what we had. As an example of the changes, here is the code to produce a phase portrait for the inverted pendulum example in FBS2e, Figure 5.4: from math import pi import control as ct def invpend_update(t, x, u, params): m, l, b, g = params['m'], params['l'], params['b'], params['g'] return [x[1], -b/m * x[1] + (g * l / m) * np.sin(x[0])] invpend = ct.nlsys( invpend_update, states=2, inputs=0, name='invpend', params={'m': 1, 'l': 1, 'b': 0.2, 'g': 1}) ct.phase_plane_plot(invpend, [-2*pi - 1, 2*pi + 1, -2, 2], 10) Before: After: The changes are not major, but if you look at the number and spacing of arrows along the streamlines, you'll see examples of the changes there (some streamlines have no arrows in "before" plot). Also note the change in location of the title.

Open Graph Description: This PR provides better I/O systems support for making phase plots: Allow I/O systems with no inputs and/or no outputs to be specified. This is useful when generating phase plots, where you might...

X Description: This PR provides better I/O systems support for making phase plots: Allow I/O systems with no inputs and/or no outputs to be specified. This is useful when generating phase plots, where you might...

Opengraph URL: https://github.com/python-control/python-control/pull/1001

X: @github

direct link

Domain: github.com

Links:

Viewport: width=device-width