无约束多维优化牛顿法vb编程

1、无约束多维优化牛顿法vb编程阻尼牛顿法程序框图转出Vb编程运行后的界面Vb程序PrivateSubCommand1_Click()DimmlAsDoubleDimm2AsDoubleDimm3AsDoubleDimm4AsDoubleDimm5AsDoubleDimm6AsDoubleDimxxAsDoubleDimfxxAsDoubleDimyp1AsDoubleDimyp2AsDoubleDimx1AsDoubleDimx2AsDoubleDimeAsDoubleDime1AsDoubleDimhAsDoubleDimp0AsDoubleDimp1AsDoubleDimp2AsDouble

2、Dimp3AsDoubleDimyp3AsDoubleDimtidu(1,0)AsDoubleDimhaise(1,1)AsDoubleDimhaiseni(1,1)AsDoubleDimsiAsDoubleDims2AsDoubleDimxxlAsDoubleDimxx2AsDoubleDimxx3AsDoublem6=Val(lnputBox(请输入常数项m6,输入框,0)xi=Val(lnputBox(请输入迭代初始点xi,输入框,0)x2=Val(lnputBox(请输入迭代初始点x2,输入框,0)mi=Val(lnputBox(请输入m2=Val(lnputBox(请输入x1A2的系

3、数mi,输入框,0)x2A2的系数m2,输入框,0)m3=Val(lnputBox(请输入m4=Val(lnputBox(请输入x1*x2的系数m3,输入框,0)xi的系数m4,输入框,0)m5=Val(lnputBox(请输入x2的系数m5,输入框,0)&=Val(lnputBox(请输入迭代收敛精度,输入框,0)&i=Val(lnputBox(请输入一维搜索的收敛精度i,输入框,0)h=Val(lnputBox(请输入一维搜索的初始步长h,输入框,0)p0=Val(lnputBox(请输入一维搜索的初始点p0,输入框,0)Labeli.Caption=目标函数：f(X)=&mi&”xiA2

4、+&m2&”乂2人2+&m3&”xix2+&m4&xi+&m5&x2+&m6Label2.Caption=迭代初始点xi=&xi&”和x2=&x2&vbLf&”迭代收敛精度&Label3.Caption=一维搜索的初始点p0=&p0&vbLf&”一维搜索的初始步长h=&vbLf&”一维搜索的收敛精度&i=&itidu(0,0)=2*mi*xi+m3*x2+m4tidu(i,0)=2*m2*x2+m3*xi+m5haise(0,0)=2*mihaise(0,i)=m3haise(i,0)=m3haise(i,i)=2*m2haiseni(0,0)=2*m2/(4*mi*m2-m3a2)haise

5、ni(0,i)=-m3/(4*mi*m2-m3a2)haiseni(i,0)=-m3/(4*mi*m2-m3a2)haiseni(i,i)=2*mi/(4*mi*m2-m3a2)DoWhile(tidu(0,0)a2+(tidu(i,0)a2&si=-(haiseni(0,0)*tidu(0,0)+haiseni(0,i)*tidu(i,0)s2=-(haiseni(i,0)*tidu(0,0)+haiseni(i,i)*tidu(i,0)pi=p0 xxi=xi+pi*si:xx2=x2+pi*s2ypi=mi*xxia2+m2*xx2a2+m3*xxi*xx2+m4*xxi+m5*xx2+

6、m6p2=p0+hxx1=x1+p2*s1:xx2=x2+p2*s2yp2=ml*xx1A2+m2*xx2A2+m3*xx1*xx2+m4*xx1+m5*xx2+m6Ifyp1yp2Thenh=2*hElseh=-0.25*hp3=p1:yp3=yp1p1=p2:yp1=yp2p2=p3:yp2=yp3EndIfp3=p0+hxx1=x1+p3*s1:xx2=x2+p3*s2yp3=ml*xx1a2+m2*xx2a2+m3*xx1*xx2+m4*xx1+m5*xx2+m6DoWhileyp2=yp3h=2*hpl=p2:yp1=yp2p2=p3:yp2=yp3p3=p0+hxx1=x1+p3*

7、s1:xx2=x2+p3*s2yp3=ml*xx1a2+m2*xx2a2+m3*xx1*xx2+m4*xx1+m5*xx2+m6LoopIfh0Thena=p1:ypa=yp1b=p3:ypb=yp3Elsea=p3:ypa=yp3b=p1:ypb=yp1EndIfpl=a+0.382*(b-a)p2=a+0.618*(b-a)xx1=x1+pl*s1:xx2=x2+pl*s2yp1=ml*xx1a2+m2*xx2a2+m3*xx1*xx2+m4*xx1+m5*xx2+m6xx1=x1+p2*s1:xx2=x2+p2*s2yp2=ml*xx1a2+m2*xx2a2+m3*xx1*xx2+m4*

8、xx1+m5*xx2+m6DoWhileAbs(b-a)1Ifyp1yp2Thenn0=0b=p2p2=p1:yp2=yp1Elsen0=1a=plpl=p2:yp1=yp2EndIfIfnO=0Thenpl=a+0.382*(b-a)xx1=x1+pl*s1:xx2=x2+pl*s2yp1=ml*xx1A2+m2*xx2A2+m3*xx1*xx2+m4*xx1+m5*xx2+m6Elsep2=a+0.618*(b-a)xx1=x1+p2*s1:xx2=x2+p2*s2yp2=ml*xx1a2+m2*xx2a2+m3*xx1*xx2+m4*xx1+m5*xx2+m6EndIfLoopxx=(a+b)/2x1=x1+xx*six2=x2+xx*s2fxx=mi*x1a2+m2*x2a2+m3*x1*x2+m4*x1+m5*x2+m6tidu(O,0)=2*mi*x1+m3*x2