python怎么安装sympy库_python中sympy库求常微分方程的用法-程序员宅基地
技术标签： python怎么安装sympy库
问题1：
程序，如下
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绘制函数图像
程序，如下
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)
2.利用python的Sympy库求解微分方程的解 y=y(x)，并尝试利用matplotlib绘制函数图像
程序，如下
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常微分方程内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
智能推荐
java常见面试题（160道）_java面试题-程序员宅基地
文章浏览阅读5.4w次，点赞89次，收藏746次。java常见面试题_java面试题
LeetCode刷题总结（C语言版）_leetcode c语言-程序员宅基地
文章浏览阅读5.4k次，点赞6次，收藏73次。编程总结每每刷完一道题后，其思想和精妙之处没有地方记录，本篇博客用以记录刷题过程中的遇到的算法和技巧001）两数之和给定一个整数数组 nums 和一个目标值 target，请你在该数组中找出和为目标值的两个整数。你可以假设每种输入只会对应一个答案。但是，你不能重复利用这个数组中同样的元素。给定 nums = [2, 7, 11, 15], target = 9因为 nums[0] ..._leetcode c语言
小程序开发者工具正常显示，但是真机调试和真机中安卓加载正常ios加载首页失败，首页的请求返回204_苹果 sec-fetch-dest-程序员宅基地
文章浏览阅读125次。检查请求头中的’sec-fetch-dest’: ‘document’ ，是否进行了特殊处理（node层）_苹果 sec-fetch-dest
ansible 批量安装zabbix-agent-程序员宅基地
文章浏览阅读321次。服务器初始化（这是在建立在新的服务器基础上做的初始化）关闭防火墙、selinux，添加epel常用源，安装常用工具、添加普通用户并禁止root1、服务器批量初始化[root@fwd ansible]# cat init.yml 系统初始化脚本---- hosts: all tasks: - name: disable selinux、firew..._ansible批量安装zabbix-agent
java日志系统--log4j配置解析过程，源码分析_log4j 源码分析 读取配置-程序员宅基地
文章浏览阅读1.4w次，点赞3次，收藏2次。Logger.getLogger(Test.class);从getLogger开始，就启动了log4j的整个工作流程，通过调用LogManager获取logger实例return LogManager.getLogger(clazz.getName());LogManager类里面有个静态块static{}，【初始化重要信息】【root logger】，做一些配置，其中url = Loader.ge_log4j 源码分析 读取配置
心灵震撼《一个8岁女孩的遗书》看完能有几人不哭…-程序员宅基地
文章浏览阅读533次。无奈的父亲有一个美丽的小女孩，她的名字叫余艳，她有一双亮晶晶的大眼睛她有一颗透明的童心.她是一个孤儿，她在这个世界上只活了8年，她留在这个世界上最后的一句话是“我来过，我很乖”她希望死在秋天，纤瘦的身体就像一朵花自然开谢的过程.在遍地黄花堆积，落叶空中旋舞的时候，她会看见横空远行的雁儿们.她自愿放弃治疗，把全世界华人捐给她的54万分成了7份，把生命当成希望的蛋糕分给了7个正徘徊在生死线上的小
随便推点
EV/HEV中的牵引逆变器驱动优化-程序员宅基地
文章浏览阅读1.6k次，点赞42次，收藏35次。什么是牵引逆变器？从本质上讲，牵引逆变器是电动汽车动力系统中的一个子系统，它从电池中获取高电压，并将其转换为交流电压——因此被称为逆变器——并基本上为电机供电。它控制电机速度和扭矩，直接影响效率和可靠性，这正成为牵引逆变器设计的设计挑战。此图片来源于网络如今的电动汽车至少有一个牵引逆变器。有些型号实际上不止一个。一个在前轴上，一个在后轴上。甚至一些高端车型实际上每个车轮都有一个牵引逆变器。因此，效率和可靠性非常重要。所以，从逆变器和电机控制的市场趋势来看——从技术趋势来看，我们看到了功率水平的提高。
Ubuntu之apt命令_ubuntu18.04 atp命令使用技巧-程序员宅基地
文章浏览阅读134次。简介apt-cache和apt-get是apt包的管理工具，他们根据/etc/apt/sources.list里的软件源地址列表搜索目标软件、并通过维护本地软件包列表来安装和卸载软件。查看本机是否安装软件：whereis package_name 或者which package_name1.搜索软件sudo apt-cache search pa..._ubuntu18.04 atp命令使用技巧
查询Dynamics 365的Audit History_dynamics 审核历史记录如何查询-程序员宅基地
文章浏览阅读150次。【代码】查询Dynamics 365的Audit History。_dynamics 审核历史记录如何查询
python yield函数的用法-程序员宅基地
文章浏览阅读1.3w次，点赞15次，收藏66次。什么是yield函数?yield函数是python里面的关键字,带有yield的函数相当于一个生成器generator.当你使用一个yield的时候，对应的函数就是一个生成器在python里面类似于return函数,他们主要的区别就是:遇到return会直接返回值,不会执行接下来的语句.但是yield并不是,在本次迭代返回之后,yield函数在下一次迭代时，从上一次迭代遇到的yield后面的代码(下一行)开始执行下面是案例分析:案例一:def gen_generator(): yiel_yield函数
【QT笔记】QFile读文件问题_qfileread后指针会移动吗-程序员宅基地
文章浏览阅读917次。如果不用seek(0)的话，默认是自己会把读取文件的指针后移的，不用手动后移；_qfileread后指针会移动吗
dw8051基本测试示例_dw8051 part1-程序员宅基地
文章浏览阅读2.5k次。整理了网上一份简单的dw8051测试示例，共享到云盘：http://pan.baidu.com/s/1bnu9lZT1.目录如下：---dut ---rtl：DW8051的core文件 ---model：ROM和RAM的model文件---testbench ---rtl.f：filelist文件 ---test_top.v：仿真的top_dw8051 part1