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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 【求助】最小二乘法求参数
71/1返回列表
查看: 2340  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-04-21 18:59:53
已阅   回复此楼   关注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人) 回复此楼
板凳要做十年冷文章不发一个字
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拟合过一个函数，和你的一样，同样的程序，每次拟合出来的参数总是不一样的，而且变化还挺大，这让我很是郁闷，怎么都找不到解决的办法。后来，没有办法，我就改用了别的软件。后来无意中发现，如果精度提高的话，就会解决这个问题了。
赞一下(2人) 回复此楼
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的数值差距不大，而 ...

很好的想法，有点预测校正的意思。
赞一下(1人) 回复此楼
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 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 23:45 by vZfe6xYu34yj
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:33 by vZfe6xYu34yj
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
关闭
请填处理意见
关闭