Форум программистов, компьютерный форум CyberForum.ru

интерполяция Лагранжа! - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 13, средняя оценка - 4.69
DimWhite
0 / 0 / 0
Регистрация: 26.03.2010
Сообщений: 14
16.05.2011, 20:27     интерполяция Лагранжа! #1
в n точках значение интерполяционного многочлена должно совпадать со значением функции, а в промежуточных точках должно быть приблизительно равно. в чем может быть ошибка?, для отрезка [0,1] n=8 в промежуточных точках очень разные значения получаем для функции и интерп-го многочлена Лагранжа.
код:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
#include <iostream>
#include <math.h>
#include <iomanip>
#include <conio.h>
using namespace std;
 
int n;
double *x= new double [n+1];
double *y= new double [n+1];
 
double func(double x)
{
double s=0;
double a=1.0/10;
    for(int k=10;k>=1;k--)
        {
        s+=a;
        a*=-(k*(19-k+1)*4*x*x)/((20-2*k+1)*(20-2*k+2));
        }
return 10*s;
}
double Ln(double arg)
{
    double l=0.0, pr;
    for(int i=0;i<=n;i++)
        {pr=1.0;
        for(int k=0;k<=n;k++)
        if (k!=i) pr*=(arg-x[k])/(x[i]-x[k]);
        l+=y[i]*pr;}
return l;
}
 
int main()
{
double a,b,h;
cout<<"a=";cin>>a;
cout<<"b=";cin>>b;
cout<<"n=";cin>>n;
 
h=(b-a)/n;
for(int i=0;i<=n;i++)
    {
    x[i]=a;
    y[i]=func(a);
    a+=h;
    }
cout<<setprecision(10);
cout<<setw(15)<<left<<"x"<<setw(20)<<"f(x)"<<setw(20)<<"Ln(x)"<<endl;
h*=0.5;
for(int i=0;i<n-1;i++)
{
    cout<<setw(15)<<left<<x[i]<<setw(20)<<y[i]<<setw(20)<<Ln(x[i])<<" - Ln(x["<<i<<"])"<<endl;
    cout<<setw(15)<<left<<x[i]+h<<setw(20)<<func(x[i]+h)<<setw(20)<<Ln(x[i]+h)<<endl;
}
getch();    
return 0;
}
 Комментарий модератора 
Используйте теги форматирования кода!
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
16.05.2011, 20:27     интерполяция Лагранжа!
Посмотрите здесь:

C++ Интерполяция Лагранжа
Обратная интерполяция Лагранжа C++
Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция) C++
Интерполяция полиномом Лагранжа C++
Оцените программу: Интерполяция через полином Лагранжа и через сокращённую формулу Ньютона C++
C++ Интерполяция Лагранжа
Интерполяция полиномом Лагранжа C++
C++ Интерполяция методом Ньютона и Лагранжа

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 02:29. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru