Python求解微分方程

一、引言
微分方程表示未知函数、未知函数的导数与自变量之间的关系的方程,叫做微分方程。微分方程种类很多,具体分类可参考以下博主的文章:
https://blog.csdn.net/air_729/article/details/139411996
微分方程的解又分成通解和特解,在工程中大多数微分方程是很难得到通解的,因此出现了数值分析或者计算方法这门学科,通过一次次迭代得到方程的某一个或某几个特解,本文中举例出针对Python对某一简单微信方程的特解做一个简单介绍,大家了解此功能即可。
二、实例分析
以以下非线性微分方程为例:
dy/dt = -2t * y + t^2
y(0) = 1。

from scipy.integrate import solve_ivp
import numpy as np
import matplotlib.pyplot as plt
def dydt(t, y):
    return -2 * t * y + t ** 2
y0 = [1]
tspan = (0, 5)
sol = solve_ivp(dydt, tspan, y0, method='RK45')
plt.plot(sol.t, sol.y[0], label='Solution')
plt.xlabel('t')
plt.ylabel('y(t)')
plt.legend()
plt.grid(True)
plt.show()

在这里插入图片描述

  • 4
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值