| 查看: 3581 | 回复: 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) |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
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
30