惩罚函数的内点法.doc_第1页
惩罚函数的内点法.doc_第2页
惩罚函数的内点法.doc_第3页
惩罚函数的内点法.doc_第4页
惩罚函数的内点法.doc_第5页
已阅读5页，还剩4页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

2013-2014 (1) 专业课程实践论文 内点法 一、算法理论内点法总是从可行域的内点出发，并保持在可行域内进行搜索，因此这种方法适用于只有不等式约束条件的问题 内点法据图计算步骤：1 给定初，允许误差，初始参数缩小系数；2 以为初始点，求解问题 Min S.t. 3若则停，得近似解；否则令回2.二、算法框图求满足 从出发，求改变约束极值方法输出结果停三、算法程序clcm=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f0=zeros(1,50);syms x1 x2 e;m(1)=1;c=0.2;a(1)=2;b(1)=-3;f=x12+x22-e*(1/(2*x1+x2-2)+1/(1-x1); f0(1)=15; fx1=diff(f,x1); fx2=diff(f,x2);fx1x1=diff(fx1,x1);fx1x2=diff(fx1,x2);fx2x1=diff(fx2,x1);fx2x2=diff(fx2,x2);for k=1:100x1=a(k);x2=b(k);e=m(k);for n=1:100f1=subs(fx1);f2=subs(fx2);f11=subs(fx1x1);f12=subs(fx1x2);f21=subs(fx2x1);f22=subs(fx2x2);if(double(sqrt(f12+f22)=0.002)a(k+1)=double(x1);b(k+1)=double(x2);f0(k+1)=double(subs(f); break; else X=x1 x2-inv(f11 f12;f21 f22)*f1 f2; x1=X(1,1);x2=X(2,1);endendif(double(sqrt(a(k+1)-a(k)2+(b(k+1)-b(k)2)=0.001)&(double(abs(f0(k+1)-f0(k)/f0(k)=0.001) a(k+1) b(k+1) k f0(k+1)break;elsem(k+1)=c*m(k);endend4、 算法实现例1. 利用内点法求解 s.t 解：改变算法中f=x12+x22-e*(1/(2*x1+x2-2)+1/(1-x1); 回车完成结果复制粘贴代码，回车出现结果例2. 利用内点法求解 解：改变算法中f=x12-x22-e*(1/(x1+x2-3)+1/(-x1); 回车完成结果例3. 利用内点法求解 解：改变算法中f=2*x12-x22-e*(1/(x1+4*x2)+1/(x2-x1); 回车完成结果例4. 利用
人人文库> 全部分类> 应用文书 > 研究报告

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论