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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 想利用最小二乘法对参数进行非线性拟合。哪种方式快捷哪？
182/2返回列表上一页12
查看: 2235  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

walk1997

金虫 (著名写手)
引用回帖:
8楼: Originally posted by loopool5358 at 2013-06-22 08:58:22
另外，需要说明的是，这几个待拟合参数的都有物理意义，就是都得大于零了。

加了这个约束进行拟合 效果也不是很理想
画出来的图看着还可以 但是参数实际的标准误差很大
而且拟合结果对第6，7个数据点比较敏感
感觉还是提供的数据有些问题
不知道是否如此？
赞一下 回复此楼
11楼2013-06-22 10:35:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loopool5358

银虫 (正式写手)
引用回帖:
10楼: Originally posted by walk1997 at 2013-06-22 10:07:00
另外 你看下 第六行的数据
这数据很突然变大 马上又下降...不觉得奇怪么
呵呵 拟合了下 总感觉有问题。。。。...

非常感谢您的热心
再者，模型和数据都是从文献中获得，只是想验证一下拟合过程。
我觉得第六个和第七个点的数据差别不算大啊。
想问您，1）用什么方法进行拟合的。
2）即使拟合结果不好，能否把结果告知哪。
赞一下 回复此楼
坚持！
12楼2013-06-22 11:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
用的mathematica 比较偷懒 建议看下拟合参数是不是还有其他限制(比如是不是只要某个局部最小的拟合)
Clear["Global`*"];
SetDirectory[NotebookDirectory[]];
test = Import["IV.txt", "Data"];
f = 38.92;
S = 1;
re = NonlinearModelFit[SetAccuracy[test, 17],
   Imax*( 1 - Exp[-f*\[Eta]]) /(K1*Exp[-(1 - \[Alpha]) f*\[Eta]] +
       K2*Exp[-f*\[Eta]] + KM/S + 1), {Imax, K1, K2,
    KM, \[Alpha]}, {\[Eta]}, Method -> NMinimize];
re[{"ParameterTable"}]
re = NonlinearModelFit[
   SetAccuracy[test,
    17], {Imax*(
       1 - Exp[-f*\[Eta]]) /(K1*Exp[-(1 - \[Alpha]) f*\[Eta]] +
        K2*Exp[-f*\[Eta]] + KM/S + 1), Imax > 0, K1 > 0, K2 > 0,
    KM > 0, \[Alpha] > 0}, {Imax, K1, K2, KM, \[Alpha]}, {\[Eta]},
   Method -> NMinimize];
re[{"ParameterTable"}]
Print[Show[ListPlot[SetAccuracy[test, 17]],
   Plot[Normal[re], {\[Eta], 0, 0.35}]]];
test = Delete[test, {{5}, {6}, {7}}];
re = NonlinearModelFit[
   SetAccuracy[test,
    17], {Imax*(
       1 - Exp[-f*\[Eta]]) /(K1*Exp[-(1 - \[Alpha]) f*\[Eta]] +
        K2*Exp[-f*\[Eta]] + KM/S + 1), Imax > 0, K1 > 0, K2 > 0,
    KM > 0, \[Alpha] > 0}, {Imax, K1, K2, KM, \[Alpha]}, {\[Eta]},
   Method -> NMinimize];
re[{"ParameterTable"}]
Print[Show[ListPlot[SetAccuracy[test, 17]],
   Plot[Normal[re], {\[Eta], 0, 0.35}]]];
－------------------------------------
结果，约束在所有参数都大于0的情况
        Estimate        Standard Error        t-Statistic        P-Value
Imax        2.99146        17.0319        0.175639        0.864466
K1        18.7077        3.02763        6.179        0.000162926
K2        0.        122.524        0.        1.
KM        1.1016*10^-6        5.6898        1.93611*10^-7        1.
\[Alpha]        0.153136        0.566881        0.270138        0.793146

－-------------------------
去掉6，7点的拟合
        Estimate        Standard Error        t-Statistic        P-Value
Imax        3.34359        0.590046        5.66666        0.00129905
K1        1.71615        1.76284        0.973518        0.367893
K2        22.5157        0.229166        98.2506        7.49176*10^-11
KM        0.107593        0.204174        0.526967        0.617122
\[Alpha]        0.579124        0.159921        3.62132        0.0110797
赞一下 回复此楼
13楼2013-06-22 12:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
btw: 你的原始数据
0.0008421052631579        0.0098314606741573
0.0538947368421053        0.5997191011235961
0.0842105263157895        1.3960674157303401
0.1204210526315790        2.1039325842696601
0.15410526315789499        2.70365168539326017
0.20042105263157900        3.0379213483146099
0.2526315789473680        2.8609550561797801
0.2686315789473680        2.9789325842696601
0.2863157894736840        2.9789325842696601
0.2981052631578950        2.9887640449438200
0.3149473684210530        2.9985955056179798
0.3351578947368420        2.9985955056179798
0.3553684210526320        2.9985955056179798
0.3747368421052630        2.9887640449438200

－-------------
第6行是3.03 第7行是2.86  和前后对比下 这结果不很奇怪么
不是实验误差的话 这区域附近应该有拐点 可以好好更密集点测量看看
赞一下(1人) 回复此楼
14楼2013-06-22 12:09:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loopool5358

银虫 (正式写手)
引用回帖:
14楼: Originally posted by walk1997 at 2013-06-22 12:09:21
btw: 你的原始数据
0.0008421052631579        0.0098314606741573
0.0538947368421053        0.5997191011235961
0.0842105263157895        1.3960674157303401
0.1204210526315790        2.1039325842696601
0.15410526315789499        2.70 ...

非常感谢您啊
原文献说的是：minimization was done using non-linear least squares estimation。minimization was performed using the quasi-Newton as implemented in Mathcad。
我的理解是，采用非线性最小二乘法获得参数拟合，在拟合过程中采用近似牛顿方法。
不过，您列出来的值，和文献中的值相差比较大。也可能是，我是从文献中取值的，所以有些误差在里面的。
我也下载您说的软件看看。
赞一下 回复此楼
坚持！
15楼2013-06-22 14:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icecream8508

金虫 (著名写手)
引用回帖:
14楼: Originally posted by walk1997 at 2013-06-22 12:09:21
btw: 你的原始数据
0.0008421052631579        0.0098314606741573
0.0538947368421053        0.5997191011235961
0.0842105263157895        1.3960674157303401
0.1204210526315790        2.1039325842696601
0.15410526315789499        2.70 ...

还有，我只是想做个带约束的，非线性最小二乘的拟合。您觉得我用哪个数学软件比较好？
赞一下 回复此楼
16楼2013-06-22 14:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
引用回帖:
16楼: Originally posted by icecream8508 at 2013-06-22 14:25:18
还有，我只是想做个带约束的，非线性最小二乘的拟合。您觉得我用哪个数学软件比较好？...

我个人是喜欢用mathematica 因为有时候有要解析计算的需要
不是特别追求速度或大型数值计算或应用特定包等的话
我觉得做一般的工作 mathematica足以胜任 内部提供的功能很有效率和强大  上手也比较快
上面的原始数据是直接下载你的附件的 可能最后几位有点小差异 但是整体性质肯定不变 会不会你上传的数据不对？
赞一下 回复此楼
17楼2013-06-22 14:49:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
引用回帖:
16楼: Originally posted by icecream8508 at 2013-06-22 14:25:18
还有，我只是想做个带约束的，非线性最小二乘的拟合。您觉得我用哪个数学软件比较好？...

原文献看来是用MathCAD 你也可以尝试下 应该是更简单的工具软件
赞一下 回复此楼
18楼2013-06-22 14:50:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 loopool5358 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +11 不知道叫什么！ 2026-03-15 13/650 2026-03-16 15:10 by 了了了了。。
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 4/200 2026-03-16 13:53 by laoshidan
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 304求调剂（085602一志愿985） +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
信息提示
关闭
请填处理意见
关闭