[原创]哪位高手帮一下！关于优化工具中的fmincon函数

haode123 · 发表于 2006-5-20 16:11

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我在解一个带约束的非线性最小化问题时，总是出现很多行如下警告， Warning: One or more output arguments not assigned during call to 'feval'. > In D:\MATLAB6p1\toolbox\optim\private\nlconst.m at line 269 In D:\MATLAB6p1\toolbox\optim\fmincon.m at line 458 In D:\MATLAB6p1\work\double_canshu_zheng.m at line 16 Warning: Reference to uninitialized variable nceqtmp in D:\MATLAB6p1\toolbox\optim\private\nlconst at line 271. > In D:\MATLAB6p1\toolbox\optim\private\nlconst.m at line 271 In D:\MATLAB6p1\toolbox\optim\fmincon.m at line 458 In D:\MATLAB6p1\work\double_canshu_zheng.m at line 16 有人建议我用大型课题优化，如果用大型课题要编写目标函数的梯度函数和hess矩阵。 我看了很多书都没提具体怎么求梯度函数和hess矩阵，那位高手帮忙解决一下。 越详细越好，先谢了！！！！

happy · 发表于 2006-5-20 16:33

把代码贴一下

haode123 · 发表于 2006-5-20 16:49

原函数程序如下： function [f,g]=fun_zheng(x) a=x(1); b=x(2); c=x(3); d=x(4); e=x(5); [qz,qf]=weibull_plot_double('E:\data\pd\hunhe\fengzhi.txt',0); qz=qz(find(qz>0)); f=-sum(log(x(5)*(x(1)/x(2))*(qz./x(2)).^(x(1)-1).*exp(-(qz./x(2)).^x(1))))-sum(log((1-x(5))*(x(3)/x(4))*(qz./x(4)).^(x(3)-1).*exp(-(qz./x(4)).^x(3)))); 主程序如下： function [b1,a1,b2,a2,p]=double_canshu_zheng(file,b10,b20,p0) [qz,fz,qf,ff]=weibull_plot_double(file,0); a10=qz(find(fz>=0.632*p0)); a20=qz(find(fz>=0.632*(p0-1)+p0)); a10=a10(1); a20=a20(1); x0=[b10,a10,b20,a20,p0]; A=[]; b=[]; Aeq=[]; beq=[]; lb=[0.01,0.01,0.01,0.01,0.01]; ub=[inf,inf,inf,inf,0.99]; %options=optimset('GradObj','on'); %options=optimset('largescale','off','display','iter'); [x,fval]=fmincon(@fun_zheng,x0,A,b,Aeq,beq,lb,ub,@mycon_zheng,options); b1=x(1) a1=x(2) b2=x(3) a2=x(4) p=x(5)

happy · 发表于 2006-5-20 16:51

E:\data\pd\hunhe\fengzhi.txt file,b10,b20,p0 没有给出来

haode123 · 发表于 2006-5-20 16:52

麻烦happy教授给解决一下，不胜感激！

happy · 发表于 2006-5-20 16:53

先把参数给一下

haode123 · 发表于 2006-5-20 16:55

function [f]=fun_zheng(x) [qz,qf]=weibull_plot_double('E:\data\pd\hunhe\fengzhi.txt',0); qz=qz(find(qz>0)); f=-sum(log(x(5)*(x(1)/x(2))*(qz./x(2)).^(x(1)-1).*exp(-(qz./x(2)).^x(1))))-sum(log((1-x(5))*(x(3)/x(4))*(qz./x(4)).^(x(3)-1).*exp(-(qz./x(4)).^x(3))));

haode123 · 发表于 2006-5-20 17:01

目标函数中的qz参数是一列试验数据.

haode123 · 发表于 2006-5-20 17:07

E:\data\pd\hunhe\fengzhi.txt 见文件 file为文件路径file=‘E:\data\pd\hunhe\fengzhi.txt’ b10，b20要求正数就可以 p0为在0到1范围内的值，在上述中可取0.1

[此贴子已经被作者于2006-5-20 17:09:44编辑过]

YooYo · 发表于 2006-5-20 19:19

你FMINCON函数不对,非线形约束的M文件写错了.

haode123 · 发表于 2006-5-20 20:06

哪里错了请指点

haode123 · 发表于 2006-5-20 20:46

哪里错了能指点一下吗？

lovewft · 发表于 2006-5-25 19:23

fmincon函数，对目标函数与约束条件有什么具体要求？？？比如说连续，可微等

账号		自动登录	找回密码
密码			我要加入

[原创]哪位高手帮一下！关于优化工具中的fmincon函数

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（haode123）[原创]哪位高手帮一下！关于优化工...

回复：（happy）回复：（haode123）[原创]哪位高手帮一...

回复：（haode123）[原创]哪位高手帮一下！关于优化工...

回复：（haode123）[原创]哪位高手帮一下！关于优化工...

传错了，原函数应该见如下：

回复：（haode123）传错了，原函数应该见如下：

回复：（happy）回复：（haode123）[原创]哪位高手帮一...

回复：（YooYo）你FMINCON函数不对,非线形约束的M文件...