So far I'm having problem with adopting your transfer function code to my current data (see above uploaded .mat file).
This portion:
% Now go get the transfer function:
options = optimset('MaxFunEvals',1000000,'TolFun',.0000001,'Display','on');
FM0=[ .1 1 0.15 0.23 ]
for n=1: 5
[FM,error] =...