Matlab分形维数与可钻性级值

本文主要是介绍Matlab分形维数与可钻性级值，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

function [A] = Brokenfractalmethod(r,Mr1,Mr2,Mr3,Mr4)
%牙轮钻头破碎岩石后，产生许多大小不一的碎块，收集所有碎块，应用分级筛网称其不同尺寸下的累计质量分布；
%然后用质量统计法求得岩石碎块尺度分布的分形维数，即破碎分形维数；
%实验表明，分形维数与可钻性级值关系密切，因此可用分形维数来表达岩石可钻性；
%本例是求分形维数以及上返岩屑筛分回归相关性系数；
%  例： A=xlsread('sampleBrokenfractalmethod.xlsx');
%       r=A(1,2:end); 不同筛网尺寸；
% Mr1=A(2,2:end);     样本深度在3270米处不同筛网尺寸下可通过质量百分比
% Mr2=A(3,2:end);     3280米
% Mr3=A(4,2:end);     3290米
% Mr4=A(5,2:end);     3300米%%   确定不同筛网尺寸及在此条件下各样本深度通过的质量百分比
A=xlsread('sampleBrokenfractalmethod.xlsx');
r=A(1,2:end);
H=A(2:end,1)';     %样本井深
Mr1=A(2,2:end);
Mr2=A(3,2:end);
Mr3=A(4,2:end);
Mr4=A(5,2:end);%%   取双对数坐标系,并求出各拟合参数
x=log(r);
y1=log(Mr1);
y2=log(Mr2);
y3=log(Mr3);
y4=log(Mr4);
p1=polyfit(x,y1,1);        %拟合函数参数值
p2=polyfit(x,y2,1);
p3=polyfit(x,y3,1);
p4=polyfit(x,y4,1);

这篇关于Matlab分形维数与可钻性级值的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---