1

我试图生成具有 11 个系数的 10 阶多项式的值。我也在尝试生成它的导数。我写了如下所示的三个函数。此代码生成多项式的值。a1 到 a10 是系数。

float polynm( float a0,float a1,float a2,float a3,float a4,float a5,float a6,float a7,float a8,float a9,float a10,float x)
     {
          float poly = a0 + a1*x + a2*pow(x,2)+a3*pow(x,3)+a4*pow(x,4)+a5*pow(x,5)+a6*pow(x,6)+a7*pow(x,7)+a8*pow(x,8)+a9*pow(x,9)+a10*pow(x,10);
             return poly;
             }

此代码生成它调用函数 deri 的多项式的导数的值

 float polynm_der(float a0,float a1,float a2,float a3,float a4,float a5,float a6,float a7,float a8,float a9,float a10,float x)
    {  float der = a1 + a2*deri(x,2)+a3*deri(x,3)+a4*deri(x,4)+a5*deri(x,5)+a6*deri(x,6)+a7*deri(x,7)+a8*deri(x,8)+a9*deri(x,9)+a10*deri(x,10);
       return der;
       }
deri is below
float deri(float x,int n)
   {   
         float term_der = n*pow(x,n-1);
           return term_der;
           }

多项式的代码效率低下。如果我想生成一个 100 阶多项式，那将变得不可能。有没有办法我可以递归地生成多项式及其导数以避免笨拙的代码。

4

2 回答 2

6

一种解决方案是接受一组系数及其长度：

float poly(int x, float coefficients[], int order)
{
    int idx;
    float total;

    for (idx = 0; idx < order; idx++)
        total += coefficients[idx] * pow(x, idx);
    return total;
}

递归解决方案会很漂亮，但这不是 Lisp。无论如何，类似的方法可以用于衍生品。请记住，在 C 语言中，函数的数组参数会变成指针，因此您不能使用很酷的东西sizeof来获取它们的长度。

编辑：作为对评论的回应，您可以在coefficients构造数组时强制执行您的要求。或者，如果您不负责该代码，则可以将其粘贴在函数中（hackishly），如下所示：

if (coefficients[0] == 0 || coefficients[1] == 0 || coefficients[order-1] == 0)
    assert(0);
于 2011-03-22T17:44:48.773 回答
0

您可以重写函数以获取 x 值、系数数组和长度。

于 2011-03-22T17:45:53.423 回答