24小时热门版块排行榜    

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

agent99

禁虫 (正式写手)

物化程序猿

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

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

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

禁虫 (正式写手)

物化程序猿

引用回帖:
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 果然有我 2026-04-10 7/350 2026-04-11 04:29 by 璃茉一定上岸
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +6 Ezra_Zhang 2026-04-09 6/300 2026-04-10 16:23 by 高维春
[考研] 314求调剂 +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[考研] 一志愿西交机械专硕求调剂 +8 求上岸的小王 2026-04-10 8/400 2026-04-10 15:09 by hemengdong
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[考研] 372分材料与化工(085600)英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 考研调剂 +5 美丽的youth_ 2026-04-04 6/300 2026-04-06 06:57 by houyaoxu
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
信息提示
请填处理意见