畫的是《重力學與固體潮》那本書的斜階梯的重力異常

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;%傾角%采樣區間% x=(-5000:50:5000); y=(-5000:50:5000); %z=-100;%測點高度 G=6.67e-11; [x2,y2]=meshgrid(x,y); %生成網線節點矩陣 if direction==1gravity=G*xigma*(pi*(h2-h1)+2*h2.*atan(((x2-x0)+h2*cot(alpha))/h2)-2*h1.*atan(((x2-x0)+h1*cot(alpha))/h1)+(x2-x0)*sin(alpha)*sin(alpha).*log(((h2+(x2-x0)*sin(alpha)*cos(alpha)).^2+(x2-x0).^2*sin(alpha)^4)./(((h1+(x2-x0)*sin(alpha)*cos(alpha)).^2+(x2-x0).^2*sin(alpha)^4)))-2.*(x2-x0)*sin(alpha)*cos(alpha).*atan(((x2-x0)*(h2-h1)*sin(alpha)*sin(alpha))./((x2-x0).^2*sin(alpha)^2+(h1+h2).*(x2-x0)*sin(alpha)*cos(alpha)+h1*h2))); end if direction==2gravity=G*xigma*(pi*(h2-h1)+2*h2.*atan(((x0-x2)+h2*cot(alpha))/h2)-2*h1.*atan(((x0-x2)+h1*cot(alpha))/h1)+(x0-x2)*sin(alpha)*sin(alpha).*log(((h2+(x0-x2)*sin(alpha)*cos(alpha)).^2+(x0-x2).^2*sin(alpha)^4)./(((h1+(x0-x2)*sin(alpha)*cos(alpha)).^2+(x0-x2).^2*sin(alpha)^4)))-2.*(x0-x2)*sin(alpha)*cos(alpha).*atan(((x0-x2)*(h2-h1)*sin(alpha)*sin(alpha))./((x0-x2).^2*sin(alpha)^2+(h1+h2).*(x0-x2)*sin(alpha)*cos(alpha)+h1*h2))); end

總結

以上是生活随笔為你收集整理的[matlab]斜阶梯式长条引起的重力异常计算的全部內容，希望文章能夠幫你解決所遇到的問題。

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