向量的SymPy替换符号
，指的是在使用SymPy库进行符号计算时，对向量中的符号进行替换的操作。
SymPy是一个Python库，用于进行符号计算和数学表达式的操作。在使用SymPy进行向量计算时，我们可以定义向量中的符号，并进行各种数学运算。而替换符号则是指将向量中的符号替换为具体的数值或表达式。
替换符号在向量计算中具有重要的作用。通过替换符号，我们可以将符号表达式转化为具体的数值，进而进行数值计算和分析。同时，替换符号也可以用于简化符号表达式，将复杂的表达式转化为简单的形式，便于理解和使用。
对于向量的SymPy替换符号，下面是一个示例：
假设我们有一个二维向量v = [x, y]，我们可以使用SymPy库定义这个向量：
import sympy as sp
x, y = sp.symbols('x y')
v = sp.Matrix([x, y])接下来，我们可以进行一些向量运算，比如向量的加法、数乘等：
w = sp.Matrix([1, 2])
u = v + w
s = 2 * v然后，我们可以使用替换符号将符号表达式转化为具体的数值。比如，将x和y替换为具体的数值：
u_sub = u.subs({x: 3, y: 4})
s_sub = s.subs({x: 2, y: 1})最后，我们可以打印替换后的结果：
print(u_sub) # [4, 6]
print(s_sub) # [4, 2]在这个示例中，我们使用了SymPy库进行向量的定义和运算，并通过替换符号将符号表达式转化为具体的数值。这样，我们就可以进行向量计算和分析。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云产品链接地址：https://cloud.tencent.com/product
- 腾讯云云服务器CVM产品：https://cloud.tencent.com/product/cvm
- 腾讯云云原生容器服务TKE产品：https://cloud.tencent.com/product/tke
- 腾讯云云数据库MySQL产品：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云人工智能服务产品：https://cloud.tencent.com/product/ai
- 腾讯云物联网产品：https://cloud.tencent.com/product/iotexplorer
- 腾讯云移动开发产品：https://cloud.tencent.com/product/mobile
- 腾讯云对象存储COS产品：https://cloud.tencent.com/product/cos
- 腾讯云区块链服务产品：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙产品：https://cloud.tencent.com/product/mu
相关·内容
1回答
我不能抨击包含衍生物的表达式。下面是解释这个问题的最低限度代码：
import sympy
x = sympy.symbols('x', real=True)
a = sympy.Function('a')(x)
b = a**2
r = sympy.diff(b, x)
sympy.lambdify((a), b) # works
r.evalf(subs={diff(a, x): 1, a: 2}) # works
b_x = sympy.lambdify((diff(a, x), a), r) # throws an error
最后一行代码引发以
浏览 0提问于2018-03-24得票数 1
回答已采纳
如何将虚拟变量(sympy.core.symbol.Dummy)转换为Sympy中的常规符号？
例如，假设我们想要在一个矩阵的核中找到满足某个方程f(x_1，x_2) = 0的所有向量(x_1，x_2)。
我会把它分成两个步骤。首先：
from sympy import Matrix
M = Matrix( [[1,0],
[0,0] ])
zeros = Matrix([[0],
[0]])
sol = M.gauss_jordan_solve(zeros)[0]
第二: solve f(sol) = 0。
但我不知道如何告诉Sym
浏览 16提问于2018-08-09得票数 0
回答已采纳
我想对lambdify的exp进行抨击，但是当我试图在sympy.Symbol上评估函数时，我遇到了一些有趣的问题。这
import sympy
t = sympy.Symbol('t')
f = sympy.lambdify(t, t**2)
f(t) # no problem
很好，但是这个
t = sympy.Symbol('t')
f = sympy.lambdify(t, sympy.exp(t))
f(t)
给出
AttributeError: 'Symbol' object has no attribute 'exp&#
浏览 0提问于2017-09-20得票数 2
回答已采纳
1回答
我试图用SymPy和SymPy的衍射包在流形上定义一个方程。由于这个方程是由用户输入的，所以它作为字符串发送到程序中，因此在代码中的流形定义之前就定义了它。为了进行有意义的计算，我试图用定义在流形上的符号替换“符号化”符号。
这是一个由用户提供的函数。
H = sympify('m*a - f')
用户还将流形的坐标作为字符串提供。
# Variables defined as symbols (non-diffgeom)
a = Symbol('a')
f = Symbol('f')
之后的一切都是自动化的。
from sympy.diffge
浏览 3提问于2016-06-30得票数 2
回答已采纳
1回答
我想用下面的代码计算函数的导数。
import sympy
pi = sympy.symbols("pi")
class H(sympy.Function):
nargs = 1
def fdiff(self, argindex=1):
x = self.args[0]
return - sympy.functions.exp(-sympy.Pow(x, 2) / 2) / sympy.sqrt(2 * pi)
def G(a):
return (
(a + 1) * H(1 / sympy.sqr
浏览 2提问于2020-02-13得票数 0
回答已采纳
3回答
我正在使用pretty_print命令来使用渐近来编写上标。这需要与一些变量(代数(x，y等))一起导入为了让pretty_print命令起作用。如果不从sympy.abc导入变量(代数)，则需要将字母放在引号''中，在引号中pretty_print命令不起作用。
例如，如果我在屏幕上打印x²，我需要导入x，以便在程序中使用它。如下所示：
import sympy
from sympy import pretty_text as exp
from sympy.abc import x
equation= x**2
exp(equation)
这将在屏幕上打印x²。
浏览 0提问于2013-02-08得票数 1
2回答
SymPy -将数字视为符号
、、、、
如何在SymPy中将数字视为符号？
例如，如果我使用符号执行因式分解，则会得到：
from sympy import factor
factor('a*c*d + a*c*e + a*c*f + b*c*d + b*c*e + b*c*f')
c*(a + b)*(d + e + f)
当我在表达式中使用数字时，我也希望有同样的行为。而不是
from sympy import factor
factor('2006*c*d + 2006*c*e + 2006*c*f + 2007*c*d + 2007*c*e + 2007*c*f')
4013*c*(d + e
浏览 4提问于2018-11-29得票数 1
回答已采纳
2回答
我需要使用渐近和符号来计算向量的范数。代码摘录如下： from sympy import *
x = Symbol('x')
sb = [2,1]
func = sympy.exp(-(sympy.sqrt((x.norm() + (x-sb).norm())**2 - sb**2)/(2)))
func_prime = func.diff(x)
ff = lambdify(x, func_prime, 'numpy')
f = -1*ff(np.array(r)) 上面的实现给出了错误：AttributeError: 'Symbol' obj
浏览 257提问于2020-09-09得票数 1
我正在尝试添加向量，然后做一个点积。但是，我不能做乘积，因为添加的向量的类型是错误的。 a = N.x*2
b = N.y*2
vec = a + b
print("Vector type", type(vec)) 输出：Vector type <class 'sympy.core.add.Add'> 预期输出：<class 'sympy.vector.vector.VectorAdd'>
浏览 86提问于2020-05-03得票数 0
回答已采纳
我想在向量模块中计算内积。如果可能的话，请让我知道如何做到这一点，而不用i或j。
from sympy import *
from sympy.vector import CoordSys3D
var('ABx ABy ACx ACy')
N =CoordSys3D('N')
AB=ABx*N.i+ABy*N.j
AC=ACx*N.i+ACy*N.j
BCx2y2=(AC-AB).components[N.i]**2+(AC-AB).components[N.j]**2
B2DtC2=AB.components[N.i]**2+AB.components[N.j
浏览 4提问于2022-06-13得票数 0
回答已采纳
2回答
我试图使用‘numpy.linalg.解题’为Ax=0找到一个非零解，但是这个包只给出了一些零向量的解。
因此，我尝试按以下方式搜索相关代码：
import sympy
import numpy as np
A = np.array([[0.8, 0.1, 0.1],
[0.7, 0.2, 0.1],
[0.6, 0.3, 0.1]]) - np.eye(3)
dim = A.shape[0]
x = sympy.symbols([f"x{idx}" for idx in range(dim)])
# Creat
浏览 20提问于2022-10-13得票数 1
1回答
如何对SymPy中的函数符号执行泰勒展开？
例如
from sympy import *
ode = f(x).diff(x, 2) - sin(f(x))
我们希望将微分方程线性化，方法如下
ode.series(f, 0, 1)
获得类似的东西
f(x).diff(x, 2) - f(x)
但我想不出如何在SymPy中做到这一点。在Maxima中，我可以定义一个依赖项，例如
depends(y, x);
ode: diff(y, x, 2) - sin(y);
taylor(ode, y, 0, 8);
它会导致
'diff(y,x,2) - y + y^3/6 - y^5/120
浏览 6提问于2016-11-18得票数 4
回答已采纳
我有一个包含dot乘积的lambdifyd渐近函数，例如，
import numpy as np
import sympy
class dot(sympy.Function):
pass
x = sympy.Symbol('x')
a = sympy.Matrix([1, 1, 1])
f = dot(x, a)
ff = sympy.lambdify((x), f, modules='numpy')
x = np.random.rand(3)
print(ff(x)) # okay
(奇怪的是，定制的dot声明是有效的。不知道为什么，但别
浏览 0提问于2016-07-08得票数 3
1回答
我试图使用SymPy来计算一些表达式，但有时单独使用.subs()函数，或者后面跟着.evalf()，都不会给出任何浮点数或整数值。
##--IMPORT
from mpmath import *
from sympy import *
from sympy.parsing.sympy_parser import parse_expr
x = symbols("x")
_fOfX = "e**x + x + 1"
if __name__== "__main__":
#Function to sympy expression
浏览 2提问于2017-05-27得票数 0
回答已采纳
我想用sympy来求解一个非线性方程组。不幸的是，这是我的计划
import sympy
x = sympy.Symbol('x')
f,g,h= sympy.symbols("f g h", cls=sympy.Function)
#here I do some computations to arrive at the equations below
eq1 = 2*f(x)*g(x)-h(x)
eq2 = f(x)**2+h(x)*g(x)
system = [eq1,eq2]
unknowns = [f(x),h(x)]
我想为f(x)和h(x)解决这个
浏览 0提问于2018-04-04得票数 1
回答已采纳
1回答
我正在使用ipython控制台和sympy编写一个小计算器。设置是 import sympy
x, y, z = symbols('x y z')
f = (x+2)**2
g = lambda y: y+3 如果能看到哪些符号已经在渐变中定义了，那就太好了。如何获取所有症状符号的列表？即包含x、y、z和f但不包含g的列表，因为g是普通的python lambda。 我正在寻找类似vars()或globals的东西，但只包含sympy相关的定义。
浏览 12提问于2020-03-31得票数 0
回答已采纳
2回答
我试着用这样的线来表示
str1="a^0_0"
ns={}
ns['a^0_0']=Symbol('a^0_0')
pprint(sympify(str1,locals=ns))
但是我得到了以下错误
Traceback (most recent call last):
File "cuaterniones_basic.py", line 114, in <module>
pprint(sympify(str1,locals=ns))
File "/usr/local/lib/python2.
浏览 1提问于2018-01-26得票数 3
回答已采纳
我希望在Symbols中扩展SymPy类，以便可以添加一个布尔属性。我可以用一个符号来完成这个任务(请参阅我的问题，以及其他人的问题)。实现这一目标的代码如下所示：
from sympy.core.symbol import Symbol
class State(Symbol):
def __init__(self, name, boolean_attr):
self.boolean_attr = boolean_attr
super(State, self).__init__(name)
但是，这个解决方案的问题是，当我定义一个多项式或涉及多个Stat
浏览 1提问于2019-08-02得票数 0
回答已采纳
2回答
假设在一个交感表达式中有许多非交换符号，比如
a, c = sympy.symbols('a c', commutative=False)
b = sympy.Symbol('b')
expr = a * c + b * c
使表达式中的所有符号交换的首选方法是什么，以便，例如，sympy.simplify(allcommutative(expr)) = c * (a + b)
在中，有人说如果不替换符号，就无法改变符号的交换性，但是也许有一种简单的方法可以在这样的表达式中用块来改变所有的符号？
浏览 5提问于2018-01-16得票数 2
回答已采纳
1回答
在下面的代码中，我定义了实数符号来获取渐近极限，然后使用subs来获得数值结果。 #To take the limit with unknown values, I have to define real symbols
from sympy import *
x = symbols("x")
H = symbols("H")
p = symbols("p")
#-----
# Position_1
#-----
dict = {H: 5,p:3,x:1}
### Some operations
H = Symbol("H"
浏览 37提问于2020-04-03得票数 0
回答已采纳
我很难在一起使用交感和兰达。我想声明一个函数(func_mx)，它接受参数S。
然后，我想将它包含在一个符号表达式X中。
下面是我要执行的代码：
import sympy as sp
S = sp.Symbol('S')
A = sp.Symbol('A')
mux_m = sp.Symbol('mux_m')
KX = sp.Symbol('KX')
func_mux = sp.Lambda(S, (mux_m * S) / (KX + S))
X = sp.sympify("(func_mux(S) - D) *
浏览 1提问于2018-03-16得票数 1
回答已采纳
2回答
对于我正在工作的一个项目，我需要一个函数对wrt cos(theta)的导数，但当使用Sympy v1.5.1时，会得到一个错误消息，指出非符号不能用作导数。这在Sympy v1.3之前是没有问题的，但是后来的版本给出了这个错误。
>>> l=1
>>> theta = symbols('theta')
>>> eq=diff((cos(theta)**2-1)**l,cos(theta),l)
Traceback (most recent call last):
File "<string>",
浏览 1提问于2020-03-01得票数 0
我要丢番图转换(t_0 --> n)
from sympy import *
var('x y t_0 n')
def myDiophantine(myf):
myf_diox = list(diophantine(myf))[0][0]
myf_dioy = list(diophantine(myf))[0][1]
print("#1#",type(myf_diox),myf_diox)
myf_diox = myf_diox.subs({t_0:n})
print("#2#",type(myf_d
浏览 4提问于2022-03-31得票数 1
回答已采纳
1回答
当矩阵和特征向量是时间相关时，我正在解一个特征值问题。矩阵维数为8x8，为hermitian矩阵。时间相关矩阵的形式如下：
import sympy as sp
t, lbd = sp.symbols(r't,\lambda', real=True)
Had = ...
print(repr(Had))
Matrix([[2*t,0, 0, 0, 0, 0, 0,0],
[ 0,-2*t, 2*t*(1 - t), 0, 0, 0,0,0],
[0, 2*t*(1 - t),0,0, 2 - 2*t, 0,0,0],
[0,0,0,0, 0, 2 - 2*t, 0,0],
浏览 3提问于2022-04-17得票数 0
回答已采纳
2回答
我希望能够获取一个索引对象，用值替换它的索引，但仍然记住它的原始索引。
from sympy.tensor import IndexedBase, Idx
C = IndexedBase("C")
i,j = Idx("i"), Idx("j")
expr = C[i,j]
expr = expr.subs({i:1, j:2})
可以从expr获得原始索引(i，j)吗？因为表达式是不可变的，所以我对原始表达式所做的任何修改都不会显示在新表达式中。
我认为这里的一般问题是如何将附加信息附加到Sympy符号。
浏览 3提问于2015-02-20得票数 0
我更新了我的Python3和sympy，突然我的代码不能工作了。我尝试使用以下代码： import sympy as sym
x, y, z= sym.symbols('x[0], x[1], x[2]')
theta = sym.atan2(y,x)
t = theta + pi
r_1 = 0.0124*sym.cos(t)+0.0414*sym.sin(2*t)
f = (r_1)*sym.cos(t)
diff_f = sym.diff(f,theta) 当我运行它时，我收到以下错误： Can't ca
浏览 374提问于2020-09-25得票数 0
1回答
我正在尝试推导一个离散分布的CDF，它的PMF由下式给出： P(d)=log_10(1+1/d) for d=1,...,9 下面是我派生CDF的代码： import sympy
#Find the CDF
d=sympy.Symbol('d')
pmf=sympy.log((1+1/d),10)
CDF=sympy.Sum(pmf,(d,1,d))
print(CDF)
#compute the CDF
CDF=CDF.doit()
print(CDF)
#simplify the expression
CDF=CDF.simplify()
print(CDF) 这不求值
浏览 37提问于2020-06-14得票数 0
回答已采纳
1回答
用SYMPY绘制微分方程解
、、、、
from sympy import *
import sympy as sp
from sympy.abc import *
import numpy as np
import matplotlib.pyplot as plt
x=sp.Function('x')
t=symbols('t')
w=int(input("enter the frequency in hz"))
eq= sp. Eq(x(t).diff(t,2)+w**2*x(t)-f*sp .cos(w*t),0)
c= sp. dsolve(eq,x(t),ics={x(0)
浏览 13提问于2022-07-05得票数 -2
1回答
如何计算Sympy CoordSys3D表达式上的点？
使用Sympy物理模块，此过程与预期的一样工作
from sympy import *
from sympy.physics.vector import ReferenceFrame, CoordinateSym, divergence
R = ReferenceFrame('R')
field = R[0]*R[1]*R[2]*(R.x+R.y+R.z)
div = divergence(vect=field, frame=R)
div_func = lambdify([R[0], R[1], R[2]], div,
浏览 1提问于2018-07-23得票数 3
回答已采纳
我正在尝试使用渐近解一个方程，但是由于某种原因，我得到了以下错误： “‘Symbol”对象没有属性“pi” 代码如下所示： solveset(Eq(parse_expr("sin(math.pi/6)+sin(a)"),parse_expr("1")),a,domain=S.Reals) 有人知道在这种情况下该怎么做吗？ 我的import语句如下： import re
from sympy import *
from sympy.parsing.sympy_parser import parse_expr
import math
浏览 20提问于2019-05-12得票数 0
1回答
我现在正在学习SymPy，我想知道是否有一种方法可以用一般形式中的一个变量来区分一个函数。 考虑这个例子： 这里有一个向量，让我们把它的分量写成x_1，x_2，x_3。这样的向量的长度r将是r= sqrt(x_1^2 + x_2^2 + x_3^2)。 我想在x_i上区分这个向量，在这里我没有指定i，它应该给出类似x_i/r^3的结果。 在SymPy中可以做到这一点吗？ 抱歉，缺少方程式渲染..
浏览 25提问于2021-06-29得票数 3
回答已采纳
如何解决多目标问题？ import sympy as sym
from sympy import lambdify
x = sym.Symbol('x')
n = sym.Symbol('n')
f = sym.sin(n*x)
derivative_f = f.diff(x)
derivative_f = lambdify(x, derivative_f)
x = float(input('x:'))
print(derivative_f(x))
print(derivative_f) 如果我输入2，那么预期的结果应该是2*cos(2*
浏览 113提问于2021-08-02得票数 0
回答已采纳
1回答
我正在尝试绘制一个方程式，这个方程式是一个solve块的结果，下面是我的代码和下面的错误消息：
%pylab inline
from sympy import init_printing;init_printing()
from sympy import *
d,vf,a,vi,t,x,h,g,theta=symbols('d vf a vi t x h g theta')
equations=[Eq(sin(theta),(0.5*g*t**2+h)/(vi*t)),Eq(cos(theta),x/(vi*t))]
ans=solve(equations,[h,t],dict
浏览 2提问于2013-07-24得票数 0
我是Python的新手，并且不断收到以下错误
..., line 27, in <module>
eq=(p**2+2)/p/sqrt(p**2+4)
AttributeError: sqrt
我尝试添加math.sqrt或numpy.sqrt，但这两个都不起作用。有人知道我哪里错了吗？
我的代码是：
from numpy import *
from matplotlib import *
from pylab import *
from sympy import Symbol
from sympy.solvers import solve
p=Symbol('
浏览 1提问于2014-07-18得票数 0
2回答
我正在尝试一种对矩阵有同情的SVD：
[3, 2, 2]
[2, 3, -2]
但是，恢复的矩阵与原始矩阵不匹配。
这是我的密码：
import sympy
A = sympy.Matrix([
[3, 2, 2],
[2, 3, -2]
])
A1 = A@A.transpose()
[U, D1] = A1.diagonalize(normalize=True)
A2 = A.transpose()@A
[V, D2] = A2.diagonalize(normalize=True)
V_T = V.transpose()
S = sympy.sqrt(D2).doit()
浏览 9提问于2020-11-04得票数 0
回答已采纳
1回答
这个Sympy代码与我预期的一样工作：
>>> x = sp.Symbol("x")
>>> y = sp.Symbol("y")
>>> z = sp.Symbol("z")
>>> (x+y+z).evalf(subs={"x":1, "y":2, "z":3})
6.0
但是，如果我使用实值符号，则表达式不会简化：
>>> x = sp.Symbol("x", real=True)
>
浏览 6提问于2022-07-29得票数 2
回答已采纳
如何用初始变量的新定义替换函数中的变量？
import sympy as sy
g, a, x = sy.symbols("g a x")
g = 5*a*x
此时在解释器中输入g打印：
5*a*x。
现在我想用α代替打印的a。为什么下列内容不起作用？
a = sy.symbols("alpha")
g = g.subs(a,a)
我清楚地知道，我可以利用以下几点：
import sympy as sy
g, a, x = sy.symbols("g a x")
g = 5*a*x
alpha = sy.symbols("alpha
浏览 1提问于2018-09-17得票数 1
回答已采纳
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
回答已采纳
2回答
我需要解一个方程，但是蜘蛛返回了这个错误：
import sympy as sym
import sympy as sympy
Re=100
Epsilon=0.00000075
D=0.01
f=symbols('f')
eq=(((1/(-2* sympy.log((Epsilon/D/ 3.7)+( 2.51 / (Re*(f**0.5))), 10)))**2 -f),0)
sym.solve((eq), (f))
print(f)
NotImplementedError: could not solve -4*f*log(20270270270270
浏览 1提问于2022-03-17得票数 0
1回答
我是sympy的新手，不知道如何关联符号和打印这种关系。
import sympy
U, R, I = sympy.symbols('U, R, I')
U = R * I
sympy.simplify(R * I)
我想打印U而不是R * I。
在这里，我正在覆盖U，我正在寻找一种方法来定义符号之间的关系，并打印该关系。这个是可能的吗？
感谢您的帮助:)
浏览 3提问于2018-11-16得票数 1
1回答
有没有办法将SymPy中的笛卡儿表达式转换为极性表达式？我有以下表达：
1/平方吨(x^2+y^2)
然而，我似乎无法让SymPy承认这是极坐标的1/r。我尝试使用“subs”命令，这两个选项都是(我以asp的形式导入了渐近，并在前面定义了所有符号)：
expr = 1/sp.sqrt(x**2+y**2)
expr.subs((x,y),(r*cos(theta),r*sin(theta))
expr.subs((sp.sqrt(x**2+y**2),sp.atan(y/x)),(r,theta))
但在这两种情况下，我只是简单地收到原来的费用。
在SymPy中，是否有方法将笛卡尔表达式转换
浏览 0提问于2018-04-21得票数 3
回答已采纳
1回答
在Sympy中，当符号是实的(real=True)时，它不能用subs()代替。即使real是false，它也不起作用。每当在real调用中指定symbols()时，替换subs()就无法工作。否则，它的作用是：
import sympy
X = sympy.symbols("X")
print('i', X.subs( {'X': 1.01} )) # works
X = sympy.symbols("X", real=True) # `real=` Seems to be the cause of the proble
浏览 7提问于2022-05-22得票数 0
回答已采纳
3回答
基本上我写了这段代码..
import sympy
number = 5.62527e-11;
f = sympy.Symbol('f');
answer = sympy.solve(Eq(f,number),f);
print(answer);
不出所料，答案是"5.62527e-11“
但是，当我将"5.62527e-11“更改为"5.62527e-67”并运行相同的代码时，...i会得到0.0作为输出，这很尴尬...it对"e-11“有效，但对"e-67”无效。有人能解释一下为什么会发生这种情况吗..and有人能告诉我这个问题的解决
浏览 9提问于2014-09-23得票数 0
我有一个方程1_x+2_y=40。我想为你解决这个问题。让它形成y=.)然后用x代替x，得到y的数字。我目前的半手工解决方案如下：
import matplotlib.pyplot as plt
import numpy as np
from sympy import symbols, solve
m=40
pa=1
pb=2
x, y = symbols('x, y')
eq = pa*x+pb*y-m
y= solve(eq, y)
print(y)
print(type(y))
print(type(y[0]))
# then watch on screen wh
浏览 2提问于2021-09-25得票数 0
回答已采纳
1回答
当v(t)不是平方时，这段代码只适用于求解微分方程v_equation。当我平方它时，它返回错误PolynomialDivisionFailed。有没有其他方法可以用Sympy来做这件事，或者我应该找一个不同的python包来做这些类型的计算。
from sympy import *
from matplotlib import pyplot as plt
import numpy as np
m = float(raw_input('Mass:\n> '))
g = 9.8
k = float(raw_input('Drag Coefficient:\n>
浏览 12提问于2016-08-22得票数 2
回答已采纳
我对使用Sympy真的很陌生，我正在尝试解决以下方程系统： a*cos(theta1) + b*cos(theta2) + c*cos(theta3) + d*cos(theta4) =0 a*sin(theta1) + b*sin(theta2) + c*sin(theta3) + d*sin(theta4) =0 其中theta1是输入变量，它会随时间变化(类似于theta1 = 2*t)，而除theta2和theta3之外的所有其他变量都是常量。 我希望Sympy根据输入变量为我提供theta2和theta3的表达式。 我尝试的解决方案如下： from sympy import *
浏览 37提问于2021-04-29得票数 1
回答已采纳
1回答
我正在复习一些线性代数，并研究Python中的一些实现。我正在处理一个与寻找矩阵A的特征向量有关的问题。
A = [[ 1, 2,-2],
[-2, 5,-2],
[-6, 6,-3]]
当我用手解决这个问题时，我得到了特征值3和-3，其中3的多重性为2，我的特征向量是[[1/3], [1/3], [1]], [[1], [1], [0]], [[-1], [0], [1]]。
尝试在NumPy中实现：
import numpy as np
A = [[ 1, 2, -2],
[-2, 5, -2],
[-6, 6, -3]]
np.linalg
浏览 0提问于2018-02-12得票数 1
2回答
评估尽可能多的参数(SymPy)
、、、、
我在SymPy中有这样的代码：
from sympy import *
par_amplitude, par_const_force = symbols('Delta_U f')
equation = par_amplitude * par_const_force
print(equation.evalf(subs={'f':3, 'Delta_U':2}))
输出显然是6。
当我只能访问一个参数时，我的问题就会出现，如下所示：
print(equation.evalf(subs={'Delta_U':2}))
我得到的是Delt
浏览 27提问于2017-08-11得票数 0
1回答
我用Sympy来计算一个复杂势的导数。我把x，y，D0，e，C，k，d，b定义为Sympy符号。然后，我继续做以下定义：
import sympy as sm
x, y, D0, e, C, k, d, b = sm.symbols("x, y, D0, e, C, k, d, b")
phi = sm.atan2(y,x)
D = D0 + e*D0*sm.sin(k*phi)
rho = sm.sqrt(x**2 + y**2)
U = (2**b)/((2*D-d)**b)*(D - rho)**b
符号"U“代表2D势。
现在，当我将此与x区分时，使用：
sm
浏览 2提问于2016-04-19得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
