interpolacja wielomianowa

Question

Czesc mam prosty kod ale gdzies jest blad.

#include <iostream>
#include <math.h>
using namespace std;
/*
x0=1.84,	y0=1.78
x1=4.25,	y1=-2.72
x2=8.74,	y2=1.72
*/
int main()
{	
	const int n = 3;
	float a,b,c;
	float x[n]={1.84,4.25,8.74};
	float y[n]={1.78,-2.72,1.72};
	
	a = y[0]/((x[0]-x[1])*(x[0]-x[2])) + y[1]/((x[1]-x[0])*(x[1]-x[2])) + y[2]/((x[2]-x[0])*(x[2]-x[1]));
	b = -y[0]*(x[1]+x[2])/((x[0]-x[1])*(x[0]-x[2]))-y[1]*(x[0]+x[2])/((x[1]-x[0])*(x[1]-x[2]))-y[2]*(x[0]+x[1])/((x[2]-x[0])*(x[2]-x[1]));
	c = y[0]*x[1]*x[2]/((x[0]-x[1])*(x[0]-x[2]))+y[1]*x[0]*x[2]/((x[1]-x[0])*(x[1]-x[2]))+y[2]*x[0]*x[1]/((x[2]-x[0])*(x[2]-x[1]));
	
	cout << "a=" << a << "\n";
	cout << "b=" << b << "\n";
	cout << "c=" << c << "\n";
	
	
	return 0;
}

a i b liczy dobrze ale c podaje zly wynik. Ktos widzi blad?

Prawidlowe c=-0.500437,

pokazuje c= 8.45258

 

Comments

Tak dla jasnosci tu wrzucam algorytm 

Answer 1

A wiec tak:

1. Laptop (kod w C++): c = 8,452579915

2. Kalkulator online: c = c = 8,452579915

3. Kalkulator graficzny casio fx-9860GII: c = c = 8,452579915

4. zwykly kalkulator: c = 8,452579915

5. recznie: c = 8,452579915

Albo ja i te wszytskie urzadzenia sie myla, albo algorytm tzn wzor na c jest bledny. Poklocilbym sie, ze ktorys ze znakow w tym wzorze jest bledny. Przy zastosowaniu tego wzoru:

c = y[0]*x[1]*x[2]/((x[0]-x[1])*(x[0]-x[2])) -
            y[1]*x[0]*x[2]/((x[1]-x[0])*(x[1]-x[2])) -
            y[2]*x[0]*x[1]/((x[2]-x[0])*(x[2]-x[1]));

wszytsko sie zgadza, tzn c=-0.500437, wiec jezeli masz wynik od wykladowcy czy kogos tam, to zglosc mu blad we wzorze

Comments

Spróbuj może użyć typ double. float może mieć zbyt małą precyzję.

---

Nie, tu ei chodzi o precyzje, a o sam fakt, ze wyniki mega sie roznica. Czy uzyje double czy float, to w tym momencie nie gra roli, bo formluka jest bledna, skoro sa az tak rozne wyniki i dopiero po jej zmianie dostaje sie prawidlowy. Watpie, zeby laptop + kalkulator + obliczenia reczne (wyniki te same) roznily sie przez precyzje.