Python编程使用matplotlib绘制动态圆锥曲线示例

 更新时间:2021年10月19日 11:56:44   作者:微小冷  
这篇文章主要介绍了Python使用matplotlib绘制动态的圆锥曲线示例实现代码,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步

作为让高中生心脏骤停的四个字,对于高考之后的人来说可谓刻骨铭心,所以定义不再赘述,直接撸图,其标准方程分别为

在这里插入图片描述

在Python中,绘制动图需要用到matplotlib中的animation包,其调用方法以及接下来要用到的参数为

ani = animation.FuncAnimation(fig, func, frames, interval)

其中fig为绘图窗口,func为绘图函数,其返回值为图像,frames为迭代参数,如果为整型的话,其迭代参数则为range(frames)

椭圆

为了绘图方便,椭圆的参数方程为

在这里插入图片描述

在这里插入图片描述

代码为:

# 这三个包在后面的程序中不再复述
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
a,b,c = 5,3,4
fig = plt.figure(figsize=(12,9))
ax = fig.add_subplot(autoscale_on=False, 
    xlim=(-a,a),ylim=(-b,b))
ax.grid()
line, = ax.plot([],[],'o-',lw=2)
trace, = ax.plot([],[],'-', lw=1)
theta_text = ax.text(0.02,0.85,'',transform=ax.transAxes)
textTemplate = '''theta = %.1f°\n
lenL = %.1f, lenR = %.1f\n
lenL+lenR = %.1f'''
xs,ys = [], []
def animate(i):
    if(i==0):
        xs.clear()
        ys.clear()
    theta = i*0.04
    x = a*np.cos(theta)
    y = b*np.sin(theta)
    xs.append(x)
    ys.append(y)
    line.set_data([-c,x,c], [0,y,0])
    trace.set_data(xs,ys)
    lenL = np.sqrt((x+c)**2+y**2)
    lenR = np.sqrt((x-c)**2+y**2)
    theta_text.set_text(textTemplate % 
        (180*theta/np.pi, lenL, lenR, lenL+lenR))
    return line, trace, theta_text
ani = animation.FuncAnimation(fig, animate, 157, 
    interval=5, blit=True)
ani.save("ellipse.gif")
plt.show()

双曲线

双曲线的参数方程为

在这里插入图片描述

设 a = 4 , b = 3 , c = 5 则代码如下

a,b,c = 4,3,5
fig = plt.figure(figsize=(12,9))
ax = fig.add_subplot(autoscale_on=False, 
    xlim=(-c,16),ylim=(-12,12))
ax.grid()
line, = ax.plot([],[],'o-',lw=2)
trace, = ax.plot([],[],'-', lw=1)
theta_text = ax.text(0.01,0.85,'',
    transform=ax.transAxes)
textTemplate = '''t = %.1f\n
lenL = %.1f, lenR = %.1f\n
lenL-lenR = %.1f'''
xs,ys = [],[]
def animate(t):
    if(t==-3):
        xs.clear()
        ys.clear()
    x = a*np.cosh(t)
    y = b*np.sinh(t)
    xs.append(x)
    ys.append(y)
    line.set_data([-c,x,c], [0,y,0])
    trace.set_data(xs,ys)
    lenL = np.sqrt((x+c)**2+y**2)
    lenR = np.sqrt((x-c)**2+y**2)
    theta_text.set_text(textTemplate % 
        (t, lenL, lenL, lenL-lenR))
    return line, trace, theta_text
frames = np.arange(-3,3,0.05)
ani = animation.FuncAnimation(fig, animate, 
    frames, interval=5, blit=True)
ani.save("hyperbola.gif")

plt.show()

在这里插入图片描述

抛物线

在这里插入图片描述

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
a,b,c = 4,3,5
p = 1
fig = plt.figure(figsize=(12,9))
ax = fig.add_subplot(autoscale_on=False, 
    xlim=(-0.6,4.5),ylim=(-3,3))
ax.grid()
ax.plot([-p/2,-p/2],[-5,5],'-',lw=2)
line, = ax.plot([],[],'o-',lw=2)
trace, = ax.plot([],[],'-', lw=1)
theta_text = ax.text(0.05,0.85,'',
    transform=ax.transAxes)
textTemplate = '''y = %.1f\n
lenL = %.1f, lenF = %.1f\n
lenL-lenF = %.1f'''
xs,ys = [],[]
def animate(y):
    if(y==-3):
        xs.clear()
        ys.clear()
    x = y**2/p/2
    xs.append(x)
    ys.append(y)
    line.set_data([-p,x,p/2], [y,y,0])
    trace.set_data(xs,ys)
    lenL = x+p/2
    lenF = np.sqrt((x-p/2)**2+y**2)
    theta_text.set_text(textTemplate % 
        (y, lenL, lenF, lenL-lenF))
    return line, trace, theta_text
frames = np.arange(-3,3,0.1)
ani = animation.FuncAnimation(fig, animate, 
    frames, interval=5, blit=True)
ani.save("parabola.gif")
plt.show()

在这里插入图片描述

极坐标方程

圆锥曲线在极坐标系下有相同的表达式,即

在这里插入图片描述

matplotlib中,极坐标图像需要通过projection='polar'来标识,其代码为

p = 2
fig = plt.figure(figsize=(12,9))
ax = fig.add_subplot(autoscale_on=False, projection='polar')
ax.set_rlim(0,8)
trace, = ax.plot([],[],'-', lw=1)
theta_text = ax.text(0.05,0.95,'',transform=ax.transAxes)
textTemplate = 'e = %.1f\n'
theta = np.arange(-3.1,3.2,0.1)
def animate(e):
    rho = p/(1-e*np.cos(theta))
    trace.set_data(theta,rho)
    theta_text.set_text(textTemplate % e)
    return trace, theta_text
frames = np.arange(-2,2,0.1)
ani = animation.FuncAnimation(fig, animate, 
    frames, interval=100, blit=True)
ani.save("polar.gif")
plt.show()

在这里插入图片描述

以上就是Python使用matplotlib绘制动态的圆锥曲线示例的详细内容,更多关于matplotlib绘制动态圆锥曲线的资料请关注脚本之家其它相关文章!

相关文章

  • python中的字符串切割 maxsplit

    python中的字符串切割 maxsplit

    这篇文章主要介绍了python中的字符串切割 maxsplit,具有很好的参考价值,希望对大家有所帮助。如有错误或未考虑完全的地方,望不吝赐教
    2022-12-12
  • 简介Python设计模式中的代理模式与模板方法模式编程

    简介Python设计模式中的代理模式与模板方法模式编程

    这篇文章主要介绍了Python设计模式中的代理模式与模板方法模式编程,文中举了两个简单的代码片段来说明,需要的朋友可以参考下
    2016-02-02
  • Python实现外星人去哪了小游戏详细代码

    Python实现外星人去哪了小游戏详细代码

    今天为大家带来一款小游戏,名叫外星人去哪了,用Python语言实现完成,代码简洁易懂,感兴趣的小伙伴快来看看吧
    2022-03-03
  • Python多层嵌套list的递归处理方法(推荐)

    Python多层嵌套list的递归处理方法(推荐)

    下面小编就为大家带来一篇Python多层嵌套list的递归处理方法(推荐)。小编觉得挺不错的,现在就分享给大家,也给大家做个参考。一起跟随小编过来看看吧
    2016-06-06
  • Python实现解析ini配置文件的示例详解

    Python实现解析ini配置文件的示例详解

    在开发过程中,配置文件是少不了的,而且配置文件是有专门的格式的,比如:ini, yaml, toml 等等。而对于 Python 而言,也都有相应的库来解析相应格式的文件,下面我们来看看 ini 文件要如何解析
    2022-09-09
  • Python用类实现扑克牌发牌的示例代码

    Python用类实现扑克牌发牌的示例代码

    这篇文章主要介绍了Python用类实现扑克牌发牌的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-06-06
  • python列出目录下指定文件与子目录的方法

    python列出目录下指定文件与子目录的方法

    这篇文章主要介绍了python列出目录下指定文件与子目录的方法,涉及Python使用os模块与glob操作目录与文件的技巧,需要的朋友可以参考下
    2015-07-07
  • python3 http.client 网络请求方式

    python3 http.client 网络请求方式

    这篇文章主要介绍了python3 http.client 网络请求方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
    2023-09-09
  • Python调用JavaScript代码的方法

    Python调用JavaScript代码的方法

    这篇文章主要介绍了Python调用JavaScript代码的方法,帮助大家在不同场景里,采用的最佳调用方式,提高程序的性能,感兴趣的朋友可以了解下
    2020-10-10
  • Pytes正确的配置使用日志功能

    Pytes正确的配置使用日志功能

    在pytest自动化测试中,如果只是简单的从应用的角度来说,完全可以不去了解pytest中的显示信息的部分以及原理,可以通过使用推荐的pytest.ini配置,从而可以做到相对来说比较通用的日志配置,这篇文章主要介绍了Pytes如何正确的配置使用日志功能,需要的朋友可以参考下
    2022-12-12

最新评论