matlab에서 trapz을 사용하여 피크 아래의 면적을 계산하는 방법은 무엇입니까?

Question

나는 matlab에 초보자이며 피크 아래에 계산 영역 문제가 있습니다. 다음은 봉우리를 찾기 위해 사용하는 코드입니다.

%% Peak flow analysis 2 (Figure 3) 
flowtable = finalCSVnew(:,[1,7:8]);           % create table containing DateAndTime, Durchflusslm, and SummeaktuellerTagm data 

peakflowEvent = flowtable{:,2} > aa ;          % determine the threshold of flow(m3/h) for peakflowEvent 

% use false as logical vector to determine transition. With function diff, 
% transitions from false (0) to true (1) will be 1 and transitions from true 
% to false will be -1. This will be 1 at the start of a dry period and -1 after the end 
peakTransitions = diff([false; peakflowEvent; false]);        
eventStarts = find(peakTransitions == 1); 
eventEnds = find(peakTransitions == -1) -1; 

% define the peak flow of each event through the flow data (peakflow) and 
% the time when peak flow is happened (peakflowtime) 
[peakflow, peakflowlocrel] = arrayfun(@(s, e) max(flowtable.Durchflusslm(s:e)), eventStarts, eventEnds); 
peakflowlocabs = peakflowlocrel + eventStarts - 1; 
peakflowtime=flowtable.DateAndTime(peakflowlocabs); 

% create result table containing start and end time for peak flow event, the duration 
% between start and end time, and peak flow 
peakflowanalysis2 = table(flowtable.DateAndTime(eventStarts), flowtable.DateAndTime(eventEnds), ... 
       flowtable.DateAndTime(eventEnds) - flowtable.DateAndTime(eventStarts), ... 
       peakflow, peakflowtime, ... 
       'VariableNames', {'Start', 'End', 'Duration','PeakFlow','PeakFlowTime'}); 

numPeakflow2 = height(peakflowanalysis2);         % calculate the number of max flow 

% plot flow and peak flow 
figure(3) 
plot(flowtable.DateAndTime,flowtable.Durchflusslm,flowtable.DateAndTime(peakflowlocabs),peakflow,'v','MarkerFaceColor','red',... 
    'MarkerSize',5) 
xlabel('Date and Time');             % define Date and Time as x-axis 
ylabel('Flow [m3/h]');              % define Flow as y-axis 
title('Peak Flow Events (2)');            % define the title of the plot 
legend('Flow','Peak Flow','Location','Northeast','Orientation','Vertical') 
grid on                  % show grid on plot 
datacursormode on               % enable to display data value interactively in the plot 

% clear temporary variables 
clearvars peakflowEvent peakTransitions eventStarts eventEnds peakflowlocrel peakflowlocabs peakflow peakflowtime 

두 개의 테이블과 한 개의 그림을 생성합니다. 피크를 포함하고 각 이벤트를 나타내는 그림입니다. 가 이 아래 표가 피크가 필터링 된 데이터를 포함

나는 아래의 면적을 계산하는 방법을 알고 싶어 피크의 수는 표는 원래의 데이터를 포함 이벤트의 수를 = 각 이벤트의 피크 (11 개 이벤트)이며 모두를 합한 것이 아닙니다. 나는 trapz()에 대해 읽었지 만, 다중 봉우리에 적용하는 방법을 혼란스럽게합니다. 저 좀 도와 주 시겠어요? 당신의 도움을 주셔서 대단히 감사합니다.

Answer 1

먼저 테이블에 큰 데이터 세트를 저장하지 마십시오. 데이터 세트가 매우 느립니다. 저는 많은 수 문학적 데이터를 가지고 작업하며 항상 데이터를위한 숫자 배열과 timeseries를위한 숫자 배열 (datenum)을 사용합니다.

은 홍수 각 피크 이상의 피크 사용 trapz 및 루프 아래의 영역을 효율적으로 활용하려면 다음, 당신이 당신의 시작을 변환하고 flowdata의 인덱스에 종료 할 필요가

time_start=datenum(Peakflowanalysis.Start,'dd.mm.yyyy HH:MM:SS'); 
time_end=datenum(Peakflowanalysis.Ende,'dd.mm.yyyy HH:MM:SS'); 
time_array=datenum(flowtable.DateandTime,'dd.mm.yyyy HH:MM:SS'); 

for i=1:height(peakflow) 
ind1=find(time_array==time_start) 
ind2=find(time_array==time_end) 
peak_sum(i)=trapz(flowdata(Start(ind1):Ende(ind2)));% trapz method 
peak_sum(i,2)=sum(flowdata(Start(ind1):Ende(ind2)));%sum method 
end 

전에.

또한 합계 방법을 사용해보세요. 데이터의 시간 해상도에 따라 trapz 결과와 매우 유사하다는 것을 알 수 있습니다.