Commit72cb23f

committed

added warning method when falling back to frd

1 parent567b4db commit72cb23fCopy full SHA for 72cb23f

File tree

-2

lines changed

-2

lines changed

Lines changed: 5 additions & 2 deletions

Original file line number	Diff line number	Diff line change
`@@ -51,6 +51,7 @@`
`51`	`51`	`"""`
`52`	`52`
`53`	`53`	`importmath`
	`54`	`+fromwarningsimportwarn`
`54`	`55`	`importnumpyasnp`
`55`	`56`	`importscipyassp`
`56`	`57`	`from .importxferfcn`
`@@ -207,7 +208,7 @@ def fun(wdt):`
`207`	`208`
`208`	`209`	`returnz,w`
`209`	`210`
`210`		`-def_numerical_inaccuracy(sys):`
	`211`	`+def_likely_numerical_inaccuracy(sys):`
`211`	`212`	`# crude, conservative check for if`
`212`	`213`	`# num(z)num(1/z) << den(z)den(1/z) for DT systems`
`213`	`214`	`num,den,num_inv_zp,den_inv_zq,p_q,dt=_poly_z_invz(sys)`
`@@ -334,7 +335,9 @@ def stability_margins(sysdata, returnall=False, epsw=0.0, method='best'):`
`334`	`335`	`elifmethod=='best':`
`335`	`336`	`# convert to FRD if anticipated numerical issues`
`336`	`337`	`ifisinstance(sys,xferfcn.TransferFunction)andnotsys.isctime():`
`337`		`-if_numerical_inaccuracy(sys):`
	`338`	`+if_likely_numerical_inaccuracy(sys):`
	`339`	`+warn("stability_margins: Falling back to 'frd' method "`
	`340`	`+"because of chance of numerical inaccuracy in 'poly' method.")`
`338`	`341`	`omega_sys=freqplot._default_frequency_range(sys)`
`339`	`342`	`omega_sys=omega_sys[omega_sys<np.pi/sys.dt]`
`340`	`343`	`sys=frdata.FRD(sys,omega_sys,smooth=True)`

Comments

(0)