def test_lyap(self):
A = matrix([[-1, 1], [-1, 0]])
Q = matrix([[1, 0], [0, 1]])
X = lyap(A, Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * A.T + Q, zeros((2, 2)))
A = matrix([[1, 2], [-3, -4]])
Q = matrix([[3, 1], [1, 1]])
X = lyap(A, Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * A.T + Q, zeros((2, 2)))
def test_lyap(self):
A = matrix([[-1, 1],[-1, 0]])
Q = matrix([[1,0],[0,1]])
X = lyap(A,Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * A.T + Q, zeros((2,2)))
A = matrix([[1, 2],[-3, -4]])
Q = matrix([[3, 1],[1, 1]])
X = lyap(A,Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * A.T + Q, zeros((2,2)))
def test_lyap_g(self):
A = array([[-1, 2], [-3, -4]])
Q = array([[3, 1], [1, 1]])
E = array([[1, 2], [2, 1]])
X = lyap(A, Q, None, E)
# print("The solution obtained is ", X)
assert_array_almost_equal(A @ X @ E.T + E @ X @ A.T + Q,
zeros((2,2)))
# Make sure that trying to solve with SciPy generates an error
with pytest.raises(ControlArgument, match="'scipy' not valid"):
X = lyap(A, Q, None, E, method='scipy')
def test_lyap(self):
A = array([[-1, 1], [-1, 0]])
Q = array([[1, 0], [0, 1]])
X = lyap(A, Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A.dot(X) + X.dot(A.T) + Q, zeros((2, 2)))
A = array([[1, 2], [-3, -4]])
Q = array([[3, 1], [1, 1]])
X = lyap(A, Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A.dot(X) + X.dot(A.T) + Q, zeros((2, 2)))
def test_lyap(self):
A = array([[-1, 1],[-1, 0]])
Q = array([[1,0],[0,1]])
X = lyap(A,Q)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,X)+dot(X,A.T)+Q,zeros((2,2)))
A = array([[1, 2],[-3, -4]])
Q = array([[3, 1],[1, 1]])
X = lyap(A,Q)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,X)+dot(X,A.T)+Q,zeros((2,2)))
def test_lyap_sylvester(self):
A = 5
B = matrix([[4, 3], [4, 3]])
C = matrix([2, 1])
X = lyap(A, B, C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * B + C, zeros((1, 2)))
A = matrix([[2, 1], [1, 2]])
B = matrix([[1, 2], [0.5, 0.1]])
C = matrix([[1, 0], [0, 1]])
X = lyap(A, B, C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * B + C, zeros((2, 2)))
def test_lyap_sylvester(self):
A = 5
B = matrix([[4, 3], [4, 3]])
C = matrix([2, 1])
X = lyap(A,B,C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * B + C, zeros((1,2)))
A = matrix([[2,1],[1,2]])
B = matrix([[1,2],[0.5,0.1]])
C = matrix([[1,0],[0,1]])
X = lyap(A,B,C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X * B + C, zeros((2,2)))
def test_lyap_sylvester(self):
A = 5
B = array([[4, 3], [4, 3]])
C = array([2, 1])
X = lyap(A, B, C)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A, X) + dot(X, B) + C, zeros((1, 2)))
A = array([[2, 1], [1, 2]])
B = array([[1, 2], [0.5, 0.1]])
C = array([[1, 0], [0, 1]])
X = lyap(A, B, C)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A, X) + dot(X, B) + C, zeros((2, 2)))
def test_lyap_sylvester(self):
A = 5
B = array([[4, 3], [4, 3]])
C = array([2, 1])
X = lyap(A,B,C)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,X)+dot(X,B)+C,zeros((1,2)))
A = array([[2,1],[1,2]])
B = array([[1,2],[0.5,0.1]])
C = array([[1,0],[0,1]])
X = lyap(A,B,C)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,X)+dot(X,B)+C,zeros((2,2)))
def test_lyap_g(self):
A = matrix([[-1, 2], [-3, -4]])
Q = matrix([[3, 1], [1, 1]])
E = matrix([[1, 2], [2, 1]])
X = lyap(A, Q, None, E)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X * E.T + E * X * A.T + Q, zeros((2, 2)))
def test_lyap_g(self):
A = matrix([[-1, 2],[-3, -4]])
Q = matrix([[3, 1],[1, 1]])
E = matrix([[1,2],[2,1]])
X = lyap(A,Q,None,E)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X * E.T + E * X * A.T + Q, zeros((2,2)))
def test_lyap_g(self):
A = array([[-1, 2], [-3, -4]])
Q = array([[3, 1], [1, 1]])
E = array([[1, 2], [2, 1]])
X = lyap(A, Q, None, E)
# print("The solution obtained is ", X)
assert_array_almost_equal(
A.dot(X).dot(E.T) + E.dot(X).dot(A.T) + Q, zeros((2, 2)))
def test_lyap_g(self):
A = array([[-1, 2], [-3, -4]])
Q = array([[3, 1], [1, 1]])
E = array([[1, 2], [2, 1]])
X = lyap(A, Q, None, E)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,dot(X,E.T)) + dot(E,dot(X,A.T)) + Q, \
zeros((2,2)))
def test_lyap_g(self):
A = array([[-1, 2],[-3, -4]])
Q = array([[3, 1],[1, 1]])
E = array([[1,2],[2,1]])
X = lyap(A,Q,None,E)
# print "The solution obtained is ", X
assert_array_almost_equal(dot(A,dot(X,E.T)) + dot(E,dot(X,A.T)) + Q, \
zeros((2,2)))
def test_lyap(self):
A = array([[-1, 1], [-1, 0]])
Q = array([[1, 0], [0, 1]])
X = lyap(A, Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A @ X + X @ A.T + Q, zeros((2,2)))
A = array([[1, 2], [-3, -4]])
Q = array([[3, 1], [1, 1]])
X = lyap(A,Q)
# print("The solution obtained is ", X)
assert_array_almost_equal(A @ X + X @ A.T + Q, zeros((2,2)))
# Compare methods
if slycot_check():
X_scipy = lyap(A, Q, method='scipy')
X_slycot = lyap(A, Q, method='slycot')
assert_array_almost_equal(X_scipy, X_slycot)
def test_lyap_sylvester(self):
A = 5
B = array([[4, 3], [4, 3]])
C = array([2, 1])
X = lyap(A, B, C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A * X + X @ B + C, zeros((1,2)))
A = array([[2, 1], [1, 2]])
B = array([[1, 2], [0.5, 0.1]])
C = array([[1, 0], [0, 1]])
X = lyap(A, B, C)
# print("The solution obtained is ", X)
assert_array_almost_equal(A @ X + X @ B + C, zeros((2,2)))
# Compare methods
if slycot_check():
X_scipy = lyap(A, B, C, method='scipy')
X_slycot = lyap(A, B, C, method='slycot')
assert_array_almost_equal(X_scipy, X_slycot)
