24小时热门版块排行榜    

查看: 1832  |  回复: 6

tmw1987

铁杆木虫 (正式写手)

[交流] 【求助】最小二乘法求参数 已有3人参与

如下式
y=A*exp(-t/B)
上式中,t为时间,y可通过实验测得数据,利用【t,y】系列数值及最小二乘法求参数A与B。我是用matlab中lsqcurefit函数拟合求解。将实验数值y分成n段求解,得到每段中A的数值差距不大,而B的波动是非常大的。

    试图寻找其中的原因和解决办法,我初步考虑可能是参数B处指数函数分母中,当指数函数数值波动一点时,相应的B值会产生很大的波动,不知这种理解是否恰当?

    所以我尝试着先将参数A确定(取平均),然后再将【A,t,y】反带回上式求解参数B,B的波动小了很多。这种处理方法可取吗?
   
    谢谢诸位大侠!

[ Last edited by javeey on 2010-4-21 at 19:10 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

★ ★ ★
小木虫(金币+0.5):恭喜抢沙发,给个红包
tmw1987(金币+10):可否帮忙进一步进行解答!!! 2010-04-21 20:39
Doctorcbw(金币+2):谢谢 2010-06-17 19:07:00
这个问题本身很简单,推导一下就可以求出A,B的表达式。令Y=ln(y),a=ln(A),b=-1/B,则有Y=a+bt,标准的最小二乘,自己套公式或者简单推导一下也可以。

关于你的“参数A确定(取平均),然后再将【A,t,y】反带回上式求解参数B,B的波动小了很多”我只能说可能是巧合。。
板凳要做十年冷文章不发一个字
2楼2010-04-21 20:20:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tmw1987

铁杆木虫 (正式写手)

引用回帖:
Originally posted by bluesine at 2010-04-21 20:20:39:
这个问题本身很简单,推导一下就可以求出A,B的表达式。令Y=ln(y),a=ln(A),b=-1/B,则有Y=a+bt,标准的最小二乘,自己套公式或者简单推导一下也可以。

关于你的“参数A确定(取平均),然后再将【A,t,y】反带回 ...

谢谢版主的亲情解答!

上贴是我的问题的简单化,我实际中用到的公式如下图所示,其中k,a为参数,h为常数,T是因变量,t为自变量。利用最小二乘拟合是k的波动小,而a的波动离奇,可否帮忙诊断一下!

3楼2010-04-21 20:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tmw1987

铁杆木虫 (正式写手)

引用回帖:
Originally posted by bluesine at 2010-04-21 20:20:39:
这个问题本身很简单,推导一下就可以求出A,B的表达式。令Y=ln(y),a=ln(A),b=-1/B,则有Y=a+bt,标准的最小二乘,自己套公式或者简单推导一下也可以。

关于你的“参数A确定(取平均),然后再将【A,t,y】反带回 ...

4楼2010-04-21 20:38:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wlxydgq

铁虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
Doctorcbw(金币+1):谢谢 2010-06-17 09:06:11
我说一下我的看法和经历,并以不一定是正确的,但是楼主可以试试看。
我原来也用matlab拟合过一个函数,和你的一样,同样的程序,每次拟合出来的参数总是不一样的,而且变化还挺大,这让我很是郁闷,怎么都找不到解决的办法。后来,没有办法,我就改用了别的软件。后来无意中发现,如果精度提高的话,就会解决这个问题了。
5楼2010-06-17 08:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by tmw1987 at 2010-04-21 12:59:53:
如下式
y=A*exp(-t/B)
上式中,t为时间,y可通过实验测得数据,利用【t,y】系列数值及最小二乘法求参数A与B。我是用matlab中lsqcurefit函数拟合求解。将实验数值y分成n段求解,得到每段中A的数值差距不大,而 ...

很好的想法,有点预测校正的意思。
6楼2010-06-17 16:17:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tmw1987

铁杆木虫 (正式写手)

引用回帖:
Originally posted by wlxydgq at 2010-06-17 08:54:51:
我说一下我的看法和经历,并以不一定是正确的,但是楼主可以试试看。
我原来也用matlab拟合过一个函数,和你的一样,同样的程序,每次拟合出来的参数总是不一样的,而且变化还挺大,这让我很是郁闷,怎么都找不到 ...

敢问大侠用的是什么软件吗?
7楼2010-06-17 18:55:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tmw1987 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-16 01:21 by Xttdmn
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 求调剂,药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 328,0703考生求调剂,一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
信息提示
请填处理意见