对于k=1:4
对于j=1:3
对于我= 1:41
dTdt(k,j,i,:)=diff(T(k,j,i,:,2))/diff(T(k,j,i,:,1));
%[val,Idx] = max(dtdt(i,1,:))
idx = find(dtdt(k,j,i,:)> 0,1,“第一”)
DT(k,j,i)= t(k,j,i,Idx,1)
T(k,j,i,idx,1)
终止
终止
终止