Witam, dostałem za zadanie zaimplementować interpolację Newtona... Tak jak nie mam problemu z liczeniem tego "na kartce" tak nie mogę znaleźć jakiejś zależności by to zaimplementować. Dodam, że muszę obejść się bez tablicy dwuwymiarowej. Do liczenia ilorazów różnicowych stworzyć tablicę o wielkości n*(n-1)/2? Czy zapisywać to jakoś do zmiennej?
Z góry dziękuję za jakieś nakierowanie mnie, nie oczekuję gotowego kodu.
Wyjaśnienie: http://www.wypcio.republika.pl/Newton.htm albo to https://www.google.pl/url?sa=t&rct=j&q=&esrc=s&source=web&cd=6&cad=rja&uact=8&ved=0ahUKEwjvlu_ntaXXAhVHqKQKHczQASgQFghEMAU&url=http%3A%2F%2Ffluid.itcmp.pwr.wroc.pl%2F~znmp%2Fdydaktyka%2Fmetnum%2Finterpolacja_newton.pdf&usg=AOvVaw2JZAzUMOJKKYfrUN6CP4hp kod: http://ganeshtiwaridotcomdotnp.blogspot.com/2009/12/c-c-code-for-newtons-interpolation.html albo to: http://www.algorytm.org/procedury-numeryczne/metoda-newtona/newton-c.html
Napisałem coś takiego:
C++ #include <iostream>
using namespace std;
void interpol(double *, double *, int, int);
double *liczIloraz(double *y, double *x, int n);
int main()
{
int n;
cout << "Podaj liczbe wezlow: " << endl;
cin >> n;
double arg;
double *x = new double[n];
double *y = new double[n];
for(int i = 0; i < n; i++)
{
cout << "x[" << i << "] = ";
cin >> x[i];
cout << "y[" << i << "] = ";
cin >> y[i];
}
cout << "Podaj punkt: ";
cin >> arg;
double* b = liczIloraz(y, x, n);
interpol(b, x, n, arg);
return 0;
}
double *liczIloraz(double *y, double *x, int n)
{
int i, j;
double z;
for (i = 0; i < n; j++)
{
z = y[0];
for (j = 0 ; j > n - i; j++)
{
y[j] = (y[j] - y[j - 1]) / (x[i + j] - x[j]);
}
y[n - i] = z;
}
return y;
}
void interpol(double *p, double *x, int n, int arg)
{
double w;
double suma = 0;
int j, i;
for (i = n - 1; i >= 0; i--)
{
w = 1;
for (j = 0; j < i; j++)
w*= (arg - x[j]);
w *= p[j];
suma += w;
}
cout << suma << endl;
}
No i program stoi po linijce
cin >> arg;
Nie robi dosłownie nic, a ja nie mam w ogóle pomysłu dlaczego..
Sprawdź dokładnie 37 linijkę, w takiej postaci ta pętla nigdy się nie wykona.
Jak nie wiesz co się dzieje to albo sobie wypisuj jakie wartości przyjmują zmienne, np:
for (i = 0; i < n; j++)
{
cout << "zmienna i: " << i << endl;
z = y[0];
for (j = 0 ; j > n - i; j++)
{
cout << "zmienna j: " << j << endl;
y[j] = (y[j] - y[j - 1]) / (x[i + j] - x[j]);
}
y[n - i] = z;
}
albo używaj debugera do śledzenia zmiennych.
Zarejestruj się i dołącz do największej społeczności programistów w Polsce.
Otrzymaj wsparcie, dziel się wiedzą i rozwijaj swoje umiejętności z najlepszymi.