def testbdalg(self): # get two SISO transfer functions h1 = TransferFunction([1], [1, 2, 2]) h2 = TransferFunction([1], [0.1, 1]) omega = np.logspace(-1, 2, 10) f1 = FRD(h1, omega) f2 = FRD(h2, omega) np.testing.assert_array_almost_equal( (bdalg.series(f1, f2)).frequency_response([0.1, 1.0, 10])[0], (bdalg.series(h1, h2)).frequency_response([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.parallel(f1, f2)).frequency_response([0.1, 1.0, 10])[0], (bdalg.parallel(h1, h2)).frequency_response([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.feedback(f1, f2)).frequency_response([0.1, 1.0, 10])[0], (bdalg.feedback(h1, h2)).frequency_response([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.negate(f1)).frequency_response([0.1, 1.0, 10])[0], (bdalg.negate(h1)).frequency_response([0.1, 1.0, 10])[0])
def testbdalg(self): # get two SISO transfer functions h1 = TransferFunction([1], [1, 2, 2]) h2 = TransferFunction([1], [0.1, 1]) omega = np.logspace(-1, 2, 10) f1 = FRD(h1, omega) f2 = FRD(h2, omega) np.testing.assert_array_almost_equal( (bdalg.series(f1, f2)).freqresp([0.1, 1.0, 10])[0], (bdalg.series(h1, h2)).freqresp([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.parallel(f1, f2)).freqresp([0.1, 1.0, 10])[0], (bdalg.parallel(h1, h2)).freqresp([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.feedback(f1, f2)).freqresp([0.1, 1.0, 10])[0], (bdalg.feedback(h1, h2)).freqresp([0.1, 1.0, 10])[0]) np.testing.assert_array_almost_equal( (bdalg.negate(f1)).freqresp([0.1, 1.0, 10])[0], (bdalg.negate(h1)).freqresp([0.1, 1.0, 10])[0])