qq_37486958的博客

地球重力学需要我们计算立方体引起的重力异常，公式见《重力学与固体潮》。这个程序取的z方向是竖直向下的，也就是说地面向下为正，地面向上为负%定义一个立方体函数%function [gravity]=draw_square(a,b,c,x0,y0,H,ph,z)%长方体模型参数说明%%a=2000;%长%b=200;%宽%c=100;%高%%质心坐标x0,y0,z0 %H=1000立方体深...

画的是《重力学与固体潮》那本书的斜阶梯的重力异常function [gravity]=draw_gravity_dixing(h1,h2,alpha,xigma,x0,z,direction)%h1=100;%上表面%h2=200;%下表面%x0=0;%移动原点位置，向右为正%xigma=3*10^3;%剩余密度pi=3.1415926535%alpha=alpha*pi;%倾角...