I fit the MC1 -> MCL TF using `vectfit4.m` (from mDV). The wrapper file is **mDV/extra/C1/ fitMC12MCL.m** .
Plotted here are the data (RED), the fit (BLUE), and the residual __x10__ (GREEN).
For the magnitude plot, residual is defined as ------ res = 1 - fit / data
For the phase plot the residual is defined as ------- res = phase(data)-phase(fit)
You can see that the agreement is very good. The phase match is better than 5 deg everywhere below 10 Hz.
This TF is so smooth that we could have probably done without using this, but its good to excercise the method anyway. |