Question

我最近开始学习python，我课程中的问题之一是对割线方法进行编程。我在重新分配循环中的变量时遇到麻烦。

无论我如何编写while语句，我都会不断得到除以零的错误。我假设这是因为我的代码使循环成为当前的x_1 = x_0并进行计算 f（x_1）-f（x_0），而不是将旧的x_1用于x_0。我的尝试：

G = 6.6741*10**-11
r_e = 6371000
r_m = 1737100
M_e = 5.9722*10**24
M_m = 7.3420*10**22
R = 3.8440*10**8
w =  2.6617*10**-6

def f(x):
    return (G*M_e)/x**2 - (G*M_m)/(R-x)**2 - w**2 * x

x_0=2*10**8
x_1=2.2*10**8

i=1

while i<=10 or 0.99 < x_1/x_0 <1.01:
    x_1=float(x_1-f(x_1)*(x_1-x_0)/(f(x_1)-f(x_0)))
    i = i+1
    x_0=x_1

print (x_1)

错误：

runfile('//myfiles/vj284/dos/python/Coursework q1.py', 
wdir='//myfiles/vj284/dos/python')
Traceback (most recent call last):

  File "<ipython-input-2-b33827a1b929>", line 1, in <module>
    runfile('//myfiles/vj284/dos/python/Coursework q1.py', 
wdir='//myfiles/vj284/dos/python')

   File "C:\Program Files\Anaconda3\lib\site- 
packages\spyder\utils\site\sitecustomize.py", line 705, in runfile
    execfile(filename, namespace)

  File "C:\Program Files\Anaconda3\lib\site- 
packages\spyder\utils\site\sitecustomize.py", line 102, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)

  File "//myfiles/vj284/dos/python/Coursework q1.py", line 20, in <module>
    x_1=float(x_1-f(x_1)*(x_1-x_0)/(f(x_1)-f(x_0)))

ZeroDivisionError: float division by zero

Answer 1

我认为这是因为我的代码使循环成为当前的x_1=x_0和计算f(x_1)-f(x_0)，而不是将旧的x_1用于x_0。

不。这是因为循环x_0 = x_1中的最后一行。为什么覆盖x_0？另外，y的用途是什么？

编辑：由于添加了文本版本，因此更容易修复代码。只需使用一个临时变量即可保存“新”结果。另外，将or更改为and并添加条件以检查收敛性（如果收敛的方法x_0与x_1相同，则代码将因除法而崩溃减零）：

i = 1
while i <= 10 and 0.99 < x_1/x_0 <1.01 and abs((x_1 - x_0) / x_1) > 1e-15: :
    tmp = x_1 - f(x_1) * (x_1 - x_0) / (f(x_1) - f(x_0))
    i = i + 1
    x_0 = x_1
    x_1 = tmp

在循环条件下为什么有and 0.99 < x_1/x_0 <1.01并不明显。我认为您应该删除它：

i = 1
while i <= 10 and abs(x_1 - x_0) > 1e-15 * abs(x_1):
    tmp = x_1 - f(x_1) * (x_1 - x_0) / (f(x_1) - f(x_0))
    i = i + 1
    x_0 = x_1
    x_1 = tmp

最后一次编辑：为了避免使用临时变量，我将重写循环，如下所示：

for _ in range(10):
    if abs(x_1 - x_0) < 1e-15 * abs(x_1):  
        break  
    x_0, x_1 = x_1, x_1 - f(x_1) * (x_1 - x_0) / (f(x_1) - f(x_0))

此外，可能您应该检查函数值而不是其参数的紧密度：通常，由于舍入错误，f(x_1)可能等于f(x_0)，即使x_1 != x_0也是如此。或者，您可能要检查f(x_1) <= eps指示已找到根。

python中割线方法的麻烦写循环

1 个答案: