鱼C论坛

 找回密码
 立即注册
查看: 826|回复: 5

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

[复制链接]
发表于 2018-4-24 13:18:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
本帖最后由 Diligence_Pytho 于 2018-4-24 19:48 编辑

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
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

 楼主| 发表于 2018-5-8 20:17:33 | 显示全部楼层
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)
分享正确代码
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2018-4-24 13:27:46 | 显示全部楼层
报错信息?
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2018-4-24 13:56:42 | 显示全部楼层
初步看了下,26行if语句错误 if f(λ[k])-f(u[k])=0:   判定应该用 ==   

建议自己先排查下,根据调试信息进行修改
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2018-4-24 14:57:56 | 显示全部楼层

这个问题,哪里都有错,BUG太多,建议运行一下
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2018-4-24 14:59:35 | 显示全部楼层
thexiosi 发表于 2018-4-24 13:56
初步看了下,26行if语句错误 if f(λ[k])-f(u[k])=0:   判定应该用 ==   

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

这个程序很多错误。
乘号有错,你说的哪里也有错。
还有其他的一些错,我也弄不出来。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-3-29 08:52

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表