24小时热门版块排行榜    

查看: 2866  |  回复: 11

fanghongmei

金虫 (初入文坛)

[求助] origin自定义拟合函数一直拟合不了曲线 已有1人参与

各位朋友,我想利用Henderson–Hasselbalch equation: pka=pH−log[(FImax − FI)/(FI − FImin)] 求算出pka的值,利用origin,函数可以编译成功,可是结果却是迭代次数为0,不收敛。我上传一下我的数据,麻烦各位朋友帮我想想办法,再此先谢谢大家了
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Book1.xlsx
  • 2015-03-23 17:13:46, 8.66 K

» 猜你喜欢

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

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:20:35
根据你问题的描述,这不是拟合。
pka可以用表示出来,直接计算不就可以了吗。pH FI数据都有的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-03-23 22:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanghongmei

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by 月只蓝 at 2015-03-23 22:27:10
根据你问题的描述,这不是拟合。
pka可以用表示出来,直接计算不就可以了吗。pH FI数据都有的。

不好意思,可能我表达的不是很清楚,FImax, FImin并不是确定的值,只能把它看做上面公式中的参数
3楼2015-03-24 09:21:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fanghongmei: 金币+10, ★★★很有帮助 2015-03-24 18:46:19
fegg7502: 金币+1, 鼓励交流 2015-06-15 11:20:51
origin做可能初值不好找,下面用 1stopt软件做。
代码:
CODE:
Variable pH,FI;
Function   pH=pka+log( (FImax-FI)/(FI-FImin)   ) ;
Data;
//pH FI
0.83        949.49652
1        945.09521
1.19        937.23456
1.39        919.08655
1.6        906.92108
2.01        863.60083
2.19        789.7934
2.4        633.59326
2.6        525.37335
2.8        386.98236
3        272.30063
3.21        189.74347
3.4        137.99492
3.6        108.07579
3.8        75.06548
4.1        42.93212
4.5        22.65061
5.1        6.69859
6.09        2.2117
7.31        3.94633

计算结果:
均方差(RMSE): 0.425891955483537
残差平方和(SSE): 3.62767915491182
相关系数(R): 0.968682222761761
相关系数之平方(R^2): 0.938345248694667
决定系数(DC): 0.934879153976819
卡方系数(Chi-Square): 0.312507477743333
F统计(F-Statistic): 273.948302748938

参数        最佳估算
----------        -------------
pka                2.84041133362774
fimax              958.492522526198
fimin              1.91435094772762

====== 结果输出 =====

No        实测值ph        计算值ph
1        0.83        0.8178440
2        1        0.9928340
3        1.19        1.1969727
4        1.39        1.4735225
5        1.6        1.5961688
6        2.01        1.8822902
7        2.19        2.1710646
8        2.4        2.5516637
9        2.6        2.7581361
10        2.8        3.0118979
11        3        3.2448723
12        3.21        3.4524330
13        3.4        3.6206925
14        3.6        3.7440763
15        3.8        3.9223610
16        4.1        4.1891263
17        4.5        4.4948834
18        5.1        5.1391414
19        6.09        6.3477302
20        7.31        5.5122891

该软件获取链接:http://muchong.com/bbs/viewthread.php?tid=4586018
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2015-03-24 13:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanghongmei

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by 月只蓝 at 2015-03-24 13:40:42
origin做可能初值不好找,下面用 1stopt软件做。
代码:

Variable pH,FI;
Function   pH=pka+log( (FImax-FI)/(FI-FImin)   ) ;
Data;
//pH FI
0.83        949.49652
1        945.09521
1.19        937.23456
1.39        919.086 ...

谢谢您的热心帮忙,1stop还没有用过,我先研究研究,非常的感谢!
5楼2015-03-24 18:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanghongmei

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by 月只蓝 at 2015-03-24 13:40:42
origin做可能初值不好找,下面用 1stopt软件做。
代码:

Variable pH,FI;
Function   pH=pka+log( (FImax-FI)/(FI-FImin)   ) ;
Data;
//pH FI
0.83        949.49652
1        945.09521
1.19        937.23456
1.39        919.086 ...

朋友,还想问一下,上述链接中的1stopt软件的主界面中并没有可以输入代码的地方,而且打开之后就关不掉了,是我电脑的问题吗?
6楼2015-03-24 19:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by fanghongmei at 2015-03-24 18:47:59
谢谢您的热心帮忙,1stop还没有用过,我先研究研究,非常的感谢!...

管理员身份运行即可

[ 发自小木虫客户端 ]
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2015-03-24 23:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanghongmei

金虫 (初入文坛)

引用回帖:
7楼: Originally posted by 月只蓝 at 2015-03-24 23:13:37
管理员身份运行即可
...

你好,我用你给的代码从新算了一下,可是总是出现附件图片里的情况,这是为什么呀?
origin自定义拟合函数一直拟合不了曲线
QQ截图20150328091447.png

8楼2015-03-28 09:19:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
fanghongmei: 金币+5, ★★★很有帮助, 谢谢您的耐心解答,真的非常感谢! 2015-03-28 19:53:46
引用回帖:
8楼: Originally posted by fanghongmei at 2015-03-28 09:19:41
你好,我用你给的代码从新算了一下,可是总是出现附件图片里的情况,这是为什么呀?

QQ截图20150328091447.png
...

复制4楼的代码时,请保证代码的完整性。
我对4楼的代码进行了测试,可以运行。
origin自定义拟合函数一直拟合不了曲线-1
附图1.png


origin自定义拟合函数一直拟合不了曲线-2
计算结果.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
9楼2015-03-28 09:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cy20131120

金虫 (初入文坛)

你好!请问下,求PKA怎么又一个FI出来?FI是什么?不确定度还是什么啊?
10楼2015-06-14 21:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanghongmei 的主题更新
信息提示
请填处理意见