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

Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция) - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Табуляция функций http://www.cyberforum.ru/cpp-beginners/thread595651.html
Составить программу, которая распечатывает заданное количество значений функций sin(x)*sqrt(x)и соответствующие значение аргумента.Таблица располагается вертикально - сверху вниз. Значение функции распечатываются с точностью 3 знака после десятичной точки, значение аргумента - с точностью 2 знака после десятичной точки. Функция табулируется на интервале двух периодов. Добавлено через 24 минуты...
C++ В массиве хранится информация о росте 35 человек.Определить сколько человек имеют самый большой рост В массиве хранится информация о росте 35 человек.Определить сколько человек имеют самый большой рост http://www.cyberforum.ru/cpp-beginners/thread595647.html
Дан массив все его элементы:увеличить в 2 раза C++
Дан массив все его элементы:увеличить в 2 раза
C++ Каждый элемент массива, больший 10, заменить на его квадратный корень
Дан массив вещественных чисел.Каждый элемент,больший 10,заменить на его квадратный корень.
C++ Найти наибольшее общее кратное двух заданных натуральных чисел http://www.cyberforum.ru/cpp-beginners/thread595636.html
Найти наибольшее общее кратное двух заданных натуральных чисел.
C++ Дано натуральное число.Определить какая цифра встречается в нем чаще 0 или 9 Дано натуральное число.Определить какая цифра встречается в нем чаще 0 или 9 подробнее

Показать сообщение отдельно
xDEAD
13 / 13 / 2
Регистрация: 10.02.2010
Сообщений: 28
04.06.2012, 08:57     Интерполяция заданной функции второй формулой Ньютона (обратная интерполяция)
Программа накапливает ошибку, да и первые значения совсем уже "приближённые". Явно что-то упустил. Посмотрите со стороны.
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
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
 
double *v;
 
double *Create_Vector(int N);
 
void node(int N, int a, int b);
 
double f(double x);
 
double Polinom(int N, double x);
 
double F3(int i, int n);
 
void Print(int N);
 
void clean(int N);
 
int main()
{
    setlocale(LC_ALL,"RUSSIAN");
    int a,b,N;
    double x,S=0;
    cout << "Функция: y(x) = pow(2,x)" << "\n";
    cout << "Введите границы интервала" << "\n" << "a: "; cin >> a; cout << "b: "; cin >> b;
    cout << "Введите количество узлов: "; cin >> N; cout << "\n";
    v=Create_Vector(N);
    node(N,a,b);
    Print(N);
    cout<<"Введите x: "; cin >> x;
    cout<<"Значение f(x): " << f(x) << "\n";
    cout<<"Значение P(x): " << Polinom(N,x) << "\n";
    clean(N);
    system("pause");
    return 0;
}
 
double *Create_Vector(int N) {
        return (double *)calloc(N,sizeof(double));
}
 
void node(int N, int a, int b) {
    double c=(double)(b-a)/(N-1);
    for(int i=0;i<N;i++)
    {
        v[i]=a+c*i;
    }
}
 
double Polinom(int n, double X) {
    double t=(X-v[n-1])/(v[1]-v[0]);
    double sum=v[n-1];
    double f=1;
    double temp=t;
    for(int i=1; i<n; i++) {
        sum=sum+((temp)/f*i)*F3(i, n-i);
        f=f*i;
        temp= temp*(t+i);
    }
    return sum;
}
 
double F3(int i, int n) {
    return pow(f(v[n])-f(v[n-1]), i);
}
 
double f(double x) {
    return pow(2,x);
}
 
void Print(int N)
{
    cout << "x[i]           ";
    cout << "f(x[i])            ";
    cout << "P(x[i])            \n\n";
    for(int i=0;i<N;i++)
    {
        cout << "x[" << i << "] = " << v[i] << "     ";
        cout << "f(x[" << i << "]) = " << f(v[i]);
        cout << "     P(x[" << i << "]) = " << Polinom(N, v[i]) << "\n";         
    }
    cout << "\n";
}
 
void clean(int N)
{
    free(v);
    v=NULL;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 07:06. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru