python在matlab中有没有对应的digits()和vpa()函数？
在Matlab中，有digits()和vpa()函数用于控制数字的显示精度和进行任意精度的计算。然而，在Python中的标准库中并没有直接对应的digits()和vpa()函数。
在Python中，可以使用第三方库sympy来实现类似的功能。sympy是一个符号计算库，提供了高精度计算和符号计算的功能。
要控制数字的显示精度，可以使用sympy库中的N()函数。N()函数可以将一个表达式或数值转换为指定精度的浮点数。例如，要将一个数值转换为10位精度的浮点数，可以使用以下代码：
from sympy import N
x = 1/3
result = N(x, 10)
print(result)要进行任意精度的计算，可以使用sympy库中的Rational()函数和evalf()方法。Rational()函数可以创建一个有理数对象，而evalf()方法可以将一个表达式求值为指定精度的浮点数。例如，要计算1/3的100位精度的结果，可以使用以下代码：
from sympy import Rational
x = Rational(1, 3)
result = x.evalf(100)
print(result)需要注意的是，sympy库中的高精度计算可能会比较耗时，因此在实际应用中需要根据具体情况进行权衡和优化。
关于腾讯云相关产品和产品介绍链接地址，由于要求不能提及具体的云计算品牌商，无法给出相关链接。但是，腾讯云提供了云服务器、云数据库、云存储等一系列云计算服务，可以根据具体需求选择适合的产品。
相关·内容
我发现自己想要使用作为拼图包的一部分提供的许可，但却对它们应该做的事情感到困惑：
交感表达式和sympy.mpmath表达式到底有什么区别？
如果我想要一个符号表达式的taylor近似，而不了解mp数包在做什么，我可以做以下工作：
#Imports
import sympy
import sympy.parsing
import sympy.parsing.sympy_parser
import Library_TaylorApproximation
#Create a sympy expression to approximate
ExampleStringExpression = &#
浏览 5提问于2015-11-12得票数 2
我正在用python.I做一些科学的计算，我在解线性方程时遇到了麻烦，其中一些系数很大~ E13，有些系数很小~E-69，它给了我一个不准确的解。
代码中的公式非常简单，就像rate[0]=rate[1]=...=rate[6] rate[7]=0.I使用slove来获得解一样，但是rate[0]!=rate[1]=rate[2]=...虽然大多数是正确的，但是物理意义是完全错误的，即unacceptable.rate~rate6必须相等。
我尝试了一些提高准确性的方法。
#convert float to symbols
kf_ = [symbols(str(k)) for k in kf_
浏览 0提问于2018-10-21得票数 0
回答已采纳
我首先定义了一些微分方程：
%% Definitions
% Constants
syms L R J Ke p
% Input
syms ud uq m
% Output
syms id iq ome theta
% Derivations
syms did diq dome dtheta
%% Equations
did=(ud/L)-(R/L)*id+ome*iq;
diq=(uq/L)-(R/L)*iq-ome*id-(Ke/L)*ome;
dome = (p/J)*((3/2)*p*Ke*iq-m);
dtheta = ome;
我现在正在计算R和L。输入和输出变量来自simulin
浏览 1提问于2015-03-19得票数 0
回答已采纳
我用python做了一个有点复杂的积分，我需要得到一个数值，例如2003708.58843。现在，我的python程序给出了一些非常奇怪的结果，其中包含了一些“超级”函数等等。如果有人能在这方面帮助我，我将非常感激。
目前的方案如下：
import sympy as sp
x = sp.symbols('x')
f1 = 5000 + ((-(-7000+x)*(-3000+x))**(1/2))
f2 = 5000 - ((-(-7000+x)*(-3000+x))**(1/2))
minimum1 = 3000
maximum1 = (500/3)*(35-2*((61
浏览 3提问于2020-06-04得票数 1
回答已采纳
我正在尝试集成函数I，它包含Legendre多项式leg_f
import math
import numpy as np
from mpmath import *
from sympy import *
from scipy.special import legendre
n = 3
c = lambda h_c,z,R : (z**2+h_c**2)**0.5
c_supp = lambda h_c,z,R : (z**2+h_c**2)**(-n)
x = lambda h_x,z,R : -4*R*(R-h_x)/c(h_x,z,R)**2
leg_f = lambda
浏览 1提问于2017-12-16得票数 0
回答已采纳
有很多方法可以近似。其中最有效的算法之一是使用以下部分和:众所周知，再匹配是相似的，收敛速度非常快。这个系列是由Chudnovsky兄弟于1988年发现的，并使用了Ramanujan首次发现的方法。此算法的版本用于对每一位Pi数字的所有破纪录计算.编写一个函数estimate_pi_chudnovsky，它使用这个公式来逼近等电点。该函数应该使用SymPy，以一个非负整数作为输入，并以至少1000位精度返回作为SymPy浮点数的估计值。如果输入不是非负整数，则函数应该引发ValueError。
my_answer：
def estimate_pi_chudnovsky(n):
if n
浏览 8提问于2022-05-02得票数 0
我正在通过Euler项目工作，以进一步提高我的编程技能。在重新查看问题3的代码后，我遇到了一个有趣的问题。下面是我的代码：
# prime numbers are only divisible by unity and themselves
# (1 is not considered a prime number by convention)
def isprime(n):
'''check if integer n is a prime'''
# make sure n is a positive integer
n
浏览 4提问于2013-10-26得票数 4
1回答
我使用SymPy作为包含许多指数的函数的代码生成器。因此，对于数值稳定性来说，不计算指数的自变量是很重要的。我想防止这种情况发生： >>> import sympy as sp
>>> x, y = sp.symbols('x y')
>>> expr = sp.exp(5.*x - 10.)
>>> print(expr)
4.53999297624849e-5*exp(5.0*x) 因为它可能导致数值上不准确的结果。 我可以防止指数求值，如下所示： >>> expr = sp.exp
浏览 18提问于2020-04-17得票数 1
回答已采纳
我无意中发现了这件事，因为它的行为不像我预料的那样。s.evalf(subs={t: 0})到底是干什么用的？
它似乎既不是s.evalf().subs({t: 0})的“捷径”，也不是s.evalf().subs({t: 0})作为
import sympy
omega, t = sympy.symbols("omega, t")
s = sympy.pi * sympy.cos(omega*t)
# both result in 3.14159265358979
print(s.evalf().subs({t: 0}))
print(s.subs({t: 0}).eva
浏览 0提问于2018-05-20得票数 4
回答已采纳
我意识到这是一个有点模糊的问题，所以我会尽我所能。如果我有一个像n = 0.783325849821429这样的数字，有没有办法显示更多的小数位？这是用format long生成的。问题是，我正在用正弦进行误差估计计算，有一个点，估计非常接近正弦的实际值，以至于MATLAB将误差计算为0，即使估计不可能是完全准确的。总会有误差，只是误差太小，MATLAB无法识别。有没有办法让MATLAB考虑更多的小数位？
浏览 0提问于2013-02-19得票数 3
1回答
我计划使用SymPy来评估一些链导数，以便手动检查我的工作。不幸的是，第一个等式不会编译。有人知道这是怎么回事吗？
from sympy import symbols, diff
from mpmath import sinh, atanh, sqrt, power
tau, epsilon = symbols("t e")
sigma = sinh(epsilon * atanh(epsilon * tau / sqrt(1 + power(tau,2))))
加薪：
-----------------------------------------------------
浏览 3提问于2018-02-28得票数 1
回答已采纳
1回答
我用这个.m文件用MATLAB计算pi
function calpi(n)
S = 0;
for i=1:n
if mod(i,2) == 0
S=S-1/(2*i-1)*(4*((1/5)^(2*i-1))-(1/239)^(2*i-1));
else
S=S+1/(2*i-1)*(4*((1/5)^(2*i-1))-(1/239)^(2*i-1));
end
end
S = 4*S;
S=vpa(S,50)
当n <= 8的时候，一切都好。但
浏览 2提问于2015-04-19得票数 1
回答已采纳
1回答
所以我试图创建一个脚本来计算从1到1000的所有奇数的乘积(使用MATLAB)。程序运行但产品不正确：%Program是用来计算从1到1000的所有奇数的乘积
% declare variable ‘product’ as zero
product = 0.;
% initialize counter, ‘n’, to 1000
n = 1000;
for i = 1:2:n
product = product + i;
end
fprintf('The product of all the odd numbers from 1 to %d is %d\n', n
浏览 3提问于2014-11-05得票数 1
我作下列进口：
from sympy.matrices import Matrix as sy_matrix
import sympy.mpmath as sy_mp
创建这样的矩阵：
sysMat = sy_matrix([[0.0]*sz1]*sz2)
resVec = sy_matrix([[0.0]]*sz2)
.Populate然后使用python复数(例如1.0+1.0j)，然后调用：
coeffVec = sy_mp.qr_solve(sysMat, resVec)
然而，我得到了以下例外：
File "..\RatSMat.py", line
浏览 1提问于2016-04-23得票数 2
回答已采纳
1回答
我试图计算一些使用很高幂指数的积分。一个例子是：
(-exp(-(x+sqrt(p)).^2)+exp(-(x-sqrt(p)).^2)).^2 ...
./( exp(-(x+sqrt(p)).^2)+exp(-(x-sqrt(p)).^2)) ...
/ (2*sqrt(pi))
其中p是常数(1000是一个典型值)，我需要x=- inf，inf的积分。如果我使用integral函数进行数值集成，结果将得到NaN。如果我将集成的限制设置为-20,20和低p (<100)，我可以避免这种情况，但理想情况下，我需要完整的范围。
我还尝试过设置syms x并使用int和vpa，但在本例中v
浏览 3提问于2014-11-10得票数 1
回答已采纳
我有下面的代码来提取对称二阶张量的平方根。
from sympy import symbols, Matrix, mpmath
import numpy as np
F11, F12, F13, F21, F22, F23, F31, F32, F33 = symbols('F11, F12, F13, F21, F22, F23, F31, F32, F33', real=True)
F = np.array([[F11, F12, F13], [F21, F22, F23], [F31, F32, F33]])
B = F.dot(F.T)
mpmath.sqrtm(Ma
浏览 1提问于2014-11-25得票数 0
我不知道如何将表达式传递给lambda参数。我认为这是sympy.core.power.pow类型的类型问题。
代码如下：
import sympy
def trapezoid(f, a, b, n):
h = float(b - a) / n
result = 0.5 * f(a) + 0.5 * f(b)
for i in range(1, n):
result += f(a + i * h)
result *= h
return result
if __name__ == "__main__":
x =
浏览 22提问于2020-07-30得票数 0
回答已采纳
我得到了一个SymPy矩阵M
In [1]: from sympy import *
In [2]: M = 1/10**6 * Matrix(([1, 10, 100], [1000, 10000, 100000]))
In [3]: M
Out[3]:
Matrix([
[1.0e-6, 1.0e-5, 0.0001],
[ 0.001, 0.01, 0.1]])
我想将输出四舍五入为3小数位，如下所示：
In [3]: M
Out[3]:
Matrix([
[ 0.000, 0.000, 0.000],
[ 0.001, 0.010, 0.100]])
在一般的Python
浏览 2提问于2018-01-28得票数 5
回答已采纳
代码是非常清楚的：
>>> from sympy.utilities.lambdify import implemented_function
>>> f = implemented_function('f', lambda x: x ** 2)
>>> g = implemented_function('g', lambda x: 2 * x)
>>> print(f( 2 ).evalf())
4.00000000000000
>>> print( g(2) .ev
浏览 2提问于2017-01-24得票数 2
回答已采纳
3回答
我正在使用Python内置解决方案和其他一些外部库(如SymPy )研究不同的舍入方法，在这样做的过程中，我偶然发现了一些需要帮助的案例，以帮助我理解其背后的原因。
前例1：
print(round(1.0065,3))
输出: 1.006
在第一种情况下，使用Python内置的舍入函数，输出是1.006而不是1.007，我可以理解这不是一个错误，因为Python舍入到了最接近的偶数，这就是所谓的银行家舍入。
这就是为什么我从一开始就开始寻找另一种方法来控制舍入行为。通过快速搜索，我找到了decimal.Decimal模块，它可以轻松地处理十进制值，并有效地使用quantize()进行舍入，如
浏览 18提问于2018-02-20得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
