0.618法求F(X)=x+20/x

题目:已知目标函数F(X)=x+20/x,取初始区间[a,b]=[2,10],ε=0.01。分别用0.618法和二次插值法求最优解。

public class Main {
    
	public static void main(String[] args) {
		//初始数据
		double a=2.0;
        double b=10.0;
        double miu=0.01;
        f_618 h=new f_618(a,b,miu);
        h.SuanFa();
	}
	
}
class f_618  {
    double a;
    double b;
    double miu;
    double x1;
    double x2;
    
	public f_618(double a,double b,double miu){
        this.a=a;
        this.b=b;
        this.miu=miu;
        this.x1=0.0;
        this.x2=0.0;
    }
    public void SuanFa() {
        int k=0;
        fun f=new fun();
        while(true) {
            k++;
            x1=a+0.382*(b-a);
            x2=a+0.618*(b-a);
            System.out.println("x1:"+x1+"x2:"+x2);
            if((b-a)<miu) break;
            else {
                if(f.F(x1)>f.F(x2)) {
                    System.out.println("左边f(x1):"+f.F(x1)+"f(x2):"+f.F(x2));
                    a=x1;
                    x1=x2;
                    x2=a+0.618*(b-a);
                }
                else {
                    System.out.println("右边f(x1):"+f.F(x1)+"f(x2):"+f.F(x2));
                    b=x2;
                    x2=x1;
                    x1=a+0.382*(b-a);
                }
            }
            System.out.println("迭代了"+k+"次,区间[a,b]=["+a+","+b+"]"+'\n'+"f(a)="+f.F(a)+'\n'+"f(b)="+f.F(b));
        }
        double f1=f.F(a);
        double f2=f.F(b);
        System.out.println("迭代了"+k+"次,最优区间[a,b]=["+a+","+b+"]"+'\n'+"f(a)="+f1+"f(b)="+f2);
    }
}
class fun  {
	public double F(double x) {
		//目标函数
	    x=x+(20.0/x);
        return x;
    }
}

答案
左边f(x1):8.944272678017402f(x2):8.944272608908868
迭代了14次,区间[a,b]=[4.470283049462319,4.479765675814997]
f(a)=8.944272678017402
f(b)=8.944284904571628
x1:4.473905412729042x2:4.476143312548274
迭代了15次,最优区间[a,b]=[4.470283049462319,4.479765675814997]
f(a)=8.944272678017402f(b)=8.944284904571628

参考:https://blog.csdn.net/weixin_49858823/article/details/118915862

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值