
Offset removal from a sinusoidal wave
8 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Mohammad Zulqarnain
el 18 de Feb. de 2021
Comentada: Star Strider
el 19 de Feb. de 2021
Hello Everyone,
I have an input waveform as shown in the figure in blue. It is a sinusoidal wave but due to different offests during positive and negative half cycles, it looks distorted. Basically, there is a big jump whenever signal changes direction from positive to negative and vice versa. So, I want to remove offsets from individual half cycles. I want an output which looks as red in the figure, which is just shifting of top and bottom parts of the waveform, taking zero as the baseline, excluding the big jumps. I have attached the input csv file as well. If anyone has any idea, how to approach it, please let me know, I will be thankful.

0 comentarios
Respuesta aceptada
Star Strider
el 18 de Feb. de 2021
Try this:
D1 = readmatrix('Output.csv');
x = 1:numel(D1);
D1HL = min(D1(D1>0));
D1HI = (D1>D1HL); % Logical Index
D1LL = max(D1(D1<0));
D1LI = (D1<D1LL); % Logical Index
minH = zeros(size(D1))+min(D1(D1HI));
maxL = zeros(size(D1))+max(D1(D1LI));
D1C = zeros(size(D1));
D1C(D1HI) = D1(D1HI)-minH(D1HI);
D1C(D1LI) = D1(D1LI)-maxL(D1LI);
figure
plot(x, D1C)
grid
xlim([0 100]) % Delete Or Change To See Larger (Or Entire) Vector
producing:

I plotted the points as well as the connecting lines to demonstrate the continuity. Change the line style to '-' to plot it without the points.
.
4 comentarios
Más respuestas (1)
darova
el 18 de Feb. de 2021
I have an idea
clc,clear
x = linspace(0,20,100);
y = sin(x);
ind = -0.5<y & y<0.5;
x1 = x(~ind);
y1 = y(~ind);
y1 = y1 - sign(y1)*0.3;
plot(x,y)
hold on
plot(x1,y1,'.r')
hold off

Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!