文档介绍：实验设计
连杆设计优化问题的解决：
问题：
设计一曲柄摇杆机构要求曲柄L1从A0转到Ai=A0+90时，摇杆L3的转角 再次恢复到初始值B0,已知道Bi的运动规律为：Bi = B0+2(Ai-A0)2/(3〃),且 已知Ll=l,机架
for a=a0:pi/18:(a0+)以pi/18为步长测试a不同值的情况
b=bO+2*(a-aO)A2/(3*pi)测试b不同值的情况
r=sqrt(26-10*cos(a)) r 为 A0 所对的边
m=acos((rA2+x(2)A2-x( 1)人2)/(2*r*x(2)))
n=acos((r A2+24)/(10*r))
bb=pi-m-n
i=(bb-b)A2
f=f+i
end循环结束
： function [c,ceq]=mycon(x)
c=[x( l)A2+x(2)A2-*x( 1 )*x(2)-16;
36-x(l)A2-x(2)A2-*x( 1 )*x(2)];
ceq=[]
最后在MATLAB命令窗口调用主程序:
format long ;
x0=[6,3]；
A=[-l-l;l-l;-l 1];
b=[-6;4;4];
lb=zeros(2,l);
零阵
以长数据格式显示
给xO赋初值
给A赋值
创建lb变量为两行一列的
options=optimset('largescale',,off);
[x,fval]=fmincon( ***@myfun,xO, A,b, [], [] ,1b, [], ***@mycon,options)调用优化函数
显示结果如下所示:
f =
**********
ceq =
[]
Optimization terminated: Magnitude of directional derivative in search direction less than 2* and maximum constraint violation is less than .
Active inequalities (to within = le-006): lower upper ineqlin ineqnonlin
fval =
**********
如下是相关截屏
约束函数的M文件
|function [c,ceq]=aycon(x)
c=[x(l)"2+x(2)"2-«x(l)*x(2)-16,
36-x(l)*2-x(2) *2-« (l)«x(2)] ceq=[]
我的电脑|
iAILAB
网上郃居
快捷方式
Editor - F:\,
■ Etftor - F:\,matiab\w...
Current Directory F
L... • X
囱廿殖扇回日
All Fites
I Fie Type
障
M-file