24小时热门版块排行榜    

查看: 455  |  回复: 1

btx97

金虫 (小有名气)

[求助] 求教多项式的稳定性

请教各位,有没有碰到过在数值计算中由于多项式导致最后结果误差比较大的情况。 你们是怎么解决的,能否告诉下。
我的原问题是一个多项式分式f(x)/g(x)的求导, 多项式f(x)的系数要比多项式g(x)的系数大很多, x存在误差。 我发现在个别x点处结果误差比较大,希望哪位高手帮忙指导下。

PS: 我在fortran里,用(系数数列)*(各项x数列)的方法算多项式的值。

[ Last edited by btx97 on 2012-9-25 at 10:07 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)

在数值方面的一本书上找到这么一段话。
We assume that you know enough never to evaluate a polynomial this way:
p=c[0]+c[1]*x+c[2]*x*x+c[3]*x*x*x+c[4]*x*x*x*x;
or (even worse!),
p=c[0]+c[1]*x+c[2]*pow(x,2.0)+c[3]*pow(x,3.0)+c[4]*pow(x,4.0);
Come the (computer) revolution, all persons found guilty of such criminal be-
havior will be summarily executed, and their programs won’t be!
鄙视一下自己。
2楼2012-09-25 17:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 btx97 的主题更新
信息提示
请填处理意见