Двовимірні інтерполяційні многочлени та ланцюгові дроби
Аналіз подвійної різниці для функції двох змінних. Інтерполяційний многочлен у формі Ньютона для функції двох змінних та інтерполяційний многочлен Лагранжа у даному випадку. Двовимірні інтерполяційні ланцюгові дроби та їх обчислення в різних випадках.
| Рубрика | Математика |
| Вид | курсовая работа |
| Язык | украинский |
| Дата добавления | 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