@@ -160,6 +160,7 @@ def _default_gains(num, den, xlim, ylim):
160
160
mymat_xl = mymat
161
161
singular_points = np .concatenate ((num .roots , den .roots ), axis = 0 )
162
162
important_points = np .concatenate ((singular_points , real_break ), axis = 0 )
163
+ important_points = np .concatenate ((singular_points , np .zeros (2 )), axis = 0 )
163
164
mymat_xl = np .append (mymat_xl , important_points )
164
165
165
166
if xlim is None :
@@ -323,20 +324,14 @@ def sgrid_func(fig):
323
324
ylocator = ax .get_yaxis ().get_major_locator ()
324
325
xlocator = ax .get_yaxis ().get_major_locator ()
325
326
326
- long_xaxis = xlocator ()[- 1 ] - xlocator ()[0 ]
327
- long_yaxis = ylocator ()[- 1 ] - ylocator ()[0 ]
328
-
329
327
angules = np .arange (- 90 , 80 , 15 )* np .pi / 180
330
328
331
- # radial lines
329
+ # zeta-constant lines
332
330
y_over_x = np .tan (angules [1 ::])* ylocator ()[- 1 ]/ xlocator ()[- 1 ]
333
-
334
331
ylim = ax .get_ylim ()
335
332
ytext_pos_lim = ylim [1 ]- (ylim [1 ]- ylim [0 ])* 0.03
336
-
337
333
xlim = ax .get_xlim ()
338
334
xtext_pos_lim = xlim [0 ]+ (xlim [1 ]- xlim [0 ])* 0.0
339
-
340
335
index = 0
341
336
zeta = np .sin (np .pi / 2 - angules [1 ::])
342
337
@@ -363,6 +358,7 @@ def sgrid_func(fig):
363
358
xtext_pos = - xtext_pos
364
359
ax .annotate (an , textcoords = 'data' , xy = [xtext_pos , ytext_pos ], fontsize = 8 )
365
360
index += 1
361
+ ax .plot ([0 , 0 ], [ylim [0 ], ylim [1 ]], color = 'gray' , linestyle = 'dashed' , linewidth = 0.5 )
366
362
367
363
angules = np .linspace (- 90 , 90 , 20 )* np .pi / 180
368
364
for xt in xlocator ():
0 commit comments