中国DOS联盟

'非对称3-4-5加速逼近
dim a,b,h,i
a=0
b=3141592.6535893
h=1
k=2
do while abs(a-b)>0.0000001
      do while a<b
         a=a+k*h:h=k*h
i=i+1
      loop
      if a>b then
      k=1:h=0.5*h:a=a-h
i=i+1
      end if
      do while a>b
         a=a-k*h:h=k*h
i=i+1
      loop
      if a<b then
      k=1:h=0.5*h:a=a+h
i=i+1
      end if
loop
msgbox(a&"逼近次数:"&i)

'y=x快速折半逼近
'只要单调，且有根，直接算出根
a=2000.5:b=211.9456
h=1 '初始步长为1
i=0

Do While Abs(a-b)>0.001
if (a-b)*(a+Sgn(b-a)*h)>0 then
   h=2*h
else
   h=h/2
end if
i=i+1
a=a+Sgn(b-a)*h
'msgbox a
Loop
MsgBox(a&"逼近次数:"&i)

'y=x快速折半逼近
'只要单调，且有根，直接算出根
a=0.1
b=3.141515926
h=1 '初始步长为1
i=0
fa=log(a)  '************log
faa=fa+1
Do While Abs(fa-b)>0.0000001
if (fa-b)*(faa-b)>0 then
   h=2*h
else
   h=h/3
end if
i=i+1
a=a+Sgn(b-fa)*h
faa=fa
fa=log(a) '*************log
'msgbox(a&"--"&fa)
'msgbox fa
Loop
MsgBox(a&"逼近次数:"&i)

'优化对分法，
'单调递减函数
h=0.5
k=2
p=1.2
a=0
i=0
tmp=a
if 0.5^(a)<>p then
a=a+h
   if 0.5^(a)>p then
      do while 0.5^(a)>p
        tmp=a:a=a+h:h=k*h '递减，+h
        i=i+1
      loop
   else
      do while 0.5^(a) <p
         tmp=a:a=a-h:h=k*h  '递减，+h
         i=i+1
      loop
   end if
end if

b=tmp

do while abs(a-b)>0.0001
    c=(a+b)/2
       i=i+1
    if (0.5^(c)-p)*(0.5^(a)-p)>0 then
        a=c
    else
        b=c
    end if
loop
c=(a+b)/2
MsgBox(c&"逼近次数:"&i)