matlab画水平圆柱体的重力异常、水平导数、垂向导数、垂向二阶导数

yeahamen

已于 2023-06-27 10:56:12 修改

阅读量241

点赞数

文章标签： java python 开发语言

于 2023-05-20 21:23:55 首次发布

本文链接：https://blog.csdn.net/amen_yanwang/article/details/130786079

版权

G =  6.67.*10.^(-11)
labota = 200%kg/m；线剩余密度为正
D = 10%m
x = linspace(-100,100,200)

%画的曲线是什么都在title里
figure(1)
delta_g = (2*G*labota*D)./((x.^2 + D.^2))
plot(x,delta_g,'b', 'Linewidth', 2, 'MarkerSize', 8)
xlabel('x')
ylabel('delta-g')
text(D,(2*G*labota*D)./((D.^2 + D.^2)),'\leftarrow X1/2,D')
text(1.414*D,(2*G*labota*D)./(((1.412*D).^2 + D.^2)),'\leftarrow X1/3,1414D')
title('I类 delta-g of cylinder')

figure(2)
delta_g = (-4*G*labota*D*x)./((x.^2 + D.^2).^(2))
plot(x,delta_g,'b', 'Linewidth', 2, 'MarkerSize', 8)
xlabel('x')
ylabel('H-gradient of delta-g')
title('II类 delta-g(H-gradient) of cylinder')

figure(3)
delta_g = (2*G*labota*(D.^2-x.^2))./((x.^2 + D.^2).^(2))
plot(x,delta_g,'b','Linewidth', 2, 'MarkerSize', 8)
hold on;
y = 0*x
plot(x,y,'g', 'Linewidth', 1)
xlabel('x')
ylabel('V-gradient of delta-g')
title('III类 delta-g(V-gradient) of cylinder')

figure(4)
delta_g = (4*G*labota*D*(D.^2-3*x.^2))./((x.^2 + D.^2).^(4))
plot(x,delta_g,'b','Linewidth', 2, 'MarkerSize', 8)
hold on;
y = 0*x
plot(x,y,'g', 'Linewidth', 1)
xlabel('x')
ylabel('VV-gradient of delta-g')
title('III类 delta-g(VV-gradient) of cylinder')