R语言计算地球固体潮汐数据

最新推荐文章于 2024-08-07 16:30:29 发布
最新推荐文章于 2024-08-07 16:30:29 发布
阅读量931 收藏 1
点赞数
文章标签： r语言 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/haidejintou2013/article/details/121254371
版权
本文记录使用R语言结合Fortran编写的ETERNA 3.4库来计算地球固体潮汐的过程。首先通过Anaconda安装环境，接着设置RStudio并配置包镜像，最后安装并运行earthtide库进行计算。
摘要由CSDN通过智能技术生成

R语言计算地球固体潮汐数据

利用R语言基于Fortan ETERNA 3.4计算地球固体潮汐

本文章算是一个笔记，给出利用jkennel在github上发行的earthtide库实现对于地球固体潮汐的计算，步骤稍显粗略。

  1. 首先安装Anaconda环境，详细步骤可以参考【Anaconda教程01】怎么安装Anaconda3
    Anaconda控制面板
  2. 安装Rstudio环境【Anaconda教程01】怎么安装Anaconda3
    在这里插入图片描述
  3. 配置Rstudio的包镜像在这里插入图片描述

Tools–>Global options–>Package–>CRAN mirror
修改至合肥

  1. 安装earthtide库

最低0.47元/天 解锁文章
價裝科研
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Tide2_固体潮_理论固体潮计算_
10-01
标题中的"Tide2_固体潮_理论固体潮计算_"表明这是一个关于固体潮的计算程序，具体是用MATLAB编程语言实现的。固体潮是地球表面上的一种周期性位移现象，主要由月球和太阳的引力引起。在海洋学、地球物理学和天文学中...
第三章、TIDE
weixin_44836552的博客
09-09 1234
针对TIBBO提供的整合开发工具TIDE（Tibbo Integrated DevelopmentEnvironment）进行介绍，教导读者如何操作TIDE平台，内容包含TIDE开发环境介绍、功能介绍、相关使用注意事项，例如如何建立项目、执行项目与除错方法等等，让读者可以对最重要的开发环境具备有基本的操作能力。 TIDE（Tibbo Integrated Development Environment）是Tibbo提供的整合式程序开发平台，与市面上常见的IDE平台具备相同主要功能，包括建立项目、撰写程序、执
【目标检测】 IterDet模型简介、TIDE评价指标介绍
weixin_47716911的博客
10-10 2632
来源 | 极链AI云（性价比最高的共享GPU算力平台，双十活动进行中 10.9-10.11，新人注册可领取198元大礼包，价值100小时的免费使用时长，领取地址：https://cloud.videojj.com/） 作者 | 小仇极链AI云技术原创奖励计划 原文地址 | [目标检测] IterDet（Iterative Scheme for Object Detection in Crowded Environments）模型简介 [目标检测-评价指标] TIDE（A General Tool..
R_ggplot2地理信息可视化_史上最全（一）
R语言中文社区
01-18 1万+
作者：李誉辉  四川大学在读研究生1.简介此次教程是关于使用ggplot2包进行地理信息可视化的，笔者将从头到尾，条分缕析的梳理ggplot2地图可视化的知识点。让读者对...
固体潮 matlab,基于MATLAB改进设计的固体潮改正程序及其应用
weixin_34816557的博客
03-18 473
科技的不断进步和发展,新的方法、软件都通算,得到每一个重力测点的固体潮改正值,并将其与过计算机不断的实现,MATLAB　以独特的优势广泛一定区域使用同一个固体潮改正经纬度计算所得到的运用于各个领域,尤其是数据的计算方面具有计算速每一点的固体潮改正值进行比较,结果采用每一点的度快,编程简单,受到广大计算机爱好者的青睐。重经纬度计算该点的固体潮改正值比一定区域使同一个力勘探作为一种物理勘探方法,越来越...
【Rust日报】2020-07-20 boids算法, tide, popol, Calcite
Rust语言学习交流
07-20 469
文章类Rust 模块系统的超清晰解释Rust 的模块系统可能对于新人来说有些困惑,这篇文章通过几个循序渐进的小例子,让你快速且清晰的对 Rust 模块系统有一个基本的认识.http://...
长度与固体密度测量实验报告带数据
12-07
长度与固体密度测量实验报告带数据 本实验报告旨在介绍长度与固体密度测量实验的过程和原理，包括游标卡尺、螺旋测微器和读数显微镜的使用和原理。 一、实验目的： 1. 了解游标卡尺、螺旋测微器的构造，掌握它们...
fem-matlab.rar_FEM_fem matlab_matlab 力学_计算固体力学
07-14
在固体力学问题中，这些参数直接影响到计算结果的准确性。例如，材料属性定义了物体的弹性模量和泊松比，边界条件则规定了物体在某些边缘或节点上的位移或应力约束。 通过这个压缩包，学习者可以深入了解如何使用...
固体物理重点计算题详细解答.pdf
06-30
在这个文档中，我们主要探讨了固体物理中关于晶体结构的一些重要计算问题。 首先，文档涉及的是晶体的正格子和倒格子之间的关系。在固体物理中，正格子是指实际的晶格结构，而倒格子是正格子的傅立叶变换，用于描述...
通过模拟退火优化空间样本_R语言_代码_下载
06-21
R语言是一种广泛应用于统计计算和图形绘制的编程语言，具有丰富的数据处理和分析库。在空间样本配置优化中，可能涉及到地理信息系统（GIS）和空间统计等领域。模拟退火算法与R语言的结合，使得在处理复杂的地理空间...
R语言绘制河流图
weifanbio的博客
06-15 1979
ggplot2绘制Streamgraph # install.packages("ggstream") 安装包 library(ggstream) library(ggplot2) colors <- c("#6181BD4E","#F348004E","#64A10E4E","#9300264E","#464E044E") #绘图 ggplot(blockbusters, aes(x = year, y = box_office, fill = genre)) + geom_stream()
利用R语言绘制世界航班路线图
R语言中文社区
03-06 4776
作者简介Introductiontaoyan：伪码农，R语言爱好者，爱开源。个人博客: https://ytlogos.github.io/公众号：生信大讲堂往期回顾R语言可视化学习笔记之相关矩阵可视化包ggcorrplotR语言学习笔记之相关性矩阵分析及其可视化ggplot2学习笔记系列之利用ggplot2绘制误差棒及显著性标记ggplot2学习笔记系列之主题（theme）设置用circlize
R各个包里面的数据集列表
四纪启悟 · 数据科学
07-07 1万+
Package  Item  Title  csv  doc datasets AirPassengers Monthly Airline Passenger Numbers 1949-1960 CSV DOC datasets BJsales Sales Data with Leading Indicator CSV DOC
R语言气象相关分析方法
weixin_39634034的博客
02-03 2960
R语言气象相关分析方法（包括二进制文件R语言读取） ######读取热带印度洋海温场资料################# read.filename &lt;- file("D:\\99\\NCEP_IOSST_30y_Wt.dat", "rb") bindata &lt;- readBin(read.filename, numeric(),16*7*30,size = 4) IOSST&lt;...
LOD的计算
Code for Coffee
02-06 1万+
在分析化学中，检测限（Limit of detection，LOD）通常与分析灵敏度相混淆（sensitivity） LOD指的是一种物质能被检测的最低含量，可以通过规定的置信水平（99%）与不存在该物质的样本进行区分 评估检测限需要：1. 空白样本的均值 2. 空白样本的SD 3. 线性的斜率（sensitivity） 4. 自定义的置信度 举个例子，某测试方法对于某样品的线性如下： y = kx + b 其中，y是被测得的信号值 k是灵敏度sensitivity x是需要被测量的浓度 b
潮汐计算php源码,科学网—潮汐再科普与再计算 - 岳东晓的博文
weixin_29530803的博客
03-28 546
我之前写过几篇潮汐的科普，并且计算了潮汐高度随角度的变化。但这个计算忽略了一项。因此，我重新计算一下。顺便再从概念上讲讲。古代中国人就知道潮汐主要是由月球引起的。牛顿在发现其运动定律及万有引力之后，对月球和太阳引起的潮汐进行了计算，但据说有错误(由于牛顿的《自然哲学的数学原理》一书只用了简单的几何和代数，其计算非常繁琐，现在很难看清)。在电影《Interstellar》中，浅水出现非常高的潮汐，应...
潮汐计算php源码,潮汐选股公式
weixin_35215045的博客
03-28 643
潮汐选股公式发布时间：2019-05-26选股公式源码C1:= AMO>40000000 AND(H -MAX(C,O))>(MIN(C,O)-L)*2 ANDH90,REF(EXIST(V>=REF(V,1)*2 AND C>REF(C,1),12),1),1) ;C2:=IF(BARSCOUNT(C)>55,MA(C,10)>MA(C,21),1)ANDIF...
R语言:使用REmap绘制超炫酷的地图
热门推荐
王亨的博客
11-10 5万+
REmap这个包是通过调用百度地图API的一个程序包，与Google的ggmap包相比，不用翻墙就可以使用。其函数主要有：remap()，remapB()、remapC()、remapH()四个。该程序包目前要github网站上，因此我们需要从github上面下载安装。 install.packages("devtools")#如果安装了就不需要这步 library(devtools) inst
【Python】requests的response.text 和 urllib.request 的 response.read()的区别
最新发布
Meggie35的博客
08-07 481
两种获取网络响应方式的不同
c语言实现地球固体潮理论值计算的源代码
05-15
地球固体潮理论值计算的源代码是一个比较复杂的程序，需要涉及到许多数学公式和计算方法。以下是一个可能的C语言实现地球固体潮理论值计算的程序： #include <stdio.h> #include <math.h> #define PI 3.14159265358979323846 // 定义圆周率常量 double double_factorial(int n) { double result = 1.0; while (n > 1) { result *= n; n -= 2; } return result; } double binomial_coefficient(int n, int k) { if (k > n) { return 0.0; } double result = 1.0; for (int i = 0; i < k; i++) { result *= (double)(n - i) / (double)(k - i); } return result; } double legendre_polynomial(double x, int n) { if (n == 0) { return 1.0; } else if (n == 1) { return x; } else { double Pn_minus_2 = 1.0; double Pn_minus_1 = x; double Pn = 0.0; for (int i = 2; i <= n; i++) { Pn = ((2.0 * i - 1.0) * x * Pn_minus_1 - (i - 1) * Pn_minus_2) / i; Pn_minus_2 = Pn_minus_1; Pn_minus_1 = Pn; } return Pn; } } double solid_tide_amplitude( double M, double omega, double phi, int l, int m, double G, double R, double k, double alpha) { double a = pow(R, l) / pow((double)l, 2.0 * l + 1.0); double b = double_factorial(2 * l - 1) * (1 - k) / pow(k, l); double c = 3.0 * pow(alpha, l + 1.0) * pow(G, l + 2.0) * M; double d = pow(omega, 2.0) / pow(G, l) / double_factorial(l); double e = binomial_coefficient(2 * l, l) * pow(sin(phi), 2.0 * l); double f = legendre_polynomial(cos(phi), l); double g = cos((double)m * phi); return a * b * c * d * e * f * g; } int main() { double M = 5.97e24; // 地球质量 double omega = 7.292e-5; // 地球自转角速度 double phi = 30.0 / 180.0 * PI; // 纬度 int l = 2; // 应力函数次数 int m = 0; // 应力函数方向 double G = 6.6743e-11; // 万有引力常数 double R = 6371e3; // 地球半径 double k = 0.3; // 岩石材料的泊松比 double alpha = 0.8; // 计算中的参数 double amplitude = solid_tide_amplitude(M, omega, phi, l, m, G, R, k, alpha); printf("地球固体潮理论值的振幅为%e m\n", amplitude); return 0; } 这个程序中包含了一些常用的数学函数，如双重阶乘函数、二项式系数、勒让德多项式等。其中，solid_tide_amplitude()函数是计算地球固体潮理论值振幅的核心函数，它根据地球的各种属性和当前位置的纬度、应力函数次数和方向等参数，计算出相应的地球固体潮理论值振幅。最后，main()函数中调用这个核心函数，并输出结果。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值