We use a new method of matrix decomposition for
r
-circulant matrix to get the determinants of
A
n
=
Circ
r
F
1
,
F
2
,
…
,
F
n
and
B
n
=
Circ
r
L
1
,
L
2
,
…
,
L
n
, where
F
n
is the Fibonacci numbers and
L
n
is the Lucas numbers. Based on these determinants and the nonsingular conditions, inverse matrices are derived. The expressions of the determinants and inverse matrices are represented by Fibonacci and Lucas Numbers. In this study, the formulas of determinants and inverse matrices are much simpler and concise for programming and reduce the computational time.