...0.23 ]
for n=1: 5
[FM,error] = lsqcurvefit('splane2_fit',FM0 ,freq ,heave_gain,[],[],options);
heavesys =tf([FM(1:2 )],[FM(3:4) 1 ])
FM0=FM+ .001*rand(1);
end
I added above code to this live script:
[h_az_tf,n]=tfestimate(act_az,heave_az,winsize,overlap,nfft,Fs)...