finding several values close to the set point

Hello!
I would like to calculate the period in a curve using the mean. However, my attempts to date have not been successful, and I think there is a simple way to do it.
I put random data - eventually I will have a lot of different data and this is to be the method for determining it.
I managed to find the closest value, but unfortunately 1, and I need a few of them to determine the curve's symmetry axis.
Edit: It doesn't search for peaks, because the data on which I work has a few of them on the main peaks, and filtration changes the shape of the plot very much.
Thank you for every advice.
A=[1.5 2.3 3.7 4.6 5.3 4.2 3.8 2.7 1.3 2.6 3.9 4.8 6.1 5.3 4.3 3.6 2.8 1.6 2.4 3.6 4.7 5.4 4.0 2.8 1.5 2.7 3.8 4.7];
A=A';
s = size(A);
t = 1:1:s(1);
plot (t,A)
a = mean(A);
Edit2:
In the attachment I add the real data I work on. From each cell I choose 1 column as X and 3 columns as Y in plot (X, Y).

 채택된 답변

Ameer Hamza
Ameer Hamza 2020년 3월 18일

1 개 추천

The following code works by detecting the peak of each cycle and then calculate the average distance between each peak
A=[1.5 2.3 3.7 4.6 5.3 4.2 3.8 2.7 1.3 2.6 3.9 4.8 6.1 5.3 4.3 3.6 2.8 1.6 2.4 3.6 4.7 5.4 4.0 2.8 1.5 2.7 3.8 4.7];
A=A';
t = 1:1:size(A,1);
plot(t,A);
[~, index_peaks] = findpeaks(A);
t_peaks = t(index_peaks);
period = mean(diff(t_peaks));

댓글 수: 8

Caroline
Caroline 2020년 3월 18일
편집: Caroline 2020년 3월 18일
Unfortunately, the real data has a few peaks on the main peak ... oh I won't work :(
In that case, please attach a sample of real data. It will be easy to suggest a working solution using your actual data.
Caroline
Caroline 2020년 3월 18일
Okay, I added in the attachment.
Check this code
load('test.mat');
a = sila{1};
s = a(:,3);
s = detrend(s);
s_ = s-mean(s);
[~, idx] = findpeaks(s, 'MinPeakProminence', max(s_));
period = mean(diff(idx));
t = 1:numel(s);
plot(t, s, '-', t(idx), s(idx), '+');
It reads the third column and then find the peak of each cycle. It then finds the average difference between peaks to calculate the period.
Caroline
Caroline 2020년 3월 18일
Thank you very much, this is the first way I was thinking but I couldn't deal with these peaks, Thank you very much !!
Glad to be of help.
Caroline
Caroline 2020년 3월 18일
Analyzing the data, I noticed that the solution has 1 drawback. Peaks occur in different places, therefore the periods are not equal :(.
Yes, that can be issue, but since you are taking the average of each interval so this difference is somewhat compensated.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Polar Plots에 대해 자세히 알아보기

태그

질문:

2020년 3월 18일

댓글:

2020년 3월 18일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by