python中sympy库求常微分方程的用法


Posted in Python onApril 28, 2020

问题1:

python中sympy库求常微分方程的用法

程序,如下

from sympy import *
f = symbols('f', cls=Function)
x = symbols('x')
eq = Eq(f(x).diff(x, x) - 2*f(x).diff(x) + f(x), sin(x))
print(dsolve(eq, f(x)))

结果

Eq(f(x), (C1 + C2*x)*exp(x) + cos(x)/2)

附:布置考试中两题

1.利用python的Sympy库求解微分方程的解 y=f(x),并尝试利用matplotlib绘制函数图像

python中sympy库求常微分方程的用法

程序,如下

from sympy import *
f = symbols('f', cls=Function)
x = symbols('x')
eq = Eq(f(x).diff(x,1)+f(x)+f(x)**2, 0)
print(dsolve(eq, f(x)))
C1 = symbols('C1')
eqr = -C1/(C1 - exp(x))
eqr1 = eqr.subs(x, 0)
print(solveset(eqr1 - 1, C1))
eqr2 = eqr.subs(C1, 1/2)
# 画图
import matplotlib.pyplot as plt
import numpy as np
x_1 = np.arange(-5, 5, 0.1)
y_1 = [-0.5/(0.5 - exp(x)) for x in x_1]
plt.plot(x_1, y_1)
plt.axis([-6,6,-10,10])
plt.grid()
plt.show()

结果

Eq(f(x), -C1/(C1 - exp(x)))
FiniteSet(1/2)

python中sympy库求常微分方程的用法

2.利用python的Sympy库求解微分方程的解 y=y(x),并尝试利用matplotlib绘制函数图像

python中sympy库求常微分方程的用法

程序,如下

from sympy import *
y = symbols('y', cls=Function)
x = symbols('x')
eq = Eq(y(x).diff(x,1), y(x))
print(dsolve(eq, y(x)))
C1 = symbols('C1')
eqr = C1*exp(x)
eqr1 = eqr.subs(x, 0)
print(solveset(eqr1 - 1, C1))
eqr2 = eqr.subs(C1, 1)
# 画图
import matplotlib.pyplot as plt
import numpy as np
x_1 = np.arange(-5, 5, 0.01)
y_1 = [exp(x) for x in x_1]
plt.plot(x_1, y_1, color='orange')
plt.grid()
plt.show()

结果

Eq(y(x), C1*exp(x))
FiniteSet(1)

python中sympy库求常微分方程的用法

到此这篇关于python中sympy库求常微分方程的用法的文章就介绍到这了,更多相关python sympy常微分方程内容请搜索三水点靠木以前的文章或继续浏览下面的相关文章希望大家以后多多支持三水点靠木!

Python 相关文章推荐
pyside写ui界面入门示例
Jan 22 Python
Python装饰器使用示例及实际应用例子
Mar 06 Python
Scrapy-redis爬虫分布式爬取的分析和实现
Feb 07 Python
Python操作csv文件实例详解
Jul 31 Python
Python将多份excel表格整理成一份表格
Jan 03 Python
python 通过字符串调用对象属性或方法的实例讲解
Apr 21 Python
Pandas之read_csv()读取文件跳过报错行的解决
Apr 21 Python
用python给csv里的数据排序的具体代码
Jul 17 Python
python 对一幅灰度图像进行直方图均衡化
Oct 27 Python
python3通过subprocess模块调用脚本并和脚本交互的操作
Dec 05 Python
selenium与xpath之获取指定位置的元素的实现
Jan 26 Python
Linux系统下升级pip的完整步骤
Jan 31 Python
tensorflow2.0的函数签名与图结构(推荐)
Apr 28 #Python
Python startswith()和endswith() 方法原理解析
Apr 28 #Python
Python如何将函数值赋给变量
Apr 28 #Python
Python多线程thread及模块使用实例
Apr 28 #Python
Python基于模块Paramiko实现SSHv2协议
Apr 28 #Python
Python内置函数locals和globals对比
Apr 28 #Python
使用python实现CGI环境搭建过程解析
Apr 28 #Python
You might like
基于php-fpm的配置详解
2013/06/03 PHP
PHP与Java进行通信的实现方法
2013/10/21 PHP
php动态生成函数示例
2014/03/21 PHP
php获得刚插入数据的id 的几种方法总结
2018/05/31 PHP
将string解析为json的几种方式小结
2010/11/11 Javascript
Three.js源码阅读笔记(Object3D类)
2012/12/27 Javascript
js string 转 int 注意的问题小结
2013/08/15 Javascript
jQuery实现等比例缩放大图片让大图片自适应页面布局
2013/10/16 Javascript
js导入导出excel(实例代码)
2013/11/25 Javascript
重写document.write实现无阻塞加载js广告(补充)
2014/12/12 Javascript
原生javascript实现隔行换色
2015/01/04 Javascript
JavaScript 事件入门知识
2015/04/13 Javascript
基于BootStrap Metronic开发框架经验小结【二】列表分页处理和插件JSTree的使用
2016/05/12 Javascript
为jQuery-easyui的tab组件添加右键菜单功能的简单实例
2016/10/10 Javascript
基于JS实现带动画效果的流程进度条
2018/06/01 Javascript
详解微信小程序获取当前时间及日期的方法
2019/04/28 Javascript
JavaScript 实现HTML DOM增删改查操作的常见方法详解
2020/01/04 Javascript
解决vue-router 切换tab标签关闭时缓存问题
2020/07/22 Javascript
JS删除对象中某一属性案例详解
2020/09/08 Javascript
Node快速切换版本、版本回退(降级)、版本更新(升级)
2021/01/07 Javascript
[03:41]DOTA2上海特锦赛小组赛第三日recap精彩回顾
2016/02/28 DOTA
用Python编写简单的定时器的方法
2015/05/02 Python
python多进程实现进程间通信实例
2017/11/24 Python
Python基于opencv实现的简单画板功能示例
2019/03/04 Python
django中使用POST方法获取POST数据
2019/08/20 Python
Pycharm+Python+PyQt5使用详解
2019/09/25 Python
八一建军节感言
2014/02/28 职场文书
保卫科工作岗位职责
2014/03/01 职场文书
《有趣的发现》教学反思
2014/04/15 职场文书
党政领导班子群众路线对照检查材料
2014/10/26 职场文书
学生不参加考试检讨书
2015/02/19 职场文书
离婚案件答辩状
2015/05/22 职场文书
暑假生活随笔
2015/08/15 职场文书
浅谈pytorch中的dropout的概率p
2021/05/27 Python
单身狗福利?Python爬取某婚恋网征婚数据
2021/06/03 Python
SQL Server代理:理解SQL代理错误日志处理方法
2021/06/30 SQL Server