Python Sympy计算变限积分积分函数含有对数（log）时，遇到的问题。

已于 2022-04-22 00:19:01 修改
阅读量1.6k 收藏 1
点赞数 1
文章标签： python
于 2022-04-21 18:38:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/libingzhe1111/article/details/124328584
版权
在使用Sympy进行含有对数的变限积分时遇到问题，使用numpy.log会导致错误，而转换为MATLAB代码则能快速得到结果。可能是Sympy库对于此类积分的计算能力有限。
摘要由CSDN通过智能技术生成 
import sympy



Z, x, t ,H1 ,H2,H3,H4= sympy.symbols('Z,x,t,H1,H2,H3,H4')



z11 = Z ** 3 + 0.01



H1 = Z ** 2 * (- (
        -16 * 41.3399927 * Z ** 3 - 48 * 41.3399927 * sympy.log(Z) - Z * (16 * -1634.89996 - 144 * 41.3399927) - (
                -240 * 41.3399927 * Z ** 9 - 120 * 41.3399927 * Z ** 6 - 20 * 41.3399927 * Z ** 3 - 15 * 41.3399927 +
                Z ** 11 * (-60 * 190592.559 + 360 * -1634.89996 - 1620 * 41.3399927) + Z ** 10 * (
                            -120 * -1634.89996 + 1080 * 41.3399927) + Z ** 8 * (
                            -15 * 190592.559 + 90 * -1634.89996 - 405 * 41.3399927) + Z ** 7 * (
                            -24 * -1634.89996 + 216 * 41.3399927) + Z ** 4 * (-15 * -1634.89996 + 135 * 41.3399927)) / (
                        15 * Z ** 12)) / 190592.559 + (
                          -15.8415841584158 * 41.3399927 * z11 ** 1.0 - 48 * 41.3399927 *
最低0.47元/天 解锁文章
libingzhe1111
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Python 计算变上限二重积分的数值模拟进阶
u011699626的博客
06-05 536
变上限二重积分数值算法
Python 计算变上限二重积分的数值模拟基础
u011699626的博客
06-05 828
Python 计算变上限二重积分的数值模拟基础
解决使用sympy解带自然对数函数log（x*y）的一元方程时，程序长时间运行无法求解的问题。sympy解方程长时间运行不报错，但是无法求出解。长时间运行后报错：
最新发布
hrDJP的博客
06-25 199
解决sympy计算带自然对数函数log的方程时，程序长时间计算没有反应。
Python求极限求积分，编程实现
智轩软件的专栏
03-18 1833
Python求极限求积分，编程实现 在PyCharm中编程实现求解极限和积分运算 1.求极限 #!/usr/bin/env python # -*- coding: UTF-8 -*- from sympy import * n = Symbol('n') s = ((n+3)/(n+2))**n print(limit(s, n, oo)) 结果输出 ： E 2求定积分 #!/usr/bin/env python # -*- coding: UT...
python log函数_用 Python 来研究数学 --- SymPy 符号工具包介绍
weixin_39875832的博客
11-03 696
用 PYTHON 来研究数学 — SYMPY 符号工具包介绍SymPy 的简单介绍SymPy 是一个符号计算的 Python 库，完全由 Python 写成，为许多数值分析，符号计算提供了重要的工具。SymPy 的第一个版本于 2007 年开源，并且经历了十几个版本的迭代，在 2019 年已经基于修正的 BSD 许可证开源了 1.4 版本。SymPy 的开源地址和官方网站分别是：GitHub 链接...
用Python实现函数微分、积分与极限【学习笔记】
m0_46820710的博客
09-22 1224
本文使用Python进行了函数微分、积分与极限的一些基本运算，得到了正确的结果。
Python SymPy2：SymPy 极限、导数、积分、级数
Amzmks的博客
07-29 614
1. 极限 limit(f(x), x, x0) x趋向x0时f(x)的极限 limit(f(x), x, x0, '-') x从左侧趋向x0时f(x)的极限，右侧+ from sympy import * x = symbols("x") expr = limit(sin(x), x, 0) print(expr) expr = limit(1 / x, x, 0, '-') print(expr) 输出： 0 -oo 2. 导数 diff(f(x), x) 返回f'(x) .
Python 矩形积分调换积分上下限（代码可优化）
u011699626的博客
05-22 537
在我们使用Python做积分运算的时候，我们本质上计算的是函数与坐标轴x围城的面积。因此在实际操作中，即使我们交换了积分的上下限，积分值也不会发生变化，但是这和我们实际积分的性质是有出入的，因此在实际操作中，我们应该先绘图观察积分的正负，如果积分值为负，则应做想应的处理。这里我们以使用矩形积分法计算sin⁡(x)\sin\left(x\right)sin(x)来进行说明。 代码如下： import numpy as np num = 101 left_vertex = 0 right_vertex = 1
Python Sympy计算梯度、散度和旋度的实例
12-23
sympy有个vector 模块，里面提供了求解标量场、向量场的梯度、散度、旋度等计算,官方参考连接： http://docs.sympy.org/latest/modules/vector/index.html sympy中计算梯度、散度和旋度主要有两种方式： 一个是...
python应用-scipy,numpy,sympy计算微积分
12-21
在Python中进行微积分计算，我们可以利用几个强大的科学计算库，如`scipy`、`numpy`和`sympy`。这些库提供了丰富的功能，能够帮助我们方便地执行积分运算和其他数学操作。本文将介绍如何使用这三个库来计算微积分。 ...
python使用sympy不定积分入门及求解
12-22
为了解决数学中一些不定积分的求解方便自己做作业，让我们一起学习如何使用python求不定积分吧！ 《matlab版点这里》 1.安装 pip3 install sympy 建议使用anaconda，里面有大量的科学包，方便使用！ 2.使用 我会...
integrate函数python_scipy.integrate Pseudo-Voigt函数，积分变为0
weixin_39611820的博客
12-10 760
我正在编写一个脚本，使用Scipy，Numpy和Matplotlib在Python中将峰形拟合到光谱数据 . 它可以同时适应多个峰值 . 峰值曲线(现在)是Pseudo-Voigt，它是高斯(aka Normal)和Lorentzian(aka Cauchy)分布的线性组合 .我有一个选项开关，我可以让软件优化高斯和洛伦兹的贡献或将其设置为固定值(其中0 =纯高斯和1 =纯洛伦兹) . 按预期工作...
用python完成类似于matlab上的积分、求导和极限的操作
scx_qsnl的博客
05-05 1618
python要完成积分、求导和极限的操作的功能仅仅是自己的功能肯定不够，要引入其他的包，而这个包就是sympy包。 这里给出用pycharm导入包教程参考https://blog.csdn.net/scx_qsnl/article/details/89705889 而这个包使用的方法，在代码中讲解 \ from sympy import * #从sympy包里导入所有模块 x = symbol...
SymPy使用教程
热门推荐
hangyudd的博客
05-21 1万+
#本文整理自实验楼，https://www.shiyanlou.com/courses/729，本文直接按照python3格式编写，可以直接执行 #SymPy是符号数学的Python库。它的目标是成为一个全功能的计算机代数系统，同时保持代码简洁、易于理解和扩展。 #SymPy完全是用Python写的，并不需要外部的库。
096: 复习习题 求导题型 Case4：变积分限函数导数；Case5：高阶导数；Case6：分段函数求导
段智华的博客
10-21 877
096: 复习习题 求导题型 Case4：变积分限函数导数；Case5：高阶导数；Case6：分段函数求导
分布函数的变上限积分的最大值问题
03-05 3718
今天，帮别人做了个优化问题问题是一个“分布函数的变上限积分的最大值问题”，还涉及到了分布函数的逆，“Inverse of the normal cumulative distribution function”可能是我遇到过的比较难算的东西。求解思想;1.积分可分解化为求和SUM, 间距划的越小计算时间越长,间距过大精确度降低....2.变上限积分同上,但间距要统一,不能因上限改变,间距改变,解决
Sympy 简介及使用
zcx999666的博客
05-13 1万+
SymPy库常用函数 简介SymPy是一个符号计算的Python库。它的目标是成为一个全功能的计算机代数系统，同时保持代码简 洁、易于理解和扩展。它完全由Python写成，不依赖于外部库。SymPy支持符号计算、高精度计算、模式匹配、绘图、解方程、微积分、组合数学、离散 数学、几何学、概率与统计、物理学等方面的功能。(来自维基百科的描述)更多内容请查看本人个人博客：https://huiyang8...
sympy库使用：第2部分 初等运算，表达式运算（翻译 tcy）
tcy-阿春
02-13 1578
用sympy进行初等运算   =========================================================== 1.初等运算 from sympy import * x,y,z = sympy.symbols('x y z') # 1.1求对数 sympy.log(sympy.E)       # 1 自然对数 sympy.log(8,2)  ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值