CODEWITHSUNDEEP
matlabsignal-processingcurve-fitting
crack_addict
crack_addict

Reputation: 653

How to make a curve smoothing in matlab?

Gradient Image

the blue plot is a noisy plot of the original plot(red). Is there any way to approximate the blue plot to nearly red plot?

Upvotes: 9

Views: 57934

Answers (4)

sergej
sergej

Reputation: 18009

gausswin() requires the Signal Processing Toolbox

smooth() requires the Curve Fitting Toolbox

If you don't have these toolboxes, here is a simple smooth() implementation:

smooth.m:

function yy = smooth(y, span)
    yy = y;
    l = length(y);

    for i = 1 : l
        if i < span
            d = i;
        else
            d = span;
        end

        w = d - 1;
        p2 = floor(w / 2);

        if i > (l - p2)
           p2 = l - i; 
        end

        p1 = w - p2;

        yy(i) = sum(y(i - p1 : i + p2)) / d;
    end
end

result for y3 = smooth(y2, 15), using @Junuxx code:

enter image description here

Upvotes: 4

0x90
0x90

Reputation: 41022

Just to add an additional option :

Use cftool in the prompt of matlab:

enter image description here

Upvotes: 2

Junuxx
Junuxx

Reputation: 14281

Let's define a wavy function:

x = 0:.1:20;
y1 = 5*sin(x) + 2*x - x.^2 +.3*x.^3 - .2*(x-15).^4 - 10*x.^2.*cos(x./3+12).^3 + .5*(x-12).^4;

And add lots of noise:

r = randi(1000,1,201) - 500;
y2 = y1+r;

Now make a 1D Gaussian filter, normalize it and convolve it with our function:

g = gausswin(20); % <-- this value determines the width of the smoothing window
g = g/sum(g);
y3 = conv(y2, g, 'same')

Let's see the result

figure;
hold on; 
plot(y1, 'r', 'linewidth', 3); 
plot(y2, 'b'); 
plot(y3, 'g', 'linewidth', 3);

Red the original function, blue the noisy version, green the smoothed, 'recovered' function.

line graph of smoothed function

Upvotes: 14

bla
bla

Reputation: 26069

another option is to use 'smooth'. I like to use it because it is a single line function. Using the code of the previous answer by @Junuxx:

x = 0:.1:20;
y1 = 5*sin(x) + 2*x - x.^2 +.3*x.^3 - .2*(x-15).^4 - 10*x.^2.*cos(x./3+12).^3 + .5*(x-12).^4;
r = randi(1000,1,201) - 500;
y2 = y1+r;

Now apply smooth:

ys = smooth(x,y2,0.25,'rloess');
plot(x,y2,x,ys)

enter image description here

For more info:

doc smooth

Upvotes: 11

Related Questions