Подставим функцию в исходное дифференциальное уравнение .
Перенесем в левую часть:
Очевидно, что выражение в скобках дает 0. Следовательно, найденная функция удовлетворяет ДУЧП.
Проверим, удовлетворяет ли функция начальным и граничным условиям.
;
;
;
найденная функция удовлетворяет всем начальным и граничным условиям.
Теперь подберем такое количество членов ряда , чтобы равенство , удовлетворялось с точностью до 10-5 при . Результат подбора при - :
0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 | ||
0 | 0.07074 | 0.5403 | 0.87758 | 1 | 0.87758 | 0.54030 | 0.07073 | 0 | 0 | 0 | |
0 | 0.07074 | 0.5403 | 0.87758 | 1 | 0.87758 | 0.54030 | 0.07073 | 0 | 0 | 0 | |
0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
График функции График разложения в ряд Фурье при
Достоверность решения доказана.
Вывод конечно-разностных уравнений, формулировка вычислительной схемы метода конечных разностей
При решении уравнений с частными производными методом конечных разностей производные заменяются соответствующими разностями применяем .
(1)
(2)
(3)
Краевая задача для волнового уравнения:
требуется найти решения уравнения , удовлетворяющие краевым условиям:
,
,
,
, .
то есть требуется найти решение в прямоугольнике, ограниченном прямыми , , , , если заданы значения искомой функции на трех его сторонах , , и значение производной при . Покроем нашу область сеткой, образованной прямыми
и будем определять приближенные значения в узлах сетки, т. е. в точках пересечения этих прямых. Введем обозначения: . Напишем вместо волнового уравнения соответствующее ему уравнение в конечных разностях для точки . В соответствии с формулами (1), (2) и (3) получим
В итоге получим систему из уравнений, решения которой являются приближенными значениями функции в узлах сетки разбиения:
Табулирование значений искомой функции в узлах сетки, покрывающей область, в которой найдено решение. Построение графика найденной функции.
Таблица точных значений функции при следующих значениях управляющих параметров: , , , , , , :
t\x | 0.00000 | 1.00000 | 2.00000 | 3.00000 | 4.00000 | 5.00000 | 6.00000 | 7.00000 | 8.00000 | 9.00000 | 10.00000 |
0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 |
1.50000 | 0.00000 | 0.00034 | 0.01433 | 0.15939 | 0.17309 | 0.17311 | 0.17309 | 0.15939 | 0.01433 | 0.00034 | -0.00000 |
3.00000 | 0.00000 | -0.14829 | -0.14839 | -0.14847 | -0.14851 | -0.14853 | -0.14851 | -0.14847 | -0.14839 | -0.14829 | -0.00000 |
4.50000 | 0.00000 | -0.00051 | 0.00908 | 0.11616 | 0.12630 | 0.12631 | 0.12630 | 0.11616 | 0.00908 | -0.00051 | -0.00000 |
6.00000 | 0.00000 | 0.00066 | 0.00132 | 0.00198 | 0.00264 | 0.00311 | 0.00264 | 0.00198 | 0.00132 | 0.00066 | 0.00000 |
7.50000 | 0.00000 | -0.00075 | -0.00899 | -0.08887 | -0.09639 | -0.09640 | -0.09639 | -0.08887 | -0.00899 | -0.00075 | 0.00000 |
9.00000 | 0.00000 | 0.08137 | 0.08142 | 0.08145 | 0.08147 | 0.08148 | 0.08147 | 0.08145 | 0.08142 | 0.08137 | 0.00000 |
10.50000 | 0.00000 | 0.00071 | -0.00413 | -0.06268 | -0.06823 | -0.06824 | -0.06823 | -0.06268 | -0.00413 | 0.00071 | 0.00000 |
12.00000 | 0.00000 | -0.00073 | -0.00145 | -0.00218 | -0.00290 | -0.00342 | -0.00290 | -0.00218 | -0.00145 | -0.00073 | -0.00000 |
13.50000 | 0.00000 | 0.00072 | 0.00554 | 0.04953 | 0.05365 | 0.05366 | 0.05365 | 0.04953 | 0.00554 | 0.00072 | -0.00000 |
15.00000 | 0.00000 | -0.04464 | -0.04466 | -0.04467 | -0.04467 | -0.04467 | -0.04467 | -0.04467 | -0.04466 | -0.04464 | -0.00000 |
16.50000 | 0.00000 | -0.00063 | 0.00179 | 0.03380 | 0.03684 | 0.03684 | 0.03684 | 0.03380 | 0.00179 | -0.00063 | -0.00000 |
18.00000 | 0.00000 | 0.00060 | 0.00119 | 0.00179 | 0.00239 | 0.00281 | 0.00239 | 0.00179 | 0.00119 | 0.00060 | 0.00000 |
19.50000 | 0.00000 | -0.00056 | -0.00336 | -0.02759 | -0.02985 | -0.02985 | -0.02985 | -0.02759 | -0.00336 | -0.00056 | 0.00000 |
21.00000 | 0.00000 | 0.02449 | 0.02449 | 0.02449 | 0.02449 | 0.02448 | 0.02449 | 0.02449 | 0.02449 | 0.02449 | 0.00000 |
22.50000 | 0.00000 | 0.00048 | -0.00072 | -0.01822 | -0.01988 | -0.01988 | -0.01988 | -0.01822 | -0.00072 | 0.00048 | 0.00000 |
24.00000 | 0.00000 | -0.00044 | -0.00087 | -0.00131 | -0.00175 | -0.00206 | -0.00175 | -0.00131 | -0.00087 | -0.00044 | -0.00000 |
25.50000 | 0.00000 | 0.00039 | 0.00202 | 0.01536 | 0.01660 | 0.01660 | 0.01660 | 0.01536 | 0.00202 | 0.00039 | -0.00000 |
27.00000 | 0.00000 | -0.01343 | -0.01342 | -0.01342 | -0.01341 | -0.01341 | -0.01341 | -0.01342 | -0.01342 | -0.01343 | -0.00000 |
28.50000 | 0.00000 | -0.00034 | 0.00025 | 0.00981 | 0.01072 | 0.01072 | 0.01072 | 0.00981 | 0.00025 | -0.00034 | -0.00000 |
30.00000 | 0.00000 | 0.00030 | 0.00060 | 0.00090 | 0.00120 | 0.00141 | 0.00120 | 0.00090 | 0.00060 | 0.00030 | 0.00000 |
Таблица приближенных значений функции при следующих значениях управляющих параметров: , , , , , , :
t\x | 0.00000 | 1.00000 | 2.00000 | 3.00000 | 4.00000 | 5.00000 | 6.00000 | 7.00000 | 8.00000 | 9.00000 | 10.00000 |
0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 | 0.00000 |
1.50000 | 0.00000 | 0.01054 | 0.01644 | 0.18001 | 0.21407 | 0.22213 | 0.21407 | 0.18001 | 0.01644 | 0.01054 | -0.00000 |
3.00000 | 0.00000 | -0.13629 | -0.12855 | -0.11852 | -0.10851 | -0.10681 | -0.10851 | -0.11852 | -0.12855 | -0.13629 | -0.00000 |
4.50000 | 0.00000 | 0.01001 | 0.02808 | 0.14679 | 0.16673 | 0.16671 | 0.16673 | 0.14679 | 0.02808 | 0.01001 | -0.00000 |
6.00000 | 0.00000 | -0.01055 | 0.01252 | 0.03100 | 0.00500 | 0.04365 | 0.00500 | 0.03100 | 0.01252 | -0.01055 | 0.00000 |
7.50000 | 0.00000 | 0.01033 | 0.01113 | -0.05924 | -0.06639 | -0.05740 | -0.06639 | -0.05924 | 0.01113 | 0.01033 | 0.00000 |
9.00000 | 0.00000 | 0.08137 | 0.09282 | 0.11153 | 0.11150 | 0.12187 | 0.11150 | 0.11153 | 0.09282 | 0.08137 | 0.00000 |
10.50000 | 0.00000 | 0.00071 | 0.01406 | -0.03270 | -0.03877 | -0.03723 | -0.03877 | -0.03270 | 0.01406 | 0.00071 | 0.00000 |
12.00000 | 0.00000 | -0.00073 | -0.00145 | 0.03009 | 0.03259 | 0.02200 | 0.03259 | 0.03009 | -0.00145 | -0.00073 | -0.00000 |
13.50000 | 0.00000 | 0.00072 | 0.00554 | 0.06903 | 0.07345 | 0.08366 | 0.07345 | 0.06903 | 0.00554 | 0.00072 | -0.00000 |
15.00000 | 0.00000 | -0.04464 | -0.04466 | -0.02567 | -0.02477 | -0.01467 | -0.02477 | -0.02567 | -0.04466 | -0.04464 | -0.00000 |
16.50000 | 0.00000 | -0.00063 | 0.00179 | 0.05387 | 0.05690 | 0.05665 | 0.05690 | 0.05387 | 0.00179 | -0.00063 | -0.00000 |
18.00000 | 0.00000 | 0.00060 | 0.00119 | 0.02179 | 0.02245 | 0.02287 | 0.02245 | 0.02179 | 0.00119 | 0.00060 | 0.00000 |
19.50000 | 0.00000 | -0.00056 | -0.00336 | -0.01759 | -0.00785 | -0.00085 | -0.00785 | -0.01759 | -0.00336 | -0.00056 | 0.00000 |
21.00000 | 0.00000 | 0.02449 | 0.02449 | 0.01450 | 0.04749 | 0.04450 | 0.04749 | 0.01450 | 0.02449 | 0.02449 | 0.00000 |
22.50000 | 0.00000 | 0.00048 | -0.00072 | -0.01822 | -0.01988 | -0.00988 | -0.01988 | -0.01822 | -0.00072 | 0.00048 | 0.00000 |
24.00000 | 0.00000 | -0.00044 | -0.00087 | -0.00131 | -0.00175 | 0.01308 | -0.00175 | -0.00131 | -0.00087 | -0.00044 | -0.00000 |
25.50000 | 0.00000 | 0.00039 | 0.00202 | 0.01536 | 0.01660 | 0.02060 | 0.01660 | 0.01536 | 0.00202 | 0.00039 | -0.00000 |
27.00000 | 0.00000 | -0.01343 | -0.01342 | -0.01342 | -0.01341 | -0.01350 | -0.01341 | -0.01342 | -0.01342 | -0.01343 | -0.00000 |
28.50000 | 0.00000 | -0.00034 | 0.00025 | 0.00981 | 0.01072 | 0.01080 | 0.01072 | 0.00981 | 0.00025 | -0.00034 | -0.00000 |
30.00000 | 0.00000 | 0.00030 | 0.00060 | 0.00090 | 0.00120 | 0.00149 | 0.00120 | 0.00090 | 0.00060 | 0.00030 | 0.00000 |
Как видно из представленных таблиц табуляции относительная погрешность между точными и приближёнными значениями функции достигает 30%. Это говорит о том, что конечно-разностный метод дает плохую сходимость для волнового уравнения. График найденной функции при вышеуказанных параметрах в области приведен в приложении.
Дата: 2019-12-10, просмотров: 243.