基于MatLa分布式编程模型的机械优化设计算法实验指导书.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 tol eps k 100 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 各个函数的语法结构和适用解决优化问题的类型，请参阅MATLAB帮助文件。另外，在MATLAB的命令窗口键入命令 help optim 结果显示该工具箱中所有函数清单,部分函数如下图示。 2.2函数调用时需要注意的问题 使用优化工具箱时，由于优化函数要求目标函数和约束条件满足一定的格式，所以需要用户在进行模型输入时注意以下几个问题：1 .目标函数最小化优化函数fminbnd、fm