- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我正在编写一个必须解决某些递推关系的应用程序,但有些关系无法使用 sympy
的 rsolve
方法解析解决>。它只返回 None
。有没有办法强制 sympy 以数字方式解决它们?
我有这样的东西:
from sympy import *
ctx = {
"f": Function("f"),
"x": var("x",integer=True)
}
initial_conditions = {
0: 1,
1: 1,
2: 1,
3: 1
}
f = sympify("-2*f(x-1)+11*f(x-2)+12*f(x-3)-36*f(x-4) +41**(x-4)+3 -f(x)", ctx)
# calculate f(10) here without creating a closed from
# The code below will not work rsolve returns None
solve_for = sympify("f(x)", ctx)
solved = rsolve(f, solve_for, initial_conditions)
我希望有人能帮助我!
最佳答案
这是我用数字评估递归关系的解决方案。当您将递归关系指定为简化的输入时,请确保值 f(x)
不是字符串的一部分。也就是说,如果你的递推关系是
f(x) = -2*f(x-1)+11*f(x-2)+12*f(x-3)-36*f(x-4) +41**(x-4)+3
你的输入字符串应该是:
"-2*f(x-1)+11*f(x-2)+12*f(x-3)-36*f(x-4) +41**(x-4)+3"
此外,此解决方案仅限于线性递推关系,但也可扩展到其他情况。
代码所做的是遍历递归关系的语法树,并通过数值计算或查找 f(x) 的已知值来评估每个节点。
from sympy import *
import operator
ctx = {
"f": Function("f"),
"x": var("x",integer=True)
}
initial_conditions = {
0: 1,
1: 1,
2: 1,
3: 1
}
func1 = sympify("-2*f(x-1)+11*f(x-2)+12*f(x-3)-36*f(x-4) +41**(x-4)+3", ctx)
def contains_function(f):
if issubclass(type(f),Function):
return True
r = map(contains_function,f.args)
return (sum(r) != 0)
def get_numeric_value(arg):
if arg.is_number:
if arg.is_integer:
return int(arg)
else:
return float(arg)
else:
return None
def evaluate_at(f, n, initial_conditions):
if f.is_Add:
result = 0
op = operator.add
elif f.is_Mul:
result = 1
op = operator.mul
elif f.is_Function:
func_arg = f.args[0]
func_arg_val = int(func_arg.subs(func_arg.free_symbols.pop(),n))
if not func_arg_val in initial_conditions:
return None
else:
return initial_conditions[func_arg_val]
else:
return None
for arg in f.args:
if arg.is_number:
result= op(result, get_numeric_value(arg))
elif contains_function(arg):
r = evaluate_at(arg,n,initial_conditions)
if r:
result=op(result,r)
else:
return None
else:
result =op(result,get_numeric_value(arg.subs(arg.free_symbols.pop(),n)))
return result
known_values = dict(initial_conditions)
for n in range(4,11):
known_values[n] = evaluate_at(func1,n,known_values)
关于python - sympy 在不使用 rsolve 的情况下计算第 n 次重复,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/53748897/
我是 Java 新手,这是我的代码, if( a.name == b.name && a.displayname == b.displayname && a.linknam
在下面的场景中,我有一个 bool 值。根据结果,我调用完全相同的函数,唯一的区别是参数的数量。 var myBoolean = ... if (myBoolean) { retrieve
我是一名研究 C++ 的 C 开发人员: 我是否正确理解如果我抛出异常然后堆栈将展开直到找到第一个异常处理程序?是否可以在不展开的情况下在任何 throw 上打开调试器(即不离开声明它的范围或任何更高
在修复庞大代码库中的错误时,我观察到一个奇怪的情况,其中引用的动态类型从原始 Derived 类型更改为 Base 类型!我提供了最少的代码来解释问题: struct Base { // some
我正在尝试用 C# 扩展给定的代码,但由于缺乏编程经验,我有点陷入困境。 使用 Visual Studio 社区,我尝试通过控制台读出 CPU 核心温度。该代码使用开关/外壳来查找传感器的特定名称(即
这可能是一个哲学问题。 假设您正在向页面发出 AJAX 请求(这是使用 Prototype): new Ajax.Request('target.asp', { method:"post", pa
我有以下 HTML 代码,我无法在所有浏览器中正常工作: 我试图在移动到
我对 Swift 很陌生。我如何从 addPin 函数中检索注释并能够在我的 addLocation 操作 (buttonPressed) 中使用它。我正在尝试使用压力触摸在 map 上添加图钉,在两
我设置了一个详细 View ,我是否有几个 Nib 文件根据在 Root View Controller 的表中选择的项目来加载。 我发现,对于 Nibs 的类,永远不会调用 viewDidUnloa
我需要动态访问 json 文件并使用以下代码。在本例中,“bpicsel”和“temp”是变量。最终结果类似于“data[0].extit1” var title="data["+bpicsel+"]
我需要使用第三方 WCF 服务。我已经在我的证书存储中配置了所需的证书,但是在调用 WCF 服务时出现以下异常。 向 https://XXXX.com/AHSharedServices/Custome
在几个 SO 答案(1、2)中,建议如果存在冲突则不应触发 INSERT 触发器,ON CONFLICT DO NOTHING 在触发语句中。也许我理解错了,但在我的实验中似乎并非如此。 这是我的 S
如果进行修改,则会给出org.hibernate.NonUniqueObjectException。在我的 BidderBO 类(class)中 @Override @Transactional(pr
我使用 indexOf() 方法来精细地查找数组中的对象。 直到此刻我查了一些资料,发现代码应该无法正常工作。 我在reducer中尝试了上面的代码,它成功了 let tmp = state.find
假设我有以下表格: CREATE TABLE Game ( GameID INT UNSIGNED NOT NULL, GameType TINYINT UNSIGNED NOT NU
代码: Alamofire.request(URL(string: imageUrl)!).downloadProgress(closure: { (progress) in
我是一名优秀的程序员,十分优秀!