tl / dr:我有一个numpy边界/初始值问题,想看看我是否正确地接近这个问题。我很擅长numpy。我正在介绍问题的简化版本。
我为a
和b
的整数值定义了2个函数t
和x
,我正在尝试计算正x
和t
(说到N
)。我想找出用numpy做这个的最好方法。
我有t=0
和x=0
的边界值,a(t,x)
仅取决于a(t-1,x-1)
和b(t-1,x-1)
而b(t,x)
取决于很多值具有较小a
,t
的{{1}}。这就是它“简单”的原因。我们有
x
和a=1
的 t=0
。
x=0
代表b=0.1
,t=0
代表b=1
。在x=0
,我们得到x=t=0
。
在内部,b=0.1
。
现在困难的部分。 a(t,x) = a(t-1,x-1) - b(t-1,x-1)
对于某些函数b(t,x) = a(t-1,x-1) S(t, t-1) + a(t-2,x-2) S(t,t-2) + ...
,S(t,y)
的总和等于f(a(t-1,1)) + f(a(t-1,2)) + ... + f(a(t-1,y))
(如果您需要特定的内容,则可以认为它只是f
)。
所以我的计划基本上是这样做的:
a + a**2
我的问题:这是一种合理的方法 - 我可以删除任何这些循环,还是应该完全采用不同的方法?
奖金问题:numpy newb是否可能出现错误以进行编码?