孤竹行者

金蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時
蟲號: 825902
注冊: 2009-08-11
性別: GG
專業(yè): 能源化工

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

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

回復此樓

» 本帖附件資源列表

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

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

» 收錄本帖的淘帖專輯推薦

納米

» 猜你喜歡

能源材料化學課題組招收碩士研究生8-10名已經(jīng)有6人回復
296求調劑已經(jīng)有5人回復
0703化學調劑，六級已過，有科研經(jīng)歷已經(jīng)有8人回復
268求調劑已經(jīng)有8人回復
材料，紡織，生物（0856、0710），化學招生啦已經(jīng)有7人回復
301求調劑已經(jīng)有5人回復
281求調劑（0805）已經(jīng)有7人回復
材料與化工求調劑已經(jīng)有6人回復
299求調劑已經(jīng)有4人回復
277調劑已經(jīng)有6人回復

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

關于matlab微分、及非線性擬合問題已經(jīng)有5人回復
matlab非線性擬合計算求助已經(jīng)有10人回復
matlab數(shù)據(jù)擬合問題求助已經(jīng)有7人回復
誠懇求助如何做IC50擬合曲線已經(jīng)有7人回復
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復
急切求助用matlab曲線擬合色散方程已經(jīng)有10人回復
求助非線性擬合問題已經(jīng)有5人回復
【求助】Matlab非線性最小二乘擬合活度系數(shù)模型（Willson、NRTL、UNIQUAC）已經(jīng)有17人回復
【求助】急求關于matlab做非線性擬合的函數(shù) 已經(jīng)有8人回復
【求助】如何使得matlab擬合的曲線強制經(jīng)過一個點已經(jīng)有8人回復
【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】已經(jīng)有15人回復
【求助】matlab 遺傳算法工具箱繪圖命令是什么？已經(jīng)有6人回復
【求助】求教matlab解非線性方程組已經(jīng)有9人回復
【求助】matlab有公式的曲線擬合已經(jīng)有7人回復
【求助】如何把origin中曲線的擬合線延長？已經(jīng)有7人回復

是你的終究是你的

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

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

孤竹行者

金蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時
蟲號: 825902
注冊: 2009-08-11
性別: GG
專業(yè): 能源化工

自己回復一下，千萬別沉啊！

贊一下

回復此樓

是你的終究是你的

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

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

littlepoint

木蟲 (職業(yè)作家)

應助: 9 (幼兒園)
金幣: 12070.6
散金: 8005
紅花: 12
帖子: 3095
在線: 334.6小時
蟲號: 1065217
注冊: 2010-07-27
專業(yè): 工程熱力學

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

孤竹行者

金蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時
蟲號: 825902
注冊: 2009-08-11
性別: GG
專業(yè): 能源化工

引用回帖:

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

謝謝你的回復，不過我主要不是研究這方面的，而且我要急著用這個程序，一時間自己也沒有時間和精力再從頭學習愛，先希望高手們能幫我寫出這個程序來，在下不勝感激！

贊一下

回復此樓

是你的終究是你的

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

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

孤竹行者

金蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時
蟲號: 825902
注冊: 2009-08-11
性別: GG
專業(yè): 能源化工

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

贊一下

回復此樓

是你的終究是你的

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

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

sujiguang

鐵蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 292.3
散金: 78
帖子: 239
在線: 52.5小時
蟲號: 1130581
注冊: 2010-10-24
專業(yè): 藥物化學

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)

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

贊一下(1人)

回復此樓

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

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

孤竹行者

金蟲 (小有名氣)

應助: 3 (幼兒園)
金幣: 1103.9
散金: 5
帖子: 137
在線: 52.7小時
蟲號: 825902
注冊: 2009-08-11
性別: GG
專業(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 ...

我這有一個自己編的程序，初步計算了一下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;
有問題再多多交流！

贊一下(1人)

回復此樓

是你的終究是你的

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

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

sujiguang

鐵蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 292.3
散金: 78
帖子: 239
在線: 52.5小時
蟲號: 1130581
注冊: 2010-10-24
專業(yè): 藥物化學

引用回帖:

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

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

贊一下

回復此樓

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

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

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

[考研] 341求調劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 326求調劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 332求調劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 085600材料與化工求調劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 材料與化工專碩調劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 材料工程專碩274一志愿211求調劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[考研] 考研調劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 藥學383 求調劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 286求調劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 0854控制工程 359求調劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 283求調劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[基金申請] NSFC申報書里申請人簡歷中代表性論著還需要在申報書最后的附件里面再上傳一遍嗎 20+5	NSFC2026我來了 2026-03-10	14/700	2026-03-15 23:53 by 不負韶華的虎
[考研] 297一志愿上交085600求調劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 330求調劑 +3	?醬給調劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 279求調劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 285 求調劑資源與環(huán)境一志愿北京化工大學 +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 304求調劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 化工學碩306求調劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu
[考研] 270求調劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu