新手求助装载机举升机构优化问题，大虾前辈们帮我！！

zjg2000

那位前辈能帮我编一下这个程序：关于装载机举升机构优化问题
目标函数：举升机构的传力比要求最大：UF2=R9*R10*sin((UG3+UG-UG5)/180*pi)/(HI4*(H4+XA));
其中R9、R10、（杆长）UG3、UG5（角度）是需要优化的变量。
其中初始给定值：R9=165；R10=70;UG3=0;UG5=12;<BR>UG=UG4：0.5：UG6；
给定常数及关系式：H1=115;
H2=300;
H3=133;
Y1=210;
Y3=27;
Xj=80;
R11=80.75;
C1=sqrt(Y3^2+H3^2);
D1=atan(Y3/H3);
L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);
H4=H3-10;<BR>UG4=(acos((Y1-Y3)/L))*180/pi;
A1=H2-Y1+C1*sin(D1+pi/4);
A2=Xj+R11+H1-C1*cos(D1+pi/4);
A3=atan(A1/A2);
UG6=(A3+pi/2)*180/pi;
HI4=sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5));
XA=L*sin (UG/180*pi);
约束条件：
50-R9*R10*sin(UG3+UG6-UG5)/sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG6-UG5))〈0
900-sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5))〈0
400-(sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG6-UG5))-sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5)))〈0
可是用fmincon等函数编写，我写的程序老是有错误，实在没办法，着急啊，请各位大虾前辈帮帮我。

[ 本帖最后由 xinyuxf 于 2007-3-14 11:33 编辑 ]

frogfish

<P>什么错误贴一下，我的现在的机器里没装matlab </P>

zjg2000

我慢慢的调试好了，是关于我毕业设计的，不过老师要求我们用基于遗传算法的工具箱进行优化，可是我不知道这个工具箱 ，你老人家能不能指教一下

zjg2000

<P>我的目标函数是：<BR>function f=Mubiaohanshu(x,UG,H4,XA)<BR>H1=115;<BR>H2=300;<BR>H3=133;<BR>Y1=210;<BR>Y3=27;<BR>Xj=80;<BR>R11=80.75;<BR>C1=sqrt(Y3^2+H3^2);<BR>D1=atan(Y3/H3);<BR>L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);<BR>UG4=(acos((Y1-Y3)/L))*180/pi; <BR>A1=H2-Y1+C1*sin(D1+pi/4);<BR>A2=Xj+R11+H1-C1*cos(D1+pi/4);<BR>A3=atan(A1/A2);<BR>UG6=(A3+pi/2)*180/pi;<BR>UG=(UG4:0.5:UG6);<BR>H4=H3-10;<BR>XA=285.*sin (UG./180.*pi);<BR>f=x(1)*x(2)*sin((x(3)+UG-x(4))/180*pi)/((sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4)))*(H4+XA)))<BR>end<BR>约束条件是：function [C,Ceq]=Yueshutiaojian(x)<BR>H1=115;<BR>H2=300;<BR>H3=133;<BR>Y1=210;<BR>Y3=27;<BR>Xj=80;<BR>R11=80.75;<BR>C1=sqrt(Y3^2+H3^2);<BR>D1=atan(Y3/H3);<BR>L=sqrt((Xj+R11+H1-C1*cos(D1+pi/4))^2+(H2-Y1+C1*sin(D1+pi/4))^2);<BR>UG4=(acos((Y1-Y3)/L))*180/pi; <BR>A1=H2-Y1+C1*sin(D1+pi/4);<BR>A2=Xj+R11+H1-C1*cos(D1+pi/4);<BR>A3=atan(A1/A2);<BR>UG6=(A3+pi/2)*180/pi;<BR>C(1)=50-x(1)*x(2)*sin(x(3)+UG6-x(4))/sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG6-x(4)));<BR>C(2)=900-sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4)));<BR>C(3)=1400-(sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG6-x(4)))-sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(x(3)+UG4-x(4))));<BR>Ceq=[];<BR>主函数是：x0=[160 70 0 12];<BR>LBnd=[130,50,-5,7];<BR>UBnd=[195,95,5,17];<BR>[x,f]=fmincon('Mubiaohanshu',x0,[],[],[],[],LBnd,UBnd,'Yueshutiaojian')<BR>R9=x(1),R10=x(2),UG3=x(3),UG5=x(4)<BR>但是结果R9 =200.2440、R10 =100.2440、UG3 = -9.7700UG5 =1.7560，超出了我的限制范围为什么？</P>

yanpi88

你好 我也是研究类似课题的 请问：约束条件你是怎么建立的啊？我参考了王国彪的《装载机工作装置优化设计》里的，建立了一些约束条件，采用了VB编程，但是在满足约束条件的情况下，得出的优化结构数据是不对的（要么越界，要么为0）。因此，我觉得可能是我的约束条件列错了。但是有许多是无法知道的。例如：在卸料性中的CC（最小安装距），DD（最大行程）等。如果知道，麻烦赐教。我的邮箱是:yanpi88@yahoo.com.cn，万分感谢。

[ 本帖最后由 mjhzhjg 于 2007-3-18 18:58 编辑 ]

jli

可以在isight中间方便的集成，你说的详细一点，我可以做个例子给你。

mjhzhjg

楼主还不尽快提供啊，你遇到了热心人

jli

UG6=(A3+pi/2)*180/pi;
HI4=sqrt(R9^2+R10^2-2*R9*R10*cos(UG3+UG4-UG5));
XA=L*sin (UG/180*pi);

UG6和XA的值小于零，这个物理意义对吗？