内点惩罚函数法子程.docVIP

下载本文档

5
0
约2.78千字
约 3页
2018-11-24 发布于河南
举报
版权申诉

内点惩罚函数法子程.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

内点惩罚函数法子程

#include hjfgf.c#include stdlib.hvoid gradient(double x[],double g[],int n){int i;??double af,f1,f2,dltx=0.000001;??for(i=0;in;i++)??g[i]=0;??f1=objf(x);??for(i=0;in;i++)??{af=*(x+i);? ?*(x+i)=af+dltx;? ?f2=objf(x);? ?g[i]=(f2-f1)/dltx;? ?*(x+i)=af;??}}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 dfpopt(double xx[],double h0,double eps,double epsg,int n){int i,j,k;??double ae,zcc;??double *s,*x,*ay[2],*df[2],*zd[2],*zc[2],*zh[2];??s=(double *)malloc(n*sizeof (double ));??x=(double *)malloc(n*sizeof (double ));??for(i=0;i2;i++)??{ay[i]=(double*)malloc(n*sizeof(double));? ?df[i]=(double*)malloc(n*sizeof(double));? ?zd[i]=(double*)malloc(n*sizeof(double));? ?zc[i]=(double*)malloc(n*sizeof(double));? ?zh[i]=(double*)malloc(n*n*sizeof(double));??}??for(i=0;in;i++)??*(ay[0]+i)=xx[i];??L1:? ?k=0;? ?for (i=0;in;i++)? ?for (j=0;jn;j++)? ? {*(zh[0]+i*n+j)=0;? ???if(i==j)? ?? ?*(zh[0]+i*n+j)=1;? ? }??for(i=0;in;i++)??{*(x+i)=*(ay[0]+i);? ?*(ay[1]+i)=*(ay[0]+i);??}??gradient(x,df[0],n);??for(i=0;in;i++)??{*(s+i)=0;? ?*(df[1]+i)=*(df[0]+i);? ?for(j=0;jn;j++)? ?*(s+i)=*(s+i)-(*(zh[0]+i*n+j))*(*(df[0]+j));??}??L2:??oneoptim(x,s,h0,epsg,n,ay[0]);??for(i=0;in;i++)? ?*(x+i)=*(ay[0]+i);??gradient(x,df[0],n);??ae=0;??for(i=0;in;i++)? ?ae=ae+(*(df[0]+i))*(*(df[0]+i));??if(ae=eps)? ?{for(i=0;in;i++)? ? *(xx+i)=*(x+i);? ? free(s);? ? free(x);? ? for(i=0;i2;i++)? ? {free(ay[i]);? ???free(df[i]);? ???free(zd[i]);? ???free(zc[i]);? ???free(zh[i]);? ? }? ? return(objf(xx));? ?}? ?if(k==n) goto L1;? ?zcc=0;? ?for(i=0;in;i++)? ?{*(zd[0]+i)=*(ay[0]+i)-(*(ay[1]+i));? ? *(zd[1]+i)=*(df[0]+i)-(*(df[1]+i));? ? *(df[1]+i)=*(df[0]+i);? ? zcc=zcc+(*(zd[0]+i))*(*(zd[1]+i));? ?}? ?for(i=0;in;i++)? ? f