remove noise from data by polynomial fitting

Hi everyone, I want to remove noise from my data and for that I want to fit a polynomial. Can anyone plz tell me how can we remove noise in matlab?

Thanks

[Matlabuser]

I have a time series of temperature and I want to fit a polynomial of this kind T=A (log t)+B (log t)^2+ .......... so on

Answer 1

[Image Analyst]

x = log(t);
coefficients = polyfit(x, T, 1);

See full demo (with a different equation) attached.

[Matlabuser]

hi, Thank You I tried to use this command, but I am not getting anything good.It inverted my time series. I did like this:

load T;
t=[10:10:2500];
x=log(t);
coff= polyfit(x,T,1);
A1=polyval(Coeff,x);
plot(x,A1);

I am new to Matlab, So I am not able to understand this problem

[Image Analyst]

You need to get t back from x by using exp() and then plot(t, A1, 'bo-'); Or better yet, make more points like I did in my demo so you can see the shape of the equation.

[Matlabuser]

I got It. One more thing I want to know that now this plot is plot is in between values of T obtained after logarithmic fitting and time.

Thank you so much !!!

[Image Analyst]

You need to give more interpolation points, like

% Setup:
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
imtool close all;  % Close all imtool figures if you have the Image Processing Toolbox.
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
t=[10:100:2500];
% T=A (log t)+B (log t)^2
A = 5;
B = 3;
T = A * log(t) + B * log(t).^2
% Add noise so it's not so perfect.
T = T + 50*rand(1, length(T));
% Plot training data
plot(t, T, 'bo', 'LineWidth', 2);
grid on;
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0.1 1 .75]);
hold on;
%--------------------------------------------------------------
% Find polynomial
x=log(t); % Transform t so it's a lienar least squares situation.
% Now T = A*x + B*x^2 and we can fit a second order polynomial.
coefficients = polyfit(x,T,2)
% Interpolate additional points for a total of 300.
extra_x = linspace(x(1), x(end), 300);
fittedData = polyval(coefficients, extra_x);
% But we want to plot against t, not x so we need to recover t from extra_x.
extra_t = exp(extra_x);
plot(extra_t, fittedData, 'r*-', 'LineWidth', 2);
xlabel('t', 'FontSize', fontSize);
ylabel('T', 'FontSize', fontSize);
title('T vs. t', 'FontSize', fontSize);
legend('Training Data', 'Fit');

Answer 2

[Torsten]

http://www.mathworks.com/matlabcentral/fileexchange/30299-savitzky-golay-smooth-differentiation-filters-and-filter-application

Best wishes

Torsten.

Answer 3

[Matt J]

See the polyfit command .

[Matlabuser]

I already tried with command, I want to fit polynomial manually

Answer 4

[Matlabuser]

Can anyone please tell me how can I fit logarithmic polynomial to time series

Thanks in advance