零购物网站怎么建设,北京网站建设公司哪家最好,中国公信网企业查询,成都中风险地区地球重力学需要我们计算立方体引起的重力异常#xff0c;公式见《重力学与固体潮》。 这个程序取的z方向是竖直向下的#xff0c;也就是说地面向下为正#xff0c;地面向上为负
%定义一个立方体函数%
function [gravity]draw_square(a,b,c,x0,y0,H,ph,z)
%长方体模型参数说…地球重力学需要我们计算立方体引起的重力异常公式见《重力学与固体潮》。 这个程序取的z方向是竖直向下的也就是说地面向下为正地面向上为负
%定义一个立方体函数%
function [gravity]draw_square(a,b,c,x0,y0,H,ph,z)
%长方体模型参数说明%
%a2000;%长%b200;%宽%c100;%高%
%质心坐标x0,y0,z0 %H1000立方体深埋深度;
%质心埋深H
%ph2*10^3;%剩余密度
%z是测点的z坐标比如如果在地面上就取0%
%采样区间%
x(-40:2:40);
y(-40:2:40);%常数%
G6.67e-11;%计算异常%
[x1,y1]meshgrid(x,y); %生成计算用的网格线
r(x0-x1).^2(y0-y1).^2(H-z).^2;
gravity-G*ph.*(a.*log(rb)b.*log(ra)-c.*atan(a*b./(r*c)))*10^5;%单位mGal%画图%
%矩阵要和前面的那个矩阵取一样的%
x(-40:2:40);
y(-40:2:40);[x1,y1]meshgrid(x,y);%生成画图用的矩阵
figure(1)%图1
mesh(x1,y1,gravity);%三维
colorbar;
xlabel(x);
ylabel(y);
title(立方体异常);
