python怎么安装sympy库_python中sympy库求常微分方程的用法-程序员宅基地

技术标签: python怎么安装sympy库  

问题1:

2020428114209551.png?2020328114221

程序,如下

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绘制函数图像

2020428114254244.png?202032811434

程序,如下

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)

2020042811414117.png

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

2020428114339336.png?2020328114350

程序,如下

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)

2020042811414118.png

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

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_39847556/article/details/110381601

智能推荐

ARM/X86工业级数据采集 (DAQ) 与控制产品解决方案_daq 设计方案-程序员宅基地

文章浏览阅读1.6k次。I/O设备,包括信号调理模块、嵌入式PCI/PCIE卡、便携式USB模块、DAQ嵌入式计算机、模块化DAQ系统,以及DAQNavi/SDK软件开发包和DAQNavi/MCM设备状态监测软件。工业I/O产品适用于各种工业自动化应用,从机器自动化控制、测试测量到设备状态监测。_daq 设计方案

ESP8266、ESP8266EX、ESP32、ESP-x、ESP8266-01、ESP8266-01S、ESP8266-12F、ATK-ESP8266的区别_esp8266ex和esp826601s_esp8266ex和esp8266的区别-程序员宅基地

文章浏览阅读526次,点赞5次,收藏20次。你的支持,我的动力;祝各位前程似锦,offer不断,步步高升!!!你的支持,我的动力;祝各位前程似锦,offer不断,步步高升!!!更多资料点击此处获qu!!_esp8266ex和esp8266的区别

微信小程序canvas层级过高,遮住弹窗解决方法_小程序中echarts因为小程序原生的canvas层级太高,而导致弹窗这类dom元素无法遮挡-程序员宅基地

文章浏览阅读2w次,点赞8次,收藏19次。页面canvas图层遮住弹出框,导致弹出框无法点击且样式存在问题因为canvas组件是由客户端创建的原生组件,原生组件层级是最高的,所以页面中的其他组件无论设置 z-index为多少都无法盖在原生组件上。想试着用cover-view组件,能覆盖canvas,map,video,camera等原生组件。但是也有一定的限制。比如:cover-view中只支持嵌套cover-view、cove..._小程序中echarts因为小程序原生的canvas层级太高,而导致弹窗这类dom元素无法遮挡

基于java实现的OPT算法_java opt 是什么-程序员宅基地

文章浏览阅读2k次,点赞2次,收藏6次。1966年,Belady提出最佳页面替换算法(OPTimal replacement,OPT)。是操作系统存储管理中的一种全局页面替换策略 。当要调入一页而必须淘汰旧页时,应该淘汰以后不再访问的页,或距最长时间后要访问的页面。它所产生的缺页数最少,然而,却需要预测程序的页面引用串,这是无法预知的,不可能对程序的运行过程做出精确的断言,不过此理论算法可用作衡量各种具体算法的标准。例子:OPT ..._java opt 是什么

linux下的tar压缩及解压缩命令详解_tar -czvf-程序员宅基地

文章浏览阅读10w+次,点赞6次,收藏30次。tar工具 tar [OPTION]...(1) 创建归档 tar -cpvf /PATH/TO/SOMEFILE.tar FILE... 解释:tar -c(c表示创建打包的意思)pvf 创建一个.tar后缀的打包文件;后面加多个文件,或者是文件夹。 p 选项 ;显示打包过程,语法:tar cvf etc.tar /app/etc(2..._tar -czvf

随便推点

小新pro13睡眠后无法唤醒_小钢炮升级版,联想小新Pro 14 2021上手体验-程序员宅基地

文章浏览阅读2.2k次。当我从盒子里拿出联想小新Pro 14 2021的95W适配器时,我发现这款笔记本的性能可能比预估的还要猛,即是搭载了MX450,这种规格的适配器在轻薄本中相当少见,这意味着联想小新Pro 14 2021的CPU和GPU的性能释放要超乎想象。(村村宠粉福利文末送上)联想小新Pro 14 2021我手上这台联想小新Pro 14 2021的核心硬件规格为:i5-1135G7、MX450、16GB双通道L..._唤醒设置小钢炮:

MATLAB 教程_matlab人工智能程序-程序员宅基地

文章浏览阅读459次。1. MATLAB语言及其应用: https://wenku.baidu.com/view/6452c659a55177232f60ddccda38376baf1fe0fb.html2. MATLAB作图教程全:https://wenku.baidu.com/view/1da5a81dc281e53a5802fff0.html3. MATLAB图像处理教程:https://wenku.ba..._matlab人工智能程序

【脚本编程】Shell脚本零基础入门_shell编程入门-程序员宅基地

文章浏览阅读4.8k次,点赞28次,收藏147次。Shell脚本零基础入门_shell编程入门

win10环境安装虚拟机并配置静态IP(CENTOS7镜像)_windows虚拟机分配静态ip-程序员宅基地

文章浏览阅读936次。win10环境安装虚拟机并配置静态IP(CENTOS7镜像)_windows虚拟机分配静态ip

数学家排行榜:高斯和黎曼谁才是近现代最伟大的数学家?-程序员宅基地

文章浏览阅读3.6k次,点赞2次,收藏6次。第一:牛顿,高斯,欧拉,阿基米德第二:柯西,庞加莱,康托尔,凯莱,哈密尔顿,黎曼,爱森斯坦,帕斯卡第三:伽罗瓦,阿贝尔,希尔伯特,克莱因,狄里克雷,莫比乌斯,莱布尼茨,笛卡尔第四:诺特,雷尔曼,欧几里得,勒让德,拉普拉斯,拉格朗日,克罗内克,雅可比,罗巴切夫斯基第五:傅里叶,量子论奠基人(狄拉克、薛定谔、爱因斯坦)http://k.sina.com.cn/arti..._高斯黎曼到底谁第一

Google的Native Client微笑着抽了HTML5一记漂亮的耳光_native client 与html5的关系-程序员宅基地

文章浏览阅读3.8k次。日前,Google在加州总部的一次会议上首次公开演示了Native Client项目(简称“NaCl”),Google已经花了三年多时间研发此项目并表示未来的浏览器中将整合该技术,同时宣布Native Client即将起飞,开发者们准备好了吗?Native Client究竟是啥东西?Google居然撇开HTML5释放如此妖怪,摘录相关介绍中的一段话:Native Client 简单来说就是一套让网_native client 与html5的关系

推荐文章

热门文章

相关标签