Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 13, средняя оценка - 4.69
DimWhite
0 / 0 / 0
Регистрация: 26.03.2010
Сообщений: 14
#1

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

16.05.2011, 20:27. Просмотров 1630. Ответов 0
Метки нет (Все метки)

в 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;
}
 Комментарий модератора 
Используйте теги форматирования кода!
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
16.05.2011, 20:27
Здравствуйте! Я подобрал для вас темы с ответами на вопрос интерполяция Лагранжа! (C++):

Интерполяция Лагранжа - C++
Доброго времени суток! У меня такая просьба, можно ли упростить данную программу? Т. е. сделать меньше строк в программе, убрать...

Интерполяция Лагранжа - C++
Добрый день! Помогите пожалуйста найти ошибку в программе, битый день делаю, я новичок, но как мне кажется очень логично написал код....

Интерполяция методом Лагранжа - C++
Не могу исправить ошибку в программе вычисляющей интерполяцию методом Лагранжа ошибка в 25,26 строках #include &lt;iostream&gt; #include...

Интерполяция полиномом Лагранжа - C++
Добрый день Я все пытаюсь разобраться с полиномом Лагранжа и столкнулся с очень интересным вопросом, по крайней мере, для себя: ...

Интерполяция полиномом Лагранжа - C++
Функция у = f(x) задана в табличной форме: x=0; y=1.763; x=0.2; y=1.917; x=0.4; y=2.143; а) cоставить алгоритм и программу...

Обратная интерполяция Лагранжа - C++
Помогите пожалуйста нужно найти все значения х по заданному у=1.5 float x = {1,2,3,4,5}; float y = {1.35,-4,0.97,-2,7.2};

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
16.05.2011, 20:27
Привет! Вот еще темы с ответами:

Интерполяция методом Ньютона и Лагранжа - C++
Задано табличное представление функции: x y 3 4 7 10 11 22 15 26 19 23 Сравнить методы интерполяции Лагранжа и Ньютона...

Оцените программу: Интерполяция через полином Лагранжа и через сокращённую формулу Ньютона - C++
Выкладываю программу интерполяции: -интерполирование через полином Лагранжа; -интерполирование через сокращённую формулу Ньютона. ...

Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция) - C++
Программа накапливает ошибку, да и первые значения совсем уже &quot;приближённые&quot;. Явно что-то упустил. Посмотрите со стороны. #include...

Интерполяция изображения c++ - C++
Как сделать интерполяцию изображения если есть например диапазон от 30 до 110. Его нужно растянуть на диапазон от 0 до 255. Допустим берем...


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

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru