24小时热门版块排行榜    

查看: 3011  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

orez

铁虫 (初入文坛)

[交流] 如何用最小二乘法拟合y=A1*exp(-x/t1)+A2*exp(-x/t2)+B,该函数的反函数如何求解? 已有3人参与

函数如下,带有2项e指数, 该函数是Origin中的拟合函数之一,用途涉及x射线的衰减领域,怎样用最小二乘法进行曲线拟合呢,最好能够编程实现的。反函数又如何求解,感觉很难求解下去,求教求教!!
如何用最小二乘法拟合y=A1*exp(-x/t1)+A2*exp(-x/t2)+B,该函数的反函数如何求解?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

咖啡渣_95

铁虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
提供个思路
关于拟合,Origin、Matlab、Python都有对应的库,直接调用就行了。实在找不到,可以用线性拟合,把e(-x/t1)\e(-x/t2)看作两个新变量就行了。就不多说了。

关于反函数求解,解析解不知道怎么整,提供一个数值求解的算法(利用函数的单调的递减特性),步骤如下(假设y取值为y0):
1、给x设定一个初值,比如1,这种情况下计算y(x),有三种情况:
    a.如果y(x)=y0,则1就是y0的反解。
    b.如果y(x)>y0,则让x增大
    c.如果y(x)<y0,则让x缩小
     注:由于不可能得到绝对准确的解,实际上a判断条件可以写为:abs( y(x)-y0 ) <= error 或者 abs( y(x)-y0 )/(y(x)+y0) < rsd, error/rsd是能接受的计算误差/相对误差,一般可以设得很小,比如1E-10。

2、x变化后,再重新计算y(x),判断属于a、b、c那种情况,如此重复,直到满足情况a。

为了保证计算速度,最好让x按下面规则增大缩小:
1)如果一直是情况b,就一直让x翻倍
2)如果一直是情况c,就一直让x减半
3)如果出现过情况b,也出现过情况c,最近一次出现b情况对应的x值设为x_min、最近一次出现c情况对应的x值设为x_max,x=sqrt(x_min*x_max)
6楼2020-07-10 13:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

独孤神宇

版主 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
Origin内置函数,直接按照拟合界面的提示操作,不需要另外编写程序。

编程可以考虑使用MATLAB拟合。
数值计算
2楼2020-06-14 09:18:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

orez

铁虫 (初入文坛)

引用回帖:
2楼: originally posted by 独孤神宇 at 2020-06-14 09:18:46
origin内置函数,直接按照拟合界面的提示操作,不需要另外编写程序。

编程可以考虑使用matlab拟合。

现在最主要纠结反函数解不出来。拟合主要想写成可执行程序,方便移植。拟合的话按解析方法,求系数比较困难,用迭代法应该可以实现,就是不知道发散如何。
3楼2020-06-14 11:50:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by orez at 2020-06-14 11:50:59
现在最主要纠结反函数解不出来。拟合主要想写成可执行程序,方便移植。拟合的话按解析方法,求系数比较困难,用迭代法应该可以实现,就是不知道发散如何。...

matlab非线性拟合函数可以实现。

至于反函数,实际上相当于隐函数x求解。

发自小木虫Android客户端
数值计算
4楼2020-06-14 11:55:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 一志愿郑大070303,338分,求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[文学芳草园] 伙伴们,祝我生日快乐吧 +12 myrtle 2026-03-10 17/850 2026-03-13 22:00 by angelyueyi
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +3 自由的_飞翔 2026-03-13 4/200 2026-03-13 10:47 by 冠c哥
[考博] 2026博士申请 +7 起泡酒 2026-03-08 7/350 2026-03-13 09:56 by mumu真真
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 材料与化工求调剂 +11 与冬清宁 2026-03-07 14/700 2026-03-10 14:23 by Linda Hu
[考研] 化工284求调剂 +12 小黑想变白 2026-03-07 12/600 2026-03-10 09:17 by Wahxp
信息提示
请填处理意见