_SHI_05_数学实验-代数方程与最优化问题_V3.1.pptx

_SHI_05_数学实验-代数方程与最优化问题_V3.1.pptx

  1. 1、本文档共90页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数 学 实 验 第五讲:代数方程与最优化问题;主要内容; 基于MATLAB 的代数方程(组)求解;1.代数方程的求解;1.1 代数方程的图解法 1.1.1 一元方程的图解法;1.1.2 二元方程的图解法;1.2 多项式型方程的准解析解法;【例4】;【例5】;【例6】;【例7】;1.3 一般非线性方程数值解;【例8】;【例9】;EqSol.m;【例10】; 基于MATLAB优化工具箱 的优化计算;优化工具箱功能概述;一、MATLAB优化工具箱中常用的优化功能函数 求解线性规划问题的主要函数是linprog。 求解二次规划问题的主要函数是quadprog。 求解无约束非线性规划问题的主要函数是fminbnd、fminunc和fminsearch。 求解约束非线性规划问题的主要函数是fgoalattain和fminimax。;一、MATLAB优化工具箱中常用的优化功能函数; 使用优化函数或优化工具箱中其它优化函数时, 输入变量见下表:;输入参数中可以用options,用于所有函数,其中??括有以下参数。 (1)Display:结果显示方式,off不显示,iter显示每次迭代的信息,final为最终结果,notify只有当求解不收敛的时候才显示结果。 (2)MaxFunEvals:允许函数计算的最大次数,取值为正整数。 (3)MaxIter:允许迭代的最大次数,正整数。 (4)TolFun:函数值(计算结果)精度,正整数。 (5)TolX:自变量的精度,正整数。 而且可以用函数optimset创建和修改。;优化函数的输出变量:;二、一般步骤;2.无约束最优化问题求解;2.1 解析解法和图解法;【例11】;2.2 基于 MATLAB 的数值解法;函数 fminbnd;函数fminbnd;函数fminsearch;函数fminsearch;【例12】;EqSol.m;2.3 全局最优解与局部最优解;EqSol.m; 1.使用格式: [x,fval,exitflag,output,grad,hessian]=fminunc(@fun,x0,options,P1,P2…);[3] fminunc为中型优化算法的步长一维搜索提供了两种算法, 由options中参数LineSearchType控制: LineSearchType=’quadcubic’(缺省值),混合的二次和三次多项式插值; LineSearchType=’cubicpoly’,三次多项式插值;管道截面积:;;;【例14】;EqSol.m;有约束最优化问题 的计算机求解;3.1 约束条件与可行解区域;【例15】;3.2 线性规划问题的数学模型及计算机求解;【例16】;【例17】;线性规划应用示例;3.确定约束条件:;3.3 二次型规划的求解;二次规划问题的数学模型;【例18】;;3.4 一般非线性规划问题的求解;;【例19】;EqSol.m;【例20】;EqSol.m;fgoalattain;多目标优化问题的数学模型;;; 设计带式输送机传动装置上的普通V带传动。已知电机额定功率P=4kW,转速n1=1440r/min,传动比i=3,采用A型V带,每天工作不超过10小时。要求传动结构紧凑(带的根数尽量少,带轮直径和中心距尽量小)。 解:(1)建立优化设计的数学模型 ①设计变量: V带传动的独立设计变量是小带轮直径dd1和带的基准长度Ld X=[dd1,Ld]T=[x1,x2]T ②目标函数: 小带轮直径: 中心距: 带的根数:; ③约束条件:含性能约束和边界约束;函数fgoalattain应用示例(续);(3)编制优化设计的M文件;(3)编制优化设计的M文件 %V带传动多目标优化设计的约束函数文件 function[g,ceq]=VDCD_3mb_YS(x) i=3;n1=1440; %已知条件:传动比,转速 g(1)=100-x(1); %小带轮直径>=Ddmin g(2)=pi*x(1)*n1/6e4-25 %带速范围V<=Vmax a1=x(2)/4-pi*x(1)*(i+1)/8; a2=x(1)^2*(i-1)^2/8; a=a1+sqrt(a1^2-a2); g(3)=120-180*(1-

文档评论(0)

***** + 关注
实名认证
内容提供者

我是自由职业者,从事文档的创作工作。

1亿VIP精品文档

相关文档