运行软件:MATLAB
% function A=datacut(data,step)
% 在智能钻井研究中,需要经使用到测井数据。测井数据在深度上,如何提取将是一个很重的事情。将测井数据换算成任意整数步长的数据,步长内的数据取平均值。
% data为数据的测井数据,step为截取深度数据的步长,A为返回值
% nepu liw 2017年10月4日
% Example: A=datacut()
%-------------------------------------------------------------------
data=importdata('z_log01.txt','\t');
step=1;
deepmax=max(data(:,1));
if step<1
disp('输入数据小于1')
elseif step>deepmax
disp('输入数据大于行数')
else
% fprintf('data1=%4.1f\n',data(1:10,1));
data(:,1)=data(:,1)/step;
% fprintf('data2=%4.2f\n',data(1:10,1));
%-------------------------------------------------------------------
deep=floor(data(:,1)); %向下取整
deep1=unique(deep);
[m_deep,~]=size(deep1);
for I=1:m_deep
B=find(deep==deep1(I));
C&