資料紹介

【東京大学】【優】黄金分割法、２分割法のアルゴリズム

資料の原本内容 ( この資料を購入すると、テキストデータがみえます。 )

数理計画と最適化第３回レポート　　　　　　　　　　　　　　　　
(A-1)
・最急降下法
b = -[-27 ;-42 ;32 ;-23]
c = 163
k = 0
x = input(' Enter matrix x ')
s = -A*x + b
while norm(s)>0.0001
f = x'*A*x/2 - x'*b + c
s = -A*x + b
a = (s'*b -s'*A*x)/(s'*A*s)
x = x + a*s
k = k + 1
endwhile
printf("f = %f\n",f)
printf("k = %f\n",k)
(初期条件に対する解)
・共役勾配法
b = -[-27 ;-42 ;32 ;-23]
c = 163
k = 0
x = input(' Enter matrix x ')
s = -A*x + b
while norm(-A*x + b)>0.0001
f = x'*A*x/2 - x'*b + c
a = (b'*s -x'*A*s)/(s'*A*s)
t = x
...

コメント0件