pusis=pi-(jBDC-jBDA); end
fx=fx+(pusi-pusis)^2; end f=fx;
图4
⑵编写约束函数M文件并以文件名mycon保存在MATLAB目录下的R2011a文件夹中，如图5所示：
约束函数的MATLAB程序如下：
function [g,cep]=mycon(x) a=1; d=5;
g=[x(1)^2+x(2)^2-(d-a)^2-2*x(1)*x(2)*cos(40/180*pi) -(x(1)^2+x(2)^2-(d-a)^2)+2*x(1)*x(2)*cos(50/180*pi)
(x(1)^2+x(2)^2-(d+a)^2)-2*x(1)*x(2)*cos(130/180*pi) -(x(1)^2+x(2)^2-(d+a)^2)+2*x(1)*x(2)*cos(140/180*pi)]; cep=[];
图5
⑶在命令窗口调用优化程序，如图6所示： 调用函数的MATLAB程序如下：
function main clc; clear all; x0=[4.5;4]; a=1;
d=5;
A=[-1,-1;1,-1;-1,1]; b=[-a-d;d-a;d-a]; lb=[0;0]; ub=[];
[x,fval]=fmincon(@fun,x0,A,b,[],[],lb,ub,@mycon)
图6
⑷运行求解结果并对结果分析 运行结果如下：
Optimization terminated:first-order optimality measure less than options.TolFun and maximum constraint violation is less than options.TolCon.
Active inequalities (to within options.TolCon =1e-006): low upper ineqlin ineqnonlin 2
x=
2.5698 5.4218 fval= 0.0688
?x1??2.5698??l2? 当x????????l?时，即连杆的长度的最优解为 x5.4218??3??2??l2?2.5698 l3?5.4218
此时目标函数返回的最优值fval=0.0688。
3、结论分析
基于MATLAB对四杆机构进行优化设计，达到了预期的目标。MATLAB优化可以很好的实现机械优化设计的算法过程，且求解方便快捷，程序简明，效率高。比起图解法和实验法，有着较大的优越性，是解决机械优化设计问题的重要工具。
4、参考文献
[1]张永恒.工程优化设计与MATLAB实现[M].北京：清华大学出版社，2011.9 [2]李明.详解MATLAB在最优化计算中的应用[M].北京：电子工业出版社，2011.5 [3]孙恒.机械原理[M].北京：高等教育出版社，2006.5
[4]李建霞，王良才.基于MATLAB的四杆机构优化设计简介[J].几机电产品开发与创新，2010，（01）：088-02