24小时热门版块排行榜    

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

agent99

禁虫 (正式写手)

物化程序猿

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

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

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

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +17 梁富贵险中求 2026-04-04 19/950 2026-04-11 00:29 by wangjihu
[考研] 人工智能320调剂08工类还有机会吗 +11 振—TZ 2026-04-10 11/550 2026-04-10 21:51 by blankyung
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 293求调剂 +4 勇远库爱314 2026-04-06 4/200 2026-04-10 11:00 by 高维春
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +8 无情的小羊 2026-04-09 9/450 2026-04-10 09:30 by 松花缸1201
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 280求调剂 +5 兮兮夜夜 2026-04-09 8/400 2026-04-09 11:15 by 兮兮夜夜
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 化学0703-一志愿211-338分求调剂 +10 vants 2026-04-05 11/550 2026-04-08 16:02 by screening
[考研] 一志愿哈工大,初试329,求环境科学与工程调剂! +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 材料调剂 +13 汉123456 2026-04-07 14/700 2026-04-07 22:53 by 来看流星雨10
[考研] 22408 一志愿双一流人工智能300分 四六级,数据分析国奖 +4 zzfeng123 2026-04-06 6/300 2026-04-07 21:02 by zzfeng123
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
信息提示
请填处理意见