实验一MATLAB的优化工具箱介绍及无约束优化问题求解.docVIP

实验一 MATLAB的优化工具箱介绍及无约束优化问题求解 实验目的 1、熟悉Matlab的基本操作命令和函数表达； 2、了解Matlab优化工具箱的基本功能； 3、掌握优化工具箱中常用的函数命令格式； 4、熟练掌握使用fminbnd和fminunc函数求解无约束优化问题。 实验原理 1、矩阵运算基本原理； 2、无约束优化问题的求解算法，如最速下降法（共轭梯度法）算法、牛顿法算法、拟牛顿法等等； 实验环境 PC一台，Windows XP系统, Matlab软件（7.0版本以上）。 实验要求 1、练习matlab的基本操作； 2、认识matlab优化工具箱； 3、完成无约束优化问题的求解操作。 实验内容 1、Matlab简介 MATLAB作为线性系统的一种分析和仿真工具，是理工科大学生应该掌握的技术工具，它作为一种编程语言和可视化工具，可解决工程、科学计算和数学学科中许多问题。 MATLAB建立在向量、数组和矩阵的基础上，使用方便，人机界面直观，输出结果可视化。 Matlab操作窗口如下图： 2、Matlab优化工具箱 工具箱是Matlab的重要组成部分，用于扩充数值计算、符号运算功能、图形建模仿真功能、文字处理功能和相关专业领域功能，如信号处理、决策优化等等，能够用于多种学科的研究计算。目前，Matlab工具箱的种类繁多，其中优化工具箱（Optimization Toolbox）在解决最优化问题上起着极其重要的作用。主要的函数包括： 　　◆ 最小化函数 　　◆ 方程求解函数 　　◆ 最小二乘（曲线拟合）函数 　　◆ 实用函数 　　◆ 大型方法的演示函数 ◆ 中型方法的演示函数 本次试验主要讨论最小化函数中的无约束求解函数minbnd和fminunc。 3、无约束优化问题求解 （1）一元函数无约束优化问题 fminbnd常用格式如下：（1）x= fminbnd (fun, x, x2) （2）x= fminbnd (fun, x, x2 , options) （3）[x , fval]= fminbnd（...）（4）[x , fval , exitflag]= fminbnd（...）（5）[x , fval , exitflag , output]= fminbnd（...） fminbnd的算法基于黄金分割法和二次插值法，它要求目标函数必须是连续函数，并可能只给出局部最优解 （2）多元函数无约束优化问题 这里X为n维变量 Fminunc常用格式为: （1）x= fminunc（fun, X0）； （2）x= fminunc（fun, X0，options）； （3）[x，fval]= fminunc（...）； （4）[x，fval，exitflag]= fminunc（...）； （5）[x，fval，exitflag，output]= fminunc（...） X0为初始值 实验步骤 完成下面计算内容 练习1：计算函数f(x)在x=5处的值, f(x)如下: 练习2：设三个复数 a＝3＋4i，b＝1＋2i， ，计算 x=ab/c （自行查找matlab复数表达方法） 练习3：计算 的结果，其中x= 45° 练习4：随机生成一个元素在1～100间的6×5矩阵（使用随机函数rand，用help命令查阅使用方法） z(x1, x2) 表示总利润； p，q，x 分别表示甲的价格、成本、销量； p，q，x 分别表示乙的价格、成本、销量； a，b，λ, ci（i，j =1，2）是待定系数. 甲的价格 p1 = b1 - a11 x1 - a12 x2 ，b1，a11，a12 0，且 a11 a12 乙的价格p2 = b2 - a21x1- a22x2 ，b2，a21，a22 0， 且 a22 a21 甲的成本 乙的成本 若根据大量的统计数据,求出系数 b1=100,a11=1,a12=0.1, b2=280,a21=0.2,a22=2, r1=30,λ1=0.015,c1=20, r2=100,λ2=0.02,c2=30, 问题转化为无约束优化问题：求甲,乙两个牌号的产量x1，x2，使总利润最大。 1 命令历史窗口 命令窗口