孤竹行者

應(yīng)助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時(shí)
蟲(chóng)號(hào): 825902
注冊(cè): 2009-08-11
性別: GG
專(zhuān)業(yè): 能源化工

[求助] matlab非線性擬合求助

小弟是個(gè)matlab菜鳥(niǎo)，現(xiàn)在有一個(gè)方程想通過(guò)matlab進(jìn)行擬合，得到其中的參數(shù)，文獻(xiàn)中說(shuō)是用Nelder-Mead simplex method，但是自己不知道怎么編寫(xiě)程序，望哪位大俠予以幫助，提供matlab程序，不勝感激！文獻(xiàn)中的公式和實(shí)驗(yàn)數(shù)據(jù)如附件所示。

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 方程及文獻(xiàn)數(shù)據(jù).doc

2013-04-05 23:28:45, 36.5 K

是你的終究是你的

1樓 2013-04-05 23:29:53

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

孤竹行者

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時(shí)
蟲(chóng)號(hào): 825902
注冊(cè): 2009-08-11
性別: GG
專(zhuān)業(yè): 能源化工

自己回復(fù)一下，千萬(wàn)別沉啊！

贊一下

回復(fù)此樓

是你的終究是你的

2樓2013-04-07 15:03:47

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

littlepoint

木蟲(chóng) (職業(yè)作家)

應(yīng)助: 9 (幼兒園)
金幣: 12070.6
散金: 8005
紅花: 12
帖子: 3095
在線: 334.6小時(shí)
蟲(chóng)號(hào): 1065217
注冊(cè): 2010-07-27
專(zhuān)業(yè): 工程熱力學(xué)

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

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與！ 2013-04-07 23:18:21
孤竹行者: 金幣+5, ★★★很有幫助 2013-04-15 12:36:21

這種問(wèn)題自己先試著解一下，看看基本原理。這種上來(lái)就求源代碼的搞法不易于自己的提高喲。http://www.pudn.com/downloads62/sourcecode/math/detail213534.html

贊一下

回復(fù)此樓

3樓2013-04-07 15:27:05

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

孤竹行者

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時(shí)
蟲(chóng)號(hào): 825902
注冊(cè): 2009-08-11
性別: GG
專(zhuān)業(yè): 能源化工

引用回帖:

3樓: Originally posted by littlepoint at 2013-04-07 15:27:05
這種問(wèn)題自己先試著解一下，看看基本原理。這種上來(lái)就求源代碼的搞法不易于自己的提高喲。http://www.pudn.com/downloads62/sourcecode/math/detail213534.html

謝謝你的回復(fù)，不過(guò)我主要不是研究這方面的，而且我要急著用這個(gè)程序，一時(shí)間自己也沒(méi)有時(shí)間和精力再?gòu)念^學(xué)習(xí)愛(ài)，先希望高手們能幫我寫(xiě)出這個(gè)程序來(lái)，在下不勝感激！

贊一下

回復(fù)此樓

是你的終究是你的

4樓2013-04-10 08:25:17

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

孤竹行者

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時(shí)
蟲(chóng)號(hào): 825902
注冊(cè): 2009-08-11
性別: GG
專(zhuān)業(yè): 能源化工

各位大俠，幫幫忙吧！小弟拜謝了！

贊一下

回復(fù)此樓

是你的終究是你的

5樓2013-04-15 12:35:30

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

sujiguang

鐵蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 292.3
散金: 78
帖子: 239
在線: 52.5小時(shí)
蟲(chóng)號(hào): 1130581
注冊(cè): 2010-10-24
專(zhuān)業(yè): 藥物化學(xué)

NRTL模型的MATLAB程序：
function f=nrtl(g)
syms r1 r2 r3 r4 a t12 t21 R g21 g12 x2 x4 s1 s2 x1 x3 T;
T=298.15;
x1=O.999889456;
x3=O.01617526;
x2=l-x1;
x4=l-x3;
a=O.20;
R=8.314;
t21=g(1)/(R*T);
t12=g(2)/(R*T);
g21=exp(-a*t21);
g12=exp(-a*t12);
r1=exp(x2^2*(t12*(g12/(x1+x2*g21))^2+t12*g12/(x2+x1*g12)^2));
r2=exp(x1^2*(t12*(g12/(x2+x1*g12))^2+t21*g21/(x1+x2*g21)^2));
r3=exp(x4^2*(t12*(g12/(x3+x4*g21))^2+t12*g12/(x4+x3*g12)^2));
r4=exp(x3^2*(t12*(g12/(x4+x3*g12))^2+t21*g21/(x3+x2*g21)^2));
s1=x1*r1-x3*r3;
s2=x2*r2-x4*r4;
f=s1^2+s2^2
s1
s2
r1
r2
r3
r4
[x,fval]=fiminsearch(x1,x2)

不懂額，你有的話(huà)發(fā)我一份吧，我也是急用額，想要NRTL和UNIQUAC的計(jì)算活度系數(shù)的matlab程序

贊一下(1人)

回復(fù)此樓

6樓2013-07-08 11:30:39

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

孤竹行者

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時(shí)
蟲(chóng)號(hào): 825902
注冊(cè): 2009-08-11
性別: GG
專(zhuān)業(yè): 能源化工

引用回帖:

6樓: Originally posted by sujiguang at 2013-07-08 11:30:39
NRTL模型的MATLAB程序：
function f=nrtl(g)
syms r1 r2 r3 r4 a t12 t21 R g21 g12 x2 x4 s1 s2 x1 x3 T;
T=298.15;
x1=O.999889456;
x3=O.01617526;
x2=l-x1;
x4=l-x3;
a=O.20;
R=8.314;
t21=g(1)/(R*T ...

我這有一個(gè)自己編的程序，初步計(jì)算了一下wilson方程還算可以，你把里面的方程改成你的NRTL和UNIQUAC方程看看行不行！
function f = wilson1(Input)
a12=Input(1);
a13=Input(2);
a21=Input(3);
a23=Input(4);
a31=Input(5);
a32=Input(6);
b12=Input(7);
b13=Input(8);
b21=Input(9);
b23=Input(10);
b31=Input(11);
b32=Input(12);
F=0;
T=[ 298.45 302.7 308.05 313.35 317.75];
X1=[0.0031549 0.0037933 0.0048003 0.0061385 0.0074595];
X2=[0 0 0 0 0 0 0 0 0 ];

for i=1:length(X1)

t=T(i);
x1=X1(i);
x2=X2(i);
x3=1-x2;

lnr1=1-log(x1+x2*(a12+b12/t)+x3*(a13+b13/t))-x1/(x1+x2*(a12+b12/t)+x3*(a13+b13/t))-x2*(a21+b21/t)/(x1*(a21+b21/t)+x2+x3*(a23+b23/t))-x3*(a31+b31/t)/(x1*(a31+b31/t)+x2*(a32+b32/t)+x3);

A=21;
B=9050;
tcal=B/(A-lnr1-log(x1));

v=(tcal-t).^2;

F=F+v;
ff=(F/(length(X1)-1))^0.5;
end
f=ff;
有問(wèn)題再多多交流！

贊一下(1人)

回復(fù)此樓

是你的終究是你的

7樓2013-07-10 14:10:04

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

sujiguang

鐵蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 292.3
散金: 78
帖子: 239
在線: 52.5小時(shí)
蟲(chóng)號(hào): 1130581
注冊(cè): 2010-10-24
專(zhuān)業(yè): 藥物化學(xué)

引用回帖:

7樓: Originally posted by 孤竹行者 at 2013-07-10 14:10:04
我這有一個(gè)自己編的程序，初步計(jì)算了一下wilson方程還算可以，你把里面的方程改成你的NRTL和UNIQUAC方程看看行不行！
function f = wilson1(Input)
a12=Input(1);
a13=Input(2);
a21=Input(3);
a23=Input(4); ...

，還是不會(huì)，哎，我再琢磨琢磨哈！為什么都是input未定義

贊一下

回復(fù)此樓

8樓2013-07-10 19:21:17

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

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

[考研] 一志愿蘇州大學(xué)材料工程（085601）專(zhuān)碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專(zhuān)利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 268求調(diào)劑 +6	好運(yùn)連綿不絕 2026-03-12	7/350	2026-03-17 14:56 by 呦呦憂(yōu)郁
[考研] 工科材料085601 279求調(diào)劑 +3	困于星晨 2026-03-17	3/150	2026-03-17 14:08 by ms629
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4	xujiaoszu 2026-03-11	9/450	2026-03-17 10:29 by xujiaoszu
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 312求調(diào)劑 +3	陌宸希 2026-03-16	4/200	2026-03-16 15:06 by peike
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 化工學(xué)碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu
[考博] 2026年博士申請(qǐng) +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43