外点法求解优化问题.docVIP

实验报告 题目：序列二次规划求解优化问题 计算机程序（matlab环境） 班级： 数学系09级5班 组成员： 任楚 刘玉 王 优化问题实验报告： 1题目 运用序列二次规划求解优化问题 2问题叙述 用外点法来求解优化问题 Min f(x)=x^2 s.t –x^2=0 x-1=0 3方法介绍 序列二次规划： 4 程 序 代 码 Matlab 环 境 程 序 代 码 运行环境为matlab function [x,i]=ec(x0,a,b,esp) f=@(x)(x^2) g1=@(x)(-x^2) g2=@(x)(x-1) i=0; d=0; i=0:100; y=0:1; for x(i)=x0:1:x(i+1)==x(i)|| x(i+1)-x(i)==esp; f=-y; g1=grad(g1(x(i)))*d+g1(x(i)); g2=grad(g2(x(i)))*d+g2(x(i))*y; y=0:1; [fval]=fmincon(f,A,b); f=1/2*d*hessian(f(x(i))); g1=grad(g1(x(i)))*d+g1(x(i)); g2=grad(g2(x(i)))*d+g2(x(i))*(-fval); [d]=quadprog(H,c,A,b,x(0)); x(i+1)=x(i)+d; i=i+1; end x=x(i); 5 结 果 分 析 由数学知识可知上述问题的最优解经过用序列二次规划法编程求的函数f（x）在条件–x^2=0 x-1=0下，的最优解。由于程序运行错误，所以无结果。