RobustDetrend

Find the best polynome order and detrend data while preserving peak features in data.
To generate the miniature :

x = linspace(0,1,1000)';
peak = 0.008.*exp(-((x-0.3).^2)./(2*(0.05^2)))+0.01.*exp(-((x-0.7).^2)./(2*(0.025^2)));
peak = peak+polyval([-0.021,+0.008,0.021,0],x);
peak = peak-0.001+0.002.*rand(size(peak));
[polorder,~,trend] = RobustDetrend(peak,6,0.975,x);
figure; plot(x,peak,'b.',x,trend,'r--')