第12章 序列二次规划法

6.用SQP方法的Matlab程序求解下列优化问题:

(1)

%%% 目标函数 f(x) %%%%%%%%%%%
function f=f1(x)
% f=-pi*x(1)^2*x(2);
f = (x(1)-2)^4 + (x(1)-2*x(2))^2;
%%%% 目标函数 f(x) 的梯度%%%%%
function df=df1(x)
% df=[-2*pi*x(1)*x(2), -pi*x(1)^2]';
df = [4*(x(1)-2)^3+2*(x(1)-2*x(2)); -4*(x(1)-2*x(2))];%%%% 约束函数 %%%%%%
function [h,g]=cons(x)
% h=[pi*x(1)*x(2)+pi*x(1)^2-150];
h=[];
% g=[x(1);x(2)];
g = [x(1)^2 + x(2)];%% 约束函数 Jacobi矩阵%%%%
function [dh,dg]=dcons(x)
% dh=[pi*x(2)+2*pi*x(1), pi*x(1)];
dh = [];
% dg=[1 0; 0 1];
dg = [2*x(1), 1];

调用:

x0 = [0 0]';
mu0 = [ ];
lam0 = [0];
[x,mu,lam,val,k] = sqpm(x0,mu0,lam0);

结果

(2)

%%% 目标函数 f(x) %%%%%%%%%%%
function f=f1(x)
% f=-pi*x(1)^2*x(2);
% f = (x(1)-2)^4 + (x(1)-2*x(2))^2;
s = -x(1)-x(2);
f = exp(s) + x(1)^2 + 2*x(1)*x(2) + x(2)^2 + 2*x(1) + 6*x(2);
%%%% 目标函数 f(x) 的梯度%%%%%
function df=df1(x)
s = -x(1)-x(2);
df = [-exp(s)+2*x(1)+2*x(2)+2; -exp(s)+2*x(1)+2*x(2)+6];
%%%% 约束函数 %%%%%%
function [h,g]=cons(x)
% h=[pi*x(1)*x(2)+pi*x(1)^2-150];
h=[];
g = [2-x(1)-x(2),x(1),x(2)];
%% 约束函数 Jacobi矩阵%%%%
function [dh,dg]=dcons(x)
% dh=[pi*x(2)+2*pi*x(1), pi*x(1)];
dh = [];
dg = [-1 -1;1 0;0 1];

调用

x0 = [2 0]';
mu0 = [ ];
lam0 = [0 0 0];
[x,mu,lam,val,k] = sqpm(x0,mu0,lam0);

结果

《最优化方法及其Matlab程序设计》马昌凤 部分习题答案相关推荐

  1. 最优化方法及其matlab程序设计马昌凤,《马昌凤-最优化方法及其Matlab程序设计》.pdf...

    最优化方法及其 Matlab 程序设计 马昌凤 2009 年12 月 内容提要 本书较为系统地介绍了非线性最优化问题的基本理论和算法及其主要算法 的Matlab 程序设计. 主要内容包括 (精确或非精 ...

  2. matlab中sqp的算法原理_最优化方法及其Matlab程序设计 马昌凤

    内容简介<最优化方法及其Matlab程序设计>较系统地介绍了非线性最优化问题的基本理论和算法,以及主要算法的Matlab程序设计,主要内容包括(精确或非精确)线搜索技术.最速下降法与(修正 ...

  3. c++语言基础教程第三版习题答案,《C++语言基础教程》吕凤翥习题答案

    <<C++语言基础教程>吕凤翥习题答案>由会员分享,可在线阅读,更多相关<<C++语言基础教程>吕凤翥习题答案(66页珍藏版)>请在人人文库网上搜索. ...

  4. c语言二分答案最小值,C语言程序设计第2版课后习题答案(贾宗璞许合利人民邮电-第2版).doc...

    C语言程序设计第2版课后习题答案(贾宗璞许合利人民邮电-第2版) C语言习题答案 第1章 1.选择题 (1) C (2) B (3) 2.填空题 (1)main (2) main (3) 有穷性.确定 ...

  5. c语言程序设计第五版课后习题答案谭浩强第八章课后题

    c语言程序设计第五版课后习题答案谭浩强第八章课后题 输入三个字符串,按由大到小的顺序输出 要求用指针实现 输入10个整数,将其中最小的数与第一个数对换,把最大的数与最后一个属兑换.写三个函数:1.输入 ...

  6. 《C语言程序设计(第五版)》习题答案

    想看更多算法题,可以扫描上方二维码关注我微信公众号"数据结构和算法",截止到目前我已经在公众号中更新了500多道算法题,其中部分已经整理成了pdf文档,截止到目前总共有800多页( ...

  7. matlab张志涌版课后习题答案,matlab教程(张志涌)课后习题答案.doc

    matlab教程(张志涌)课后习题答案 matlab教程(张志涌)课后习题答案 导读:就爱阅读网友为您分享以下"matlab教程(张志涌)课后习题答案"资讯,希望对您有所帮助,感谢 ...

  8. c语言程序设计教程中国农业出版社答案,C语言程序设计教程杨路明课后习题答案北京邮电大学出版社.pdf...

    C语言程序设计教程杨路明课后习题答案北京邮电大学出版社 1.算法描述主要是用两种基本方法:第一是自然语言描述,第二是使用专用工具进行算法描述 2.c语言程序的结构如下: 1.c语言程序由函数组成,每个 ...

  9. java语言程序设计第二版 张思民_java语言程序设计(张思民编著)习题答案.doc

    java语言程序设计(张思民编著)习题答案.doc 1第1章JAVA概述[1]为什么说JAVA的运行与计算机硬件平台无关解答JAVA编译器能够产生一种与计算机体系结构无关的字节指令BYTECODE,只 ...

最新文章

  1. 利用Python实现十大经典排序算法(附代码流程)
  2. Java笔记14:泛型初探
  3. 我去,JS自执行匿名函数竟然有20几种写法!
  4. web前端html图片轮播,如何使用LayUI实现网页轮播图_WEB前端开发,layui,轮播图
  5. first() mysql_EF6配合MySQL或MSSQL(CodeFirst模式)配置指引
  6. 说说windows10自带浏览器Edge的好与不好
  7. Windows系统 services.msc命令详解,Windows命令行查看本地服务
  8. 判断一个整数是否是回文数
  9. silverlight安装后仍提示未安装_CAD提示“许可管理器不起作用或未正确安装”?来,我来教你方法...
  10. [转]Java 的强引用、弱引用、软引用、虚引用
  11. 捷联惯导系统学习4.2(捷联惯导误差方程)
  12. kettle案例21-将字段值设置为常量
  13. 架构师聊的四层代理和七层代理,都在聊什么?
  14. Win 10 系统下搭建 Python 编程环境,有手就行
  15. Linux挂载点的概念
  16. 【web】java后台获取checkbox、radio的值并输出
  17. 穿戴式心电信号采集系统设计(任务书+lunwen+答辩PPt+查重报告)
  18. Flask—静态资源配置
  19. PUT和POST的区别?
  20. 爬虫煎蛋网的mm图片

热门文章

  1. 三、C语言的数据类型—浮点型数据
  2. php游戏实例,php实现的简易扫雷游戏实例_PHP
  3. 中国象棋总体设计 中国象棋课设02
  4. 小程序发布上线流程_微信小程序发布流程:公司、个人如何注册小程序
  5. mapbox+threejs实现三维气泡球体专题图
  6. 计算机毕业设计php的毕业设计选题管理系统
  7. oracle 考试资料,Oracle考试复习题
  8. Android滑动浮层(滑动布局中使其中子布局一个浮动)
  9. helm部署jenkins到k8s并创建pipeline构建项目
  10. 自动写稿机器人下载,写稿机器人有哪些,写稿机器人软件下载