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;%倾角 ...