无限长水平圆柱体重力异常MATLAB
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/19 11:24:46
无限长水平圆柱体重力异常MATLAB
无限长水平圆柱体重力异常MATLAB
无限长水平圆柱体重力异常MATLAB
你也来了
于是我回忆起,有一次,在天上……
我就捕捉住它
在阳光中跳跃,在风的嬉戏和絮语中
你们把秘密保守得如此小心!
已记不清后天该是什么一个子哈哈
clear;
x=-10:0.1:10;
[x]=meshgrid(x);
r=1;
h=3;p=0.8;
p1=pi*r^2*p;
G=6.67;
gg=2*G*p1*h./(x.^2+h^2);
g=gg*10^-2;
Vxz=-4*G*p1*x.*h./(x.^2+h.^2).^2;
Vzz=2*G*p1*(h^...
全部展开
clear;
x=-10:0.1:10;
[x]=meshgrid(x);
r=1;
h=3;p=0.8;
p1=pi*r^2*p;
G=6.67;
gg=2*G*p1*h./(x.^2+h^2);
g=gg*10^-2;
Vxz=-4*G*p1*x.*h./(x.^2+h.^2).^2;
Vzz=2*G*p1*(h^2-x.^2)./(x.^2+h.^2).^2;
Vzzz=4*G*p1*h.*(-3*x.^2+h.^2)./(x.^2+h.^2).^3;
figure(1);
[c,h]=contour(g,20);title('Gravity contour');
figure(2);
subplot(221);plot(x,g(),'b');title('profile anomaly');
subplot(222);plot(x,Vxz(),'b');legend('Vxz');
subplot(223);plot(x,Vzz());legend('Vzz');
subplot(224);plot(x,Vzzz());legend('Vzzz');
figure(3);
r1=3;
r2=6;
g1=2*G*pi*r1^2*p*h./(x.^2+h^2);
g2=2*G*pi*r2^2*p*h./(x.^2+h^2);
plotyy(x(11,:),g1(11,:),x(11,:),g2(11,:));
title('the gravity of the different radius');
收起