|
Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
|
|
Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ В ЧАСТНЫХ ПРОИЗВОДНЫХ. ⇐ ПредыдущаяСтр 8 из 8
ЦЕЛЬ РАБОТЫ Сформировать у студентов представление о применении дифференциальных уравнений в частных производных в различных областях физики и техники; выработать навыки решения смешанной задачи для уравнения теплопроводности методом сеток. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ 1. Изучить теоретическую часть. Выполните задания, соответствующие номеру Вашего варианта, и продемонстрируйте их преподавателю. 2. Оформите отчет по лабораторной работе, который должен содержать: · титульный лист; · исходные данные варианта; · решение задачи; · результаты решения задачи.
МЕТОДИЧЕСКИЕ РЕКОМЕНДАЦИИ Пример 16.1. Используя метод сеток, решить уравнение теплопроводности (уравнение параболического типа): Решение: 1. В командном окне программы МАТLAB наберем следующую последовательность операторов: > > h=0.1; s=1/6; > > k=h^2*s k = 0.0017 > > Nx=0.6/h+1 % число узлов координатной сетки Nx = 7.0000 > > Nt=0.01/k+1 % число шагов по времени Nt = 7.0000 > > for i=1: Nx+1 % задание пространственно-временной сетки x(i)=(i-1)*h; end; > > x x = 0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 > > for i=1: Nt+1 t(i)=(i-1)*k; end; > > t t = 0 0.0017 0.0033 0.0050 0.0067 0.0083 0.0100 % задание начальных и граничных условий > > for i=1: Nx+1 tab(1, i)=(3*x(i)*(1-x(i)))+0.12; end; > > for j=1: Nt+1 tab(j, 1)=2*(t(j)+0.06); tab(j, 7)=0.84; end; % матрица краевых условий > > tab tab = 0.1200 0.3900 0.6000 0.7500 0.8400 0.8700 0.8400 0.1233 0 0 0 0 0 0.8400 0.1267 0 0 0 0 0 0.8400 0.1300 0 0 0 0 0 0.8400 0.1333 0 0 0 0 0 0.8400 0.1367 0 0 0 0 0 0.8400 0.1400 0 0 0 0 0 0.8400 % заполнение оставшейся части матрицы > > for j=2: Nt+1 for i=2: Nx tab(j, i)=1/6*(tab(j-1, i-1)+4*tab(j-1, i)+tab(j-1, i+1)); end; end; > > tab tab = 0.1200 0.3900 0.6000 0.7500 0.8400 0.8700 0.8400 0.1233 0.3800 0.5900 0.7400 0.8300 0.8600 0.8400 0.1267 0.3722 0.5800 0.7300 0.8200 0.8517 0.8400 0.1300 0.3659 0.5704 0.7200 0.8103 0.8444 0.8400 0.1333 0.3607 0.5612 0.7101 0.8009 0.8380 0.8400 0.1367 0.3562 0.5526 0.7004 0.7920 0.8322 0.8400 0.1400 0.3524 0.5445 0.6911 0.7834 0.8268 0.8400 % визуализация численного решения (рис 16.1) > > surf(tab)
Рис. 16.1. Графическое решение уравнения теплопроводности ПРИМЕРНЫЕ ВОПРОСЫ НА ЗАЩИТЕ РАБОТЫ 1. Уравнения, какого вида называются дифференциальными уравнениями в частных производных (УЧП), и что является решением УЧП? 2. Какое УЧП называется линейным? 3. Назовите основные типы линейных ДУ с постоянными коэффициентами, и условие, в зависимости от которого происходит эта классификация. 4. Назовите и запишите примеры простейших УЧП эллиптического и параболического типов. 5. Что значит решить задачу Коши для УЧП? 6. Что является начальными и краевыми условиями для УЧП? 7. Что значит решить краевую и смешанную задачи для УЧП? 8. Какая задача носит название «Задачи Дирихле» для уравнения Лапласа, и что значит решить эту задачу? 9. Какие функции называются гармоническими в теории УЧП? 10. Какие методы численного решения УЧП? 11. В чем состоит суть метода конечных разностей (сеток)? 12. Какие точки называются внутренними, граничными I и II рода при решении ДУ методом сеток? 13. В чем заключается суть решения задачи Дирихле методом сеток? 14. По какой формуле находятся значения функции во внутренних узлах и внешних (граничных) при решении задачи Дирихле методом сеток? 15. За счет чего происходит грубое приближение искомых значений при решении задачи Дирихле методом сеток? 16. В чем заключается суть процесса усреднения Либмана, и для чего он применяется? 17. Как строится решения задачи Дирихле методом сеток? 18. Что значит решить уравнение теплопроводности методом сеток? 19. По какой формуле рассчитывается значение искомой функции во внутренних узлах сетки, при решении уравнения теплопроводности? 20. Алгоритм решения уравнения теплопроводности методом сеток. ЗАДАНИЕ Используя метод сеток, составить функцию, реализующую решение смешанной задачи для дифференциального уравнения параболического типа: Варианты заданий:
РАБОТА №18 ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ЦЕЛЬ РАБОТЫ Сформировать у студентов представление об интегральных уравнениях, выработать умения применять квадратурные методы решения интегральных уравнений Фредгольма и Вольтера, составлять и применять алгоритмы и программы для их решения, дать навыки в использовании программных средств для решения интегральных уравнений.
ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ 1. Изучить теоретическую часть. Выполните задания, соответствующие номеру Вашего варианта, и продемонстрируйте их преподавателю. 2. Оформите отчет по лабораторной работе, который должен содержать: · титульный лист; · исходные данные варианта; · решение задачи; · результаты решения задачи.
МЕТОДИЧЕСКИЕ РЕКОМЕНДАЦИИ Пример 18.1. Найти в пакете MATLAB решение интегрального уравнения Фредгольма:
Решение: 1. Создайте файл Q17.m (листинг 17.1), содержащий описание функции, возвращающей значения функции Q(t, s). Листинг 17.1. Файл Q17.m function z=Q17(t, s) z=1./sqrt(t+s.^2);
2. Создайте файл F17.m (листинг 17.2), содержащий описание функции, возвращающей значения функции f(t). Листинг 18.2. Файл F17.m function z=F17(t) z=sqrt(t+1)-sqrt(t+9)+t;
3. Создайте файл Solve_Q17.m (листинг 17.3), содержащей описание функции, возвращающей решение интегрального уравнения. Листинг 18.3. Файл Solve_Q17.m function [X, Y]=Solve_Q17(a1, b1, N, Lambda) % задание временной сетки h=(b1-a1)/(N-1); i=1: N; t(i)=a1+h*(i-1); s=t; %задание коэффициентов квадратурной формулы методом трапеций A(1)=0.5; m=2: N-1; A(m)=1; A(N)=0.5; %вычисление значений функции Q(t, s) в узлах сетки for i=1: N for j=1: N q(i, j)=Q17(t(i), s(j)); end; end; %вычисление значений функции f(t)в узлах временной сетки F=F17(t); for i=1: N for j=1: N if i==j M(i, j)=1-Lambda*A(i)*q(i, j)*h; else M(i, j)=-Lambda*A(i)*q(i, j)*h; end; end; end; %нахождение решения интегрального уравнения X=t; Y=M^-1*F';
4. Выполнить следующую последовательность команд: > > a1=1; > > b1=2; > > N=300; > > Lambda=1; > > [X, Y]=Solve_Q17(a1, b1, N, Lambda); > > plot(X, Y) Получим следующий график:
Рис. 17.1. Численное решение интегрального уравнения Фредгольма
Пример 18.2. Найти в пакете MATLAB решение интегрального уравнения Вольтерра:
Решение: 1. Создайте файл Q17_2.m (листинг 17.4), содержащий описание функции, возвращающей значения подынтегральной функции. Листинг 18.4. Файл Q17_2.m function z=Q17_2(t, s) z=t*sin(t*s.^3).^2;
2. Создайте файл F17_2.m (листинг 17.5), содержащий описание функции, возвращающей значения функции f(t). Листинг 18.5. Файл F17_2.m function z=F17_2(t) z=t.^2-1/4*tan(t.^4);
3. Создайте файл Solve2_Q17.m (листинг 17.6), содержащей описание функции, возвращающей решение интегрального уравнения. Листинг 18.6. Файл Solve2_Q17.m function [T, Y]=Solve2_Q17(t1, t2, N) % задание временной сетки h=(t2-t1)/(N-1); i=1: N; t(i)=t1+h*(i-1); s=t; %задание коэффициентов квадратурной формулы методом трапеций A(1)=0.5; m=2: N-1; A(m)=1; A(N)=0.5; %вычисление значений функции Q(t, s) в узлах сетки for i=1: N for j=1: N q(i, j)=Q17_2(t(i), s(j)); end; end; %вычисление значений функции f(t)в узлах временной сетки F=F17_2(t); %вычисление решения интегрального уравнения x(1)=F(1)./(1-A(1)*q(1, 1)); for m=2: N S=0; for j=1: m-1 S=S+h.*A(j).*q(m, j).*x(j); end; x(m)=F(m)+S./(1-h.*A(m).*q(m, m)); end; T=t; Y=x;
4. Выполнить следующую последовательность команд: > > t1=0; > > t2=5; > > N=300; > > [X Y]=Solve2_Q17(t1, t2, N); > > plot(X, Y) Получим следующий график:
Рис. 18.2. Численное решение интегрального уравнения Вольтерра
ПРИМЕРНЫЕ ВОПРОСЫ НА ЗАЩИТЕ РАБОТЫ 1. Что называется интегральным уравнением? 2. Назовите основные типы интегральных уравнений. 3. Уравнение, какого вида называют интегральным уравнением второго рода Фредгольма? 4. Уравнение, какого вида называют интегральным уравнением второго рода Вольтерра? 5. В чем состоит отличие интегрального уравнения Фредгольма от интегрального уравнения Вольтерры? 6. В чем состоит суть квадратурного метода решения интегральных уравнений Фредгольма? 7. В чем состоит суть квадратурного метода решения интегральных уравнений Вольтерры?
ЗАДАНИЕ 1. Построить каркас приближенного решения интегрального уравнения Фредгольма 2. Применяя квадратурную формулу прямоугольников на отрезке [a2; b2] с шагом h2, найти каракас приближенного решения интегрального уравнения Вольтерра
Исходные параметры для уравнений Фредгольма и Волтерры.
Исходные параметры для заданий.
ЛИТЕРАТУРА
1. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. – М.: Наука, 1987. 2. Вержбицкий В.М. Основы численных методов: Учебник для вузов. – 2-е изд.– М.: Высш.шк., 2005. 3. Гутер Р.С., Овчинский Б.В. Элементы численного анализа и математической обработки результатов опыта. – М.: Изд-во «Наука», 1970. 4. Данко П.Е. Высшая математика в упражнениях и задачах. В 2 ч. Учебное пособие для вузов/ П.Е. Данко, А.Г. Попов, Т.Я. Кожевникова. - М., 2005. 5. Поршнев С.В. Вычислительная математика. Курс лекций. – СПб.: БХВ-Петербург, 2004. 6. Поршнев С.В., Беленкова И.В. Численные методы на базе Mathcad. – СПб.: БХВ-Петербург, 2005. Содержание Предисловие...................................................................................................... 3 ЛАБОРАТОРНАЯ РАБОРА №1...................................................................... 4 ЛАБОРАТОРНАЯ РАБОРА №2..................................................................... 12 ЛАБОРАТОРНАЯ РАБОТА №3................................................................... 17 ЛАБОРАТОРНАЯ РАБОТА №4................................................................... 20 ЛАБОРАТОРНАЯ РАБОТА №5.................................................................... 28 ЛАБОРАТОРНАЯ РАБОТА №6................................................................... 38 ЛАБОРАТОРНАЯ РАБОТА №7................................................................... 43 ЛАБОРАТОРНАЯ РАБОТА №8.................................................................... 49 ЛАБОРАТОРНАЯ РАБОТА №9................................................................... 51 ЛАБОРАТОРНАЯ РАБОТА №10................................................................. 55 ЛАБОРАТОРНАЯ РАБОТА №11................................................................. 58 ЛАБОРАТОРНАЯ РАБОТА №12................................................................. 63 ЛАБОРАТОРНАЯ РАБОТА №13.................................................................. 64 ЛАБОРАТОРНАЯ РАБОТА №14................................................................. 72 ЛАБОРАТОРНАЯ РАБОТА №15................................................................. 77 ЛАБОРАТОРНАЯ РАБОТА №16................................................................. 80 ЛАБОРАТОРНАЯ РАБОТА №17.................................................................. 84 лиТЕРАТУРА…………………………………………………………………90
Учебное издание
ГОНЧАРЕНКО ГЕОРГИЙ ГРИГОРЬЕВИЧ
ЛАБОРАТОРНЫЙ ПРАКТИКУМ
Методические указания по выполнению лабораторно-практических работ
Компьютерная верстка Подписано в печать.Формат 60х84/16 Гарнитура Times. Печать оперативная. Бумага потребительская. Усл. печ. л.. Тираж 100 экз. Заказ №.
Московский авиационный институт (государственный технический университет) «МАИ» Отпечатано в типографии МАИ
|
Последнее изменение этой страницы: 2017-03-14; Просмотров: 990; Нарушение авторского права страницы
при заданных начальных условиях: u(x; 0)=3x(1-x)+0.12; u(0; t)=2(t+0.06); u(0.6; t)=0.84, где х 
[0; 0.6]. Решение найти при h=0, 1 для t 
.
и u(0.6; t)= 
, где х 
.
.
на сетке точек ti отрезка [a1; b1] c шагом h1, пользуясь какой-либо квадратурной формулой. На основе полученного каркаса записать приближенное решение в виде непрерывной функции (используя интерполяционные формулы) и с ее помощью вычислить приближенные значения x(c1) и x(d1).
. Представить полученное дискретное решение интерполяционным многочленом третьей степени, построенным по первым четырем узлам заданной сетки, и вычислить приближенно x(c2) и x(d2).
