最优化程序设计 内点法_第1页
最优化程序设计 内点法_第2页
最优化程序设计 内点法_第3页
最优化程序设计 内点法_第4页
最优化程序设计 内点法_第5页
全文预览已结束

下载本文档

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

文档简介

1、北方民族大学运筹学内点法课程设计学院：信息与计算科学学院班级：08级统计（2 ）班姓名： 吕慧兵学号：20080546指导老师：王蜻设计日期：2011.6.12011.6.8一）内点的应用（1）惩罚函数定义于可行域内,序列迭代点在可行域内不断趋于约 束边界上的最优点.（2）只适合求解具有不等式约束的优化问题例子min X12 + % 22st 2 x + x - 2 2-七+1 0二）实验目的通过内点法的学习让我们掌握利用罚函数解决线性规划为解决相应 问题的一种思路与策略。实验思路与步骤【实验编译程序】clcm=zeros(1,50);a=zeros(1,50);b=zeros(1,50);f

2、0=zeros(1,50);syms x1 x2 e;m(1)=1;c=0.2;a(1)=2;b(1)=-3;f=x1 2+x2 2-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:100 x1=a(k);x2=b(k);e=m(k);for n=1:100f1=subs(fx1);f2=subs(fx2);f11=subs(fx1x1);f12

3、=subs(fx1x2);f21=subs(fx2x1);f22=subs(fx2x2);if(double(sqrt(f2+f22)=0.002)a(k+1)=double(x1);b(k+1)=double(x2);f0(k+1)=double(sub s(f);break;elseX=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)kf0(k+1)break;elsem(k+1)=c*m(k);endend四）【实验结论】9. OQCOWE =-7. h =34ans Ch 9&99利用Matlab作出上述结果，最优解为x*（，1

人人文库> 全部分类> 图纸下载 > 毕业设计

温馨提示

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

最新文档

评论

0/150

提交评论