24小时热门版块排行榜    

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

agent99

禁虫 (正式写手)

物化程序猿

[求助] 稳定的高阶多项式拟合

我在implement一个新的量子化学方法,里面需要高阶多项式拟合(原作者用的是20阶)。他用的是mathematica,但我们老板让我写成程序,这样可以整合到我们实验室的library里方便调用。因为计算量很大(除了这个拟合以外还有很多大计算量的步骤)所以我用的fortran。现在其他部分没有问题,就是20阶多项式拟合总是不稳定。请问高手如何解决这个问题?定重谢!

我尝试过的解决办法:
1.采用32位(相当于C的long double)精度,但是并没有彻底解决问题,某些时候还是会不稳定。
2.用Gaussian-Hermite法取代拉格朗日方法,对于某些简单的情况有帮助,但有些情况还是不行。
回复此楼
理论与计算化学,化学物理,生物物理
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agent99

禁虫 (正式写手)

物化程序猿

引用回帖:
2楼: Originally posted by chembetsey at 2013-02-26 13:25:14
多项式拟合是个线性问题,可用多元线性回归方法解决,稳定的方法是奇异值分解SVD。

谢谢,我试过SVD的,但是因为我要求的精度很高,SVD貌似做不到。比如那些高次项,两三个数量级的误差都会造成最后结果偏差很大(拟合得到的系数最后要乘在一系列复杂的函数上)
理论与计算化学,化学物理,生物物理
4楼2013-02-28 03:43:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

chembetsey

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-02-28 00:20:13
agent99: 金币+5, 有帮助 2013-02-28 04:04:12
多项式拟合是个线性问题,可用多元线性回归方法解决,稳定的方法是奇异值分解SVD。
2楼2013-02-27 09:25:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-02-28 00:20:28
agent99: 金币+5, 有帮助 2013-02-28 04:04:25
这是个算法问题,与FORTRAN语言本身无关。
3楼2013-02-27 20:12:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

agent99

禁虫 (正式写手)

物化程序猿

引用回帖:
3楼: Originally posted by nono2009 at 2013-02-27 00:12:20
这是个算法问题,与FORTRAN语言本身无关。

请问有什么好的算法吗?我试过SVD,但是因为我要求的精度很高,SVD貌似做不到。比如那些高次项,两三个数量级的误差都会造成最后结果偏差很大(拟合得到的系数最后要乘在一系列复杂的函数上)
理论与计算化学,化学物理,生物物理
5楼2013-02-28 03:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 284求调剂 +11 archer.. 2026-04-10 12/600 2026-04-10 18:57 by HPUCZ
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考研] 本9一志愿2 0854低分专硕286求调剂 +10 芒种111 2026-04-04 10/500 2026-04-10 12:31 by luosha500
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +16 路痴小琪 2026-04-05 16/800 2026-04-10 08:08 by kangsm
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 291分调剂 +3 上岸小莹加油 2026-04-09 4/200 2026-04-09 17:28 by 1753564080
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 一志愿沪9,求生物学调剂,326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
信息提示
请填处理意见