应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab非线性拟合求助
81/1返回列表
查看: 2072  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

孤竹行者

金虫 (小有名气)

[求助] matlab非线性拟合求助

小弟是个matlab菜鸟，现在有一个方程想通过matlab进行拟合，得到其中的参数，文献中说是用Nelder-Mead simplex method，但是自己不知道怎么编写程序，望哪位大侠予以帮助，提供matlab程序，不胜感激！文献中的公式和实验数据如附件所示。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 方程及文献数据.doc
    • 2013-04-05 23:28:45, 36.5 K

» 收录本帖的淘帖专辑推荐

纳米

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
是你的终究是你的
1楼 2013-04-05 23:29:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孤竹行者

金虫 (小有名气)
自己回复一下，千万别沉啊！
赞一下 回复此楼
是你的终究是你的
2楼2013-04-07 15:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littlepoint

木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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发消息 送TA红花 TA的回帖

孤竹行者

金虫 (小有名气)
引用回帖:
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发消息 送TA红花 TA的回帖

孤竹行者

金虫 (小有名气)
各位大侠，帮帮忙吧！小弟拜谢了！
赞一下 回复此楼
是你的终究是你的
5楼2013-04-15 12:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sujiguang

铁虫 (小有名气)
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)


不懂额，你有的话发我一份吧，我也是急用额，想要NRTL和UNIQUAC的计算活度系数的matlab程序
赞一下(1人) 回复此楼
6楼2013-07-08 11:30:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孤竹行者

金虫 (小有名气)
引用回帖:
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发消息 送TA红花 TA的回帖

sujiguang

铁虫 (小有名气)
引用回帖:
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发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孤竹行者 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭