% 一维直流电测深反演
clear
% close all
clc
% 电性参数
nlayer = 3
res(1) = 10
res(2) = 1
res(3) = 30
% res(4) = 8
thick(1) = 5
thick(2) = 5
thick(3) = 100
% thick(4) = 50
% AB/2的数据,实际上AB/2 = r,单位是m
r(1) = 1
r(2) = 2
r(3) = 3
r(4) = 4.5
r(5) = 6
r(6) = 9
r(7) = 12
r(8) = 15
r(9) = 20
r(10) = 30
r(11) = 45
r(12) = 60
r(13) = 90
r(14) = 120
r(15) = 150
r(16) = 200
r(17) = 300
r(18) = 450
r(19) = 600
r(20) = 900
r(21) = 1200
r(22) = 2000
% 一维直流电测深正演子程序
nr = length(r)
[Rs] = sdc1dfwd(nr,r,nlayer,res,thick)
% 画一维直流电测深视电阻率曲线类型图
figure(1)
loglog(r,Rs,'r->')
grid on
xlabel('AB/2 (m)')
ylabel('\rho_s (\Omega-m)')
axis equal
% 保存计算数据
r = r(:)
Rs = Rs(:)
clear a
a = [r Rs]
dlmwrite('smod_0001.txt',a)
% b = [res thick]
% dlmwrite('mod_qq_0066.txt',b)