greetrunners

金蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 487.8
散金: 544
紅花: 6
帖子: 1134
在線: 190.1小時(shí)
蟲號(hào): 1632555
注冊(cè): 2012-02-21
性別: GG
專業(yè): 分離過程

[求助] MATLAB, 代碼錯(cuò)誤，請(qǐng)大俠幫忙看看。

[t12 t21 t13 t31 t23 t32]=lsqnonlin(@lsqfun3,1,1,1,1,1,1)
Undefined function or variable 'a2'.

Error in lsqfun3 (line 8)
k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21+t31))./((a1+a2.*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2);

Error in lsqnonlin (line 197)
initVals.F = feval(funfcn{3},xCurrent,varargin{:});

Caused by:
Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot continue.

為了減少輸出的參數(shù)，編寫以下函數(shù)，出現(xiàn)了上述問題。
function obj=lsqfun3(t12,t21,t13,t31,t23,t32)
%t12,t21,t13,t31,t23,t32，是所需要關(guān)聯(lián)的NRTL方程參數(shù)；
%a1,a2,a3,是組分1,2,3，液相平衡數(shù)據(jù)；
%y1,y2,y3,是組分1,2,3，氣相平衡數(shù)據(jù)；
%ps1，ps2，ps3，是組分1,2,3的飽和蒸汽壓；
%p==101.35;
p=101.35;
k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21+t31))./((a1+a2.*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2);
k2=(a2.^2*t12*exp(-0.3*t12)+a2.*a3*exp(-0.3*t12)*exp(-0.3*t32)*(t12-t32))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2);
k3=(a3.^2*t13*exp(-0.3*t13)+a2.*a3*exp(-0.3*t13)*exp(-0.3*t23)*(t13-t23))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2);
k4=(a1.^2*t21*exp(-0.3*t21)+a1.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21-t31))./((a1+a2*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2);
k5=(a1.^2*t12*exp(-0.3*t12)^2+a3.^2*t32*exp(-0.3*t32)^2+a1.*a3*exp(-0.3*t12)*exp(-0.3*t32)*(t12+t32))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2);
k6=(a3.^2*t23*exp(-0.3*t23)+a1.*a3*exp(-0.3*t13)*exp(-0.3*t23)*(t23-t13))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2);
k7=(a1.^2*t31*exp(-0.3*t31)+a1.*a2*exp(-0.3*t21)*exp(-0.3*t31)*(t31-t21))./((a1+a2*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2);
k8=(a2.^2*t32*exp(-0.3*t32)^2+a1.*a2*exp(-0.3*t12)*exp(-0.3*t32)*(t32-t12))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2);
k9=(a1.^2*t13*exp(-0.3*t13)^2+a2.^2*t23*exp(-0.3*t23)^2+a1.*a2*exp(-0.3*t13)*exp(-0.3*t23)*(t13+t23))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2);
gama1=exp(k1+k2+k3);
gama2=exp(k4+k5+k6);
gama3=exp(k7+k8+k9);
y1c=ps1.*gama1.*a1./p;
y2c=ps2*gama2*a2./p;
y3c=ps3*gama3*a3./p;
obj=sum(sum((y1c-y1).^2+(y2c-y2).^2+(y3c-y3).^2));
>> y1=[0.8907
0.8924
0.733
0.7731
0.698
0.7778
0.7622
0.7796
0.7934];
y2=[0.1082
0.1067
0.6626
0.5623
0.732
0.4326
0.3978
0.3766
0.3473];
y3=[0
0
0
0
0
0
0
0
0];
a1=[0.8918
0.8933
0.3374
0.4377
0.268
0.5674
0.6022
0.6234
0.6527];
a2=[0.10723
0.10649
0.62613
0.53391
0.70928
0.42518
0.39054
0.37007
0.34364];
a3=[0.00097
0.00021
0.03647
0.02839
0.02272
0.00742
0.00726
0.00653
0.00366];

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

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

請(qǐng)大俠幫忙看一下matlab利用lsqcurvefit曲線擬合問題已經(jīng)有5人回復(fù)
請(qǐng)大俠幫忙看一下我這個(gè)matlab小程序結(jié)果為什么不返回矩陣已經(jīng)有7人回復(fù)
MATLAB關(guān)于lsqcurvefit程序幾個(gè)錯(cuò)誤，求大神指導(dǎo)一下已經(jīng)有6人回復(fù)
matlab中關(guān)于非線性擬合，用nlinfit.程序總出錯(cuò)。謝謝各位大俠已經(jīng)有3人回復(fù)
用matlab 的 fminu 函數(shù)做非線性最優(yōu)化出現(xiàn)錯(cuò)誤，求高手幫忙已經(jīng)有4人回復(fù)
幫忙給看看這個(gè)matlab優(yōu)化函數(shù) 問題已經(jīng)有8人回復(fù)
幫忙修改Matlab中出現(xiàn)的錯(cuò)誤已經(jīng)有4人回復(fù)
MATLAB運(yùn)行錯(cuò)誤已經(jīng)有5人回復(fù)
采用matlab求解微分方程遇到問題，請(qǐng)幫忙分析一下錯(cuò)誤-重謝已經(jīng)有3人回復(fù)
下載了MATLAB2012可是安裝不了。大俠幫幫忙謝謝已經(jīng)有13人回復(fù)
怎么用matlab畫極坐標(biāo)三維圖？各位大俠幫幫忙��！已經(jīng)有5人回復(fù)
求助matlab一程序運(yùn)行的問題，幫忙看看哪里不對(duì) 已經(jīng)有4人回復(fù)
求大俠幫忙解決一個(gè)B類錯(cuò)誤��！謝謝已經(jīng)有5人回復(fù)
最優(yōu)化計(jì)算求助�。�1stopt或者M(jìn)atlab）已經(jīng)有16人回復(fù)
matlab 問題求助，是否計(jì)算結(jié)果錯(cuò)誤，如何改正？已經(jīng)有3人回復(fù)
matlab程序錯(cuò)誤已經(jīng)有5人回復(fù)
用matlab求解方程出問題，請(qǐng)幫忙看看已經(jīng)有3人回復(fù)
幫忙看看這個(gè)平面圖用matlab怎么畫已經(jīng)有19人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
matlab程序運(yùn)行總是出現(xiàn)下面的錯(cuò)誤，請(qǐng)各位高手指點(diǎn)一下，謝謝了！急急急已經(jīng)有5人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】MATLAB錯(cuò)誤已經(jīng)有6人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)

努力！努力！努力！

1樓 2013-01-30 10:48:24

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

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲號(hào): 696514
注冊(cè): 2009-02-05
專業(yè): 計(jì)算機(jī)軟件

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
greetrunners: 金幣+10, ★★★很有幫助, 謝謝 2013-01-31 09:39:29
xzhdty: 金幣+1, 謝謝驃騎將軍 2013-01-31 11:08:38

function obj=lsqfun3(t12,t21,t13,t31,t23,t32)
%t12,t21,t13,t31,t23,t32，是所需要關(guān)聯(lián)的NRTL方程參數(shù)；
%a1,a2,a3,是組分1,2,3，液相平衡數(shù)據(jù)；
%y1,y2,y3,是組分1,2,3，氣相平衡數(shù)據(jù)；
%ps1，ps2，ps3，是組分1,2,3的飽和蒸汽壓；
%p==101.35;
p=101.35;
k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)

===============
函數(shù)中a2,a3沒定義啊,怎么就開始用了

贊一下

回復(fù)此樓

matlab/VB/python/c++/Java寫程序請(qǐng)發(fā)QQ郵件:790404545@qq.com

2樓2013-01-30 17:17:00

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

greetrunners

金蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 487.8
散金: 544
紅花: 6
帖子: 1134
在線: 190.1小時(shí)
蟲號(hào): 1632555
注冊(cè): 2012-02-21
性別: GG
專業(yè): 分離過程

引用回帖:

2樓: Originally posted by libralibra at 2013-01-30 17:17:00
function obj=lsqfun3(t12,t21,t13,t31,t23,t32)
%t12,t21,t13,t31,t23,t32，是所需要關(guān)聯(lián)的NRTL方程參數(shù)；
%a1,a2,a3,是組分1,2,3，液相平衡數(shù)據(jù)；
%y1,y2,y3,是組分1,2,3，氣相平衡數(shù)據(jù)；
%ps1，ps2，ps3，是 ...

應(yīng)該如何定義？

回復(fù)此樓

努力！努力！努力！

3樓2013-01-30 19:08:17

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

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲號(hào): 696514
注冊(cè): 2009-02-05
專業(yè): 計(jì)算機(jī)軟件

引用回帖:

3樓: Originally posted by greetrunners at 2013-01-30 19:08:17
應(yīng)該如何定義？...

這...你自己知道吧?
起碼函數(shù)的定義應(yīng)該是

function ret_value = foo(input_value)
param1 = ...
param2 = ...
% then use param1 and param2 to calculate the others
...
end

你得先定義參數(shù),然后使用,你的函數(shù)前面幾行看著寫的挺清楚,可是那是注釋,你沒有定義a2,a3....幾乎所有要用的參數(shù)

贊一下

回復(fù)此樓

matlab/VB/python/c++/Java寫程序請(qǐng)發(fā)QQ郵件:790404545@qq.com

4樓2013-01-30 20:43:55

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

greetrunners

金蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 487.8
散金: 544
紅花: 6
帖子: 1134
在線: 190.1小時(shí)
蟲號(hào): 1632555
注冊(cè): 2012-02-21
性別: GG
專業(yè): 分離過程

引用回帖:

可以幫我修改一下嗎？

贊一下

回復(fù)此樓

努力！努力！努力！

5樓2013-01-31 12:14:29

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

greetrunners

金蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 487.8
散金: 544
紅花: 6
帖子: 1134
在線: 190.1小時(shí)
蟲號(hào): 1632555
注冊(cè): 2012-02-21
性別: GG
專業(yè): 分離過程

引用回帖:

就是NRTL方程參數(shù)關(guān)聯(lián)，k1,k2,---k9,是方程拆開的部分，

贊一下

回復(fù)此樓

努力！努力！努力！

6樓2013-01-31 12:15:27

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

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

[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 321求調(diào)劑 +3	何潤采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 314求調(diào)劑 +8	無懈可擊的巨人 2026-03-12	8/400	2026-03-18 14:50 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級(jí)證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike