clear,clc,clf
rho1=100;
h1=10;a=0.011;
k=[-0.9,-0.7,-0.5,-0.3,0,0.3,0.5,0.7,0.9];
l=2.1:0.2:15;
for x=1:length(k)
rho2(x)=(k(x)+1)/(1-k(x))*rho1;
for y=1:length(l)
if l(y)<=h1
s=0;
for z=1:10000
s=s+k(x)^z*log((2*z*h1+l(y)+sqrt(a^2+(2*z*h1+l(y))^2))/((2*z*h1-l(y)+sqrt(a^2+(2*z*h1-l(y))^2))));
end
R(x,y)=rho1/(2*pi*l(y))*(log((l(y)+sqrt(a^2+l(y)^2))/a)+s);
end
if l(y)>h1
s=0;
for z=1:10000
s=s+k(x)^z*log((2*z*h1+l(y)+sqrt(a^2+(2*z*h1+h1)^2))/((2*z*h1-h1+sqrt(a^2+(2*z*h1+h1)^2))));
end
R(x,y)=rho1*rho2(x)/(2*pi*(rho2(x)*h1+rho1*(l(y)-h1)))*(log((l(y)+sqrt(a^2+l(y)^2))/a)+s);
end
end
end
for i=1:length(k)
plot(l,R(i,:));
hold on
end
新建文件夹.rar_土壤电阻率计算_接地电阻_接地计算
版权申诉
5星 · 超过95%的资源 166 浏览量
2022-07-13
21:28:12
上传
评论 1
收藏 1KB RAR 举报
JaniceLu
- 粉丝: 79
- 资源: 1万+