Commitdaf2536

authored

Merge pull request#447 from sawyerbfuller/fix-dt-impulse-size

fix impulse size for discrete-time impulse response

2 parents37d1121 +0b85e3b commitdaf2536Copy full SHA for daf2536

File tree

+10

-1

lines changed

+10

-1

lines changed

Lines changed: 9 additions & 0 deletions

Original file line number	Diff line number	Diff line change
`@@ -218,6 +218,15 @@ def test_impulse_response(self):`
`218`	`218`	`np.testing.assert_array_almost_equal(`
`219`	`219`	`yy,np.vstack((youttrue,np.zeros_like(youttrue))),decimal=4)`
`220`	`220`
	`221`	`+deftest_discrete_time_impulse(self):`
	`222`	`+# discrete time impulse sampled version should match cont time`
	`223`	`+dt=0.1`
	`224`	`+t=np.arange(0,3,dt)`
	`225`	`+sys=self.siso_tf1`
	`226`	`+sysdt=sys.sample(dt,'impulse')`
	`227`	`+np.testing.assert_array_almost_equal(impulse_response(sys,t)[1],`
	`228`	`+impulse_response(sysdt,t)[1])`
	`229`	`+`
`221`	`230`	`deftest_initial_response(self):`
`222`	`231`	`# Test SISO system`
`223`	`232`	`sys=self.siso_ss1`

Lines changed: 1 addition & 1 deletion

Original file line number	Diff line number	Diff line change
`@@ -822,7 +822,7 @@ def impulse_response(sys, T=None, X0=0., input=0, output=None, T_num=None,`
`822`	`822`	`new_X0=B+X0`
`823`	`823`	`else:`
`824`	`824`	`new_X0=X0`
`825`		`-U[0]=1.`
	`825`	`+U[0]=1./sys.dt# unit area impulse`
`826`	`826`
`827`	`827`	`T,yout,_xout=forced_response(sys,T,U,new_X0,transpose=transpose,`
`828`	`828`	`squeeze=squeeze)`

Comments

(0)