版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

muyouyou

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 234.8
散金: 5
帖子: 22
在線: 21.4小時(shí)
蟲(chóng)號(hào): 3250631
注冊(cè): 2014-06-03
專業(yè): 交通工程

[求助] matlab優(yōu)化求解中遇到問(wèn)題，望大神指點(diǎn)。。。已有1人參與

在matlab中
目標(biāo)函數(shù)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')
運(yùn)行結(jié)果為：??? 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，總是出現(xiàn)這種錯(cuò)誤，不知是怎么回事。望高手指點(diǎn)。。。

回復(fù)此樓

» 猜你喜歡

311求調(diào)劑已經(jīng)有4人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)
354求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有3人回復(fù)
274求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有15人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有8人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有10人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有8人回復(fù)
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷已經(jīng)有10人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求大神用matlab求下圖陰影部分的面積和周長(zhǎng)�。� 已經(jīng)有4人回復(fù)
求matlab大神指點(diǎn)，關(guān)于logistic擬合已經(jīng)有13人回復(fù)
現(xiàn)在用MATLAB求解一個(gè)循環(huán)后畫(huà)出曲線遇到問(wèn)題，求大神指點(diǎn)��！已經(jīng)有7人回復(fù)
matlab最優(yōu)化問(wèn)題求解中，復(fù)雜的目標(biāo)函數(shù)和約束條件無(wú)法直接表示的問(wèn)題已經(jīng)有3人回復(fù)
matlab求解，求大神已經(jīng)有1人回復(fù)
matlab求解非線性方程組報(bào)錯(cuò)，請(qǐng)各位大神指點(diǎn)！已經(jīng)有8人回復(fù)
matlab畫(huà)圖問(wèn)題�。�！急！希望會(huì)的大神來(lái)看看�。� 已經(jīng)有11人回復(fù)
求大神幫忙matlab擬合函數(shù)求參數(shù) 已經(jīng)有9人回復(fù)
菜鳥(niǎo)求助，老師給了一個(gè)分子讓優(yōu)化，可總是錯(cuò)誤，求大神指點(diǎn) 已經(jīng)有22人回復(fù)
求大神用matlab給我處理一組數(shù)據(jù)，有人沒(méi)有��？已經(jīng)有15人回復(fù)
Matlab 哪位大神解決此問(wèn)題重金懸賞已經(jīng)有4人回復(fù)
求助：遺傳算法的matlab編程求解資源優(yōu)化問(wèn)題已經(jīng)有7人回復(fù)
MATLAB微分方程參數(shù)擬合問(wèn)題，求大神已經(jīng)有7人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

約束條件是非線性的。
對(duì)于非線性約束情況，fmincon按附圖1中格式調(diào)用。
附圖1中正是一個(gè)例子。

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2014-11-05 21:21:38
約束條件是非線性的。
對(duì)于非線性約束情況，fmincon按附圖1中格式調(diào)用。
附圖1中正是一個(gè)例子。

重新上傳附圖1。

matlab優(yōu)化求解中遇到問(wèn)題，望大神指點(diǎn)。。。

附圖1.png

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

muyouyou

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 234.8
散金: 5
帖子: 22
在線: 21.4小時(shí)
蟲(chóng)號(hào): 3250631
注冊(cè): 2014-06-03
專業(yè): 交通工程

引用回帖:

3樓: Originally posted by 月只藍(lán) 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
我按你給的步驟計(jì)算了一下，結(jié)果是??? Error: File: F:\software\matlab\work\Fmincon.m Line: 3 Column: 1
Missing variable or function..是怎么回事呢？謝謝你了。。。

贊一下

回復(fù)此樓

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

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計(jì)算結(jié)果:\n');

fprintf(' \t目標(biāo)函數(shù)值 fval = %.6e\n\n',fval);

fprintf(' \t約束條件1函數(shù)值 f1 = %.6e\n',ff(1));

fprintf(' \t約束條件2函數(shù)值 f2 = %.6e\n',ff(2));

fprintf(' \t約束條件3函數(shù)值 f3 = %.6e\n',ff(3));

fprintf(' \t約束條件4函數(shù)值 f4 = %.6e\n',ff(4));

fprintf('\n\t參數(shù) x1 = %.6f',x(1))

fprintf('\n\t參數(shù) x2 = %.6f',x(2))

fprintf('\n\t參數(shù) 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人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

5樓的代碼全部復(fù)制進(jìn)一個(gè)m文件，運(yùn)行即可。運(yùn)行成功之前，一個(gè)字母都不要改。
運(yùn)行成功后，你再看看，是否可以改改初值。
特別需要注意的是，請(qǐng)檢查約束條件中的等式約束和非等式約束是否定義正確。

贊一下

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

6樓2014-11-06 09:59:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

muyouyou

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 234.8
散金: 5
帖子: 22
在線: 21.4小時(shí)
蟲(chóng)號(hào): 3250631
注冊(cè): 2014-06-03
專業(yè): 交通工程

引用回帖:

6樓: Originally posted by 月只藍(lán) at 2014-11-06 09:59:48
5樓的代碼全部復(fù)制進(jìn)一個(gè)m文件，運(yùn)行即可。運(yùn)行成功之前，一個(gè)字母都不要改。
運(yùn)行成功后，你再看看，是否可以改改初值。
特別需要注意的是，請(qǐng)檢查約束條件中的等式約束和非等式約束是否定義正確。...

知道了，特別感謝。。。

回復(fù)此樓

7樓2014-11-11 19:14:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

muyouyou

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 234.8
散金: 5
帖子: 22
在線: 21.4小時(shí)
蟲(chóng)號(hào): 3250631
注冊(cè): 2014-06-03
專業(yè): 交通工程

引用回帖:

x1=[21,26,31,21,26,31,21,26,31,21,26,31,21,26,31,21,26,31];
x2=[133,143,153,133,143,153,143,153,133,153,133,143,143,153,133,153,133,143];
x3=[3.9,4.4,4.9,4.4,4.9,3.9,3.9,4.4,4.9,4.9,3.9,4.4,4.9,3.9,4.4,4.4,4.9,3.9];
y=[40.69803,40.24597,39.66466,39.65027,39.3111,41.64728,39.43865,39.07176,42.09825,36.28181,42.55311,41.77374,37.45602,40.11952,43.03311,37.21668,40.57048,42.8215];
如何擬合成y=f（x1，x2，x3）的多項(xiàng)式，望高手指點(diǎn)指點(diǎn)。。。

贊一下

回復(fù)此樓

8樓2014-11-11 19:28:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 muyouyou 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 321求調(diào)劑 +3	何潤(rùn)采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +3	Ymlll 2026-03-18	5/250	2026-03-18 19:32 by Ymlll
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 【0856】化學(xué)工程（085602）313 分，本科學(xué)科評(píng)估A類院�；瘜W(xué)工程與工藝，誠(chéng)求調(diào)劑 +7	小劉快快上岸 2026-03-11	8/400	2026-03-17 16:57 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655

24小時(shí)熱門(mén)版塊排行榜

[求助] matlab優(yōu)化求解中遇到問(wèn)題，望大神指點(diǎn)。。。 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] matlab優(yōu)化求解中遇到問(wèn)題，望大神指點(diǎn)。。。已有1人參與