同时具有sympy.series.gruntz.calculate_series的x和log(x) - logx变元的Sympy级数
Sympy是一个Python库，用于符号计算和数学建模。它提供了一系列功能强大的工具，可以进行符号计算、代数运算、微积分、方程求解等操作。
在Sympy中，sympy.series.gruntz.calculate_series函数用于计算给定表达式的级数展开。它接受两个参数，第一个参数是要展开的表达式，第二个参数是关于哪个变量进行展开。
对于给定的表达式log(x) - logx，我们可以使用sympy.series.gruntz.calculate_series函数来计算其级数展开。具体代码如下：
import sympy
x = sympy.symbols('x')
expr = sympy.log(x) - sympy.log(x)
series_expansion = sympy.series(expr, x)
print(series_expansion)输出结果为0，这是因为log(x) - logx在x=1附近的级数展开为0。
Sympy的级数展开功能可以在数学建模、微积分、物理学等领域中得到广泛应用。它可以帮助我们分析和理解复杂的数学表达式，并进行近似计算。
腾讯云提供了一系列云计算产品，可以满足各种不同场景的需求。其中与数学计算相关的产品包括云服务器、云函数、人工智能平台等。您可以通过访问腾讯云官方网站（https://cloud.tencent.com/）了解更多关于这些产品的详细信息。
相关·内容
5回答
我想在不使用math.h函数的情况下找到In(x)的值
我正在遵循我在数学书中找到的这个公式：
log(m) base e = 2[ (m-1/m+1) + (1/3) * (m-1/m+1)^3 + (1/5) * (m-1/m+1)^5 + ... ]
下面是我的代码：
i = 3;
logx = 0 ;
ty = (x-1)/(x+1) ;
do
{
logx = logx + ty ;
tty = ty ;
ty = (ty * ((x-1)/(x+1)) * ((x-1)/(x+1))) / i ;
i = i + 2 ;
} while(tty -
浏览 3提问于2010-07-27得票数 0
假设我有一个变量和pow/log常量的方法：
double LogX(int x){
return Math.Log(x, constantBase);
}
double PowX(int x){
return Math.Pow(constant, x);
}
我不知道在使用数学函数时，这些函数的正确时间复杂度是多少。我的概念印象是，PowX将是O(n)，因为它必须将常数n乘以1，但我知道还有其他方法来实现幂函数，并且无法找到关于我应该在那里假设什么的明确答案。如果它是常数时间，那么它是O(n)吗？同样，我也不知道如何正确地处理LogX。
如果这很重要，我将使用C#，但我希望对此
浏览 1提问于2015-12-13得票数 2
回答已采纳
以下函数在给定log(x)和log(y)的值的情况下计算log(x + y)，从而避免在x或y非常大或非常小时发生溢出或下溢：
double log_add(double logx, double logy)
{
return max(logx, logy) + log1p(exp(-fabs(logx - logy)));
}
必须有一个类似的log_sub函数来计算log(x - y)。那是什么？
更一般地，我需要在给定log(x)、log(y)和log(z)的情况下计算log(x - y - z)。从log_add和log_sub我可以分两步计算log(x - y - z)，但也
浏览 0提问于2015-04-23得票数 3
我认为Sympy在计算trig函数的时出错了。例如：
from sympy import fourier_transform, sin
from sympy.abc import x, k
print fourier_transform(sin(x), x, k)
的预期答案是
但是Sympy返回0。有时，由于fourier_transform(Heaviside(t)*cos(t),t,omega)和fourier_transform(Heaviside(t)*sin(t),t,omega)返回正确的答案，函数工作得很好。我认为Sympy可能使用Laplace变换来计算傅里叶变换。
浏览 10提问于2014-01-18得票数 9
回答已采纳
我想展开一个包含一个变量x (假设它很小)的表达式。展开将涉及在x中为多项式的项，以及在log(x)中为多项式的项。示例： (log(x/(1-x)))^2 = log(x)^2 + 2 x log(x) + x^2 log(x) + x^2 + .... 我没有找到一种方法来实现这个目标，但是有一种感觉告诉我这是可能的，因为如果我查看一下渐近/序列/gruntz.py的源代码，并寻找函数calculate_series I find def calculate_series(e, x, logx=None):
""" Calculates at least
浏览 19提问于2020-07-22得票数 0
2回答
Python中的阶乘求和
、、、
我想在python中用符号代数计算阶乘和。我可以生成的问题的最简单版本是：
from sympy.abc import j
from math import factorial
from sympy import summation
summation(factorial(j), (j, 1, 4))
我得到了以下错误：
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "sympy/core/expr.py", line 194
浏览 4提问于2015-02-18得票数 2
回答已采纳
2回答
从这个数据框开始
myDF = structure(list(Value = c(-2, -1, 0, 1, 2)), .Names = "Value", row.names = c(NA, 5L), class = "data.frame")
假设我想对myDF$Value的每一行运行此函数
getNumberInfo <- function(x) {
if(x %% 2 ==0) evenness = "Even" else evenness="Odd"
if(x > 0) positivity = "Po
浏览 0提问于2009-09-08得票数 5
回答已采纳
我正在尝试使用SymPy 功能，但遇到了属性错误。
import sympy as sy
from sympy.abc import x
sy.fps(x / (1 - sy.exp(-x))).truncate()
但我得到了
AttributeError: 'Mul' object has no attribute 'truncate'
我读到过关于NumPy的名称空间问题的文章，所以我尝试为SymPy指定指数函数，但它不起作用。我只是在计算泰勒级数展开
x / (1 - exp(-x))
以及它的各种功能，例如：
(x / (1 - exp(-x)))**
浏览 0提问于2019-11-24得票数 1
1回答
B类扩展了A类。我称A为父类，B为子类。两者都有构造函数。B在其构造函数内部调用super()。两者都有一个同名的方法。也许只是巧合或错误，两者都有一个'this.x‘变量。这样就没有办法访问父级的this.x变量了。然后，它成为孩子和父母之间可能无意中进行交流的一个点。
class A {
constructor(){
this.x = "super x!";
}
logx(){
console.log(this.x);
}
}
class B extends A{
浏览 25提问于2019-04-15得票数 1
1回答
我想使用dsolve来解决初始值问题，但是初始条件ics似乎没有效果。示例中的两个dsolve都给出了相同的结果。
from sympy import *
x = symbols('x')
f = Function('f')
eq = Eq(Derivative(f(x), x), f(x))
# No initial condition
pprint(dsolve(eq, f(x)))
# With initial condition f(0) = 1
pprint(dsolve(eq, f(x), ics={f(0):1}))
在这两种情况下，我得到了
浏览 5提问于2017-10-27得票数 3
回答已采纳
1回答
在matplotlib中显示对数分布数据的直方图(在本例中是随机对数正态分布噪声)很容易做到：
import numpy as np
import matplotlib.pyplot as plt
logx = np.log10(x)
bins = 10**np.linspace(min(logx), max(logx), 101)
plt.hist(x, bins=bins)
plt.xscale('log')
但是，当我试图通过将density=True添加到plt.hist来规范这个分布时，它会由于不同的bin宽度而扭曲分布。
我如何才能绘制一个lo
浏览 5提问于2022-06-01得票数 1
1回答
我有以下数据集：
data have;
input x10 x22 x13 x64;
cards;
20 10 30 1
;
run;
我想创建四个新列，名为log_x10、log_x22、log_x13、log_x64，它们是原始列的日志。我知道这可能是一个非常简单的数组，循环过程，但我对数组相当陌生，无法很好地理解语法。以下是我所拥有的：
data want;
set have;
array var[*] x: ;
do j=1 to dim(var);
logx[j]=log(var[j]);
end;
run;
它不会总是四个变量
浏览 3提问于2014-10-15得票数 2
回答已采纳
我已经看到了以下代码记录宏的一些变体：
(defmacro log
"for debugging, output code and code->val to stdout, returns val"
[code]
`(let [c# ~code]
(prn '~code)
(clojure.pprint/pprint c#)
c#))
但是，我还没有看到一个递归版本可以将所有的宏或函数形式包装在代码体中。有没有人有好的工作实现？
编辑：
为了清楚起见，以下面的代码为例：
(map #(+ % 10) (r
浏览 3提问于2011-02-09得票数 1
回答已采纳
2回答
我正在编写一个相当简单的函数来执行在日志空间中应用梯形规则的积分。 我想添加axis参数来实现类似于numpy.trapz函数的功能，但是我对如何正确地实现它感到有点困惑。 不可广播的函数如下所示： import numpy as np
def logtrapz(y, x):
logx = np.log(x)
dlogx = np.diff(logx)
logy = np.log(y)
dlogy = np.diff(logy)
b = dlogx + dlogy
a = np.exp(logx + logy)
dF = a[
浏览 16提问于2019-02-21得票数 1
回答已采纳
1回答
我已经尝试了R接口rSymPy到SymPy，它工作得很好。但是，我无法找到正确的语法来使用一些更复杂的特性，比如查找Taylor系列。例如，我尝试了以下几种方法：
library(rSymPy)
sympy("var('p')")
#
##### Cannot make this work ???
#
sympy("from sympy.mpmath import *")
xt <- sympy("p=taylor(exp, 0, 10)")
但是它抛出了一个错误：
Error in .jcall("RJavaToo
浏览 6提问于2014-08-19得票数 5
回答已采纳
2回答
数组的函数
、、
我有一个这样定义的数组(提取数据集的第三列)：
value=[]
value.append((p[3]))
x=np.array(value)
如果我想要获得一个包含数组x的Log10 (o一个不同的函数)的新数组，我该怎么做？我已经尝试过了：
logx=np.array(log(x))
但它给出了以下错误：
TypeError: 'numpy.ufunc' object is not subscriptable.
我哪里错了？
浏览 0提问于2013-02-19得票数 1
回答已采纳
1回答
下面是我所做的：
from sympy import *
x = symbols("x")
y = Function("y")
dsolve(diff(y(x),x) - y(x)**x)
我得到的答案(SymPy 1.0)是：
Eq(y(x), (C1 - x*(x - 1))**(1/(-x + 1)))
但这是错误的。Mathematica和Maple都不能解决这个问题。这是怎么回事？
浏览 0提问于2018-01-08得票数 5
回答已采纳
所以，我有这个密码
from __future__ import division, print_function
import sympy as sp
import numpy as np
from sympy.utilities.lambdify import *
u = np.random.uniform(4, 6, 500)
w, k = sp.symbols('w k')
f = sp.log((k - w) * sp.exp((k - w)**5))
l = sum(f.subs(dict(k=k)) for k in u)
现在我想使用l作为w的一个函数。所以我
浏览 3提问于2014-09-11得票数 2
回答已采纳
詹姆斯·斯图尔特的问题“演算早期超越”，第1.2节，第27节
其主要思想是将幂函数作为x的函数，其中
X= c(10,104,8958,11423,76184,114511)
Y= c(5、9、40、39、84、76)
教科书提供的答案是y= 2.3356 x^{0.3072}。
然而，我使用下面的R代码来拟合曲线
x = c(10, 104, 8958, 11423, 76184, 114511)
y = c(5, 9, 40, 39, 84, 76)
df <- data.frame(x, y)
fitCurve <- nls( y ~ I(a*(x^n)), data = d
浏览 2提问于2014-12-04得票数 0
回答已采纳
1回答
from sympy import exp,log
import sympy as sym
x = sym.Symbol('x')
expr = log(exp(x))
print(expr)
print(sym.simplify(expr))
diff_expr = x - log(exp(x))
if diff_expr == 0:
print('Expression was simplified!')
显然，SymPy不能像这样简化表达式。有没有任何方法可以做到这一点，而不必处理它与字符串操作？也就是这样的东西：
func_comb = [
浏览 0提问于2018-08-04得票数 1
回答已采纳
我一直在尝试将一些数据拟合到特定集合x和y的最佳拟合线上。我尝试了很多次，但都没有成功，似乎找不到一种使用yscale('log')和xscale('log')来拟合数据的方法。我得到了这个奇怪的结果，但我似乎找不到为什么它会得到这个奇怪的结果
结果。
我的代码：
#!/usr/bin/env python
# import the necessary modules
import numpy as np
import matplotlib.pyplot as plt
# Generate x and y values which the curve wil
浏览 3提问于2015-03-24得票数 2
2回答
当需要从子级联到父级联方法调用时，组合是继承的有用替代方案，请参见child extended class method calls its super version but that still only sees child data 但是，为了让子对象呈现与父类型兼容的接口，必须实现大量可能都具有相同形式的存根方法，即它们只是将调用传递给前一个父对象(现在是一个组件)。 我的问题是，有没有可能写一个通用的方法？当没有调用任何其他方法时，将调用一个catch-all方法。然后，catch-all将仅将调用中继到父组件。因此，存根只需写入一次。变体可以用来对多重继承等进行排序。 如下所示
浏览 6提问于2019-04-15得票数 0
1回答
为了我的学习目的，我正在研究蟒蛇的交感库。我想计算下一个极限：
因此，我有代码：
import sympy as sp
x = sp.Symbol('x')
f = sp.sign(sp.cos(sp.pi/x))
sp.limit(f, x, 2, '-')
但下一个错误是：
TypeError Traceback (most recent call last)
<ipython-input-31-56b1fd6b589d> in <module>()
1 x
浏览 0提问于2018-11-06得票数 0
回答已采纳
我正在绘制一个函数的曲线，它在p点是正切的，我想独立地管理曲线和切线的xlim。在此代码中，切线半长应为1：
from sympy import init_printing, symbols, N, plot
from sympy import diff
from sympy import log, cos, atan
init_printing()
x = symbols('x')
# Plot a tangent at point (p_x, p_y), of length l
def plot_line(p_x, p_y, x, a, l):
# Comput
浏览 4提问于2022-09-11得票数 0
回答已采纳
我试图用SymPy计算依赖于三角函数sinc ()的函数的taylor级数，为了简化我的问题，我们可以假定我需要的函数是：
f(x1, x2) = sinc(x1) * sinc(x2)
我的问题是，在导入sympy.mpmath时，经常会遇到以下错误：
不能从.创建强积金。
我试过使用近似或 (数字1)，但它们似乎都失败了。例如，对于后面的选项，行：
(sinc(x)*sinc(y)).series(x,0,3).removeO().series(y,0,3).removeO()
返回：
不能从x创建强积金
我还尝试将函数定义为表达式和lambda函数。但似乎什么都起不到作用
浏览 0提问于2014-05-22得票数 3
回答已采纳
2回答
对不起，如果这是简单的，但我不能让它工作，甚至我已经阅读了许多类似的主题数小时。我不知道还能找到什么。
我想调用回调函数，但是调用对象中的回调函数，而不仅仅是全局(?)中的函数。作用域。
var something = function (x) {
this.x = x;
};
something.prototype.alertx = function() {
alert(this.x);
};
something.prototype.logx = function() {
console.log(this.x);
};
something.prototype.m
浏览 0提问于2015-07-20得票数 0
回答已采纳
我试着用Sympy来找出这个序列的极限
((-1)**(n))*(2**(-n))
使用PyCharm调试器和IPython控制台。
python模块是：
import numpy as np
from plotly.graph_objs import Scatter, Figure
from sympy import Symbol, srepr, limit_seq, oo
from sympy.abc import n
import plotly.plotly as py
py.sign_in('****', '**********
浏览 3提问于2018-01-21得票数 3
1回答
我正在使用IPython笔记本。我有一堆相同的情节。我需要使用不同的轴参数来显示这些图，即对于以下四种情况：
Cases x-axis y-axis
1 non-log non-log
2 non-log log
3 log non-log
4 log log
有没有这样简单的方法：
#many lines of code for generating bunch of plots
plt.show()
#figure shown with non-log axis
ax.set_yscale('log'
浏览 1提问于2017-05-10得票数 0
回答已采纳
1回答
使用熊猫数据图中的绘图功能，我试图得到一个正确的对数x轴:样例代码：
import pandas as pd
import matplotlib
import matplotlib.pyplot as plt
fig,ax = plt.subplots()
df = pd.DataFrame({'Freq':[63,125,250,500],'A':[1,2,3,4]})
ax.set_xscale('log')
ax.set_xticks(df['Freq'])
ax.set_xticklabels(df['Freq
浏览 2提问于2020-03-31得票数 0
1回答
我试图通过$log ()和$delegate扩展AngularJS提供的AngularJS，以便捕获对$log和$broadcast的所有调用。这样，其他控制器就可以捕获$log消息并独立地更新网站的其他部分：
var app = angular.module('testApp', [])
.config(function($provide) {
$provide.decorator('$log', function($delegate, $rootScope, logX) {
return logX($delegate);
})
浏览 3提问于2014-10-03得票数 1
回答已采纳
1回答
我试图实现一种算法，该算法涉及两个变量"x“和"y”中的多项式，但有时我需要将它们解释为单变量多项式(即，将x作为常数)，例如，为了使用函数 (扩展的欧几里得算法)。是否有一种简单的方法可以将"x“解释为常量而不是变量？
我试过以下几种方法：
import sympy
x = sympy.Symbol('x', constant=True)
y = sympy.Symbol('y')
f = sympy.Poly(x*y + y**2)
g = sympy.Poly(x+y)
(s, t, gcd) = sympy.gcdex(f,g)
浏览 3提问于2016-10-16得票数 3
回答已采纳
在Sympy中是否有按部分声明函数的规范方法？我试过了
import sympy
import sympy.functions.special.delta_functions as special
sympy.init_printing()
x = sympy.symbols('x', real=True)
V = x*x * (special.Heaviside(x + 1) - special.Heaviside(x - 1)) \
+ (1 + 2*sympy.log(x)) * special.Heaviside(x - 1) \
+ (1 + 2*s
浏览 0提问于2018-08-27得票数 1
回答已采纳
我有一个关于TypeScript进口的新手问题。我试图创建一个类，它在静态变量中保存一些数据，并且数据是在getInstance()方法中延迟初始化的。
myStaticClass.ts：
class MyData {
x = 1;
}
export class MyStaticClass {
private static data: MyData;
static getInstance() {
if (MyStaticClass.data == null) {
console.log('data is null, initialising'
浏览 4提问于2019-10-22得票数 1
回答已采纳
1回答
我正在使用Python的为sin x创建泰勒系列
以下是sin(x)的Taylor级数。参考
然后是如何编写python代码来创建Taylor系列表达式
代码：
import sympy as sym
import math
x = sym.symbols('x')
# Technique 1: Compute the 5 first terms in a Taylor series for sin(x) around x = 0.
T_sin5 = 1 - (x**3/math.factorial(3)) + (x**5/math.factorial(5)) - (x*
浏览 4提问于2020-10-25得票数 1
回答已采纳
1回答
我知道在sympy中向量场的散度是可能的。例如：
>>> from sympy.physics.vector import ReferenceFrame
>>> from sympy.physics.vector import divergence
>>> R = ReferenceFrame('R')
>>> vfield = 3*R[0]**2*R.x+2*R[1]*R.y
>>> divergence(vfield,R)
6*R_x + 2
但是，张量场的散度是可能的吗？我在文档里没
浏览 3提问于2017-08-09得票数 0
1回答
正域上函数的反演
、、
所以，我试图取一个函数的反函数，然后我会把它区分开来。我只对正实领域的解决方案感兴趣。这里有一个子问题，就是如何将指数作为一些例子来处理：
如果我输入x^2 = u，我希望它给我u^(1/2)。
如果我给它u=logx，我希望它给我指数的u。等等
有什么简单的方法吗？问题是，它返回了太多的解决方案，有没有一种方法只是放弃消极的解决方案？
from sympy import *
x, b, a, u, t, dt, dW = symbols('x b a u t dt dW', real = True)
utility = Eq(x**2, u)
invutility = sol
浏览 2提问于2019-06-06得票数 0
回答已采纳
1回答
我如何才能将像log(exp(exp(x)))这样的表达式简化为exp(x)呢？它似乎适用于更简单的表达式，如exp(log(x)) => x。这是一个很小的例子，展示了我迄今所做的尝试：
import sympy
from sympy import exp, log
x = sympy.symbols('x')
a = exp(log(x))
print a
# Gives `x` automatically, no call to simplify needed
b = log(exp(exp(x)))
print sympy.simplify(b), sym
浏览 3提问于2015-05-18得票数 2
回答已采纳
1回答
在博士中，我们可以读到：
manualintegrate模块具有返回所使用的步骤的函数(有关更多信息，请参见模块docstring )。
但是打电话给help(sympy.integrals.manualintegrate)我们得到：
Help on function manualintegrate in module sympy.integrals.manualintegrate:
manualintegrate(f, var)
manualintegrate(f, var)
Compute indefinite integral of a single variable usi
浏览 0提问于2018-02-26得票数 6
回答已采纳
当交互在lm中指定时，R默认包括主要效果，没有任何选项来抑制它们。这通常是适当和方便的，但在某些情况下(在估计器、比率LHS变量等)，这是不合适的。
我有这样的代码，它将日志转换变量与响应变量相匹配，独立于数据子集中。
下面是一个愚蠢但可重复的例子：
id = as.factor(c(1,2,2,3,3,3,4,4,4,4,5,5,5,5,6,7,7,8,8,8,9,9,9,9,10))
x = rexp(length(id))
y = rnorm(length(id))
logx = log(x)
data = data.frame(id,y,logx)
for (i in data$id
浏览 6提问于2013-05-03得票数 0
回答已采纳
1回答
假设我有以下SAS数据集
data have;
input id$ x1 x2 x3;
datalines;
1 10 5 15
2 15 10 7
3 12 10 9
4 14 15 10
;
run;
我想计算几何平均值，几何变异系数(公式是100*(exp(ASD^2)-1)^0.5，ASD是对数转换数据的算术SD )。每个变量x1、x2和x3。如何执行这些操作？由于我在SAS是全新的，我将感谢您的支持。
浏览 12提问于2022-07-20得票数 1
回答已采纳
我的df是： df<- data.frame(logy, logx1, logx2)
dput(head(df, 20))
structure(list(logy = c(2.86483129264695, 2.71422609892467,
6.0029648649718, 6.23407114218406, 3.10441832707604, 3.31883951518659,
2.74899119270203, 3.33693389469922, 3.08234859652005, 2.86894009277142,
3.14037873461243, 6.1199962
浏览 15提问于2019-06-14得票数 1
我为这个基本问题道歉，但无论出于什么原因，我真的被困住了。我希望从y= a*x^b得到'a‘和'b’幂曲线的输出值。
x y log10(x) log10(y)
7 240 0.84509804 2.380211242
45 610 1.653212514 2.785329835
14 340 1.146128036 2.531478917
30 500 1.477121255 2.698970004
24 450
浏览 4提问于2013-10-15得票数 3
回答已采纳
1回答
我想象征性地确定无穷级数的收敛值。使用Mathematica和Sum命令，我得到了一个结果(我不知道是否正确)，我试图用Python重现它，但它无法详细说明它。
我的代码是：
import sympy as sym
from sympy.abc import k, x, d,f
from sympy import Sum, oo
result = sym.Sum(x**(2*(k-1))*sym.exp(sym.I*(2*k-1)*d), (k, 1, oo)).doit()
有什么问题吗？有没有更有效的方法来执行这个操作？
谢谢
浏览 0提问于2021-06-01得票数 0
import numpy as np
import matplotlib.pyplot as plt
x_ = np.array([1.0, 2.0, 4.0, 8.0, 16.0, 33.0])
y_ = np.array([0.4, 0.55, 0.62, 0.72, 0.76, 0.8])
我想拟合一个y = x ** a形式的回归模型，其中估计了a。其目的是推断出(更大)更大的、看不见的x值。注:我需要强迫我的模型通过(0, 0)。
不确定这是否是最好的方法，但到目前为止，我已经尝试将其转换为找到最适合log(y) / log(x)的常量。我有：
x_ = np.array([1.0,
浏览 0提问于2016-06-28得票数 3
回答已采纳
1回答
Sympy返回log而不是ln
、、、、
我有一个等式：
import sympy as sp
x = sp.Symbol(‘x’, real = True)
fx = sp.log(x,3)
sp.diff(fx, x)
Sympy返回：
1/(x*log(3))
Sympy应该返回：
1/(x*ln(3))
为什么Sympy返回日志函数而不是自然日志函数？
浏览 6提问于2020-01-03得票数 3
回答已采纳
1回答
我试图使用fminunc函数进行凸优化。但是，在我的例子中，我用logx作为梯度。如果我的目标函数是F，那么梯度将是
dF/dx = (dF/dlogx) * (1/x)
= > dF/dlogx = (dF/dx) * x
所以
logx_new = logx_old + learning_rate * x * (dF/logx)
x_new = exp(logx_new)
如何在fminunc中实现这一点？
浏览 3提问于2012-05-14得票数 3
1回答
我在交感网站上找到了lambdify的文档：
尝试使用复数的例子似乎分崩离析：
SympyExpression = sympy.parsing.sympy_parser.parse_expr('0.2*exp(1.6*I*pi*x)*log(x - 1.5)')
print "\nSympyExpression.subs(sympy.Symbol('x'), 0.0):"
print SympyExpression.subs(sympy.Symbol('x'), 0.0)
PythonFunction = sympy.lambd
浏览 2提问于2016-06-01得票数 2
回答已采纳
1回答
我有一首颂歌：
x‘= -x + f(x)
看起来很简单，但是x是100维的。
X= x1，.，x100
此外，
fi(x) = ln(xi)/(ln(x1)+...+ln(x100))
其中I介于1到100之间，f(x) = f1(x)，.，f100(x)
在MATLAB的网站上，它说我应该首先创建一个函数：
但我该怎么做呢？我有100个变量，所有的变量都是通过这个高度非线性的函数耦合的。
任何帮助都是非常感谢的！
浏览 3提问于2017-02-04得票数 0
回答已采纳
1回答
我有这个初始的图。
这是它的代码
x = subhalos['SubhaloVelDisp']
y = (subhalos['SubhaloBHMass'] * 1e10 / 0.704) # in units of M_sol h^-1
logx = np.log(x)
logy = np.log(y)
plt.plot(logx, logy, '.')
plt.ylabel('$\log(M_{BH}$ / M$_{\odot}$)' )
plt.xlabel(' $\log(\sigma_{1D}$
浏览 0提问于2016-04-11得票数 0
为什么不在SymPy中工作？
我试过这个：
from sympy import *
var('x y')
print(simplify(log(x*y)))
print(expand(log(x*y)))
print(collect(log(x*y),x))
print(solve(log(x*y),x))
# log(x*y)
# log(x*y)
# log(x*y)
# [1/y]
浏览 0提问于2018-07-30得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
