@@ -424,7 +424,6 @@ def gen_zero_centered_series(val_min, val_max, period):
424
424
# Nyquist plot
425
425
#
426
426
427
-
428
427
def nyquist_plot (syslist , omega = None , Plot = True ,
429
428
labelFreq = 0 , hl = 0.1 , hw = 0.1 , color = None , * args , ** kwargs ):
430
429
"""
@@ -499,11 +498,11 @@ def nyquist_plot(syslist, omega=None, Plot=True,
499
498
c = p [0 ].get_color ()
500
499
ax = plt .gca ()
501
500
# Plot arrow to indicate Nyquist encirclement orientation
502
- ax .arrow (x [0 ], y [0 ], ((x [1 ]- x [0 ])/ 2 ), ((y [1 ]- y [0 ])/ 2 ), fc = c , ec = c ,
503
- head_width = hw , head_length = hl , color = color )
501
+ ax .arrow (x [0 ], y [0 ], ((x [1 ]- x [0 ])/ 2 ), ((y [1 ]- y [0 ])/ 2 ),
502
+ fc = c , ec = c , head_width = hw , head_length = hl ,
503
+ color = color )
504
504
505
- plt .plot (x , - y , '-' , color = color , * args , ** kwargs )
506
- # plt.figure()
505
+ plt .plot (x , - y , '-' , color = c , * args , ** kwargs )
507
506
ax .arrow (
508
507
x [- 1 ], - y [- 1 ], ((x [- 1 ]- x [- 2 ])/ 2 ), ((y [- 1 ]- y [- 2 ])/ 2 ),
509
508
fc = c , ec = c , head_width = hw , head_length = hl , color = color )
@@ -548,8 +547,6 @@ def nyquist_plot(syslist, omega=None, Plot=True,
548
547
#
549
548
550
549
# TODO: think about how (and whether) to handle lists of systems
551
-
552
-
553
550
def gangof4_plot (P , C , omega = None , ** kwargs ):
554
551
"""Plot the "Gang of 4" transfer functions for a system
555
552
Generates a 2x2 plot showing the "Gang of 4" sensitivity functions
@@ -650,8 +647,6 @@ def gangof4_plot(P, C, omega=None, **kwargs):
650
647
#
651
648
652
649
# Compute reasonable defaults for axes
653
-
654
-
655
650
def default_frequency_range (syslist , Hz = None , number_of_samples = None ,
656
651
feature_periphery_decades = None ):
657
652
"""Compute a reasonable default frequency range for frequency
@@ -772,7 +767,6 @@ def default_frequency_range(syslist, Hz=None, number_of_samples=None,
772
767
# KLD 5/23/11: Two functions to create nice looking labels
773
768
#
774
769
775
-
776
770
def get_pow1000 (num ):
777
771
"""Determine exponent for which significand of a number is within the
778
772
range [1, 1000).
@@ -824,4 +818,4 @@ def find_nearest_omega(omega_list, omega):
824
818
# Function aliases
825
819
bode = bode_plot
826
820
nyquist = nyquist_plot
827
- gangof4 = gangof4_plot
821
+ gangof4 = gangof4_plot
0 commit comments