绘制大西洋中脊剖面轮廓#abs求绝对值#append

本文介绍如何利用数学公式描述大西洋中脊剖面的轮廓变化,给出具体参数D=2.5km, C=0.11km/Myr, u=50km/Myr,通过Python实现海底轮廓线的绘制,展示了x坐标与深度、海底减薄速率的关系。
摘要由CSDN通过智能技术生成

 大西洋中脊横截剖面轮廓形状如上右图红色曲线所示。对于某处剖面,以中脊对称轴向上方向为 y 轴、以海平面为 x 轴建立直角坐标系,则轮廓线上某点 A 的横坐标 x 表示 A 离对称轴的距离(单位为千米),而纵坐标 y 表示 A 点的海拔(单位为千米,一般为负数),x 和 y 存在如下关系:

  其中,∣x∣ 表示求 x 的绝对值,D 为对称轴处离海面的距离(单位为 km)、C 为给定的常量、u 为海底减薄速率(单位为 km/Myr),假设 D=2.5km、C=0.11、u=50km/Myr,本关任务是绘制这种情况下的剖面轮廓线。

import matplotlib.pyplot as plt
import math
########## Begin ##########
D,C,U=2.5,0.11,50
X=[]
Y=[]
d,delta=1200,200
for x in range(-d,d+1,delta):
    X.append(x)
    y=-(D+C*(10*abs(x)/U)**0.5)
    Y.append(y)
########## End ##########
#以下代码用于打印结果
for i in range(len(X)):
    print('%5s\t%.2f' % (int(X[i]),Y[i]))  
#以下代码用于保存图片
plt.plot(x,y,'r-')
plt.axis([-d,d,-4.2,0])
plt.savefig( 'src/step7/student/pic.png' )
plt.close()
  • 2
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

省心省事省时

你的鼓励是我最大的动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值