鲍威尔法Word格式文档下载.docx

鲍威尔法Word格式文档下载.docx

《鲍威尔法Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《鲍威尔法Word格式文档下载.docx（7页珍藏版）》请在冰豆网上搜索。

for（i=0;

i<

N;

i++）

z[i]=x[i]+b*s[i];

f[1]=gs（z）;

if（f[1]<

f[0]）

{do{t=2*t;

b=b+t;

f[0]=f[1];

if（f[0]>

f[1]）

a=b-t;

elsebreak;

}

while（f[0]>

f[1]）;

else

{t=0-t;

do

{a=a+t;

f[1]=f[0];

i++）z[i]=x[i]+a*s[i];

f[0]=gs（z）;

if（f[1]>

{b=a-t;

t=2*t;

while（f[1]>

f[0]）

}

do

{j=j+1;

A[0]=b-q*（b-a）;

z[i]=x[i]+A[0]*s[i];

A[1]=a+q*（b-a）;

z[i]=x[i]+A[1]*s[i];

do{if（f[0]>

{a=A[0];

A[0]=A[1];

j=j+1;

if（j>

50）

break;

{b=A[1];

A[1]=A[0];

if（（f[1]-f[0]）/f[1]>

0）

e=（f[1]-f[0]）/f[1];

elsee=0-（f[1]-f[0]）/f[1];

}

while（e>

c）

if（（A[1]-A[0]）/A[0]>

e=（A[1]-A[0]）/A[0];

elsee=0-（A[1]-A[0]）/A[0];

if（e<

=d）break;

b=A[1];

d）

if（f[0]<

{for（i=0;

return（f[0]）;

return（f[1]）;

voidmain（）

{floatxo[N],ss[N][N],z[N],c,d,t,fmin,xooo[N],xoo[N],fi[N],s[N],x[N],xe[N],f[N],fio,dm,dml,df,dx,dxl;

inti,j,k,ks,m;

{printf（"

请输入初始点x[%d]="

i）;

scanf（"

%f"

&

xo[i]）;

printf（"

请输入迭代精度Ef,Ex="

）;

%f,%f"

c,&

d）;

请输入步长T="

t）;

k=0;

ks=0;

xooo[i]=xo[i];

for（j=0;

j<

j++）

{if（i==j）ss[i][j]=1;

elsess[i][j]=0;

xoo[i]=xooo[i];

k=k+1;

{

{for（j=0;

s[j]=ss[i][j];

fmin=ywyh（xooo,t,s,z）;

j++）x[j]=z[j];

fi[i]=fmin;

xooo[j]=x[j];

{s[i]=x[i]-xoo[i];

xe[i]=2*x[i]-xoo[i];

f[0]=gs（xoo）;

f[1]=fmin;

f[2]=gs（xe）;

fio=f[0];

dm=fio-fi[0];

m=0;

N-1;

{dml=fi[i]-fi[i+1];

if（dm<

dml）

{dm=dml;

m=i+1;

else;

}if（（f[0]+f[2]-2*f[1]）*（f[0]-f[1]-dm）*（f[0]-f[1]-dm）<

0.5*dm*（f[0]-f[2]）*（f[0]-f[2]）&

&

f[2]<

f[0]）

{fmin=ywyh（xooo,t,s,z）;

f

or（j=0;

x[j]=z[j];

{if（i<

m）

ss[i][j]=ss[i][j];

ss[i][j]=ss[i+1][j];

ss[N-1][i]=s[i];

ks=ks+1;

{if（fmin>

f[2]）

i++）

x[i]=xe[i];

f[1]=f[2];

fmin=f[2];

xooo[i]=x[i];

if（（f[0]-f[1]）/f[1]>

df=（f[0]-f[1]）/f[1];

elsedf=0-（f[0]-f[1]）/f[1];

if（xoo[0]-x[0]>

0）dx=xoo[0]-x[0];

else

dx=x[0]-xoo[0];

for（i=1;

{if（xoo[i]-x[i]>

dxl=xoo[i]-x[i];

elsedxl=x[i]-xoo[i];

if（dx<

dxl）dx=dxl;

while（dx>

c||df>

z[i]=x[i];

\n最优点坐标x[%i]=%f\n"

i,z[i]）;

\n最优函数值f（x）=%f\n"

fmin）;

printf（"

\n搜索累计次数k=%d\n"

k）;

\n替换方向累计次数ks=%d\n"

ks）;

}