fortran 黄金分割法求极小值编程.docxVIP

programhjfgimplicit noneinteger n ,hreal x10,x20,x1,x2,y1,y2,a,b,c1,c2,d1,d2,jwrite(*,*)运算次数为n write(*,*)初始点为(2,2). write(*,*)初始探索区间为(1.0，7.0) write(*,*)初始一维探索方向为(0.7071,0.7071) write(*,*)初始迭代精度为0.1x10=2.0x20=2.0 d1=0.7071 d2=0.7071 j=0.1 a=1.0 b=7.0do n=0c1=a+0.382*(b-a)x1=x10+d1*c1x2=x20+d2*c1 y1=x1*x1+x2*x2-8*x1-12*x2+52 c2=a+0.618*(b-a)x1=x10+d1*c2x2=x20+d2*c2 y2=x1*x1+x2*x2-8*x1-12*x2+52do while (ABS(a-b)j) n=n+1c1=a+0.382*(b-a)x1=x10+d1*c1x2=x20+d2*c1 y1=x1*x1+x2*x2-8*x1-12*x2+52 c2=a+0.618*(b-a)x1=x10+d1*c2x2=x20+d2*c2 y2=x1*x1+x2*x2-8*x1-12*x2+52if(y1y2) then a=c1 c1=c2 y1=y2 c2=a+0.618*(b-a)x1=x10+d1*c2x2=x20+d2*c2 y2=x1*x1+x2*x2-8*x1-12*x2+52elseif(y1y2) then b=c2 c2=c1 y2=y1 c1=a+0.382*(b-a)x1=x10+d1*c1x2=x20+d2*c1 y1=x1*x1+x2*x2-8*x1-12*x2+52elseif(y1==y2) then a=c1 b=c2endif write(*,(3X,a3,I3,//6(3X,a,F12.6,5X//)))n=,n,x1=,x1,y1=,y1,x2=,x2,y2=,y2,a=,a,b=,benddo write(*,*)运算结果为： write(*,(3X,a3,I3,//6(3X,a,F12.6,5X//)))n=,n,x1=,x1,y1=,y1,x2=,x2,y2=,y2,a=,a,b=,b10 write(*,*)改变初始点请输入1write(*,*)改变一维探索方向请输入2write(*,*)改变迭代精度请输入3write(*,*)改变探索区间请输入4write(*,*)停止请输入5read(*,*) hif (h==1) then write(*,*) 请输入两个实数:read(*,*) x10,x20 d1=0.7071 d2=0.7071 j=0.1 a=1.0 b=7.0elseif (h==2) then write(*,*) 请输入两个实数:read(*,*) d1,d2x10=2.0x20=2.0 j=0.1 a=1.0 b=7.0elseif (h==3) then write(*,*) 请输入一个实数:read(*,*) jx10=2.0x20=2.0 d1=0.7071 d2=0.7071 a=1.0 b=7.0elseif(h==4) then write(*,*)请输入两个实数read(*,*) a,bx10=2.0x20=2.0 d1=0.7071 d2=0.7071 j=0.1elseif (h==5) thenexitelse write(*,*)请重新输入goto 10endifenddoendprogramecczimplicit nonereal x1,x2,x3,xp,f1,f2,f3,fp,j,c1,c2integern,ij=0.4x1=1.0x3=7.0x2=0.5*(x1+x3)f1=(x1-3)**2f2=(x2-3)**2f3=(x3-3)**2write(*,*)初始精度为0.4write(*,*)初始区间为(1,7)don=0write(*,*)jdo while (ABS(x1-x3)j)c1=(f3-f1)/(x3-x1)c2=((f2-f1)/(x2-x1)-c1)/(x2-x3)xp=0.5*(x1+x3-c1/c2)fp=(xp-3)**2 n=n+1if (xp=x2) thenif(fp=f2) thenx3=xp