Interpolyatsion kubatur formulalar.
Integral ostidagi funksiyani 2 o’lchovli interpolyatsion ko’phad bilan almashtiramiz
Agar
Li x , y
ko’phadlarni quyidagicha
Li x j , y j
aniqlab olsak, u holda
1,
0,
agar i j ,
agar i j ,
i, j
1, N
N
L x, y
i 1
f xi , yi
Li x , y
( 2.25)
ko’phad x j , y j
nuqtada
f x j , y j
qiymatni qabul qiladi. Integral ostidagi
funksiyani (2.25) bilan almashtiramiz:
f
N
x , y dxdy L x , y dxdy Ai
f xi , yi ,
Bu yerda
Ai Li x , y dxdy
i 1
bo’lib, uni murakkab bo’lmagan sohalar uchun hisoblash qiyin emas.
Faraz qilaylik, soha to’g’ri to’rtburchak bo’lsin: a x b ,
Integrallash to’ri sifatida
c y d .
b a d c
xi a ih, y j c jk i 0, m ; j 0, n , h , k
m n
to’g’ri chiziqlarning kesishishlaridan hosil bo’lgan nuqtalar to’plamini olamiz, u holda quyidagi interpolyatsion formulaga ega bo’lamiz:
f x , y f
xi , y j
t s .
i 0
j 0
t 0
xt
s 0
y j y s
t i s j
Buni to’g;ri to’rtburchak bo’ylab integrallasak,
b d m n
f
x, y dxdy
Aij f
xi , y j
a c
hosil bo’ladi, bu yerda
i 0
b m x x
j 0
d n
y y
A
t dx
s dy
ij x x y y
a t 0 i t c s 0 j s
yoki
Aij
b a d
t i s j
c I i , m 1 I j , n 1
o’rinishda yozish mumkin, formulasining koeffitsiyentlaridir.
I i , m 1
va I j , n 1
lar esa Nyuton-Kotes
Kubatur formulanini Simpson formulasi orqali hisoblash
f ( x y )
1
( x y ) 2
a 4
b 5
c 0
d 1
5 1
x0 a
x1
( a b ) 2
x2 b
Z
4
f ( x y ) dx dy
0
y 0 c
y 1
( c d ) 2
y 2 d
Z
Z dastur chiqargan natija
T Simpson formulasining kubatur formulada chiqargan natijasi
R1 xatoligi
h
I
( b a) 2
h k 9
k
( d c) 2
P f x0 y 0 f x2 y 0 f x0 y 2 f x2 y 2 4 f x1 y 0 f x0 y 1 f x2 y 1 f x1 y 2 16 f x1 y 1
T I P
R1 Z T T
R1
bob bo’yicha qisqacha xulosa
Dissertasiyani ikkinchi bobida interpolyatsion kubatur formulalar kurib chiqilgan, ularni xatoliklari tahlil qilinib, effektivligi ko’rib chiqilgan va interpolyatsion kubatur formular uchun algoritm va dastur tuzilib misollarda qo’llanilgan.
BOB. SOBOLEV FAZOSIDA DAVRIY BO'LMAGAN FUNKSIYALAR UCHUN OPTIMAL INTERPOLYATSION FORMULALAR QURISH. 3.1. Sobolev fazosi va unda kubatur formulalar
Oxirgi vaqtlarda integrallarni taqribiy hisoblash uchun kubatur formulalar qurishda sfera sirtida, sferik garmonikalar nazariyasidan foydalaniladi. Bunda integral ostidagi funksiya sferik garmonikalar bo’yicha qatorlarga yoyiladi. Bu bobda esa asosan oddiy va vaznli kubatur formulalarni qarab chiqamiz. Quyidagi oddiy kubatur formulani qarab chiqamiz
f d
S
N
C f ( ( ) ) ,
1
(3.1)
2
L(m ) ( S ) fazoda va uning xatolik funksionalini normasini hisoblaymiz.
2
L(m ) ( S ) fazo berilgan bo’lsin
x
| x |
(
1 , 2
,...,
n 1 ) .
D f ( )
| | 1
(1 ,...., n 1 ) S , ( 1 , 2 , , n 1 ),
| | 1 2 .... n 1 è ! 1 ! n 1 !;
2
sferada berilgan va m tartibli umumlashgan hosilalari kvadrati bilan jamlanuvchi funksiyalar fazosi, funksiya normasi quyidagicha kiritilgan.
k ,
a
( n ,k )
|| f
| Lm ( S ) ||2
2 k m k n 2 m
(3.2)
2
k 1 1
2
L(m ) ( S ) fazoda aslida Sobolev ta’rifi bo’yicha norma quyidagicha kiritilgan.
f L2
S
S m
D f
!
ds
, (3.3)
Shubday qilib
L2 ( S )
- S birlik sferada berilgan va kvadrati bilan jamlanuvchi
funksiyalar fazosini bildiradi.
Agar funksiya
L2 ( S )
fazoga tegishli bo'lsa
f L2 ( S ),
unda uni
ortonormallangan sferik garmonikalar bo’yicha qatorga yoyish mumkim [1]
( n ,k )
f
ak , Yk ,
Yk ,
ak ,
f Yk ,
S
d , k – tartibli ko’rinishdagi ortonormallangan sferik
garmonikalar ;
n , k 2 k n 2 k n 3 !
n 2 !k !
- k – tartibli chiziqli bog’liq bo’lmagan sferik
2
Lm S
bo’lishini [1]:
fazo aniqlaydiki,
f L2
( S ),
funksiyalar uchun yarim norma
k ,
a
( n ,k )
|| f
| Lm ( S ) ||2
2 k m k n 2 m .
2
k 1 1
2
Lm ( S )
fazo elementlarining tuzilishi bo’yicha
L( m ) S
S.L.Sobolev fazosi
2
bilan ustma ust tushadi [2].
(3.1) ko’rinishdagi kubatur formulaning xatolik funksionali quyidagicha bo’ladi.
N
( ) C
1
(
( ) ) , (3.4)
bu yerda - Dirakning delta funksiyasi, s ( )
- S sohaning xarakteristik funksiyasi
( S )
1, ( x )
0,
x S x S
< (x), f(x)>=f(0)
Kubatur formulaning xatoligi deganda biz quyidagi ayirmani tushunamiz:
N , f
S
f ( )d
N
1
C f ( ( ) )
R n
N ( ) f ( )d ,
N
( ) C
1
(
( ) )
( S )
1, ( x )
0,
x S
x S , < (x), f(x)>=f(0)
2
funksionalni tashkil qiladi. Bundan
2 m n ,
L(m ) (S ) C (S )
kelib chiqadi [1].
2
Quyidagi teoremani qarab chiqamiz. Bu teorema G’.N. Salixov tomonidan kiritilgan va isbotlangan lekin isboti qisqa keltirilgan. Biz uning isbotini to’liq keltiramiz sababi keyingi olinadigan asosiy natijalarda ham biz undan foydalanamiz.
Do'stlaringiz bilan baham: |