使用matlab / octave找到最小反射点

Question

我有一个红线数据集。 如果我拍摄数据集的反射/镜像，我试图找到以黄色突出显示的最小点。

请参阅下面的示例代码/图表我试图找到一种方法来查找数据集的反射/镜像图像中以黄色突出显示的最小点（蓝色线条） ）低于反射线（黑线）。

请注意，这只是一个简单的数据集，会有大得多的数据集，大约100000 +

PS：我使用Octave 3.8.1，就像matlab

clear all,clf, clc,tic

x1=[0.;2.04;4.08;6.12;8.16;10.2;12.24;14.28;16.32;18.36]
y1=[2;2.86;4;2;1;4;5;2;7;1]
x2=[0.;2.04;4.08;6.12;8.16;10.2;12.24;14.28;16.32;18.36]
y2=abs(y1-max(y1));
data1 = y2; 
reflection_line=max(y1)/2
[pks3 idx3] = findpeaks(data1,"DoubleSided","MinPeakHeight",0.1);

line([min(x1) max(x1)], [reflection_line reflection_line]);
hold on;
plot(x1,reflection_line)
hold on;
plot(x1,y1,'-r',x2,y2,'-b')

Answer 1

我没有反映您的原始数据，而是找到原始值的局部最大值，这些值大于您给定的行。

使用findpeaks的DIY替代品： 如果值大于（或等于）其前任和后继者，则值为局部最大值。

%% Setup
x1 = [0.;2.04;4.08;6.12;8.16;10.2;12.24;14.28;16.32;18.36];
y1 = [2;2.86;4;2;1;4;5;2;7;1];
reflection_line = max(y1)/2;
%% Sort by x value
[x1, I] = sort(x1);
y1 = y1(I);
%% Compute peaks
maxima = @(y) [true; y(2:end)>=y(1:end-1)] & ... % Value larger than predecessor
              [y(1:end-1)>=y(2:end); true];      % Value larger than successor
maximaLargerThanLine =  maxima(y1) & (y1>reflection_line);
%% Plotting
plot(x1,y1);
hold on;
plot(x1(maximaLargerThanLine),y1(maximaLargerThanLine),'rx');
line([min(x1) max(x1)], [reflection_line reflection_line]);