Рішення систем диференціальних рівнянь за допомогою неявної схеми Адамса 3-го порядку

Курсовой проект - Математика и статистика

Другие курсовые по предмету Математика и статистика

чно не відчутно. Збільшення кроку сітки приведе до підвищення точності рішення, однак це збільшить і час роботи обчислювального процесу.

Задана точність досягається за мінімальну кількість ітерацій (13 ітерації).

Нижче наведений графік функцій отриманого й точного рішень:

 

Рис.5.1 Графік отриманого й точного рішення

 

Рис.5.2 Графік отриманого й точного рішення

Як видно з малюнків 5.1, 5.2, розбіжність кривих спостерігається тільки при досить великому збільшенні графіка.

Запропонована задача Коші була також вирішена в математичному пакеті Mathcad 11 двома методами: методом Рунге-Кутта 5-го порядку й методом Рунге-Кутта з непостійним кроком. Реалізація рішення системи диференціальних рівнянь в Mathcad 11 і таблиці результатів наведені нижче:

Реалізація рішення задачі Коші методом Рунге-Кутта 5-го порядку:

 

 

Таблиця 5.1 Результати рішення задачі Коші методом Рунге-Кутта 5-го порядку.

xu(x)v(x)xu(x)v(x)247,38905613,16,222,197952,024,047,53832493,126,2422,646382,044,087,69060923,146,2823,103872,064,127,84596983,166,3223,57062,084,168,00446893,186,3624,046752,14,28,16616993,26,424,532532,124,248,33113753,226,4425,028122,144,288,49943763,246,4825,533722,164,328,67113763,266,5226,049542,184,368,84630623,286,5626,575772,24,49,02501353,36,627,112642,224,449,20733083,326,6427,660352,244,489,39333133,346,6828,219132,264,529,58308913,366,7228,789192,284,569,77668043,386,7629,370772,34,69,97418243,46,829,96412,324,6410,1756743,426,8430,569412,344,6810,3812373,446,87999931,186962,364,7210,5909513,466,91999931,816982,384,7610,8049033,486,95999932,459722,44,811,0231763,56,99999933,115452,424,8411,2458593,527,03999933,784432,444,8811,4730413,547,07999934,466922,464,9211,7048113,567,11999935,16322,484,9611,9412643,587,15999935,873542,54,999999912,1824943,67,19999936,598232,525,039999912,4285973,627,23999937,337572,545,079999912,6796713,647,27999938,091842,565,119999912,9358173,667,31999938,861342,585,159999913,1971383,687,35999939,646392,65,199999913,4637383,77,39999940,44732,625,239999913,7357233,727,43999941,264392,645,279999914,0132043,747,47999942,097992,665,319999914,2962893,767,51999942,948422,685,359999914,5850933,787,55999943,816042,75,399999914,8797323,87,59999944,701182,725,439999915,1803223,827,63999945,604212,745,479999915,4869853,847,67999946,525472,765,519999915,7998433,867,71999947,465352,785,559999916,1190213,887,75999948,424212,85,599999916,4446473,97,79999949,402452,825,639999916,7768513,927,83999950,400442,845,679999917,1157653,947,87999951,41862,865,719999917,4615273,967,91999952,457322,885,759999917,8142733,987,95999853,517032,95,799999818,17414547,99999854,598152,925,839999818,5412872,945,879999818,9158462,965,919999819,2979722,985,959999819,68781635,999999820,0855373,026,039999820,4912913,046,079999820,9052433,066,119999821,3275573,086,159999821,758402

Реалізація рішення задачі Коші методом Рунге-Кутта з непостійним кроком:

 

 

Таблиця 5.2 Результати рішення задачі Коші методом Рунге-Кутта з непостійним кроком.

Xu(x)v(x)247,3890560992,24,49,0250134862,44,811,023176342,65,213,463737962,85,616,444646633620,085536693,26,424,532529813,46,829,964099443,67,236,598233483,87,644,7011834854,59814775

Як видно з отриманих таблиць результатів, точність рішення в 0.0001 при рішенні методом Рунге-Кутта з непостійним кроком досягається всього за 10 кроків, у той час, коли для досягнення цієї ж точності при рішенні методом Рунге-Кутта 5-го порядку з постійним кроком потрібно близько 100 кроків.

Порівнюючи отримані результати з результатами роботи програми Adams3.exe, доходимо висновку, що неявна схема Адамса третього порядку досить ефективна при чисельному рішенні задачі Коші (швидкість, висока точність рішення), однак по своїх характеристиках вона уступає більше зробленим методам, що застосовуються в різних математичних пакетах.

 

Висновок

 

Результатом виконання курсового проекту є готовий програмний продукт, що дозволяє вирішувати задачу Коші для системи диференціальних рівнянь за допомогою неявної схеми Адамса 3-го порядки, що демонструє можливості чисельного рішення поставленої задачі із заданим ступенем точності.

Готовий програмний продукт може знайти широке застосування при рішенні багатьох прикладних технічних програм, а зокрема, ефективне використання застосованої схеми Адамса 3-го порядки для рішення так званих твердих систем диференціальних рівнянь, для яких існує лише чисельне рішення.

Дана програма вирішує задану користувачем систему диференціальних рівнянь із зазначеною точністю за мінімальний проміжок часу. При цьому користувачеві надається можливість візуально оцінити неточність рішення, порівнюючи графіків отриманого й точного рішень.

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

 

Література

 

1. АрхангельськийО.Я.Програмування в С++ Builder 6. К., 2004

2. Каліткин М.М. Чисельні методи. К., 2003

3. СамарськийА.А., Гулін А.В. Чисельні методи. К., 2003

4. Синіцин О.К., НавроцкийА.А.Алгоритми обчислювальної математики. К., 2003

5. Синіцин О.К.Програмування алгоритмів у середовищі Builder C++. К., 2003

6. Страуструп Бьерн. Язык программирования C++. М., 2002

7. Шилд Г. Программирование на Borland C++ для профессионалов М., 1999.

 

 

Додаток 1

 

1. Блок-схема алгоритму

 

 

2. Блок-схема рішення задачі Коші неявною схемою Адамса 3-го порядки

 

3. Блок-схема алгоритму перетворення рядка у зворотний польський запис

 

4. Блок-схема обчислення функцій