地球重力學需要我們計算立方體引起的重力異常，公式見《重力學與固體潮》。
這個程序取的z方向是豎直向下的，也就是說地面向下為正，地面向上為負

%定義一個立方體函數% function [gravity]=draw_square(a,b,c,x0,y0,H,ph,z) %長方體模型參數說明% %a=2000;%長%b=200;%寬%c=100;%高% %質心坐標x0,y0,z0 %H=1000立方體深埋深度; %質心埋深H %ph=2*10^3;%剩余密度 %z是測點的z坐標，比如如果在地面上就取0% %采樣區間% x=(-40:2:40); y=(-40:2:40);%常數% G=6.67e-11;%計算異常% [x1,y1]=meshgrid(x,y); %生成計算用的網格線 r=(x0-x1).^2+(y0-y1).^2+(H-z).^2; gravity=-G*ph.*(a.*log(r+b)+b.*log(r+a)-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('立方體異常');

總結

以上是生活随笔為你收集整理的立方体引起的引力异常计算＆画图的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。