| 查看: 2106 | 回复: 6 | |||
agent99禁虫 (正式写手)
物化程序猿
|
[求助]
稳定的高阶多项式拟合
|
|
我在implement一个新的量子化学方法，里面需要高阶多项式拟合（原作者用的是20阶）。他用的是mathematica，但我们老板让我写成程序，这样可以整合到我们实验室的library里方便调用。因为计算量很大（除了这个拟合以外还有很多大计算量的步骤）所以我用的fortran。现在其他部分没有问题，就是20阶多项式拟合总是不稳定。请问高手如何解决这个问题？定重谢！ 我尝试过的解决办法： 1.采用32位（相当于C的long double）精度，但是并没有彻底解决问题，某些时候还是会不稳定。 2.用Gaussian-Hermite法取代拉格朗日方法，对于某些简单的情况有帮助，但有些情况还是不行。 |
回复此楼» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- mathematica拟合多项式
已经有7人回复
- Matlab多项式拟合
已经有6人回复
- Matlab多项式拟合中参数的不确定度及相关系数
已经有3人回复
- 关于多项式拟合的问题
已经有4人回复
- Gleeble做出加工硬化曲线，为何加工硬化率一直大于0呢？感谢高手指点
已经有11人回复
- 数据拟合多项式
已经有6人回复
- 从excel读取数据多项式拟合出错
已经有3人回复
- R软件如何做局部多项式拟合的预测
已经有10人回复
- R软件如何做多项式拟合线
已经有3人回复
- 【致H兄】数列题的万能多项式拟合求法
已经有28人回复
- 求问SPSS软件进行二次多项式的回归拟合问题。
已经有4人回复
- 【求助】怎样用chebyshev多项式拟合XRPD 背景?
已经有6人回复
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
2楼2013-02-27 09:25:14
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555210
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2013-02-27 20:12:20
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
4楼2013-02-28 03:43:46
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
5楼2013-02-28 03:47:10
chembetsey
木虫 (小有名气)
- 应助: 125 (高中生)
- 金币: 3820.1
- 散金: 5
- 红花: 9
- 帖子: 262
- 在线: 281.8小时
- 虫号: 1781694
- 注册: 2012-04-27
- 专业: 理论和计算化学
6楼2013-02-28 07:23:19
agent99
禁虫 (正式写手)
物化程序猿
- 应助: 75 (初中生)
- 贵宾: 0.001
- 金币: 3630.7
- 散金: 70
- 红花: 8
- 沙发: 1
- 帖子: 392
- 在线: 177.5小时
- 虫号: 1247746
- 注册: 2011-03-28
- 性别: GG
- 专业: 理论和计算化学
7楼2013-03-05 03:21:11