| 查看: 3611 | 回复: 2 | ||
[求助]
MATLAB用fmincon函数求最小优化问题，提示输入参数不足，求解答
|
|
如下是自己编写的一段程序，目的是想通过fmincon函数求优化求参数k；但是运行就出现问题，Sample Text function f=fop(k) S=[39.2486,18.6974,5.5453,12.0719,6.8712;21.3467,30.485,4.1298,15.0155,7.1906;53.3333,31.0388,6.4493,16.407,7.4912;55.3806,24.7649,3.9371,25.6651,7.0956;45.4191,16.4854,6.0034,23.342,7.5939;42.5722,20.3236,3.7494,21.9155,6.8254;31.3707,21.8137,4.3522,13.5309,6.9808;49.7532,28.9058,5.7443,24.999,6.6091;32.6931,14.9153,6.8971,17.3615,7.0777;26.6678,15.8494,3.4524,13.1626,7.4408;48.5182,18.8443,7.3267,28.9634,7.3066;36.088,26.5438,6.7461,26.9071,7.7217;28.8349,27.7746,5.2497,18.8,7.343;38.182,25.7221,4.8418,21.2241,6.72;23.5811,23.0133,4.9581,19.5143,7.6823]; %S样本数据（15*5） Y=[0.132;0.155;0.636;0.963;0.648;0.534;0.123;0.854;0.225;0.015;0.948;0.734;0.345;0.504;0.268];%Y响应（15*1） %k=[k(1),k(2),k(3),k(4),k(5)];%给定的系数，目前还未解决如何确定这个系数 F=ones(15,1);%采用x的0阶多项式 result = zeros(15*15,5); for i = 1:15 for j = 1:15 result((i-1)*15+j,: ) = abs(S(i,: )-S(j,: ));%所有的S行元素两两相减，得到225*5 end end result_=(result').^2;%平方,得到5*225矩阵 R=exp(reshape(-[k(1),k(2),k(3),k(4),k(5)]*result_,15,15));%得到R(k;S)15*15,对角元素为1的对称矩阵 jeta=inv(F'*inv(R)*F)*F'*inv(R)*Y;%一个数 thetaz2=1/15*(Y-F*jeta)'*inv(R)*(Y-F*jeta);%一个数 f=(abs(det(R)))^(1/15)*thetaz2;%函数表达式 option=optimset;option.LargeScale='off';option.Display='off'; k0=[1;1;1;1;1];lb=[0;0;0;0;0];%初值，下界 [k,f]=fmincon('fop',k0,[],[],[],[],lb,[],[],option) |
回复此楼» 猜你喜欢
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
申请科研项目纯靠实力能成功吗？
已经有7人回复
-
极限拉扯，如何抉择
已经有14人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
2楼2018-01-13 16:35:11
chenglei2
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 190
- 散金: 366
- 帖子: 80
- 在线: 84小时
- 虫号: 1363197
- 注册: 2011-08-08
- 性别: GG
- 专业: 金属材料的合金相、相变及
3楼2019-09-16 21:02:41
10