There has been some input matrix diagonalization in the past by Yuta and Kiwamu, but I find the automation to be not totally satisfactory.
It would be better if we could automatically fit the data to find the Suspended optic eigenfrequencies and then use that to get the matrix. So** I wrote a peak fitter to get the matrix**.
It gets the data from mafalda with NDS2, then it makes the PSDs, and then starts with some initial guesses (based on looking at the plots) and them uses fminsearch to get the peak frequencies and Q's.
Using the output of this, we can use Yuta's method and take the passive transfer functions with the free swing data (from April 30, so *we got do do it quick*) to get the input matrix.
Doing the SUS input matrix is nice for having good damping (as long as we remember to include SIDE), but my **motivation is to produce a good null stream from the 4 face sensors so that we can estimate the sensor noise**s at all times. |