电法勘探实验代码

weixin_45862237 2020-05-20 12:20:29
电阻率法数值模拟实验和激发极化法数字模拟实验报告 求matlab代码 数据越简单易懂越好
...全文
80 2 打赏 收藏 转发到动态 举报
写回复
用AI写文章
2 条回复
切换为时间正序
请发表友善的回复…
发表回复
weixin_45862237 2020-05-21
  • 打赏
  • 举报
 回复
引用 1 楼 hyff123的回复:
哈哈你也是新学的吧，这里是我的代码，探讨一下。
dy=1
ny=101
ymin=-50
dx=1
nx=101
xmin=-50
x=xmin:dx:(xmin+(nx-1)*dx)
y=ymin:dy:(ymin+(ny-1)*dy)
[X,Y]=meshgrid(x,y)
p1=10
u=0.02
r=5
h=9
%中间梯度装置
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
plot(x,ps),title('中间梯度主刨面异常')
grid on;
hold on
j=0
figure
for y=-5:1:5
j=j+1
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
q(j)=plot(x,ps)
hold on
end
grid on;
title('中间梯度随刨面位置变化')
figure
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+X.^2+Y.^2).^5/2]
pcolor(X,Y,ps);shading interp
title('中间梯度等值线图')
% 联合刨面
MN=4
AM=10
AN=14
BN=10
BM=14
dx=1
nx=101
x=xmin:dx:(xmin+(nx-1)*dx)
KK=AM*AN/MN
da=((x-2-AM).^2+h^2).^(1/2)
db=((x+2+BN).^2+h^2).^(1/2)
rm=((x-2).^2+h^2).^(1/2)
rn=((x+2).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM^2)./(da.*rm.*2)
cosb=(da.^2+rn.^2-AN^2)./(da.*rn.*2)
psA=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN^2)./(db.*rn.*2)
psB=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
figure
xx=plot(x,psA)
hold on
yy=plot(x,psB)
legend("psA","psB")
title("联合刨面低阻球体上方视电阻率区线")
%对称四极
figure
for x=0:1:100
AB=20:10:1000
MN=4
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((x-2-AM).^2+h^2).^(1/2)
db=((x+2+BN).^2+h^2).^(1/2)
rm=((x-2).^2+h^2)^(1/2)
rn=((x+2).^2+h^2)^(1/2)
cosa=(da.^2+rm^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
psAB=(psA1+psB1)./2
plot(log10(AB/2),psAB)
hold on
end
title("视电阻率测深曲线")
figure
x=-50:1:50
ab=10:4:100;
[AB,X]=meshgrid(ab,x)
MN=4
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((X-2-AM).^2+h^2).^(1/2)
db=((X+2+BN).^2+h^2).^(1/2)
rm=((X-2).^2+h^2).^(1/2)
rn=((X+2).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn.^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
psAB=(psA1+psB1)./2
pcolor(X,log10(AB./2),psAB),shading interp
set(gca,'YDir','reverse');
colorbar
xlabel('X(m)'),ylabel('AB/2(m)');
title("拟断面图")
%这是数值模拟实验
% 联合刨面
r=5
h=10
dy=1
ny=101
ymin=-50
dx=1
nx=101
xmin=-50
p1=10
p2=8
MN=1
AM=10
AN=14
BN=10
BM=14
dx=1
nx=101
n=0.3
u=p2/p1
y=ymin:dy:(ymin+(ny-1)*dy)
x=xmin:dx:(xmin+(nx-1)*dx)
[X,Y]=meshgrid(x,y)
KK=AM*AN/MN
da=((x-0.5-AM).^2+h^2).^(1/2)
db=((x+0.5+BN).^2+h^2).^(1/2)
rm=((x-0.5).^2+h^2).^(1/2)
rn=((x+0.5).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM^2)./(da.*rm.*2)
cosb=(da.^2+rn.^2-AN^2)./(da.*rn.*2)
psA=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN^2)./(db.*rn.*2)
psB=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u=p3/p1
psA2=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
nsA=1-psA./psA2
nsB=1-psB./psB2
figure
plot(x,nsA)
hold on
plot(x,nsB)
title("低阻球体上方连剖视极化率曲线")
legend("nsA","psB")
xlabel("X"),ylabel("ns")
%中间梯度装置
p1=10
p2=8
u=p2/p1
ps1=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
ns=1-ps1./ps2
figure
plot(x,ns)
title("中间梯度视极化率主刨面图")
xlabel("X"),ylabel("ns")
hold on
figure
for y=-5:1:5
u=p2/p1
ps1=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
ns=1-ps1./ps2
plot(x,ns)
hold on
end
grid on;
title('随刨面位置变化')
xlabel("X"),ylabel("ns")
figure
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+X.^2+Y.^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+Y.^2+X.^2).^5/2]
ns=1-ps1./ps2
pcolor(X,Y,ns);shading interp
colorbar
title('中间梯度等值线图')

%对称四极
figure
for x=0:1:50
u=p2/p1
AB=2:4:100
MN=1
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((x-0.5-AM).^2+h^2).^(1/2)
db=((x+0.5+BN).^2+h^2).^(1/2)
rm=((x-0.5).^2+h^2)^(1/2)
rn=((x+0.5).^2+h^2)^(1/2)
cosa=(da.^2+rm^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u2=p3/p1
psA2=p1*[1+2*(u2-1)./(2*u2+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*(u2-1)/(2*u2+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
nsA=1-psA1./psA2
nsB=1-psB1./psB2
ns=(nsA+nsB)/2
plot(log10(AB/2),ns)
hold on
end
grid on
title("对称四级视极化率")
xlabel("AB/2"),ylabel("ns")
figure
x=-50:1:50
ab=10:4:100;
[AB,X]=meshgrid(ab,x)
MN=1
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((X-0.5-AM).^2+h^2).^(1/2)
db=((X+0.5+BN).^2+h^2).^(1/2)
rm=((X-0.5).^2+h^2).^(1/2)
rn=((X+0.5).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn.^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u=p3/p1
psA2=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
nsA=1-psA1./psA2
nsB=1-psB1./psB2
ns2=(nsA+nsB)/2
pcolor(X,AB,ns2); shading interp
set(gca,'YDir','reverse');
xlabel('X(m)'),ylabel('AB/2(m)');
colorbar
title("拟断面图")
%这是激发极化法
用不同的数据运行一下。
谢谢兄弟的代码，我回去好好研究一下
hyff123 2020-05-20
  • 打赏
  • 举报
 回复
哈哈你也是新学的吧，这里是我的代码，探讨一下。
dy=1
ny=101
ymin=-50
dx=1
nx=101
xmin=-50
x=xmin:dx:(xmin+(nx-1)*dx)
y=ymin:dy:(ymin+(ny-1)*dy)
[X,Y]=meshgrid(x,y)
p1=10
u=0.02
r=5
h=9
%中间梯度装置
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
plot(x,ps),title('中间梯度主刨面异常')
grid on;
hold on
j=0
figure
for y=-5:1:5
j=j+1
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
q(j)=plot(x,ps)
hold on
end
grid on;
title('中间梯度随刨面位置变化')
figure
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+X.^2+Y.^2).^5/2]
pcolor(X,Y,ps);shading interp
title('中间梯度等值线图')
% 联合刨面
MN=4
AM=10
AN=14
BN=10
BM=14
dx=1
nx=101
x=xmin:dx:(xmin+(nx-1)*dx)
KK=AM*AN/MN
da=((x-2-AM).^2+h^2).^(1/2)
db=((x+2+BN).^2+h^2).^(1/2)
rm=((x-2).^2+h^2).^(1/2)
rn=((x+2).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM^2)./(da.*rm.*2)
cosb=(da.^2+rn.^2-AN^2)./(da.*rn.*2)
psA=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN^2)./(db.*rn.*2)
psB=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
figure
xx=plot(x,psA)
hold on
yy=plot(x,psB)
legend("psA","psB")
title("联合刨面低阻球体上方视电阻率区线")
%对称四极
figure
for x=0:1:100
AB=20:10:1000
MN=4
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((x-2-AM).^2+h^2).^(1/2)
db=((x+2+BN).^2+h^2).^(1/2)
rm=((x-2).^2+h^2)^(1/2)
rn=((x+2).^2+h^2)^(1/2)
cosa=(da.^2+rm^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
psAB=(psA1+psB1)./2
plot(log10(AB/2),psAB)
hold on
end
title("视电阻率测深曲线")
figure
x=-50:1:50
ab=10:4:100;
[AB,X]=meshgrid(ab,x)
MN=4
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((X-2-AM).^2+h^2).^(1/2)
db=((X+2+BN).^2+h^2).^(1/2)
rm=((X-2).^2+h^2).^(1/2)
rn=((X+2).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn.^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
psAB=(psA1+psB1)./2
pcolor(X,log10(AB./2),psAB),shading interp
set(gca,'YDir','reverse');
colorbar
xlabel('X(m)'),ylabel('AB/2(m)');
title("拟断面图")
%这是数值模拟实验
% 联合刨面
r=5
h=10
dy=1
ny=101
ymin=-50
dx=1
nx=101
xmin=-50
p1=10
p2=8
MN=1
AM=10
AN=14
BN=10
BM=14
dx=1
nx=101
n=0.3
u=p2/p1
y=ymin:dy:(ymin+(ny-1)*dy)
x=xmin:dx:(xmin+(nx-1)*dx)
[X,Y]=meshgrid(x,y)
KK=AM*AN/MN
da=((x-0.5-AM).^2+h^2).^(1/2)
db=((x+0.5+BN).^2+h^2).^(1/2)
rm=((x-0.5).^2+h^2).^(1/2)
rn=((x+0.5).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM^2)./(da.*rm.*2)
cosb=(da.^2+rn.^2-AN^2)./(da.*rn.*2)
psA=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN^2)./(db.*rn.*2)
psB=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u=p3/p1
psA2=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(cosa./(da.^2+rm.^2)-cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*KK*(u-1)/(2*u+1)*r^3*(-cosc./(db.^2+rm.^2)+cosd./(db.^2+rn.^2))]
nsA=1-psA./psA2
nsB=1-psB./psB2
figure
plot(x,nsA)
hold on
plot(x,nsB)
title("低阻球体上方连剖视极化率曲线")
legend("nsA","psB")
xlabel("X"),ylabel("ns")
%中间梯度装置
p1=10
p2=8
u=p2/p1
ps1=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2)./(h^2+x.^2).^5/2]
ns=1-ps1./ps2
figure
plot(x,ns)
title("中间梯度视极化率主刨面图")
xlabel("X"),ylabel("ns")
hold on
figure
for y=-5:1:5
u=p2/p1
ps1=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*x.^2+y^2)./(h^2+x.^2+y^2).^5/2]
ns=1-ps1./ps2
plot(x,ns)
hold on
end
grid on;
title('随刨面位置变化')
xlabel("X"),ylabel("ns")
figure
ps=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+X.^2+Y.^2).^5/2]
p3=p2/(1-n)
u=p3/p1
ps2=p1.*[1+2*(u-1)./(2*u+1).*r^3*(h^2-2*X.^2+Y.^2)./(h^2+Y.^2+X.^2).^5/2]
ns=1-ps1./ps2
pcolor(X,Y,ns);shading interp
colorbar
title('中间梯度等值线图')

%对称四极
figure
for x=0:1:50
u=p2/p1
AB=2:4:100
MN=1
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((x-0.5-AM).^2+h^2).^(1/2)
db=((x+0.5+BN).^2+h^2).^(1/2)
rm=((x-0.5).^2+h^2)^(1/2)
rn=((x+0.5).^2+h^2)^(1/2)
cosa=(da.^2+rm^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u2=p3/p1
psA2=p1*[1+2*(u2-1)./(2*u2+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*(u2-1)/(2*u2+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
nsA=1-psA1./psA2
nsB=1-psB1./psB2
ns=(nsA+nsB)/2
plot(log10(AB/2),ns)
hold on
end
grid on
title("对称四级视极化率")
xlabel("AB/2"),ylabel("ns")
figure
x=-50:1:50
ab=10:4:100;
[AB,X]=meshgrid(ab,x)
MN=1
AM=(AB-MN)/2
BN=(AB-MN)/2
AN=AM+MN
BM=BN+MN
K=AM.*AN/MN
da=((X-0.5-AM).^2+h^2).^(1/2)
db=((X+0.5+BN).^2+h^2).^(1/2)
rm=((X-0.5).^2+h^2).^(1/2)
rn=((X+0.5).^2+h^2).^(1/2)
cosa=(da.^2+rm.^2-AM.^2)./(2*da.*rm)
cosb=(da.^2+rn.^2-AN.^2)./(2*da.*rn)
psA1=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
cosc=(db.^2+rm.^2-BM.^2)./(db.*rm.*2)
cosd=(db.^2+rn.^2-BN.^2)./(db.*rn.*2)
psB1=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
p3=p2/(1-n)
u=p3/p1
psA2=p1*[1+2*(u-1)./(2*u+1)*r.^3*(K.*cosa./(da.^2+rm.^2)-K.*cosb./(da.^2+rn.^2))]
psB2=p1*[1+2*(u-1)/(2*u+1)*r^3*(-K.*cosc./(db.^2+rm.^2)+K.*cosd./(db.^2+rn.^2))]
nsA=1-psA1./psA2
nsB=1-psB1./psB2
ns2=(nsA+nsB)/2
pcolor(X,AB,ns2); shading interp
set(gca,'YDir','reverse');
xlabel('X(m)'),ylabel('AB/2(m)');
colorbar
title("拟断面图")
%这是激发极化法
用不同的数据运行一下。
电法勘探matlab代码-OlinAmateurRadio:奥林学院业余无线电俱乐部成员进行的软件定义的无线电实验的资料库
电法勘探matlab代码 名称： 软件无线电-奥林学院业余无线电俱乐部 作者： 斯巴什·班萨尔（Sparsh Bansal） 版本： 1.0 奥林学院业余无线电俱乐部 这是奥林学院业余无线电俱乐部成员进行的软件定义的无线电实验的存储库。 安装 设置和运行ADALM-PLUTO硬件的最简单，最快的方法是： 安装Zadig应用 install Communications Toolbox in MATLAB install Communications Toolbox Support Package for ADALM-PLUTO Radio in MATLAB install Communications Toolbox Support Package for RTL-SDR Radio in MATLAB 探索SDR的其他想法： 一世： 三次SDR ii： GNURadio同伴 iii： GQRX SDR iv： 检验 v： 石灰套房 vi： Pothos流 文献资料 我们为可能有助于某人理解该程序的每一行代码添加了注释。 贡献作品 我们使用了以下信息： 一世： MATLAB文档-Ma
二层电测深
这是成都理工大学地电学二层电测深实验源代码，需要的同学自己下载
基于粒子群优化的双频激电数据联合反演 (2013年)
双频激电法是一种非常有效的探测极化目标体的勘探方法，可以同时观测获取视电阻率和极化率数据。为了提高双频激电观测数据解译的可靠性，利用粒子群优化算法对电阻率数据和极化率进行联合反演。在对异常目标体采用旋转椭球体进行几何近似模拟的基础上，基于Core-Core散射理论实现双频激电法的快速数值模拟计算，为粒子群优化联合反演提供正演基础。在对粒子群优化算法参数进行分析设定的基础上设计了联合反演算法，并采用加入不同程度噪声的模拟双频激电数据对算法进行了实验测试。测试结果表明：粒子群优化联合反演算法能有效实现电阻率与
勘探重力实验matlab,重力场与重力勘探实验指导.ppt
重力场与重力勘探实验指导 实验一、重力仪认识操作 一、实验目的 了解掌握金属弹簧重力仪的测量原理； 学习操作金属弹簧重力仪，了解金属弹簧重力仪的工作范围、操作步骤、性能特点等，为实际野外重力测量实习奠定良好基础。 二、实验内容 金属弹簧重力仪的测量原理； 金属弹簧重力仪野外测量相关情况介绍； 金属弹簧重力仪的认识操作。 三、实验原理 四、实验报告 内容包括实验目的、实验内容、实验原理、实验结果...
视电阻率正演matlab,大地电磁测深一维正演——地电学实验报告分析.doc
大地电磁测深一维正演——地电学实验报告分析实验报告课程名称： 地电学课题名称： 大地电磁层状模型数值模拟实验专 业： 地球物理学姓 名： xx班 级： 06xxxx完成日期： 2016 年 11月26日目录实验名称……………………………………………………3实验目的…………………...
其他

249

社区成员

6,554

社区内容

发帖
与我相关
我的任务
社区描述
其他产品/厂家
社区管理员
  • 其他
加入社区
  • 近7日
  • 近30日
  • 至今

加载中

查看更多榜单
社区公告
暂无公告

试试用AI创作助手写篇文章吧

+ 用AI写文章