Двовимірні інтерполяційні многочлени та ланцюгові дроби

Аналіз подвійної різниці для функції двох змінних. Інтерполяційний многочлен у формі Ньютона для функції двох змінних та інтерполяційний многочлен Лагранжа у даному випадку. Двовимірні інтерполяційні ланцюгові дроби та їх обчислення в різних випадках.

Рубрика Математика
Вид курсовая работа
Язык украинский
Дата добавления 13.11.2017
Размер файла 126,6 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Var Nx,Ny,Cx,Cy:Integer;

X:Array[0..MaxX] Of Real;

Y:Array[0..MaxY] Of Real;

B:MyArr;

Xa,Xb,Ya,Yb:Real;

D1,D2:^MyArr;

cc,cc1:Integer;

Function Func(x,y:Real):Real;

Begin

Func:=1/(x*x+y*y+x*y);

End;

Procedure DataInput;

Var i,j:Integer;

Begin

{ Write('Input Xa : '); ReadLn(Xa);

Write('Input Xb : '); ReadLn(Xb);

Write('Input Ya : '); ReadLn(Ya);

Write('Input Yb : '); ReadLn(Yb);}

Xa:=1; Xb:=2; Ya:=1; Yb:=2;

{ Write('Input Nx : '); ReadLn(Nx);

Write('Input Ny : '); ReadLn(Ny);}

nx:=cc; ny:=cc1*2-1;

{ For i:=0 To Nx Do X[i]:=(Xa+Xb)/2+(Xb-Xa)*Cos(Pi*i/Nx)/2;

For i:=0 To Ny Do Y[i]:=(Ya+Yb)/2+(Yb-Ya)*Cos(Pi*i/Ny)/2;}

For i:=0 To Nx Do X[i]:=Xa+(Xb-Xa)*i/Nx;

For i:=0 To Ny Do Y[i]:=Ya+(Yb-Ya)*i/Ny;

End;

Procedure BuildCoefTable;

Function Xij(i,j:Integer):Real;

Begin

If i>j Then Xij:=X[i]-X[j] Else Xij:=1;

End;

Function Yij(i,j:Integer):Real;

Begin

If i>j Then Yij:=Y[i]-Y[j] Else Yij:=1;

End;

Function Teta(t,s:Integer):Integer;

Begin

If s>t Then Teta:=-1 Else Teta:=0;

End;

Function Delta(k,i,j:Integer):Real;

Begin

Delta:=Xij(i,k)*Yij(j,k)/

( D1^[i,j]+

Teta(k,j)*D1^[i,k]+

Teta(k,i)*D1^[k,j]+

Teta(k,i)*Teta(k,j)*D1^[k,k]

);

End;

Var i,j,s,k,Mx:Integer;

Begin

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

D1^[i,j]:=Func(X[i],Y[j]);

End;

k:=0;

D2^:=D1^;

If Nx>Ny Then Mx:=Nx Else Mx:=Ny;

While k<Mx+1 Do

Begin

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

If i>j Then s:=i Else s:=j;

If s=k Then B[i,j]:=D2^[i,j];

End;

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

D2^[i,j]:=Delta(k,i,j);

End;

D1^:=D2^;

k:=k+1;

End;

End;

Function Drib(xx,yy:Real):Real;

Var n:Integer;

Function GetH(m,k:Integer):Real;

Begin

If m=n+1 Then GetH:=0

Else

Begin

GetH:=(xx-X[m-1])/(B[m,k]+GetH(m+1,k));

End;

End;

Function GetL(m,k:Integer):Real;

Begin

If m=n+1 Then GetL:=0

Else

GetL:=(yy-Y[m-1])/(B[k,m]+GetL(m+1,k));

End;

Function GetG(k:Integer):Real;

Begin

If k=n+1 Then GetG:=0

Else

GetG:=(xx-X[k-1])*(yy-Y[k-1])/

(B[k,k]+GetH(k+1,k)+GetL(k+1,k)+GetG(k+1));

End;

Begin

If Nx<Ny Then n:=Nx Else n:=Ny;

Drib:=B[0,0]+GetH(1,0)+GetL(1,0)+GetG(1);

End;

Function Polinom(xx,yy:Real):Real;

Var p,q,s,s1,p1,q1:Real; i,j,k:Integer;

Begin

s:=0;

For i:=0 To Nx Do

For j:=0 To Ny Do

Begin

p:=1; q:=1;

For k:=0 To Nx Do If k<>i Then p:=p*(xx-X[k])/(X[i]-X[k]);

For k:=0 To Ny Do If k<>j Then q:=q*(yy-Y[k])/(Y[j]-Y[k]);

s1:=p*q*Func(X[i],Y[j]);

s:=s+s1;

End;

Polinom:=s;

End;

Procedure GetMaxError;

Var i,j:Integer; dx,dy,MaxErr1,p1,p2,p3,VidnErr1,MaxErr2,VidnErr2:Real; F:Text;

Begin

MaxErr1:=0; VidnErr1:=0; MaxErr2:=0; VidnErr2:=0;

dx:=(Xb-Xa)/Cx; dy:=(Yb-Ya)/Cy;

For i:=0 To Cx Do

For j:=0 To Cy Do

Begin

p1:=Func(Xa+i*dx,Ya+j*dy);

p2:=Drib(Xa+i*dx,Ya+j*dy);

p3:=Polinom(Xa+i*dx,Ya+j*dy);

If Abs(p1-p3)>MaxErr1 Then

Begin

MaxErr1:=Abs(p1-p3); VidnErr1:=Abs((p1-p3)/p1);

end;

If Abs(p1-p2)>MaxErr2 Then

Begin

MaxErr2:=Abs(p1-p2); VidnErr2:=Abs((p1-p2)/p1);

End;

End;

Assign(f,'mix.txt'); Append(f);

WriteLn(f,nx:4,ny:4,MaxErr2:19:12,VidnErr2:19:12,MaxErr1:19:12,VidnErr1:19:12);

Close(f);

End;

Begin

For cc:=1 To 10 Do For cc1:=1 To 5 Do

Begin

DataInput; cx:=33; cy:=33;

WriteLn('Nx=',nx,' Ny=',ny);

New(D1); New(D2); BuildCoefTable; Dispose(D1); Dispose(D2);

GetMaxError;

End;

WriteLn('Press <ENTER>'); ReadLn;

End.

Размещено на Allbest.ru

...

Подобные документы

  • Сутність інтерполяційних поліномів. Оцінка похибок інтерполяційних формул, їх застосування. Програма обчислення наближених значень функції у випадку, коли функція задана таблично, використовуючи інтерполяційні формули для рівновіддалених вузлів.

    курсовая работа [956,4 K], добавлен 29.04.2011

  • Функція двох змінних, методика визначення її головних параметрів. Поняття екстремуму функцій двох змінних, необхідні та достатні умови її існування. Особливості визначення екстремуму функції за деяких умов, які обмежують область зміни аргументів.

    курсовая работа [1,0 M], добавлен 22.10.2014

  • Побудова дотичної площини та нормалі до поверхні. Геометричний зміст диференціала функції двох змінних. Поняття скалярного поля, зв'язок між градієнтом і похідною в даній точці. Формула Тейлора для функції двох змінних та її локальні екстремуми.

    реферат [713,9 K], добавлен 14.05.2011

  • Будування сіткової функції. Методи прямокутників і трапецій, підвищення їх точності. Інтерполяційний многочлен Лагранжа другого степеня. Формула Сімпсона для чисельного інтегрування. Похибка формули Сімпсона. Обчислення наближеного значення інтеграла.

    презентация [99,6 K], добавлен 06.02.2014

  • Суть функції багатьох змінних, її означення і символіки. Границя і неперервність функції багатьох змінних. Визначення відкритої та замкненої області. Множина точок площини, для яких задана формула має зміст, як область визначення. Функція двох змінних.

    реферат [289,8 K], добавлен 01.05.2011

  • Сутність, особливості та історична поява чисел "пі" та "е". Доведення ірраціональності та трансцендентності чисел "пі" та "е". Методи наближеного обчислення чисел "пі" та "е" за допомогою числових рядів та розкладу в нескінченні ланцюгові дроби.

    курсовая работа [584,5 K], добавлен 18.07.2010

  • Сутність фізичного та геометричного змісту похідної, особливості його використовування у математичних задачах. Означення диференціалу, формула його обчислення. Екстремуми функцій двох змінних. Правила знаходження найбільшого і найменшого значення функції.

    презентация [262,6 K], добавлен 20.05.2015

  • Метод решения задачи, при котором коэффициенты a[i], определяются непосредственным решением системы - метод неопределенных коэффициентов. Интерполяционная формула Ньютона и ее варианты. Построение интерполяционного многочлена Лагранжа по заданной функции.

    лабораторная работа [147,4 K], добавлен 16.11.2015

  • Кінцеві різниці різних порядків. Залежність між кінцевими різницями і функціями. Дискретний і неперервний аналіз. Поняття про розділені різниці. Інтерполяційна формула Ньютона. Порівняння формул Лагранжа і Ньютона. Інтерполяція для рівновіддалених вузлів.

    контрольная работа [75,6 K], добавлен 06.02.2014

  • Загальні формули прямокутників. Похибка методу прямокутників. Площа криволінійної трапеції. Формула парабол (Сімпсона). Інтерполяційний багаточлен Лагранжа. Формула трьох восьмих. Абсолютна похибка обчислення. Наближення підінтегральної функції.

    лабораторная работа [298,1 K], добавлен 26.03.2011

  • Частинні похідні та диференційованість функції: поняття та теореми. Повний диференціал функції та його застосування до обчислення функцій і похибок. Диференціали вищих порядків. Інваріантність форми повного диференціала. Диференціювання неявної функції.

    реферат [278,8 K], добавлен 02.05.2011

  • Определение абсолютной и относительной погрешностей приближенных чисел. Оценка погрешностей результата. Интерполирование и экстраполирование данных, интерполяционный многочлен Лагранжа и Ньютона, их основные характеристики и сравнительное описание.

    лабораторная работа [74,8 K], добавлен 06.08.2013

  • Построить интерполяционный многочлен Ньютона. Начертить график и отметить на нем узлы интерполяции. Построить интерполяционный многочлен Лагранжа. Выполнить интерполяцию сплайнами третьей степени.

    лабораторная работа [70,8 K], добавлен 06.02.2004

  • Вычислительные методы линейной алгебры. Интерполяция функций. Интерполяционный многочлен Ньютона. Узлы интерполяции. Интерполяционный многочлен Лагранжа. Интерполяция сплайнами. Коэффициенты кубических сплайнов.

    лабораторная работа [70,5 K], добавлен 06.02.2004

  • Метод Гаусса, метод прогонки, нелинейное уравнение. Метод вращения Якоби. Интерполяционный многочлен Лагранжа и Ньютона. Метод наименьших квадратов, интерполяция сплайнами. Дифференцирование многочленами, метод Монте-Карло и Рунге-Кутты, краевая задача.

    курсовая работа [4,8 M], добавлен 23.05.2013

  • Поняття диференційованості функції в даній точці, основні формули. Диференціал функції однієї змінної, його застосування. Основні означення, які відносяться до функції кількох змінних. Похідна алгебраїчної суми скінченного числа диференційованих функцій.

    реферат [101,8 K], добавлен 02.11.2015

  • Коротка біографія видатного математика Б. Тейлора. Тейлорова формула із залишковим членом у формі Пеано та у Лагранжовій формі. Розвинення деяких елементарних функцій за формулою Тейлора. Формула Тейлора для многочлена та для функції однієї змінної.

    курсовая работа [547,0 K], добавлен 20.05.2015

  • Первая дробь, с которой познакомились люди в Египте. Числитель и знаменатель дроби. Правильная и неправильная дробь. Смешанное число. Приведение к общему знаменателю. Неполное частное. Целая и дробная часть. Обратные дроби. Умножение и деление дробей.

    презентация [48,9 K], добавлен 11.10.2011

  • Теоретико-методологические основы формирования математического понятия дроби на уроках математики. Процесс формирования математических понятий и методика их введения. Практическое исследование введения и формирования математического понятия дроби.

    дипломная работа [161,3 K], добавлен 23.02.2009

  • Тождества, используемые для системы Жигалкина. Многочлен Жигалкина функции одной, двух и трех переменных. Содержание теоремы. Практический пример преобразования многочлена с помощью метода цепочки и неопределенных коэффициентов. Закон полного поглощения.

    контрольная работа [95,5 K], добавлен 06.08.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.