How to find the end value of a flat minimum peak in a noisy data?

Question

Bálint Kovács 2022-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1784500-how-to-find-the-end-value-of-a-flat-minimum-peak-in-a-noisy-data

评论： Bálint Kovács 2022-8-24

angle.mat

Hi!

I have a problem finding the end value of minimum flat peaks. I have angular displacement data form a dinamometer. But this data is a tranmitted signal to an EMG system, so the y axis shows voltage, that's why the scaling is unusual. The original data had noise which I filtered (see below). The offset of the signal is varying from measurament to measurament.

I used findepeak for locating the peaks which are the end points of the tasks. I used high tresholds therfore it was easy to find these peaks. I red circled the points I'm trying to locate, that's where the task is started. However, I cannot solve the issue here. The findpeak function does not have an option for a flat peak, and islocalmin influanced by the small noise in the dataset. It is already filtered but still have small noise wich interfere with the islocalmin function. I don't know how should I filter more the dataset or how to locate the flat end with the slightly noisy data with different approach. Does anybody can help? I think it is not a big problem only for me, but please let me know if you have any suggestion.

Thanks!

I use this script:

load ('angle.mat');
%% I used this basic lowpass filter on the raw data
% Fs = 2000;  
% N  = 4;   
% Fc = 10;  
% h  = fdesign.lowpass('N,F3dB', N, Fc, Fs);
% Hd = design(h, 'butter');
% Angle = filter (Hd, Angle);
%% I found the max peak like this (These are the end points for the task)
[Maxima2,MaxIdx2] = findpeaks(Angle, 'MINPEAKHEIGHT',10, 'MINPEAKDISTANCE', 5000);
index_end=[Maxima2,MaxIdx2];
%%I don't have idea how to find the starting points.