Solve Rational Matrix Equation

Given a symmetrix positive definite matrix Q and a non-singular matrix L, Find symmetric positive definite solution X such that X = Q + L (X inv) L^T.

Solve Rational Matrix


Q <- rbind(c(2,-1,0),c(-1,2,1),c(0,1,2))
L <- rbind(c(2,-2,18),c(2,1,0),c(1,2,0))

##QR Decomposition

##LQ Decompostion

##Solve Rational Matrix
X <- sol.rationalmatrix.euqation(Q, L)

#Checking the results
X - (Q + L %*% solve(X) %*% t(L))

Try the SolveRationalMatrixEquation package in your browser

Any scripts or data that you put into this service are public.

SolveRationalMatrixEquation documentation built on May 2, 2019, 2:45 p.m.