I fit the MC1 -> MCL TF using
vectfit4.m (from mDV). The wrapper file is
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.