24小时热门版块排行榜    

查看: 727  |  回复: 4

yarlyju

金虫 (小有名气)

[求助] 帮忙看下1stopt代码

nrtl方程拟合,怎么不行?


Parameters x(1),x(2),x(3);
Variable x1 ,x2 ,r,T;
ConstStr
tao21=(x(2)-x(1))/(8.314*T),
tao12=(x(1)-x(2))/(8.314*T),
g21=exp(-x(3)*tao21),
g12=exp(-x(3)*tao12)
Function r=exp(x2^2*(tao21*g21^2/(x1+x2*g21)^2+tao12*g12^2/(x2+x1*g12)^2));
data;
0.0317        0.9683        4.7509        366.29
0.0424        0.9576        4.717        364.32
0.0863        0.9137        3.5852        360.33
0.13        0.87        2.905        358.1
0.1666        0.8334        2.5357        357.16
0.2137        0.7863        2.1719        356.07
0.293        0.707        1.7506        354.97
0.3525        0.6475        1.5273        354.5
0.395        0.605        1.4167        353.99
0.4531        0.5469        1.3164        353.4
0.506        0.494        1.2182        353.01
0.5629        0.4371        1.1619        352.64
0.6142        0.3858        1.1126        352.33
0.6395        0.3605        1.0951        352.15
0.6794        0.3206        1.0715        351.95
0.724        0.276        1.0497        351.77
0.774        0.226        1.0267        351.57
0.8436        0.1564        1.0082        351.48
0.8612        0.1388        1.0058        351.44
0.902        0.098        1.0018        351.42
0.9464        0.0536        0.9941        351.48
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
csgt0: 金币+2, 多谢应助 2012-08-03 11:48:05
yarlyju: 金币+3, ★★★★★最佳答案 2012-08-03 15:08:06
内容已删除
2楼2012-08-03 11:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yarlyju

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-08-03 11:32:37
参数和变量定义重名了,另外:
tao21=(p2-p1)/(8.314*T),
tao12=(p1-p2)/(8.314*T),

可简化为:
tao21=p1/(8.314*T),
tao12=-tao21,

否则会出现参数不唯一现象

Parameters p1,p2,p3;
Variable x1 ,x2 ...

谢谢
3楼2012-08-03 15:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yarlyju

金虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-08-03 11:32:37
参数和变量定义重名了,另外:
tao21=(p2-p1)/(8.314*T),
tao12=(p1-p2)/(8.314*T),

可简化为:
tao21=p1/(8.314*T),
tao12=-tao21,

否则会出现参数不唯一现象

Parameters p1,p2,p3;
Variable x1 ,x2 ...

x1和x(1)会被识别成一个符号吗
4楼2012-08-03 15:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应该是的。
5楼2012-08-03 15:58:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yarlyju 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见