基于MatLab分布式编程模型的机械优化设计算法实验指导书.docVIP

基于MATLAB分布式编程模型的机械优化设计算法实验指导书 一、实验所需基础知识 MATLAB是matrix laboratory两个词的组合意为矩阵工厂矩阵实验室。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。在很大程度上摆脱了传统非交互式程序设计语言如C、的编辑模式function [x,minf] = minHJ(f,a,b,isDisplay,eps) % 一维黄金分割法 %f: 目标函数的符号表达式 %a: 区间左端 %b: 区间右端 %isDisplay: 是否显示求解结果 %eps: 收敛精度 if nargin = 3 isDisplay = false; end format long; if nargin = 4 eps = 1.0e-6; end l = a + 0.382*(b-a); u = a + 0.618*(b-a); k=1; tol = b-a; while toleps k100 fl = subs(f , findsym(f), l); fu = subs(f , findsym(f), u); if fl fu a = l; l = u; u = a + 0.618*(b - a); else b = u; u = l; l = a + 0.382*(b-a); end k = k+1; tol = abs(b - a); end if k == 100000 disp(找不到最小值！); x = NaN; minf = NaN; return; end x = (a+b)/2; minf = subs(f, findsym(f),x); if isDisplay disp ----------黄金分割法求解结果----------; fprintf(迭代次数 k=%6d\n,k); fprintf(最优点 x*=%6.2f\n,x); fprintf(最优值 f(x*)=%6.2f\n,minf); end format short; 1.2 命令文件 命令文件没有输入参数也没有输出参数，只是一些MATLAB命令和函数的组合。命令文件可以操作工作空间的变量，也可以生成新的变量，但是当命令文件执行结束后新变量将保存在工作空间中，不会被自动清除。 例1.2：调用黄金分割法求解一维函数的最优解minHJ_test.m clc; syms t; f = t^3 - 3*t + 2; [x,minf]=minHJ(f,-3,5,true); 在编辑器界面按F5键运行该程序，其运行结果如下： ----------黄金分割法求解结果----------- 迭代次数 k= 34 最优点 x*= 1.00 最优值 f(x*)= 0.00 1.3 子函数 一个M文件中可以定义多个函数，其中第一个定义的函数称为主函数，主函数名必须与M文件名相同，其它函数称为子函数。子函数只能被定义它的文件中的函数调用，其它函数和文件不能调用。 1.4 局部函数 把一些文件放在private目录下，那么这些文件中定义的函数就是局部函数。局部函数只能被其父母目录中的函数或文件调用，其它函数和文件无权调用。 1.5 M文件的执行 当Matlab 遇到一个新的名称（命令、函数、变量）时，按如下次序查找这个名称： （1）当前的工作空间； （2）子函数名称； （3）局部函数名称； （4）set path 指定的路径。 要想使自己编写的函数或命令能够在命令窗口自由的使用，一个比较好的处理方法是建立一个目录，并把这个目录加到Set Path指定的目录中，将自己编写的M 文件都存放到这个目录中。 2．MATLAB优化工具箱(Optimization Toolbox) 我们除了可以运用MATLAB语言自己编写优化算法外，还可以直接调用优化工具箱来解决工程设计中的优化问题。 2.1 MATLAB优化工具箱的常用函数 MATLAB优化工具箱由一系列函数组成，主要包括： 边界约束条件下的非线性最小化 Fminbnd 求解多变量函数的最小化 Fminunc 求解无约束非线性最小化 Fminsearch 求解线性规划问题 Linprog 求解二次规划问题 Quadprog 求解有约束的非线性最小化 Fmincon 求解多目标规划的优化问题 Fgoalattain 求解最小、最大化问题 Fminimax 各个函数的语法结构和