1 кўринишидаги интегрални ҳисоблаш масаласи қўйилаяпти, бу ерда а ва b
Download 0.79 Mb.
|
6-tajriba. Aniq integralni taqribiy hisoblash
Xi= a + ih қадамида аппроксимацияловчи нолинчи тартибли полином, яъни ўзгармас катталик Y = f(xI) билан алмаштирилади (16-расм). Арифметик прогрессияни ташкил этувчи X = {xI}, i = 0,1,2, … , N сонлар кетма-кетлиги учун ўзгармас катталикларни ҳисоблаймиз:
Y0 = f(x0), Y1 = f(x1), Y2 = f(x2), … YI = f(xI), … YN = f(xN). Демак, эгри чизиқли трапеция N-1 та тўртбурчаклар билан алмаштирилади ва интегрални қиймати тўртбурчаклар юзаларининг йиғиндисига тенг бўлади: ry f(x) YS0 S1 S2 S3 Si SN-2 SN-1 X
16–расм. S0 = hY0 , S1= hY1 , S2 = hY2 , … , Si = hYi , … , SN-1 = hYN-1. (4) интегрални қийматини R хатолик билан чап тўртбурчаклар усулида ҳисоблаш алгоритмлари 17-расмда келтирилган. Б Б Б a , b , N a , b , N a , b , N h = (b-a)/N h = (b-a)/N h = (b-a)/N S = f(a) x=a , S = f(x) S = 0 i = 1, N-1 i = 1 , N-1 x = a, b - h , h xi = a + ih x = x + h S = S + f(xi) S = S + f(x) S = S + f(x) S = hS S = hS S = hS S , N S , N S , N Тамом Тамом Тамом 17-расм.
ry Y f(x) SN-1 SN-2 Si S3 S2 S1 S0 X a h h h h h h b 18 – Расм. X0 = b , X1 = b – 1h, X2 = b – 2h, … , XI = b – ih, … , XN = b – Nh = a. Y0 = f(x0), Y1 = f(x1), Y2 = f(x2), … , Yi = f(xi), … , YN-1 = f(xN-1), xi = b-ih, i = 0,1, …, N-1. Si = hf(xi) , xi = b – ih , i = 0,1, … , N-1. (5) xi = b – ih , i = 1,2, … , N-1. Интегрални қийматини R хатолик билан ўнг тўртбурчаклар усулида ҳисоблаш алгоритмлари 19-расмда келтирилган. Б Б Б a , b , N a , b , N a , b , N h = (b-a)/N h = (b-a)/N h = (b-a)/N S = f(b) x=b , S = f(x) S = 0 i = 1, N-1 i = 1 , N-1 x = b , a+h, -h xi = b - ih x = x - h S = S + f(xi) S = S + f(x) S = S + f(x) S = hS S = hS S = hS S , N S , N S , N Тамом Тамом Тамом 19 – Расм. Итегрални қийматини чап (4) ёки ўнг (5) тўртбурчаклар усулида ҳисоблаганда R хатолик нисбатан анча юқори қийматга эга бўлади. R хатоликни қийматини баҳолашни кейинги усулда кўриб чиқамиз. Ўрта қийматли тўрбурчаклар усули. чап ва ўнг тўртбурчаклар усулида ҳосил бўладиган хатоликни камайтириш мақсадида (n ни қийматини ўзгартирмаган ҳолда) интеграл остидаги f(x) функцияни интеграллаш қадамига мос [xi , xi+1] кесмани ўртасидаги қийматга мос келувчи ўзгармас катталик билан алмаштирамиз, яъниУ ҳолда битта тўртбурчакни юзаси қуйидагича ҳисобланади (20-расм) Y f(x) f(xi) Si X a xi xi+1 b 20-расм. Ўрта тўрт бурчакни юзаси қуйидагича ҳисобланади Изланаётган интегрални қиймати қуйидагича ҳисобланадиЎрта қийматли тўртбурчаклар усули билан интегрални қийматини ҳисоблаш алгоритми 21–расмда келтирилган. Б Б Б a , b , N a , b , N a , b , N h = (b-a)/N h = (b-a)/N h = (b-a)/N S = 0 x=a + h/2, S=f(x) S = 0 i = 1, N i = 1 , N x = a+h/2, b, h xi = a+h/2+(i–1)h x = x + h S = S + f(xi) S = S + f(x) S = S + f(x) S = hS S = hS S = hS S , N S , N S , N Тамом Тамом Тамом 21–расм. Ўрта қийматли тўртбурчаклар усулида хатоликни баҳолаш. Умумий хатолик R = rУ + rЯ . rУ - усулни чекланиш хатолиги; rЯ - яхлитлаш хатолиги. Усулни чекланиш хатолиги қуйидагига тенг: Rу = Jа,к - Jт,к , Jа,к - Интегрални аниқ қиймати; Jт,к - Интегрални тақрибий ҳисобланган қиймати. Хатолик Rу ни баҳолаш учун интеграл остидаги функцияни xi нуқта атрофида Тейлор қаторига ёямиз. (8) Интеграл остидаги f(x) функцияни ўрнига функцияни нуқтадаги Тейлор қаторига ёйилмасини қўйиб ҳадлаб интеграллаш натижасида қуйидаги ифода ҳосил қилинади. (9) Юқоридаги ифодани ҳадларини қийматини ҳисоблаш пайтида, ифодадаги нинг жуфт даражали ҳадларининг қиймати нолга айланиб кетади. Чексиз йиғиндини ифодаловчи (9) ифодадаги биринчи ҳад ўрта тўртбурчаклар усулини формуласидир, қолган ҳадлари эса мазкур усулни хатолигини ташкил қилади. Интеграллаш қадамининг сон қийматини жуда кичиклигини ҳисобга олсак усулни хатолигини иккинчи ҳадни қиймати билан чекланса бўлади, яъни (10) Бутун [a , b] оралиқнинг интеграллаш қадамларида [xi , xi + h] ҳосил бўлган чекланиш хатоликларини жамлаб усулнинг умумий хатолигини аниқлаймиз. (11) Ҳосил қилинган ифода, (11) формула, интегрални қийматини ўрта тўртбурчаклар усули билан ҳисоблашдаги чекланиш хатолигини назарий жиҳатдан баҳолаш формуласидир. Бу формула чекланиш хатолигини амалий жиҳатдан ҳисоблаш учун ноқулайдир. Трапеция формуласи: Трапеция формуласининг қолдиқ ҳадини қиймати қуйидагича ҳисобланади: Агар M=max деб олинса . (12) [1] адабиётда келтирилган маълумотлар асосида яхлитлаш хатолигини юкори чегарасини қуйидагича баҳолаш мумкин: (13) бу ерда - барча f(x ) ларнинг ўртача киймати ; - битта арифметик амални яхлитлаш хатолиги; t - вергулдан кейинги ишончли рақамлар сони. Симпсон формуласи: , бу ерда ; c= 1 ; i=1,3,5, ... . ; Симпсон формуласининг қолдиқ ҳадини қиймати қуйидагича ҳисобланади: , (14) бу ерда M=max ; . Демак, Симпсон формуласи аниқлиги бўйича 4 – тартибли бўлиб қолдиқ ҳадида жуда кичкина қийматли коэффициентга ( ) эга. Симпсон формуласи аниқ ечимни беради агар интеграл остидаги функцияни =0 lV – тартибли ҳосиласи нолга тенг бўлса. Акс ҳолда 2 – тартибли усуллар анча юқори аниқликда ҳисоблаши мумкин, Симпсон усулига нисбатан. Масалан, функция учун n=2 да [-1,1] интеграллаш оралиғида т\рапеция усули J=4 га тенг бўлган аниқ натижани беради, айни шу шароитда Симпсон формуласи анча кўп хатоликка йўл қўяди (натижа J= -8/3). Юқоридаги (11), (12), (13), (14) формулаларни таҳлил қилиш натижасида усулни чекланиш хатолиги h2 га, яхлитлаш хатолиги эса 1/h га боғлик, эканлиги кўриниб турибди. Демак, h ни камайиши билан хатолик rу камаяди, rЯ эса ортаборади. Назарий жиҳатдан h ни камайтира бориб интегрални хоҳлаганча аниқликда ҳисоблашга имконият бор дейиш мумкин, лекин амалда бундай бўлмайди, сабаби яхлитлаш хатолиги бунга йўл қўймайди. Шу сабабдан [а;b] оралиқда интегрални ҳисоблашда h ни хоҳлаганча кичик қийматини олиб бўлмайди, чунки ҳар бир интеграл остидаги функция учун оптимал бўлакларга бўлиш қиймати мавжуд, яъни N0 (14-расм). Аниқ интегрални тақрибий ҳисоблашга доир мисоллар Download 0.79 Mb. Do'stlaringiz bilan baham: |
Ma'lumotlar bazasi mualliflik huquqi bilan himoyalangan ©fayllar.org 2024
ma'muriyatiga murojaat qiling
ma'muriyatiga murojaat qiling