wolf_powell线搜索方法

[复制链接]

字体大小: 正常放大

ilikenba

1万主题	49 听众	2万积分

TA的每日心情

	奋斗 2024-6-23 05:14

签到天数: 1043 天

[LV.10]以坛为家III

群组: 万里江山

群组: sas讨论小组

群组: 长盛证券理财有限公司

群组: C 语言讨论组

群组: Matlab讨论组

电梯直达

1^#

发表于 2004-4-30 11:21 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

> function wolf_powell(A,b,dir,gradt,x) result(wp_result)
real,dimension(

,intent(in)::b,dir,gradt,x
real,dimension(:,

,intent(in)::A
real::a1,a2,f1,f2,a0,p,q,t,r,wp_result,gradt1(size(gradt)) ,r1
t=2.0
a0=1
p=0.1
q=0.9
a1=0
a2=100000
r=dot_product(dir,gradt)
f1=f(x,A,b)
10 f2=f(x+a0*dir,A,b)
if(f2>f1+p*r*a0)then
   a2=a0
  a0=a1+0.5*(a0-a1)/(1+(f1-f2)/(a0-a1)*r)
  goto 10
else
   gradt1=matmul(A,x+a0*dir)+b
  r1=dot_product(gradt1,dir)
  if(r1>=q*r)then
   wp_result=a0
  else
   a1=a0
f1=f2
r=r1
a0=a0+(a0-a1)*r1/(r-r1)
goto 10
  endif
endif
end function
</P>
<

>本程序由Fortran 90语言编写，在Visual Fortran 5上编译通过！本程序由沙沙提供！</P>

zan