How can i calculate FRF from random vibraition data recoreded using accelerometer?

I have output-only response data of time vs acceleration recorded using accelerometer. Input force is not known. Anyone can help me in detail, How to calculate the frf and then I want to identify the resonant frequencies, damping and mode shapes of the structure? It would be very helpful if someone can share MATLAB code regarding this above mentioned problem.