clc;clear;close all;
cd("D:\桌面\毕设\材料\摩阻扭矩\MATLAB")
%% 数据读取
trajectory_parameters = xlsread('LH11-1-B4H3.xlsx','轨道','A2:C1000'); % 读取轨道数据，依次为井深、井斜角、井斜方位角
pipe_parameters = xlsread('LH11-1-B4H3.xlsx','管柱数据','C5:G100'); % 读取套管柱数据，依次为内径，外径，接头外径，强度，长度
well_parameters = xlsread('LH11-1-B4H3.xlsx','井身结构','A2:C3'); % 读取井身结构数据，依次为1上层导管下入深度，2上层套管内径，3井眼直径
drill_parameters = xlsread('LH11-1-B4H3.xlsx','钻进数据','A2:G10000'); % 读取钻进数据，依次为1井深，2钻压，3钻头扭矩，4转速，5管柱运动速度，6排量，7钻井液密度
HK_TOR = xlsread('LH11-1-B4H3.xlsx','测量值','A2:C1000'); % 读取1井深2钩载3转盘扭矩
%%
drill_parameters = drill_parameters(1:57,:);
HK_TOR = HK_TOR(1:57,:);
% drill_parameters = drill_parameters(58:105,:);
% HK_TOR = HK_TOR(58:105,:);
tic
[alpha,k,kal,kph,tau,dk_ds,tz,nz,bz,dalpha_ds,dpha_ds]=traj_para_linear(trajectory_parameters); % 自定义函数计算插值后的计算轨迹参数
dL = 10;
[mu1, mu2, delta] = cal_mu1_and_mu2_fmin(pipe_parameters,drill_parameters, well_parameters,HK_TOR,alpha,k,kal,kph,tau,dk_ds,tz,nz,bz,dL);
toc
fprintf("摩阻系数1: %.4f\n",mu1)
fprintf("摩阻系数2: %.4f",mu2)
