版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

muyouyou

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 234.8
散金: 5
帖子: 22
在线: 21.4小时
虫号: 3250631
注册: 2014-06-03
专业: 交通工程

[求助] matlab优化求解中遇到问题，望大神指点。。。已有1人参与

在matlab中
目标函数function=myfun（x）
f=0.5*1.935*10^(-5)*x(1)^2*x(2)*(x(3)+1.5)+0.5*5.75*x(2)^3*x(3)/x(1)^4;
约束条件function[C,Ceq]=mycon(x)
c=[1.35*10^(5)/x(1)+12.76*10^(6)*x(2)^2-6.9*10^5*x(1)*x(2)/(12.56*x(1)^3*x(2)-12.56*x(1)^4);x(1)^2*10^10/(x(2)*x(3));(x(1)*x(3)+x(1))/x(2)+86.25*10^3/x(2);x(1)*x(3)+1.5*x(1)-215];
ceq=[];
主程序：
x0=[38.6 191 4.06];
lb=[25 185 3.5];
ub=[40 250 7.0];
options=optimset('LargeScale','off');
[x,fval]=fmincon(‘myfun’,x0,[],[],[],[],lb,ub,'mycon')
运行结果为：??? Error using ==> fmincon
FMINCON cannot continue because user supplied nonlinear constraint function
failed with the following error:
One or more output arguments not assigned during call to 'F:\software\matlab\work\mycon.m (mycon)'.
本人初接触matlab，总是出现这种错误，不知是怎么回事。望高手指点。。。

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
孩子确诊有中度注意力缺陷已经有6人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求大神用matlab求下图阴影部分的面积和周长！！已经有4人回复
求matlab大神指点，关于logistic拟合已经有13人回复
现在用MATLAB求解一个循环后画出曲线遇到问题，求大神指点！！已经有7人回复
matlab最优化问题求解中，复杂的目标函数和约束条件无法直接表示的问题已经有3人回复
matlab求解，求大神已经有1人回复
matlab求解非线性方程组报错，请各位大神指点！已经有8人回复
matlab画图问题！！！急！希望会的大神来看看！！已经有11人回复
求大神帮忙matlab拟合函数求参数已经有9人回复
菜鸟求助，老师给了一个分子让优化，可总是错误，求大神指点已经有22人回复
求大神用matlab给我处理一组数据，有人没有啊？已经有15人回复
Matlab 哪位大神解决此问题重金悬赏已经有4人回复
求助：遗传算法的matlab编程求解资源优化问题已经有7人回复
MATLAB微分方程参数拟合问题，求大神已经有7人回复

1楼2014-11-05 09:18:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

2楼: Originally posted by 月只蓝 at 2014-11-05 21:21:38
约束条件是非线性的。
对于非线性约束情况，fmincon按附图1中格式调用。
附图1中正是一个例子。

重新上传附图1。

附图1.png

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2014-11-05 21:22:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

约束条件是非线性的。
对于非线性约束情况，fmincon按附图1中格式调用。
附图1中正是一个例子。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2014-11-05 21:21:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

muyouyou

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 234.8
散金: 5
帖子: 22
在线: 21.4小时
虫号: 3250631
注册: 2014-06-03
专业: 交通工程

引用回帖:

3楼: Originally posted by 月只蓝 at 2014-11-05 21:22:11
重新上传附图1。

附图1.png
...

function Fmincon
>>clear all;clc
>> x0=[38.6 191 4.06];
>> lb=[25 185 3.5];
>> ub=[40 250 7.0];
>> [x,fval]=fmincon(@func,x0,[],[],[],[],lb,ub,@NlinCons)
%---------------------------------------------------------
function[C,Ceq]=NlinCons(x)
c=[1.35*10^(5)/x(1)+12.76*10^(6)*x(2)^2-6.9*10^5*x(1)*x(2)/(12.56*x(1)^3*x(2)-12.56*x(1)^4);x(1)^2*10^10/(x(2)*x(3));(x(1)*x(3)+x(1))/x(2)+86.25*10^3/x(2);x(1)*x(3)+1.5*x(1)-215];
ceq=[];
%---------------------------------------------------------
functionf=func(x)
f=0.5*1.935*10^(-5)*x(1)^2*x(2)*(x(3)+1.5)+0.5*5.75*x(2)^3*x(3)/x(1)^4
我按你给的步骤计算了一下，结果是??? Error: File: F:\software\matlab\work\Fmincon.m Line: 3 Column: 1
Missing variable or function..是怎么回事呢？谢谢你了。。。

赞一下

回复此楼

4楼2014-11-06 09:32:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

引用回帖:

4楼: Originally posted by muyouyou at 2014-11-06 09:32:10
function Fmincon
>>clear all;clc
>> x0=;
>> lb=;
>> ub=;
>> =fmincon(@func,x0,[],[],[],[],lb,ub,@NlinCons)
%------------------------------------------------------ ...

CODE:

function test_1106

clear all;clc,format long

tic,

 x0=[38.6 191 4.06];

 lb=[25 185 3.5];

 ub=[40 250 7.0];

 OPTIONS=optimset('MaxFunEvals',5000,'TolFun',1e-6,'TolX',1e-6,'Algorithm','active-set');%,'Algorithm','interior-point','sqp','active-set');

 [x,fval]=fmincon(@func,x0,[],[],[],[],lb,ub,@NlinCons,OPTIONS);

 [ff fff]=NlinCons(x);

fprintf('\n\n计算结果:\n');

fprintf(' \t目标函数值 fval = %.6e\n\n',fval);

fprintf(' \t约束条件1函数值 f1 = %.6e\n',ff(1));

fprintf(' \t约束条件2函数值 f2 = %.6e\n',ff(2));

fprintf(' \t约束条件3函数值 f3 = %.6e\n',ff(3));

fprintf(' \t约束条件4函数值 f4 = %.6e\n',ff(4));

fprintf('\n\t参数 x1 = %.6f',x(1))

fprintf('\n\t参数 x2 = %.6f',x(2))

fprintf('\n\t参数 x3 = %.6f',x(3))

%---------------------------------------------------------

function[C,Ceq]=NlinCons(x)

C=[1.35*10^(5)/x(1)+12.76*10^(6)*x(2)^2-6.9*10^5*x(1)*x(2)/(12.56*x(1)^3*x(2)-12.56*x(1)^4);x(1)^2*10^10/(x(2)*x(3));(x(1)*x(3)+x(1))/x(2)+86.25*10^3/x(2);x(1)*x(3)+1.5*x(1)-215];

Ceq=[];

%C=[];

%Ceq=[1.35*10^(5)/x(1)+12.76*10^(6)*x(2)^2-6.9*10^5*x(1)*x(2)/(12.56*x(1)^3*x(2)-12.56*x(1)^4);x(1)^2*10^10/(x(2)*x(3));(x(1)*x(3)+x(1))/x(2)+86.25*10^3/x(2);x(1)*x(3)+1.5*x(1)-215];

%---------------------------------------------------------

function f=func(x)

f=0.5*1.935*10^(-5)*x(1)^2*x(2)*(x(3)+1.5)+0.5*5.75*x(2)^3*x(3)/x(1)^4;

toc

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2014-11-06 09:57:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答