24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2973  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08600生物与医药-327 +9 18755400796 2026-04-05 9/450 2026-04-06 22:35 by 52305043001
[考研] 085600材料与化工专硕329 求调剂 +11 额cc 2026-04-06 11/550 2026-04-06 22:33 by chenzhimin
[考研] 生物学学硕求调剂:351分一志愿南京师范大学生物学专业 +6 …~、王…~ 2026-04-06 7/350 2026-04-06 18:54 by macy2011
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 0854电子信息319求调剂(接受跨专业调剂) +3 星星不眨眼喽 2026-04-05 3/150 2026-04-05 20:20 by 啵啵啵0119
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 307分材料专业求调剂 +7 Hll胡 2026-04-05 7/350 2026-04-05 18:47 by 无际的草原
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666!
[考研] 工科求调剂 +15 11ggg 2026-04-03 15/750 2026-04-05 16:24 by zzx2138
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
信息提示
请填处理意见