利用Python求解黄金分割法求解的问题

Diligence_Pytho

e=float(input('请输入收敛精度e='))
a=int(input('请输入区间左值='))
b=int(input('请输入区间右值='))
M=(5**0.5-1)/2
n=1
def f(x):
  return x**2+2*x
print("n      a    \t b    \t x1    \t x2    \t y1    \t     y2\n")

def loop():
    x1=a+(1-M)*(b-a)
    x2=a+M*(b-a)
    y1=f(x1)
    y2=f(x2)
    print('%d %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n'%(n,a,b,x1,x2,y1,y2))
    n=n+1

    while True:
        if (b-a)<e:
            x=(b+a)/2
            break
        
        if (y1-y2)<0:
            b=x2
            x2=x1
            y2=y1
            x1=a+(1-M)*(b-a)
            y2=y1
            y1=f(x1)
            print("%d\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n",n,a,b,x1,x2,y1,y2)
            n=n+1

        elif (y1-y2)>0:
            a=x1
            x1=x2
            x2=a+M*(b-a)
            y1=y2
            y2=f(x2)
            print("%d\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n",n,a,b,x1,x2,y1,y2)
            n=n+1

        else:
            a=x1
            b=x2
            loop()
            
loop()


显示错误：
Traceback (most recent call last):
  File "G:/Python3.6.3（32-bit）/HJ.py", line 47, in <module>
    loop()
  File "G:/Python3.6.3（32-bit）/HJ.py", line 11, in loop
    x1=a+(1-M)*(b-a)
UnboundLocalError: local variable 'a' referenced before assignment

题目是  
用黄金分割法求解
min λ2+2λ  -3≤λ≤5

Diligence_Pytho

e=1e-08
a=int(input('请输入区间左值='))
b=int(input('请输入区间右值='))
M=(5**0.5-1)/2
n=1
def f(x):
  return x**2+2*x
print("n      a    \t b    \t x1    \t x2    \t y1    \t     y2\n")

def loop(a,b,e,M,n):
    x1=a+(1-M)*(b-a)
    x2=a+M*(b-a)
    y1=f(x1)
    y2=f(x2)
    print('%d %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n'%(n,a,b,x1,x2,y1,y2))
    n=n+1

    while True:
        if (b-a)<e:
            x=(b+a)/2
            break
        
        if (y1-y2)<0:
            b=x2
            x2=x1
            y2=y1
            x1=a+(1-M)*(b-a)
            y2=y1
            y1=f(x1)
            print('%d %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n'%(n,a,b,x1,x2,y1,y2))
            n=n+1

        elif (y1-y2)>0:
            a=x1
            x1=x2
            x2=a+M*(b-a)
            y1=y2
            y2=f(x2)
            print('%d %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\t %.4lf\n'%(n,a,b,x1,x2,y1,y2))
            n=n+1

        else:
            a=x1
            b=x2
            loop(a,b,e,M,n)
            
loop(a,b,e,M,n)
分享正确代码

BngThea

报错信息？

thexiosi

初步看了下，26行if语句错误 if f(λ[k])-f(u[k])=0:   判定应该用 ==   

建议自己先排查下，根据调试信息进行修改

Diligence_Pytho

BngThea 发表于 2018-4-24 13:27
报错信息？

这个问题，哪里都有错，BUG太多，建议运行一下

Diligence_Pytho

thexiosi 发表于 2018-4-24 13:56
初步看了下，26行if语句错误 if f(λ[k])-f(u[k])=0:   判定应该用 ==   

建议自己先排查下，根据调试信 ...

这个程序很多错误。
乘号有错，你说的哪里也有错。
还有其他的一些错，我也弄不出来。