[理化生]第3章优化设计.ppt

MATLAB基础知识 一、变量 MATLAB不要求用户在输入变量的时候进行声明。当用户在MATLAB工作空间内输入一个新的变量时，MATLAB会自动给该变量分配适当的内存，若用户输入的变量已经存在，则MATLAB将使用新输入的变量替换原有的变量。 无论给MATLAB中的变量赋什么值，在内存中均以矩阵的形式表示。 在MATLAB中，变量名是以字母开头，后接字母、数字和下划线的字符序列，最多31个字符。另外，在MATLAB中，变量名区分字母大小写。 二、矩阵变量的建立 (1) 把矩阵元素列入方括号中； (2) 每行内的元素间用逗号或空格分隔； (3) 行与行之间用分号分隔； 三、MATLAB运算 2.关系运算 3.逻辑运算 (2)从向量、矩阵或数组中挑选出指定的行列元素。 二维图形的输出 使用函数fmincon 数学模型： X=[x1 x2]T Min F(x)=120x1+x2 Sub.to g1(X)=x1≥0 g2(X)=x2 ≥0 g3(X)=1/4x2-1 ≥0 g4(X)=7/45x1x2-1 ≥0 g5(X)==7/45x13x2-1 ≥0 g6(X)=1/321x1x22-1 ≥0 -x2/4=-1 -7x1x2/45+1=0 -7x13x2/45+1 =0 -x1x22/321+1 =0 f=inline(120*x(1)+x(2)) x0=[10;10] a=[0 -0.25] b=-1 aeq=[] beq=[] lb=[0;0] ub=[] [x,fval,exit]=fmincon(f,x0,a,b,aeq,beq,lb,ub,@mycon) function [c,ceq]=mycon(x) c=[-7*x(1)*x(2)/45+1;-7*x(1)^3*x(2)/45+1;-x(1)*x(2)^2/321+1] ceq=[] mycon函数内容： -x2/4 * * 第2章 [例1] 确定职工编制问题 某厂每日八小时的产量不低于1800件。为了进行质量控制，计划聘请两种不同水平的检验员。一级检验员的标准为：速度 25件/小时，正确率 98%，计时工资 4元/小时；二级检验员的标准为：速度 15件/小时，正确率 95%，计时工资 3元/小时。检验员每错检一次，工厂要损失2元。现有可供厂方聘请的检验员人数为一级8人和二级10人。为使总检验费用最省，该工厂应聘一级、二级检验员各多少名？ [例2] 工件加工任务分配问题 某车间有两台机床甲和乙，可用于加工三种工件。假定这两台机床的可用台时数分别为700和800，三种工件的数量分别为300、500和400，且已知用三种不同机床加工单位数量的不同工件所需的台时数和加工费用（如表2所示），问怎样分配机床的加工任务，才能既满足加工工件的要求，又使总加工费用最低？ 表2 机床加工情况表 800 8 12 11 1.3 1.2 0.5 乙 700 10 9 13 1.0 1.1 0.4 甲 工件3 工件2 工件1 工件3 工件2 工件1 可用 台时数 单位工件的加工费用 单位工作所需加工台时数 机床类型 3.已知功率P=20kw,小齿轮转速n1=1000r/min; 齿数比u=3;载荷稳定,预期寿命10年,每年300个工作日,工作时间占20%。动力机为电动机,工作中有中等振动,传动不逆转,齿轮对称布置。人工选择齿轮材料和热处理:小齿轮40Cr调质,硬度260HB;大齿轮45调质,硬度240HB;许用应力为[σH1]= 798Mpa ，[σH2]=690Mpa; [σF1]=456Mpa , [σF2]=349Mpa。设计单级圆柱齿轮传动，使齿轮体积最小。 以体积最小为目标，设计变量为法向模数mn、小齿轮齿数z1及齿宽系数φd，则数学模型为： 约束条件为： 模数约束： 1.5≤x1≤20; 齿宽系数约束条件： 0.2≤x3≤1.2; 根切约束： g1(x)=17-x2≤0; 齿面接触应力约束：g2(x)=σH2-[σ]H2≤0; 齿根弯曲应力约束：g3(x)=σF1-[σ]F2≤0; g4(x)=σF2-[σ]F2≤0。 设计变量为： X = [x1,x2, x3]T= [mn, z1, φd]T 目标函数为： (1) 在命令窗口中直接输入 a=[1 2 3;4 5 6;7 8 9] 或 a=[1,2,3;4,5,6;7,8,9] (2) 在M文件中创建矩阵 M文件是一种包含MATLAB代码的文本文件。 M文件的创建 在MATLAB中，可以用以下方式输入