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

Интерполяция функции кубическим сплайном - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Составьте программу, демонстрирующую работу с каждым из 2х классов http://www.cyberforum.ru/cpp-beginners/thread1060311.html
Необходимо разработать заголовочный файл, содержащий описание иерархии из 2х классов, включающий функции для ввода-вывода объектов и заданную функцию. Определение функций-членов класса разместить в отельном файле. П.С. составить программу, демонстрирующую работу с каждым из классов, если можно, с объяснением. Базовый класс-сфера, производный класс-конус,функция объем.
C++ Определить, является ли данная последовательность символов правильным изображением целого 8-го числа Вводится последовательность символов, которая завершается символом точка. Определить, является ли данная последовательность символов правильным изображением целого 8-го числа http://www.cyberforum.ru/cpp-beginners/thread1060306.html
C++ Найти сумму первых n элементов последовательности
Дана последовательность чисел, числа n и k. Найти 1) сумму первых n элементов последовательности 2) первый элемент последовательности который по модулю меньше или больше, (в зависимости от того, убывает последовательность или возрастает), заданного числа k. Последовательность: 2,−1,3,4,6,13,23,... Добавлено через 19 минут Каким образом задать эту последовательность?
C++ Составить блок-схему
Выполнить слияние двух рассортированных текстовых файлов в третий упорядоченный файл. Сортировок не использовать. Дополнительных массивов и файлов не использовать.
C++ Составить блок-схему http://www.cyberforum.ru/cpp-beginners/thread1060276.html
Ввести символьную строку. Вывести на экран все символы, повторяющиеся в строке более 1 раза. Новых строк не создавать.
C++ [Проверить программу] Дан файл, компоненты которого являются действительными числами. Найти произведение компонент файлов Дан файл, компоненты которого являются действительными числами. Найти произведение компонент файлов. Постоянно выдает в ответ "1". В чем заключается проблема? Спасибо за ответы. #include <iostream.h> #include <fstream.h> int main() { ifstream f("f.txt"); double v, s = 1; подробнее

Показать сообщение отдельно
katyasambuka
0 / 0 / 0
Регистрация: 29.12.2013
Сообщений: 13
01.01.2014, 19:09  [ТС]     Интерполяция функции кубическим сплайном
Вот методичка,которой пользовался при написании программы:
миниатюрыИнтерполяция функции кубическим сплайном

Вот переделанная программа (пользовался прямой y=x):
прога
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
#include <stdio.h>
#include <math.h>
#include <conio.h>
#include <malloc.h>
#include <vcl.h>
void main (void)
{system("chcp 1251");
system("cls");
clrscr();
float x[200],y[200],a,b,s,x0;
float h;
float m[200];
int i,n=0;
char tmp[10],tmp1[10];
FILE*p;
p=fopen("I:\\DAT21.TXT","r");
i=0;
fscanf(p,"%s%s",&tmp,&tmp1);
  while(!feof(p))
  {
 fscanf(p,"%f%f",&x[i],&y[i]);
/*printf("\n x=%f y=%f",x[i-1],y[i-1]);*/
  i++;
  n++;
  }
  i=1;
h=x[2]-x[1];
m[0]=(4*y[1]-y[2]-3*y[0])/(2*h);
m[n]=(3*y[n]+y[n-2]-4*y[n-1])/(2*h);
for(i=1;i<n;i++)
{
m[i]=(y[i+1]-y[i-1])/(2*h);
}
p=fopen("I:\\proizvodnie.txt","w");
for(i=0;i<n;i++)
{
fprintf(p,"m[%d]=%f\n",i,m[i]);
}
fclose(p);
x0=(x[1]-x[0])/2;
p=fopen("I:\\resultat.txt","w");
fprintf(p,"ïðîìåæóòî÷Г*Г*Гї òî÷êГ*:\tГ§Г*Г*Г·ГҐГ*ГЁГҐ ïîëèГ*îìГ* Гў ýòîé òî÷êå:");
for(i=1;i<n-2;i++)
{
s=(x[i+1]-x0)*(x[i+1]-x0)*(2*(x0-x[i])+h)/(h*h*h)*y[i]+(x0-x[i])*(x0-x[i])*(2*(x[i+1]-x0)+h)/(h*h*h)*y[i+1]+(x[i+1]-x0)*(((x[i+1]-x0)*(x0-x[i]))/(h*h)*m[i])-((x0-x[i])*(x0-x[i])*(x0-x[i-1])/(h*h)*m[i+1]);
printf("x0=%f s=%f \n",x0,s);
fprintf(p,"\nx0=%f\t\ts(x0)=%f",x0,s);
x0+=h;
}
fclose(p);
getchar();
getchar();
}


производные
m[0]=1.000000
m[1]=1.000000
m[2]=1.000000
m[3]=1.000000
m[4]=1.000000
m[5]=1.000000
m[6]=-2.500000
m[7]=-3.000000

интерполирующие значения для промежуточных точек

промежуточная точка: значение полинома в этой точке:
x0=0.500000 s(x0)=0.750000
x0=1.500000 s(x0)=1.750000
x0=2.500000 s(x0)=2.750000
x0=3.500000 s(x0)=3.750000
x0=4.500000 s(x0)=5.187500
 
Текущее время: 01:02. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru