Linear Quadratic regulation (LQR) method

I have a model in state space : x^∙=Ax+Bu Y=Cx+Du

A =

         0    1.0000         0         0
         0  -85.0455   58.8600         0
         0         0         0    1.0000
         0   28.3485  -22.8900         0

B =

         0         0
    0.5551         0
         0         0
   -0.1850    0.0006

C =

     1     0     0     0
     0     1     0     0
     0     0     1     0
     0     0     0     1

D =

     0     0
     0     0
     0     0
     0     0
 when I use Linear Quadratic regulation (LQR) method When using this method you receive an error 
Error using lqr (line 43)
In the "lqr(A,B,Q,R,N)" command, the A and Q matrices must have the same size.
I want a solution to this problem