【東京大学】【優】最急降下法、共役勾配法のアルゴリズム
数理計画と最適化 第2回レポート
黄金分割法(1)
a=input('Enter a');
b=input('Enter b');
t=(sqrt(5)-1)/2;
c=0.000001;
r=a+(1-t)*(b-a);
s=a+t*(b-a);
p=[1 0 10 -199 1000];
while b-a>c
fr=polyval(p,r);
fs=polyval(p,s);
if fr>fs
a=r;
r=s;
s=a+t*(b-a);
else
b=s;
s=r;
r=a+(1-t)*(b-a);
endif
endwhile
a
polyval(p,a)
黄金分割法(2)
a=input('Enter a');
b=input('Enter b');
t=(sqrt(5)-1)/2;
c=0.000001;
r=a+(1-t)*(b-a);
s=a+t*(b-a);
p=[5 -6 5];
while b-a>c
fr=1-1/polyval(p,r);
fs=1-1/...