matlab 黄金分割算法

function y=fungold(a)
y=a^2+2*a;


end
function goldmin
clear;
clc;
a=-3;
b=5;
tolerance=0.01;
n=0;
i=100;
for k=1:i
    if(abs(b-a)<=tolerance)
        break;
    else
        a1=b-0.618*(b-a);
        a2=a+0.618*(b-a);
        if (fungold(a1)>=fungold(a2))
            a=a1;
        else
            b=a2;
            y=fungold((b+a)/2);
        end
        n=n+1;
    end
end
n,a0=(b+a)/2,y=fungold((b+a)/2)

            


end


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值