我已经写了'r'(Radius)值的一般轨道方程的代码,并试图绘制x=r*cos(theta)和{}的图形,但是我的代码给出了相同的“x”值来表示不同值的“r”,所以我用一条直线而不是椭圆来结束enter code here。在from numpy import *

import matplotlib.pyplot as plt

import pprint

L = 9.11*10**38 #L = angular momentum

m = 3.28*10**23 #m = mass of mercury

M = 1.99*10**30 #M = mass of sun

a = 5.8*10**7 #a = semi-major axis

G = 6.674*10**-11 #G = Gravitationl constant

k = G*M*m

E = -k/(2*a) #E = energy

p = L**2/(m*k)

c = 1 + (2*E*L**2)/m*k**2

e = sqrt(-c) #e = eccentricity

def fx(x):

r = p/(1 + e*cos(x))

return r

n = 1000

phi =linspace(0,2*pi,n)

radius = zeros([n])

theta = zeros([n])

x = zeros([n])

y = zeros([n])

for i in range(0,n):

radius[i] = fx(phi[i])

theta[i] = 180*phi[i]/pi

for i in range(0,n):

x[i] = radius[i]*cos(phi[i])

for i in range(0,n):

y[i] = radius[i]*sin(phi[i])

print('r =',radius)

print('x =',x)

print('y =',y)

plt.plot(x,y)

plt.show()

这是我的数学书截图:

python行星名称输出_用python绘制行星轨道相关推荐

  1. python简单图形输出_基于 Python Matplotlib 模块的高质量图形输出

    Matplotlib 是一个由 John Hunter 等开发的,用以绘制二维图形的 Python 模块.它利用了 Python 下的数值计算模块 Numeric 及 Numarray,克隆了许多 M ...

  2. python 表格格式输出_利用python对excel中一列的时间数据更改格式操作

    问题场景:需要将下列的交期一列的数据格式更改成2019/05/10 存货编码 尺寸 数量 交期 0 K10Y0190000X B140 200 2019-05-10 00:00:00 1 K10Y01 ...

  3. python咖啡名称提取_(转)Python中文全攻略

    插入数据库时,总是空,我做了如何操作,就好了.数据是采集过来的,程序是gbk编码person_sql="insert into analyst(education_id,alys_name, ...

  4. python中隔行输出_在Python中的分隔行上打印列表元素

    我试图打印出Python路径文件夹使用: import sys print sys.path 输出是这样的: >>> print sys.path ['.', '/usr/bin', ...

  5. python时间间隔标准化输出_利用Python将时间或时间间隔转为ISO 8601格式方法示例...

    前言 大家都知道,Python自带的datetime库提供了将datetime转为ISO 8610格式的函数,但是对于时间间隔(inteval)并没有提供转换的函数,下面我们动手写一个. 下面话不多说 ...

  6. python按列输出_在python中创建漂亮的列输出

    我试图在python中创建一个很好的列列表,用于我创建的命令行管理工具. 基本上,我想要一个列表,如: [['a', 'b', 'c'], ['aaaaaaaaaa', 'b', 'c'], ['a' ...

  7. python turtle画滑稽_使用python的turtle函数绘制一个滑稽表情的方法

    Turtle库是Python语言中一个很流行的绘制图像的函数库,想象一个小乌龟,在一个横轴为x.纵轴为y的坐标系原点,(0,0)位置开始,它根据一组函数指令的控制,在这个平面坐标系中移动,从而在它爬行 ...

  8. python 标准错误输出_过程的实时标准输出/错误捕获

    我正在开发一个C应用程序,我需要运行一个外部控制台进程(例如python脚本)并实时接收脚本的输出数据.python脚本如下所示:import time while 1: print("Hi ...

  9. python常用语法和示例_使用Python中的示例进行输入和输出操作

    python常用语法和示例 A Program needs to interact with the user to accomplish the desired task; this is done ...

最新文章

  1. windows域设计best practice
  2. 【转载】DotNet 资源大全中文版
  3. JavaScript中的剪贴板(clipboardData)
  4. POJ 1321 棋盘问题(DFS 状压DP)
  5. 解锁一个新技巧,让你舒服到不想起床
  6. sql 替换text字段中的指定字符
  7. 单关闭功能_广告弹窗关闭
  8. libuv 原理_Nodejs的运行原理-libuv篇
  9. 让敏捷交付优秀的软件
  10. 为什么开源应该是云原生环境的首选
  11. XSS学习笔记(一)
  12. Observer设计模式中-委托事件-应用在消息在窗体上显示
  13. 清空缓存的命令_布隆过滤器应用——解决Redis缓存穿透问题
  14. element-ui中table 去掉表头、去掉边框线、去掉鼠标悬停背景颜色
  15. alpha-beta剪枝 个人理解
  16. ArcGis Server10.2 授权文件教程
  17. Linux环境变量PSI指什么,psi是什么单位,pSI指标应用原则
  18. 启发式算法与机器学习的区别_使用强化学习训练受启发的四足机器人
  19. 安装向日葵远程后电脑亮度和夜间模式不能使用了——故障笔记
  20. 客厅窗帘最实用的色彩组合搭配,大气稳重-窗帘十大品牌江南爱

热门文章

  1. VC++ 获取窗体句柄,并发送键盘消息(这种方法也可以打开某些应用程序)
  2. ACAD DWG to PDF Converter 9.8.2.4版本更新啦
  3. bilibili_Linux4_Linux常用命令
  4. mac apache 重启
  5. 当SAP遇上马云,智慧企业加速中国数字化转型
  6. 炒菜机器人放食材的顺序_机器人能否消灭厨师这个职业?
  7. (实战)Node.js 实现抢票小工具短信通知提醒
  8. linux nvme文件系统,Intel NVMe驱动器扇区大小不是4096的xfs文件系统的性能下降
  9. 实现表格中的el-switch 开关切换
  10. 第三十八章 短语动词