clc
clear all
close all
M = csvread('shui_A-Aref.csv', 1, 2);
N = csvread('kongA-Aref.csv', 1, 2);
baseline_1 = mean(M(1:16));
baseline_2 = mean(N(1:16));
A = length(M);
B = length(N);
j = 0;
k = 0;
a = [];
c = [];
d = [];
threshold = -1050;
将A数组的数据一个个输入
for i = 17:A
if M(i) < baseline_1
baseline_1 = M(i);
j = 0;
elseif M(i)>baseline_1 && M(i) < threshold
baseline_1 = baseline_1*0.75+0.25*M(i);
j = 0;
elseif M(i)>=threshold
j = j+1 ;
if j >= 8
a = [a;i,M(i)];
fprintf('chufa\n')
end
end
end
% 将B数组的数据一个个输入
for h = 17:B
if N(h) < baseline_2
baseline_2 = N(h);
k=0;
elseif N(h)>baseline_2 && N(h)<threshold
baseline_2 = baseline_2*0.75+0.25*N(h);
elseif N(h)>=threshold
k = k+1 ;
if k == 8
c = [c;(h-7),N(h-7)];
elseif k>8 && N(h+1)<threshold
d = [d;h,N(h)];
end
end
end
plot(N)
hold on
scatter(c(:,1),c(:,2),'*')
hold on
scatter(d(:,1),d(:,2),'x')
hold off