Wielomian interpolacyjny - bezsensowne wyniki

0

Witam. Napisałem (z użyciem biblioteki Eigen) program, który ma wyznaczyć współczynniki wielomianu interpolacyjnego. Problem w tym, że wyniki zwracane przez ten program są co najmniej dziwne. W programie korzystam z macierzy Vandermonde'a oraz metody LU. Czy możliwe jest, że zagadnienie jest aż tak źle uwarunkowane, że jedyny współczynnik, który się zgadza, ma wartość 1, natomiast pozostałe wyniki to 0 lub -0?

Kod programu:

#include <iostream>
#include <cmath>
#include <iomanip>
#include <eigen3/Eigen/Dense>
#include <eigen3/Eigen/LU>

int main()
{
    const int size=5;

    double x[]={-1.23, -1.19, -0.74, 0.11, 2.56};   //węzły
    double fx[]={1.5129, 1.4161, 0.5476, 0.0121, 6.5536};   //wartości w węzłach

    Eigen::MatrixXd X(size, size);
    Eigen::VectorXd a(size);
    Eigen::VectorXd Y(size);
    //Eigen::IOFormat format(50, 0, ", ", "\n", "[", "]");

    for(int i=0; i<size; ++i)
    {
        for(int j=0; j<size; ++j)
        {
            X(i, j)=std::pow(x[i], size-1-j);
        }

        Y(i)=fx[i];
    }
    
    Eigen::PartialPivLU<Eigen::MatrixXd> lu(X);

    a=lu.solve(Y);
    
    for(int k=0; k<size; ++k)
    {
        std::cout<<"a"<<k<<" = "<<a(size-1-k)<<std::endl;
    }
}
0

czemu size-1-j? Dla macierzy Vandermonde’a jest inaczej. Czemu odwróciłeś indeksy?

Offtopic: naucz się nadawać symbolom nazwy, które będą ci pomagać, a nie utrudniać czytanie kodu. x i X to jest prośba o kopa w ... .

0

Macierz wypełniłem zgodnie z tym, co podaje wykład. Jeśli dobrze rozumiem, w pierwszej kolumnie powinna być potęga zerowa, natomiast w ostatniej n-1?

Co do offtopa: dzięki za radę. Faktycznie, nazwy mogą być mylne. Zmienię je, gdy tylko będę miał dostęp do komputera.

1

Dobra już widzę, w czym problem:
Problem polega na tym, że NIE MA PROBLEMU :), jest dobrze!
Zwróć uwagę, że dane wejściowe dokładnie pokrywają się z funkcją f(x) = x^2 :), więc ma wyjść 0 0 1 0 0 (wartość współczynników a_i).
Niedokładności w reprezentacji tych liczb w systemie binarnym, powodują, że nie wychodzą idealne zera.

0

Ok. Dzięki wielkie. Sugerowałem się wynikami z internetowego kalkulatora oraz Mathematici, dlatego nie wpadłem na coś takiego :). Tam pozostałe współczynniki były rzędu 10^-15 i niżej.
Jeszcze z ciekawości wypróbuję metodę z ćwiczeń (wyznaczenie wielomianu, którego wartość w zerze będzie równa a_i oraz wyinterpolowanie jej).

1 użytkowników online, w tym zalogowanych: 0, gości: 1