Форум программистов [Powered by Invision Power Board]

Код


#include <iostream.h>
#include <conio.h>
#include <stdio.h>
float x[5] = {1, 2, 3, 4, 5}; 
float y[5] = {2, 3, 4, 5, 6}; 

/* Функция, вычисляющая коэффициенты Лагранжа
x - аргумент функции
n - степень многочлена (или число x-ов)
i - номер узла
*/

float L(float xp,int n,int i) {
// числитель и знаменатель 
float Chesl; 
float Znam;

Chesl = 1; Znam = 1;

int k;
// вычисление числителя
for (k = 0; k!= n; k++ ) {

if ( k == i ) continue; 
// убираем множитель x - x(i)
Chesl *= xp - x[k];
}
// вычисление знаменателя
for(k= 0; k!= n;k++) {

if (x[i] == x[k]) continue; 
// убираем, а то ноль в знаменателе

Znam *= x[i] - x[k];
}

return Chesl/Znam; 

}
int main(int argc, char* argv[])
{

// вычисляем степень полинома
int n = sizeof(y)/sizeof(float); 

// начальное значение 
float R = 0; 

// произвольная точка для проверки 
float px = 1;

// вычисляем значение интерполяционного многочлена

for (int i = 0; i != n; i++) { 

R += y[i]*L(px,n,i);
}

printf("Result : %f ",R);
    getch();

// вывод результата
return 0;
}

Код


#include <stdio.h>
#include <conio.h>

main(void){
 const  int n = 5;
 double Xi[n] = {1,2,3,4,5};
 double Yi[n] = {2,3,4,5,6};
 double X = 0.1;
 double f, LN, XXX, XX=1.;
 int  i, j, k;

 for (i=1, LN=Yi[0]; i<n; i++)
 {
  XX *= (X-Xi[i-1]);
  for (j=0,f=0; j<=i; j++)
  {
   for (k=0,XXX=1.; k<=i; k++)
   {
    if (k!=j)
     XXX *= Xi[j]-Xi[k];
   }
   f += Yi[j]/XXX;
  }
  LN += XX * f;
 }
 printf("Y[X] = %g",LN);
 getch();
 
}

Код


double Lagrange(double x){
    double F,res=0;
    for (int i=0;i<n;i++){
        F=1;
        for (int j=0;j<n;j++){
            if (i!=j) F*=(x-x[j])/(x[i]-x[j]);
        }
        res+=F*y[i];
    }
    return res;
}

Код


double Newton_f(double x){
    double res=y[0],F,den;
    int i,j,k;
    for(i=1;i<n;i++){
        F=0;
        for(j=0;j<=i;j++){//следующее слагаемое полинома
            den=1;
//считаем знаменатель разделенной разности
            for(k=0;k<=i;k++){
                if (k!=j) den*=(x[j]-x[k]);
            }
//считаем разделенную разность
            F+=y[j]/den;
        }
//домножаем разделенную разность на скобки (x-x[0])...(x-x[i-1])
        for(k=0;k<i;k++)F*=(x-x[k]);
        res+=F;//полином
    }
return res;
}

Цитата(igolka6662 @ 28.11.2010, 11:09)

А как вывести на экран сам полином?