[计算机]Powell优化算法是利用仪器测井理建立误差函数.docVIP

Powell优化算法是利用仪器测井理建立误差函数（非相关函数），借助Powell方向加速法求出非相关函数达到最小时的解，对于气，水两相流动，从预设的气，水流量初始值出发，沿不同的广向进行搜索，可求出气，水两相流动中可能最大产量。与目前常用的生产测井解释方法相比，文中提出的方法具有精度高，实用性强等优点，在测井曲线有缺陷时，仍有可能得到较好的结果。 powell.c代码如下： CODE: ＃i nclude hjfgf.c double oneoptim(double x0[],double s[],double h0,double epsg,int n,double x[]) {double *a,*b,ff; a=(double *)malloc(n*sizeof(double)); b=(double *)malloc(n*sizeof(double)); jtf(x0,h0,s,n,a,b); ff=gold(a,b,epsg,n,x); free(a); free(b); return (ff); } double powell(double p[],double h0,double eps,double epsg,int n,double x[]) {int i,j,m; double *xx[4],*ss,*s; double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d; ss=(double *)malloc(n*(n+1)*sizeof(double)); s=(double *)malloc(n*sizeof(double)); for(i=0;in;i++) {for(j=0;j=n;j++) *(ss+i*(n+1)+j)=0; *(ss+i*(n+1)+i)=1; } for(i=0;i4;i++) xx[i]=(double *)malloc(n*sizeof(double)); for(i=0;in;i++) *(xx[0]+i)=p[i]; for(;;) {for(i=0;in;i++) {*(xx[1]+i)=*(xx[0]+i); x[i]=*(xx[1]+i); } f0=f1=objf(x); dlt=-1; for(j=0;jn;j++) {for(i=0;in;i++) {*(xx[0]+i)=x[i]; *(s+i)=*(ss+i*(n+1)+j); } f=oneoptim(xx[0],s,h0,epsg,n,x); df=f0-f; if(dfdlt) {dlt=df; m=j; } } sdx=0; for(i=0;in;i++) sdx=sdx+fabs(x[i]-(*(xx[1]+i))); if(sdxeps) {free(ss); free(s); for(i=0;i4;i++) free(xx[i]); return(f); } for(i=0;in;i++) *(xx[2]+i)=x[i]; f2=f; for(i=0;in;i++) {*(xx[3]+i)=2*(*(xx[2]+i)-(*(xx[1]+i))); x[i]=*(xx[3]+i); } fx=objf(x); f3=fx; q=(f1-2*f2+f3)*(f1-f2-dlt)*(f1-f2-dlt); d=0.5*dlt*(f1-f3)*(f1-f3); if((f3f1)||(qd)) {if(f2=f3) for(i=0;in;i++) *(xx[0]+i)=*(xx[2]+i); else for(i=0;in;i++) *(xx[0]+i)=*(xx[3]+i); } else {for(i=0;in;i++) {*(ss+(i+1)*(n+1))=x[i]-(*(xx[1]+i)); *(s+i)=*(ss+(i+1)*(n+1)); } f=oneoptim(xx[0],s,h0,epsg,n,x); for(i=0;in;i++) *(xx[0]+i)=x[i]; for(j=m+1;j=n;j