昨天晚上在写educoder的实训作业的时候,发现了一些问题。
大概是这个这样的用二分法求函数的根的问题。还需要将根逐行输出。
首先二分法是通过不断寻找中的数与目标数进行比较的方法,运用于函数需要注意左右端点的正负号的问题。
这个题目分析可知,分别在(0,1),(4,5)内存在零点,而函数在0点的值大于0,在4点的值小于0,因此在二分法的左右端点的变化中需要注意。
正确代码如下
import numpy as np
E=1e-6 #表示10的负六次方
def f(x):
return 6*np.exp(x) -113*x +17
def binary_search1(i,j): #这是左边大于0的情况
while i<=j:
k=(i+j)//2
if abs(f(k))<E:
return k
if f(k)<-E:
j=k
else:
i=k #可以直接用else 因为是上面是直接进行下来的
def binary_search2(i,j): #这是右边大于0的情况
while i<=j:
k=(i+j)//2
if abs(f(k))<E:
return k
if f(k)<-E:
i=k
else:
j=k
print('%.4f'%binary_search1(0,1))
print('%.4f'%binary_search2(4,5))
但我犯了错误,在于想直接逐行输出解。
def div_wrong(m,n):
while m<n:
a=(m+n)/2
if abs(f(a)) <= E:
print('{:.4f}'.format(a)) #遇到return会立马退出函数体 所以用return写函数是最好的
if f(a)< -E:
m=a
else:
n=a
将return写为了print,结果是一直输出同一解。
所以在函数中,如果碰到了return是会直接跳出函数体的。另外要多注意函数的解的二分法求解与一般的二分法不同,在于左右端点的正负值不同,区间的左右赋值会出现不同的情况。