Question

我在for循环中使用scipy.optimize包中的bisection方法。我们的想法是使用“eps_komp”向量中每个元素（值）的二分法获得“sig”值。我编码很多：

import numpy as np
import scipy.optimize as optimize

K=300
n = 0.43
E = 210000
Rm = 700
sig_a = []
RO_K = 300
RO_n = 0.43

eps_komp =     [0.00012893048999999997,
0.018839115269999998,
0.01230539995,
0.022996934109999999,
-0.0037319012899999999,
0.023293921169999999,
0.0036927752099999997,
0.020621037629999998,
0.0063656587500000002,
0.020324050569999998,
-0.0025439530500000001,
0.018542128209999998,
0.01230539995,
0.019730076449999998,
0.0045837363899999999,
0.015275270549999997,
-0.0040288883499999999,
0.021215011749999999,
-0.0031379271699999997,
0.023590908229999999]

def eps_f(i):
    return eps_komp[i]

for j in range(len(eps_komp)):
    eps_komp_j = eps_f(j)
    if j <= len(eps_komp):
        def func(sig):
            return eps_komp_j - sig/E - (sig/RO_K)**(1/RO_n) 
        sig_a.append(optimize.bisect(func, 0, Rm))
    else:
        break

print(sig_a)

现在，如果我将eps_f（j）中“j”的值更改为0：

eps_komp_j = eps_f(0)

它有效，所以它对我手动插入的所有其他值都有效，但是如果我保持它就像在for循环中那样，“j”值不会自动改变而且我得到一个错误：

f（a）和f（b）必须有不同的符号

有谁知道问题是什么，如何解决？

此致

→

P.S。昨天我确实发布了关于这个问题的另一个话题，但我对这个问题并不十分具体，并得到负面反馈。但是，我今天确实需要解决这个问题，因此我不得不再次发布它，但是我确实设法让我的代码更进一步，然后我在之前的帖子中做了，所以它不是重新发布... < / p>

Answer 1

如果您read the docs，您会发现：

基本二分程例程，在参数a和b之间找到函数f的零。 f（a）和f（b）不能有相同的符号。慢但确定。

在您的代码中：

    def func(sig):
        return eps_komp_j - sig/Emod - (sig/RO_K)**(1/RO_n) 
    sig_a.append(optimize.bisect(func, 0, Rm))

您正在传递func(0)和func(700)。

将optimize.bisect行替换为print(func(0), func(700))，我得到以下输出：

0.00012893048999999997 -7.177181168628421
0.018839115269999998 -7.158470983848421
0.01230539995 -7.165004699168421
0.02299693411 -7.15431316500842
-0.00373190129 -7.1810420004084206
0.02329392117 -7.154016177948421
0.0036927752099999997 -7.173617323908421
0.02062103763 -7.156689061488421
0.00636565875 -7.17094444036842
0.02032405057 -7.156986048548421
-0.00254395305 -7.17985405216842
0.018542128209999998 -7.15876797090842
0.01230539995 -7.165004699168421
0.019730076449999998 -7.157580022668421
0.00458373639 -7.172726362728421
0.015275270549999997 -7.162034828568421
-0.00402888835 -7.181338987468421
0.02121501175 -7.156095087368421
-0.0031379271699999997 -7.1804480262884205
0.02359090823 -7.153719190888421

请注意具有相同符号的多对。 optimize.bisect无法处理这些问题。我不知道你尝试完成了什么，但这是错误的做法。

Python：函数在for循环中没有收到值

1 个答案: