文档介绍：机械优化实例及matlab工具箱
POWERPOINT
1
（2）抗皱稳定性
扭转切应力不得超过扭转稳定得临界切应力:
空心传动轴的扭转稳定的临界切应力为:
设计实例1:
整理得:
设计实例1ilename is a
x = 0::2*pi;
y = sin(x);
plot(x,y)
M-文件的操作
M-文件的操作
M-文件的操作
M-文件的操作
第2部分 优化计算工具
线性规划优化函数
无约束非线性优化函数
约束优化函数
MATLAB解决的线性规划问题的标准形式为：
线性规划优化函数
[函数] linprog
[格式] x = linprog(f,A,b,Aeq,beq)
x = linprog(f,A,b,Aeq,beq,lb,ub,x0)
f——是优化参数x的系数矩阵；
A——线性不等式约束系数阵
b——线性不等式约束常数向量
Aeq——线性等式约束系数阵
Beq——线性等式约束常数向量
线性规划及其优化函数
[应用举例]
求使函数 取最小值的x值，
且满足约束条件：
[代码] f = [-5; -4; -6];
A = [1 -1 1;3 2 4;3 2 0];
b = [20; 42; 30];
lb = zeros(3,1);
[x,fval] = linprog(f,A,b,[],[],lb)
[结果] x =
fval = -
无约束非线性优化函数
[函数]fminunc
[格式]
x = fminunc(fun,x0)
x = fminunc(fun,x0,options)
[x,fval] = fminunc(…)
无约束非线性优化函数
[应用举例] 求 的最小值
[代码]%
function f=myfun(x)
f=3*x(1)^2+2*x(1)*x(2)+x(2)^2
%然后调用函数 fminunc
x0=[1,1];
[x,fval]=fminunc(myfun,x0)
无约束非线性优化函数
[结果]
x =
-008 *
-
fval =
-016
约束优化函数
[函数] fmincon
[格式]
x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)
x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options)
[x,fval] = fmincon(…)
约束优化函数—盖板优化实例：
目标函数：
约束：
盖板优化实例
盖板优化实例
盖板优化实例
运行结果：
x =
fval =
盖板优化实例
前面空心轴的问题：
clear all
x0=[23,19,4];
options=optimset('largescale','off','display','iter','tolx',1e-6);
[x,fval,exitflag,output]=fmincon('myfun1',x0,[],[],[],[],[],[],'confun1',options)
function f=myfun1(x)
f=*(x(1)^2-x(2)^2)*x(3)*10e-6