Commit2c25b73

committed

Added a unit test for residualization of unstable systems

1 parentbf6c446 commit2c25b73Copy full SHA for 2c25b73

File tree

+12

-0

lines changed

+12

-0

lines changed

Lines changed: 12 additions & 0 deletions

Original file line number	Diff line number	Diff line change
`@@ -50,6 +50,18 @@ def testModredMatchDC(self):`
`50`	`50`	`np.testing.assert_array_almost_equal(rsys.C,Crtrue,decimal=3)`
`51`	`51`	`np.testing.assert_array_almost_equal(rsys.D,Drtrue,decimal=2)`
`52`	`52`
	`53`	`+deftestModredUnstable(self):`
	`54`	`+# Check if an error is thrown when an unstable system is given`
	`55`	`+A=np.matrix('4.5418, 3.3999, 5.0342, 4.3808;\`
	`56`	`+ 0.3890, 0.3599, 0.4195, 0.1760;\`
	`57`	`+ -4.2117, -3.2395, -4.6760, -4.2180;\`
	`58`	`+ 0.0052, 0.0429, 0.0155, 0.2743')`
	`59`	`+B=np.matrix('1.0, 1.0; 2.0, 2.0; 3.0, 3.0; 4.0, 4.0')`
	`60`	`+C=np.matrix('1.0, 2.0, 3.0, 4.0; 1.0, 2.0, 3.0, 4.0')`
	`61`	`+D=np.matrix('0.0, 0.0; 0.0, 0.0')`
	`62`	`+sys=ss(A,B,C,D)`
	`63`	`+np.testing.assert_raises(ValueError,modred,sys, [2,3])`
	`64`	`+`
`53`	`65`	`deftestModredTruncate(self):`
`54`	`66`	`#balanced realization computed in matlab for the transfer function:`
`55`	`67`	`# num = [1 11 45 32], den = [1 15 60 200 60]`

Comments

(0)