[matlab]斜阶梯式长条引起的重力异常计算

画的是《重力学与固体潮》那本书的斜阶梯的重力异常

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 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/489895.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

[Android]使用MVP解决技术债务(翻译)

以下内容为原创,欢迎转载,转载请注明 来自天天博客:http://www.cnblogs.com/tiantianbyconan/p/5892671.html 使用MVP解决技术债务 原文:https://medium.com/picnic-engineering/tackling-technical-debt-with-mvp-67e805ed5103#.…

Science:睡眠剥夺影响大脑思考竟是因为蛋白质罢工了!

来源:生物探索睡眠会影响我们的思维,当我们获得充足的睡眠后,大脑思维会变得清晰;而当我们睡眠不足时,大脑会变得迟钝。那么进入睡眠状态后,大脑又是如何调整以保证睡醒后脑回路清晰的呢?近日&a…

【Matlab】滤波器常用命令

在命令行中输入: designfilt然后就会弹出滤波器的种类,选定后会让你填参数。 参数填好之后点确定,会在命令行窗口生成这个滤波器,复制到编辑器里直接用就好。 %读取音频文件% [x,Fs]audioread(C:\Users\16000\Desktop\testfile.…

element ui input视图没刷新_[Selenium自动化测试实战] 如何在UI自动化测试中加入REST API的操作...

问题当我们描述一个“好的自动化测试用例”时,经常出现标准是:精确。自动化测试用例应该测试一件事,只有一件事。与测试用例无关的应用程序的某个部分中的错误不应导致测试用例失败。独立。自动化测试用例不应该受测试套件中任何其他测试用例…

语句中如何结束本循环进入下一循环_Python3基础语法(八)--控制循环 while...

一、while 简介Python 的循环有 for 和 while 两种&#xff0c;while 为条件控制循环&#xff0c;通过循环控制条件表达式控制循环结束。流程图如下&#xff1a;Python 中 while 语句的格式如下&#xff1a;while <条件表达式>:【语句块】释&#xff1a;当 while 的 <…

关于生命、宇宙和万事万物的42个终极问题

来源&#xff1a;世界科技创新论坛" 我们的宇宙是否稳定&#xff0c;黑洞熵的起源和温度是什么&#xff0c;爱因斯坦的相对论和标准场论总是有效的吗&#xff0c;时空几何中是否存在奇异的性质&#xff0c;化学、应用物理和科技的极限是什么……“在达到完全开悟的道路上&…

【matlab】画图的文字调整大小

hxlabel(x); set(h,Fontsize,14);hylabel(df(x)/dx); set(h,Fontsize,14);htitle(精确解和二阶差分对比); set(h,Fontsize,14);hlegend(精确的一阶偏导,二阶差分得到的偏导); set(h,Fontsize,14);这样就行了

python decimal_【进阶】嫌弃Python慢,试试这几个方法?

(给机器学习算法与Python学习加星标&#xff0c;提升AI技能)选自towardsdatascience&#xff0c;作者&#xff1a;Martin Heinz本文转自机器之心(nearhuman2014)本文将介绍如何提升 Python 程序的效率&#xff0c;让它们运行飞快&#xff01;计时与性能分析在开始优化之前&…

王道8套有变化吗_求求你别再套花艺设计公式了

花艺设计也有公式吗&#xff1f;确实有花艺设计只有公式吗&#xff1f;并不是无论是哪门设计学科&#xff0c;公式这种东西&#xff0c;谈多了是否有种千篇一律的感觉&#xff1f;设计风格相似的花艺师要越来越多&#xff0c;一时间竟然以为都是一个人。就和网红一样&#xff0…

对象构造函数的原型图

对象的定义其实很广泛,万物皆为对象,我们创建对象一般都是用构造函数来创建的,这里我们来说说构造函数创建对象的原型图把. 这个问题有点抽象,举个例子来说,方便一点: 这是我们构造函数,这里我们要结合一张 图来说明就更清楚了,这里我们就用一个实例p1好了,其他两个就不用了. 这…

【服务器】sbatch 提交作业脚本

在已经写好的脚本的#!bin/bash下面加上&#xff1a; #SBATCH --get-user-env #SBATCH --mail-typeend #SBATCH -J simple_module #SBATCH --nodes4 #SBATCH --ntasks-per-node24然后退出即可。 运行时使用 sbatch ./name_of_bash即可 其中 #SBATCH -J simple_module指定了提…

站在AI与神经科学交叉点上的强化学习

来源&#xff1a; 混沌巡洋舰一&#xff0c;强化学习概述让机器来决策&#xff0c;首先体现在如何模仿人类的决策。对于决策这个问题&#xff0c; 对于人类是困难的&#xff0c; 对于机器就更难。而强化学习&#xff0c; 就是一套如何学习决策的方法论。强化学习最初的体现就是…

澜起科技云计算服务器_服务器严重缺货!云应用大爆发!云计算正强势起爆(附龙头)...

催化因素&#xff1a;这两天全国上千万企业、近两亿人开启在家办公模式。阿里、华为、腾讯等各大网络办公平台纷纷告急。对云服务的需求大增也让服务器生产企业开足马力&#xff0c;春节假期里&#xff0c;山东浪潮集团就接到了1500台服务器的订单。目前&#xff0c;多家软件服…

【matlab】零相位延迟滤波器

使用designfilt命令自动生成滤波器&#xff0c; 然后my_filt‘生成的滤波器代码’&#xff0c; 最后filtfilt(my_filt,signal)

博客迁移到简书

由于博客园对markdown的支持不是很好&#xff0c;hexo访问速度比较慢&#xff0c;后续博客迁移到简书&#xff0c;简书地址 http://www.jianshu.com/users/7966a93b66ce/latest_articles 。 转载于:https://www.cnblogs.com/actionke/p/5899368.html

车险赔付率分析报告_车险有变!价格…

各位车友请注意&#xff01;《商业车险综合示范条款(2020版征求意见稿)》于近日发布向社会公开征求意见从修订版条款的内容来看大幅删减了责任免除项目扩展了保险责任在最大化让利于消费者的同时努力提升消费者体验那么&#xff0c;此次修订版有哪些具体的亮点呢&#xff1f;一…

“众声喧哗”中的VR,谁来买单?

来源&#xff1a;VR每日必看未来智能实验室是人工智能学家与科学院相关机构联合成立的人工智能&#xff0c;互联网和脑科学交叉研究机构。未来智能实验室的主要工作包括&#xff1a;建立AI智能系统智商评测体系&#xff0c;开展世界人工智能智商评测&#xff1b;开展互联网&…

.npy文件_Numpy库使用入门(六)文件的存取

ERNIE&#xff1a;BERT&#xff0c;你看到我的npy了吗&#xff0c;我记得我放在这个文件夹里的呀(」&#xff1e;&#xff1c;)」BERT&#xff1a;就你还有npy&#xff1f;我还单着呢&#xffe3;へ&#xffe3;ERNIE&#xff1a;你想什么呢&#xff1f;我指的是numpy储存数据的…

【linux】设置镜像源

首先找到默认的源 $cd /etc/apt $cp source.list source.list_default然后替换掉source.list文件为国内源。 清华源

Java易混小知识——equals方法和==的区别

一、equals方法和的区别 1.equals是String对象的方法&#xff0c;可以通过".“调用。 2. 是一个运算符。 二、常用的比较用法 1、基本数据类型比较。 equals和都比较两个数值 是否相等。相等即为true,不相等则为false。 2、引用对象的比较。 equals和都比较栈内存中的地址…