位置: IT常识 - 正文

如何用Python求解微分方程组(如何用python求解航天器追逃博弈)

编辑:rootadmin
如何用Python求解微分方程组 文章目录odeint简介示例odeint简介

推荐整理分享如何用Python求解微分方程组(如何用python求解航天器追逃博弈),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:如何用python求解一元二次方程程序,如何用python求解有约束条件的方程组,如何用python求解一元二次方程,如何用python求解一元二次方程,如何用python求解航天器追逃博弈,如何用python求解一元二次方程,如何用python求解微分,如何用python求解一元二次方程,内容如对您有帮助,希望把文章链接给更多的朋友!

scipy文档中将odeint函数和ode, comples_ode这两个类称为旧API,是scipy早期使用的微分方程求解器,但由于是Fortran实现的,尽管使用起来并不方便,但速度没得说,所以有的时候还挺推荐使用的。

其中,odeint的参数如下

scipy.integrate.odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5, printmessg=0, tfirst=False)

其中func为待求解函数;y0为初值;t为自变量列表,其他参数都有默认选项,可以不填,而且这些参数非常多,其中常用的有

args func中除了t之外的其他变量Dfun func的梯度函数,当此参数不为None时,若将col_deriv设为True,则可提升效率。full_output 如果为True,则额外返回一个参数字典ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5,printmessg 为True时打印信息。tfirst 当为False时,func的格式为func(y,t...),否则格式为func(t, y...)示例

对于常微分方程

如何用Python求解微分方程组(如何用python求解航天器追逃博弈)

θ′′(t)+bθ′(t)+csin⁡θ(t)=b=0.25;c=5θ()=π−0.1;θ′()=\theta''(t)+b\theta'(t)+c\sin\theta(t)=0\\ b=0.25;\quad c=5\\ \theta(0)=\pi-0.1;\quad \theta'(0)=0θ′′(t)+bθ′(t)+csinθ(t)=0b=0.25;c=5θ(0)=π−0.1;θ′(0)=0

将其中的二阶导数项用一个新变量替代,ω(t)=θ′(t)\omega(t)=\theta'(t)ω(t)=θ′(t),则常微分方程可拆分成微分方程组

θ′(t)=ω(t)ω′(t)=−bω(t)−csin⁡θ(t)\begin{aligned} \theta'(t)&=\omega(t)\\ \omega'(t)&=-b\omega(t)-c\sin\theta(t) \end{aligned}θ′(t)ω′(t)​=ω(t)=−bω(t)−csinθ(t)​

令y=[θ,ω]y=[\theta, \omega]y=[θ,ω],则y′=[θ′,ω′]y'=[\theta', \omega']y′=[θ′,ω′],据此可设计函数func

import numpy as npdef pend(y, t, b, c): th, om = y dydt = [om, -b*om - c*np.sin(th)] return dydt

然后调用并求解

from scipy.integrate import odeinty0 = [np.pi-0.1, 0]t = np.linspace(0, 10, 101)sol = odeint(pend, y0, t, args=(0.25, 5))

然后绘制一下结果

import matplotlib.pyplot as pltplt.plot(t, sol[:,0], label="theta")plt.plot(t, sol[:,1], label="omega")plt.legend()plt.show()

这个形状还是比较离奇的。

本文链接地址:https://www.jiuchutong.com/zhishi/297854.html 转载请保留说明!

上一篇:vue 项目中刷新当前路由(vue页面刷新时原有的数据还在吗)

下一篇:[JavaScript你真的以为这么简单吗]JavaScript高级

  • 利用博客评论推广网站的方法分享(利用博客评论推销违法吗)

    利用博客评论推广网站的方法分享(利用博客评论推销违法吗)

  • 苹果11录屏能录通话吗(苹果11录屏能录多长时间)

    苹果11录屏能录通话吗(苹果11录屏能录多长时间)

  • 为什么手机充到80就充不进去了(为什么手机充到97就充不进去了)

    为什么手机充到80就充不进去了(为什么手机充到97就充不进去了)

  • 抖音直播涂鸦礼物在哪里(抖音直播涂鸦礼物多少钱)

    抖音直播涂鸦礼物在哪里(抖音直播涂鸦礼物多少钱)

  • 拼多多步数那个怎么找不到了(拼多多步数是什么)

    拼多多步数那个怎么找不到了(拼多多步数是什么)

  • b站如何把头像挂件取消(b站如何把头像保存下来)

    b站如何把头像挂件取消(b站如何把头像保存下来)

  • 三星手机左下角出现安全模式怎么办(三星手机左下角出现安全模式)

    三星手机左下角出现安全模式怎么办(三星手机左下角出现安全模式)

  • 淘宝显示系统繁忙请稍后再试 怎么办(淘宝的系统繁忙)

    淘宝显示系统繁忙请稍后再试 怎么办(淘宝的系统繁忙)

  • qq vip能干什么(qq vip有什么作用)

    qq vip能干什么(qq vip有什么作用)

  • 拍立得开机键在哪(拍立得按开机没有反应)

    拍立得开机键在哪(拍立得按开机没有反应)

  • 淘宝不发货多久会自动取消订单(淘宝不发货多久会自动收货)

    淘宝不发货多久会自动取消订单(淘宝不发货多久会自动收货)

  • 怎么关闭自己的朋友圈不让别人看(怎么关闭自己的微信朋友圈)

    怎么关闭自己的朋友圈不让别人看(怎么关闭自己的微信朋友圈)

  • 抖音用昵称搜不到人(抖音用昵称搜不到用户)

    抖音用昵称搜不到人(抖音用昵称搜不到用户)

  • 为什么付款方式被拒(为什么付款方式显示无效)

    为什么付款方式被拒(为什么付款方式显示无效)

  • iphone4s什么时候出的(iphone4s什么时候出的上市)

    iphone4s什么时候出的(iphone4s什么时候出的上市)

  • word页眉文字怎么设置(word页眉文字怎么设置居中)

    word页眉文字怎么设置(word页眉文字怎么设置居中)

  • 手机刷机以后怎么恢复以前的系统(手机刷机以后怎样找回以前的信息)

    手机刷机以后怎么恢复以前的系统(手机刷机以后怎样找回以前的信息)

  • 剪映教程音乐识别字幕(如何剪映音乐)

    剪映教程音乐识别字幕(如何剪映音乐)

  • oppo手机隐藏图标怎么打开(oppo手机隐藏图集密码忘记了怎么办)

    oppo手机隐藏图标怎么打开(oppo手机隐藏图集密码忘记了怎么办)

  • 苹果耳机可以连安卓吗(苹果耳机可以连华为平板吗)

    苹果耳机可以连安卓吗(苹果耳机可以连华为平板吗)

  • 微信公众号星标什么意思(微信公众号星标是什么意思)

    微信公众号星标什么意思(微信公众号星标是什么意思)

  • 大麦如何开通免密支付(开通大麦vip有什么用)

    大麦如何开通免密支付(开通大麦vip有什么用)

  • 【前端灵魂脚本语言JavaScript⑤】——JS中数组的使用(前端 自动化脚本 怎么写)

    【前端灵魂脚本语言JavaScript⑤】——JS中数组的使用(前端 自动化脚本 怎么写)

  • Erphpdown 插件使用教程(erphpdown插件安装教程)

    Erphpdown 插件使用教程(erphpdown插件安装教程)

  • 价外税怎么理解
  • 普票千元版和万元版板面区别
  • 手撕定额发票不是免税怎么还报税
  • 出售本月转入处置的固定资产应交税费
  • 发票隔月报送可以么?
  • 国地税合并后有多少个税种
  • 期末调汇的会计科目
  • 个税起征点调整至5000
  • 企业如何申请科研项目
  • 外资企业所得税和内资区别
  • 应收账款进行债务转让
  • 库存现金存入银行怎么记账
  • 事业单位收个人部分社保怎么入账
  • 公司业务招待费占比
  • 税收的凭证有哪几种
  • 融资租赁购入固定资产开票吗
  • 预付油卡怎么入账
  • 房地产企业综合税率是多少
  • 小规模的成本票
  • 债券折价摊销属于借款费用吗
  • 无票收入确认收入
  • 高新企业研发费用认定条件
  • 苹果macos catalina10.15.7
  • 旅游门票怎么做分录
  • 出口业务怎么做分录
  • 资产处置损益是什么科目
  • php面向对象是什么意思
  • php socket_read
  • ajax无刷新更新数据
  • 破产资产拍卖
  • 如何直接访问ip地址
  • 非货币性资产交换
  • 年营业收入平均增长
  • 公司购买的大型安装设备什么时候转入固定资产
  • 银行同业利息 水利基金
  • 预提费用会计处理
  • 投影仪属于固定资产哪一类
  • 仙岩寺庙
  • css加小手
  • ai当前的发展
  • phpcms建站流程
  • 查账征收的纳税人能否简易注销
  • 以前年度多计提折旧
  • 印花税申报怎么更正采集处理
  • 小微企业增值税优惠政策最新2023
  • 长期借款属于什么会计分录
  • sql server 2008数据库文件
  • 预付账款指的是哪些
  • 资产负债表其他综合收益
  • 计提银行存款利息
  • 水电发票如何开
  • 网上缴税app叫什么
  • 公司举办活动的好处
  • 新开办的企业怎么做账
  • 费用报销的凭证是什么
  • 个体工商户是否要交税
  • 主营业务收入多栏式明细账怎么登
  • 房屋装修支出怎么做账
  • 速动资产包括哪些流动资产
  • 高新技术企业享受优惠时间
  • 在mysql中使用视图的限制
  • 正确使用显微镜的七个步骤
  • 目前默认系统%1
  • xp系统没有安装好,请重新运行安装程序
  • linux 命令帮助
  • squid ssh
  • win8pe下载iso
  • skyleap用不了
  • win8使用教程和技能
  • sqlservr.exe - sqlservr是什么进程 有什么作用
  • Linux VPS vsftp/proftpd FTP时间差八小时的解决方法
  • vs2008安装教程
  • Please ensure that adb is correctly located at 'D:Androidandroid-sdkplatform-toolsadb.exe' and
  • jquery filedownload
  • js与jquery的关系
  • 已申报的纳税申报表怎么修改
  • 税务谈心谈话记录2021
  • 江苏 税务局
  • 智能财税证书含金量如何
  • 严守纪律守规矩
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设