Метод Гира

Одним из одношаговых методов получим решения y1,y2,y3 задачи Коши (7.2,7.2’) в точках x1,x2,x3. В окрестности узлов x0,x1,x2,x3,x4 искомое решение y(x) приближенно заменим интерполяционным полиномом Ньютона четвертой степени:

y(x)= + y0+y01(x-x0)+y012(x-x0)(x-x1)+ y0123(x-x0)(x-x1)(x-x2)+y01234(x-x0)(x-x1)(x-x2)(x-x3),   (7.25)

где y01, y012, y0123, y01234 - разделенные разности порядков с первого по четвертый.

Левую часть уравнения (7.2), т.е. производную y’(x), приближенно найдем путем дифференцирования по x полинома (7.25):

y’(x)= + - y01+y012(2x-x0-x1)+y0123[3x2-2x(x0+x1+x2)+ x0x1+x0x2+x1x2]+ y01234[4x3-3x2(x0+x1+x2)+2x(x0x1+x0x2+x0x3+x1x2+x1x3+x2x3)- x0x1x2-x0x1x3-x0x2x3-x1x2x3];   (7.26)

Разделенные разности для равноотстоящих узлов выражаются через узловые значения аппроксимируемой функции:

y01 =(y1-y0) / h,  
y012 =(y2-2y1+y0) / (2h2), (7.27)
y0123 =(y3-3y2+3y1- y0) / (6h3),  
y01234 =(y4-4y3+6y2-4y1+y0) / (24h4).  

Полагая в (7.26) x=x4 и учитывая (7.27), получим:

y’(x4)=(3y0-16y1+36y2-48y3+25y4) / (12h). (7.28)

C другой стороны, исходное дифференциальное уравнение (7.2) при x=x4 принимает вид:

y’(x4)=f(x4,y4). (7.29)

Приравнивая правые части (7.28) и (7.29), находим:

y4=[3(4hf(x4,y4)-y0)+16y1-36y2+48y3] / 25. (7.30)

Формула (7.30) представляет собой неявную схему Гира четвертого порядка для решения задачи Коши (7.2,7.2’). Выражение (7.30) есть уравнение относительно y4, для решения которого можно применить метод простых итераций. Начальное приближение к y4 можно получить из следующих соображений. Полагая в выражении (7.26) x=x3, имеем:

y’(x3)=(-y0+6y1-18y2+10y3-y4) / (12h). (7.31)

Приравнивая правые части (7.2) при x=x3 и выражения (7.31), получим так называемую схему прогноза

y4=4hf(x3,y3)+(y0-10y3)/3-2y1+6y2, (7.32)

которую и можно использовать в качестве начального приближения для решения уравнения (7.30).








Дата добавления: 2015-02-16; просмотров: 2344;


Поиск по сайту:

При помощи поиска вы сможете найти нужную вам информацию.

Поделитесь с друзьями:

Если вам перенёс пользу информационный материал, или помог в учебе – поделитесь этим сайтом с друзьями и знакомыми.
helpiks.org - Хелпикс.Орг - 2014-2024 год. Материал сайта представляется для ознакомительного и учебного использования. | Поддержка
Генерация страницы за: 0.008 сек.