how to remove noise from curves and take their derivates

12 visualizaciones (últimos 30 días)
KIRAN noor
KIRAN noor el 16 de En. de 2023
Editada: KIRAN noor el 8 de Feb. de 2023
Image Analyst ha marcado con alerta este/a pregunta
Hello,
I have some curves which are not smooth, I have to take their derivative. therefore first requirement is remove the noise and take the derivative.
I am doing this work through curve fitting using rat35, poly9 etc. and then taking the derivative. but everytime i run the script, result changes slighty.
i am attaching the curve , their zoom version and then warnings which appeared in workspace, would you please guide me how i should handle this issue.
Regards,
Kiran

Respuesta aceptada

Bruno Luong
Bruno Luong el 16 de En. de 2023
Editada: Bruno Luong el 16 de En. de 2023
I would chose spline fitting rather than polynomial (bad) or rational.
Look at the shape of your data, I doubt any simple analytic model can fit it well.
c=readtable('100_li_4th_discharge1.xlsx');
x=table2array(c(:,1));
y=table2array(c(:,2));
% https://fr.mathworks.com/matlabcentral/fileexchange/25872-free-knot-spline-approximation
pp=BSFK(x,y,[],[],[],struct('sigma',1e-9,'KnotRemoval','none'));
ppd=ppder(pp); % derivative
xi=linspace(min(x),max(x),1024);
figure
h1=plot(x,y,'.g',xi,ppval(pp, xi),'r');
yyaxis right
h2=plot(xi,ppval(ppd,xi),'k')
legend([h1; h2],'data','fit','derivative');
  3 comentarios
Bruno Luong
Bruno Luong el 17 de En. de 2023
Editada: Bruno Luong el 17 de En. de 2023
ppd is dy/dx, if you want dx/dy it is 1./ppd.
BTW in this case you could fit x vs y instead of the opposite.
c = readmatrix('https://fr.mathworks.com/matlabcentral/answers/uploaded_files/1265640/100_li_4th_discharge1.xlsx')
Q=c(:,1);
V=c(:,2);
% https://fr.mathworks.com/matlabcentral/fileexchange/25872-free-knot-spline-approximation
pp=BSFK(V,Q,[],[],[],struct('sigma',1e-9,'KnotRemoval','none'));
ppd=ppder(pp); % derivative
yi=linspace(min(V),max(V),1024);
figure
h1=plot(V,Q,'.g',yi,ppval(pp, yi),'r');
xlabel('V')
ylabel('Q')
yyaxis right
h2=plot(yi,ppval(ppd,yi),'k')
legend([h1; h2],'data','fit','derivative');
ylabel('dQ/dV)*')
Matt J
Matt J el 17 de En. de 2023
@IQRA KIRAN If it worked, then you should Accept-click the answer.

Iniciar sesión para comentar.

Más respuestas (2)

Matt J
Matt J el 16 de En. de 2023
Editada: Matt J el 16 de En. de 2023
If you provide a StartPoint in your fit options, the warning will go away and the fit will be the same when you rerun.
However, the differences you see with each run are not meaningful. You get a slightly different result because the start point of the iterative fitting process is different each time, but in this case they all seem to be converging to essentially the same curve.
  6 comentarios
Matt J
Matt J el 16 de En. de 2023
Why would you need another way?
Matt J
Matt J el 16 de En. de 2023
Editada: Matt J el 16 de En. de 2023
Use a lower order polynomial so that the parameter estimates won't be as sensitive. Or use a non-polynomial model. Image Analyst suggested an exponential fit. A smoothing spline might also work,
Basically, you should experiment with other models to see if a stiffer model would work.

Iniciar sesión para comentar.


Image Analyst
Image Analyst el 16 de En. de 2023
You can use sgolayfilt, which is a sliding polynomial fit within the sliding window, however if I were going to take the derivatives of your curves, I'd fit them to an exponential decay (to get a formula with numerical coefficients), which has a known derivative. See my attached exponential decay demo which uses fitnlm. I have many other functions fit with fitnlm, so it you want them, just ask.

Categorías

Más información sobre Spline Postprocessing en Help Center y File Exchange.

Productos


Versión

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by