Лабораторная работа №4 Тема. Численные методы решения линейных обыкновенных дифференциальных уравнений с краевыми условиями. Метод конечных разностей

Задание. Решить краевую задачу методом конечных разностей, построив два приближения (две итерации) с шагом h и с шагом h/2.

Проанализировать полученные результаты. Варианты заданий приведены в приложении 4.

Порядок выполнения работы

1. Постройте вручную конечноразностную аппроксимацию краевой задачи (конечноразностную СЛАУ) с шагом h, заданным вариантом.

2. Используя метод конечных разностей, сформируйте в Excel систему линейных алгебраических конечно-разностных уравнений для шага h разбивки отрезка [a, b]. Запишите эту СЛАУ на рабочем листе книги Excel. Расчетная схема приведена на рис.4.1.

3. Полученную СЛАУ решите методом прогонки.

4. Проверьте правильность решения СЛАУ с помощью надстройки Excel Поиск решения.

5. Уменьшите шаг сетки в 2 раза и еще раз решите задачу. Результаты представьте в графическом виде.

6. Сравните полученные результаты. Сделайте вывод о необходимости продолжения или о прекращении счета.

Решение краевой задачи с использованием электронных таблиц Microsoft Excel.

Пример 4.1. Методом конечных разностей найти решение краевой задачи , y(1)=1, y^’ (2)=0,5 на отрезке x[1, 2] с шагом h=0,2 и с шагом h=0,1. Сравнить полученные результаты и сделать вывод о необходимости продолжения или о прекращении счета.

Расчетная схема для шага h=0,2 приведена на рис.4.1.

Полученное решение (сеточную функцию) Y{1.000, 1.245, 1.474, 1.673, 1.829, 1.930}, Х {1; 1,2; 1,4; 1,6; 1,8;2} в столбце L и B можно принять за первую итерацию (первое приближение) исходной задачи.

Рис.4.1.

Для нахождениявторой итерации сделайте сетку вдвое гуще (n=10, шаг h=0,1) и повторите приведенный выше алгоритм.

Это можно проделать на том же или на другом листе книги Excel. Решение (второе приближение) приведено на рис.4.2.

Сравните полученные приближенные решения. Для наглядности можно построить графики этих двух приближений (двух сеточных функций), рис.4.3.

Рис.4.2.

Порядок построения графиков приближенных решений краевой задачи

1. Постройте график решения задачи для разностной сетки с шагом h=0,2 (n=5).

2. Активизируйте уже построенный график и выберите команду меню Диаграмма\Добавить данные

3. В окне Новые данные укажите данные x_i, y_i для разностной сетки с шагом h/2 (n=10).

4. В окне Специальная вставка установите флажки в полях:

• новые ряды,

• категории(значение оси х) в первом столбце.

Как видно из приведенных данных, два приближенных решения краевой задачи (две сеточные функции) отличаются друг от друга не более, чем на 5%. Поэтому за приближенное решение исходной задачи принимаем вторую итерацию, т.е.

Y {1, 1.124, 1.246, 1.364, 1.478, 1.584, 1.683, 1.772, 1.849, 1.914, 1.964}

Рис.4.3.

Лабораторная работа №5 Тема. Численные методы оптимизации. Графический метод

Задание. Решить задачу линейного программирования, приведенную в приложении 5, графическим методом и на ЭВМ. Сравнить полученные решения.