greetrunners

金蟲 (著名寫手)

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

[求助] MATLAB, 代碼錯誤，請大俠幫忙看看。

[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ù)此樓

» 猜你喜歡

0703化學(xué)調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有10人回復(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ù)

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

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

至尊木蟲 (著名寫手)

驃騎將軍

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

【答案】應(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寫程序請發(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小時
蟲號: 1632555
注冊: 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

至尊木蟲 (著名寫手)

驃騎將軍

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

引用回帖:

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寫程序請發(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小時
蟲號: 1632555
注冊: 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小時
蟲號: 1632555
注冊: 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ā)表

[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +5	小材化本科 2026-03-18	5/250	2026-03-18 21:53 by 無懈可擊111
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant