| 查看: 464 | 回复: 1 | ||
[求助]
求教多项式的稳定性
|
|
请教各位，有没有碰到过在数值计算中由于多项式导致最后结果误差比较大的情况。 你们是怎么解决的，能否告诉下。 我的原问题是一个多项式分式f(x)/g(x)的求导, 多项式f(x)的系数要比多项式g(x)的系数大很多， x存在误差。 我发现在个别x点处结果误差比较大，希望哪位高手帮忙指导下。 PS: 我在fortran里，用（系数数列）*（各项x数列）的方法算多项式的值。 [ Last edited by btx97 on 2012-9-25 at 10:07 ] |
回复此楼
|
在数值方面的一本书上找到这么一段话。 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
