Remove strong Noise in signal

Question

Thiago Petersen 2019-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/457432-remove-strong-noise-in-signal

评论： Thiago Petersen 2019-4-22

SignalNoise.mat

Hello,

I am thinking if is possible to remove a very strong noise in a signal that i recorded. The idea is to remove the noise from the files and only plot the signal (like in the first five spikes). Is it possible to locate only the major spikes and so remove it? Other problem is that sometimes the minor spikes happens in the same times than major spikes (this happens in the 6th spike). Follow the data attached. I will appreciate any help. Thanks!

2 个评论
显示无隐藏无

KALYAN ACHARJYA 2019-4-20

Do you want to remove larger spike (amplitude) or clip its aplitude in certain range?

Thiago Petersen 2019-4-20

The idea is to remove all the spikes that have larger amplitude than the first 5 spikes (the minor ones)

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2019-4-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/457432-remove-strong-noise-in-signal#answer_371484

在 MATLAB Online 中打开

Try this:

% Get threshold based on first 5 points.  It should be bigger than the tallest of those, say by 10% or whatever
threshold = max(signal(1:5)) * 1.10;
% Get "good" indexes - those with signal values less than the threshold.
goodIndexes = signal < threshold;
filteredSignal = signal(goodIndexes);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Thiago Petersen 2019-4-22

在 MATLAB Online 中打开

s2.jpg

Hi! Thanks for the reply, but it solves partially my problem. The problem in your solution is that the signal (minor spikes) and the noise (major spikes) will be affected. I only want to the bigger spikes be removed, not the smaller ones.

I was thinking to make a new vector with the noise inverted and add it in the same time that the noise occurs (using the indexes) to remove it. The noise has a fixed waveform, so maybe I could use the ginput to collect the noise data (like the .jpg attached), insert it to a vector and sum to the signal data. I tried to do it, but I got stuck after I got the noise shape.

You think this could work?

I will appreciate your reply.

mga = 2.5;
dist1 = 10000; 
fs = 48828.125; % sample rate
time =[1/fs:1/fs:length(signal)/fs];
plot(time,signal)
hold on
lim1 = max(signal(1:dist1))*mga;
[t2,peakind]= findpeaks(signal,fs,'MinPeakProminence',lim1)
plot(peakind,t2,'or')
% to get the indexes of the noise 
[peakValues, indexes] = findpeaks(signal,'MinPeakProminence',lim1) 
% this is to get the noise shape
pause % so you can zoom to find a clear noise (bigger spike)
ind2 = ginput(2) % select start and end of the noise spike
ind21 = find(time<ind2(1));
ind21 = ind21(end);
ind22 = find(time<ind2(2));
ind22 = ind22(end);
s2 = signal(ind21:ind22);
ts2 = [1/fs:1/fs:length(s2)/fs];
is2 = s2*-1; %invert the noise shape 