| 查看: 1197 | 回复: 5 | ||
[求助]
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) 请指教，谢谢 |
回复此楼» 猜你喜欢
有在开发铸钢丸的吗？
已经有0人回复
-
2026国际能源化学会议将于10月23-25日在南昌召开
已经有51人回复
-
工程热物理与能源利用论文润色/翻译怎么收费?
已经有76人回复
-
2026国际能源化学会议将于10月23–25日在南昌召开
已经有23人回复
-
锂电极片模切机是做什么的，如何使用和操作好设备？
已经有0人回复
-
英国拉夫堡大学-CSC（国家留学基金）奖学金招收2027/2028学年博士
已经有0人回复
-
想配vasp计算用主机，各位同行有推荐的配置吗? 求指点
已经有3人回复
-
不同封口压力对纽扣电池循环性能影响有多大？
已经有0人回复
-
铜合金材料加工课件资源
已经有0人回复
-
一条软包电池实验线，折边机是必须配的设备吗？
已经有0人回复
-
掺杂Er3+的体系 有没有更新的测温模式？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB优化问题
已经有12人回复
- MATLAB优化问题
已经有0人回复
- matlab 优化问题
已经有7人回复
- matlab优化报错
已经有8人回复
- MATLAB编程优化问题
已经有5人回复
- matlab 装配序列优化问题的matlab代码
已经有1人回复
- matlab优化求解中遇到问题，望大神指点。。。
已经有7人回复
- matlab最优化问题求解中，复杂的目标函数和约束条件无法直接表示的问题
已经有3人回复
- MATLAB 矩阵优化优化问题
已经有0人回复
- MATLAB 矩阵优化优化问题
已经有1人回复
- MATLAB最优化
已经有4人回复
- matlab优化
已经有5人回复
- matlab 培养基组分优化问题
已经有2人回复
- 帮忙给看看这个matlab优化函数 问题
已经有8人回复
- matlab优化问题
已经有2人回复
- 关于matlab优化的问题
已经有1人回复
- matlab线性矩阵不等式求解做优化问题
已经有10人回复
- matlab用fmincon优化问题
已经有2人回复
- 关于MATLAB数学建模最优化问题
已经有2人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2016-12-10 10:42:32
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
3楼2016-12-10 11:20:34
4楼2016-12-10 15:18:03
5楼2016-12-10 15:20:20
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 209小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
6楼2016-12-11 10:46:36
20