应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 机械 » 机械设计 » MATLAB优化问题
51/1返回列表
查看: 915  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

老张666

新虫 (初入文坛)

[求助] MATLAB优化问题已有2人参与

我是初学者，请教一个问题，在做无约束优化时，zeros（1,16）这个1和16 是怎么来的啊？
我查了很多资料，都没有解释的
另外这样写目标函数，约束函数和主函数对不对？
目标函数：
function f=mubiaohanshou(x)      
OM=x(1);
ON=x(2);
F1=x(3);
OO1=0.55;
O1O2=0.55;
m1=23.67;
m2=11.7;
m3=7.9;
m4=57.9;
OG1=0.415;
OG2=0.26;
for i=4:10;%大臂转角在30°-90°范围时
    alfa(i)=(i-1)*10*pi/180;
    j=1:13;%机械手小臂转角在0°-120°
    beta(j)=(j-1)*10*pi/180;
    MN(i)=sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(0.5*pi-alfa(i)));
    Ld(i)=OG1*cos(alfa(i));
    Lx(j)=OO1*cos(alfa(i)+O1G2*sin(beta(j)));
    Lf(j)=OO1*cos(alfa(i)+O1O2*sin(beta(j)));
    Lp(i)=(x(1)*x(2)*sin(0.5*pi-alfa(i)))/MN(i);
    Md(i)=m1*Ld(i);
    Mx(j)=m2*Lx(j);
    Mk(j)=m3*Lf(j);
    Mf(j)=m4*Lf(j);
    Mp(i)=2*x(3)*Lp(i);
    ditam1(i)=Md(i)+Mx(j)+Mk(j)-Mp(i);
    ditam2(i)=Md(i)+Mx(j)+Mk(j);
    ditam3(i)=Md(i)+Mx(j)+M(j)-Mp(i);
    ditam4(i)=Md(i)+Mx(j)+M(j);
    end
for i=11:16;%大臂转角在90°-150°范围时
    alfa(i)=(i-1)*10*pi/180;
    j=1:16;%机械手小臂转角在0°-150°
    beta(j)=(j-1)*10*pi/180;
    MN(i)=sqrt(x(1)^2+x(2)^2-2*x(1)*x(2)*cos(alfa(i)-0.5*pi));
    Ld(i)=OG1*cos(pi-alfa(i));
    Lx(j)=abs(OG1*cos(pi-alfa(i))-O1G2*sin(beta(j)));
    Lf(j)=abs(OO1*cos(pi-alfa(i)-O1O2*sin(beta(j))));
    Lp(i)=(x(1)*x(2)*sin(alfa(i)-0.5*pi))/MN(i);
    Md(i)=m1*Ld(i);
    Mx(j)=m2*Lx(j);
    Mk(j)=m3*Lf(j);
    Mf(j)=m4*Lf(j);
    Mp(i)=2*x(3)*Lp(i);
   ditam1(i)=Md(i)+Mx(j)+Mk(j)-Mp(i);
   ditam2(i)=Md(i)+Mx(j)+Mk(j);
   ditam3(i)=Md(i)+Mx(j)+M(j)-Mp(i);
   ditam4(i)=Md(i)+Mx(j)+M(j);
f1=max(ditam1(i));
f2=max(ditam2(i));
f3=max(ditam3(i));
f4=max(ditam4(i));
end
f=max(f1,f3);
约束函数：
function[g,ceq]=yueshuhanshu(x)   
g(1)=0-x(1);
g(2)=x(1)-1000;
g(3)=10-x(2);
g(4)=x(2)-40;
g(5)=10-x(3);
g(6)=x(3)-500;
ceq=[];
主函数：
x0=[50;20;30];
A=[-1,0,0;1,0,0;0,-1,0;0,1,0;0,0,-1;0,0,1];
b=[0,1000,-10,40,-10,500];
Aeq=[];
beq=[];
lb=[0;10;10];
ub=[1000;40;500];
[x,fn]=fmincon(@zhulijixieshoumubiaohanshu,x0,A,b,Aeq,beq,lb,ub,@zhulijixieshouyueshuhanshu);
disp('最优解')
fprintf('  x(1)=%3.4f  \n',x(1))
fprintf('  x(2)=%3.4f  \n',x(2))
fprintf('  x(3)=%3.4f  \n',x(3))
fprintf('  f=%6.4f  \n',f)
请指教，谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2016-12-09 20:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuburfransth

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1

chendequan

铁虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1

老张666

新虫 (初入文坛)

老张666

新虫 (初入文坛)
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭