WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Pages:   || 2 | 3 |

«АЛГ ОР ИТ МЫ МЕ Т О Д О В ВЗВЕ Ш Е ННЫ Х НЕВЯЗОК ДЛЯ РЕШЕНИЯ ЛИНЕЙНЫХ ЗАДАЧ МАТЕМАТИЧЕСКОЙ ФИЗИКИ И ИХ РЕАЛИЗАЦИЯ В СИСТЕМЕ MATHCAD Учебное пособие Ульяновск 2006 УДК 519.6 (075) ББК 22.311 ...»

-- [ Страница 1 ] --

А. В. Анкилов, П. А. Вельмисов, А. С. Семёнов

АЛГ ОР ИТ МЫ МЕ Т О Д О В

ВЗВЕ Ш Е ННЫ Х НЕВЯЗОК

ДЛЯ РЕШЕНИЯ ЛИНЕЙНЫХ ЗАДАЧ

МАТЕМАТИЧЕСКОЙ ФИЗИКИ

И ИХ РЕАЛИЗАЦИЯ В СИСТЕМЕ MATHCAD

Учебное пособие

Ульяновск 2006

УДК 519.6 (075)

ББК 22.311 я7

A 67

Рецензенты: Кафедра прикладной математики Ульяновского государственного университета (зав. кафедрой доктор физико-математических наук, профессор А. А. Бутов);

Доктор физико-математических наук, проф. УлГУ В. Л. Леонтьев.

Утверждено редакционно-издательским советом университета в качестве учебного пособия.

Анкилов, А. В. Алгоритмы методов взвешенных невязок для решения линейных задач математической физики и их реализация в системе A MathCAD: учебное пособие / А. В. Анкилов, П. А. Вельмисов, А. С. Семёнов. – Ульяновск : УлГТУ, 2006. – 168 с.

ISBN 5-89146-913- Пособие содержит изложение алгоритмов численного решения некоторых линейных краевых или начально-краевых задач математической физики.

Приведены постановки лабораторных работ, с использованием математически ориентированного пакета M athCAD. Даны примеры.

Пособие предназначено для студентов, изучающих специальные курсы современных численных методов, для аспирантов и инженеров, применяющих численные методы к решению прикладных задач.

Работа выполнена на кафедре «Высшая математика» УлГТУ. Печатается в авторской редакции.

УДК 519.6 (075) ББК 22.311 я Учебное издание АНКИЛОВ Андрей Владимирович ВЕЛЬМ ИСОВ Петр Александрович СЕМ ЁНОВ Алексей Степанович Алгоритмы методов взвешенных невязок для решения линейных задач математической физики и их реализация в системе MathCAD Учебное пособие Подписано в печать 20.12.2006. Формат 60 84/16.

Бумага офсетная. Усл. печ. л. 10,00.

Тираж 200 экз. Заказ Ульяновский государственный технический университет 432027, Ульяновск, ул. Северный Венец, 32.

Типография УлГТУ, 432027, Ульяновск, ул. Северный Венец, 32.

© А. В. Анкилов, П. А. Вельмисов, А. С. Семёнов, © Оформление. УлГТУ, ISBN 5-89146-913-

ОГЛАВЛЕНИЕ

Предисловие ……………………………………………………………… Программы по лабораторным работам ….…………………..……….. 1. Математическое моделирование физических задач …………….. 1.1. Вывод уравнений одномерной теплопроводнос ти ………..….. 1.




2. Постановка начально-краевой задачи одномерной стационарной теплопроводнос ти …………………………………… 1.3. Постановка начально-краевой задачи одномерной нестационарной теплопроводности ………………………………… 1.4. Постановка краевых задач двухмерной стационарной теплопроводности …………………………………………………… 1.5. Вывод уравнений поперечных колебаний струны …………… 1.6. Вывод уравнения продольных и крутильных колебаний cтержня ………………………………………………………………. 1.7. Постановка статических краевых задач для струны и стержня 1.8. Краевые задачи в теории колебаний струн и стержней ……… 2. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка методом Галеркина ……………………………………………………. 2.1. Подс тановка задачи и алгоритм метода ………………………. 2.2. Построение систем пробных и поверочных функций ………. 2.3. Задание к лабораторной работе ………………………………. 2.4. Выполнение работы в компьютерном классе ……………….. 2.5. Порядок выполнения лабораторной работы …………………. 2.6. Тестирующий пример ………………………………………….. 2.7. Вопросы для самоконтроля ……………………………………. 3. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка вариационным методом Ритца ……………………………………… 3.1. Подс тановка задачи и алгоритм метода ……………………… 3.2. Построение систем пробных и поверочных функций ………. 3.3. Задание к лабораторной работе ……………………………….. 3.4. Выполнение работы в компьютерном классе ……………….. 3.5. Порядок выполнения лабораторной работы …………………. 4. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка интегральным методом наименьших квадратов …………………… 4.1. Подс тановка задачи и алгоритм метода ……………………… 4.2. Задание к лабораторной работе ……………………………….. 4.3. Выполнение работы в компьютерном классе ………………… 4.4. Порядок выполнения лабораторной работы …………………. 4.5. Тестирующий пример ………………………………………….. 4.6. Вопросы для самоконтроля ……………………………………. 5. Решение начально-краевой задачи для одномерного параболического уравнения методом Галеркина..…………………. 5.1. Подс тановка задачи и алгоритм метода ………………………. 5.2. О построении функции u 0 ( x, t) ………………………………… 5.3. Задание к лабораторной работе ……………………………….. 5.4. Выполнение работы в компьютерном классе ……………….. 5.5. Порядок выполнения лабораторной работы …………………. 5.6. Тестирующий пример …………………………………………. 5.7. Вопросы для самоконтроля …………………………………… 6. Решение начально-краевой задачи для одномерного гиперболического уравнения методом Галеркина..……………….. 6.1. Подс тановка задачи и алгоритм метода ……………………… 6.2. Задание к лабораторной работе …………………………….…. 6.3. Выполнение работы в компьютерном классе ……………….. 6.4. Порядок выполнения лабораторной работы …………………. 6.5. Тестирующий пример …………………………………………. 6.6. Вопросы для самоконтроля ……………………………………. 7. Решение первой краевой задачи для двухмерного эллиптического уравнения методом Галеркина..………………….. 7.1. Подс тановка задачи и алгоритм метода ………………………. 7.2. Задание к лабораторной работе ……………………………….. 7.3. Выполнение работы в компьютерном классе ……………….. 7.4. Порядок выполнения лабораторной работы ………………….. 7.5. Тестирующий пример ………………………………………….. 7.6. Вопросы для самоконтроля ……………………………………. 8. Прикладной математический пакет «MathCAD»..……………… 8.2. Основные понятия и функции ………………………………… 8.3. Операторы математического анализа ………………………… 8.4. Функции и операторы матриц ………………………………… 8.5. Создание декартовых графиков на плоскости и Библиографический список ………..…………………………………… Целый ряд современных методов, предназначенных для решения самых разнообразных задач математической физики, базируется на идеях ученых Б. В. Галеркина и В. Ритца. К этим методам относятся, например, методы взвешенных невязок и вариационные методы [1,2].





В настоящем пособии представлены возможные алгоритмы применения метода Галеркина и интегрального метода наименьших квадратов, относящихся к группе методов взвешенных невязок, и вариационного метода Ритца при численном решении краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка, метода Галеркина при численном решении линейной начально-краевой задачи для одномерного параболического, одномерного гиперболического уравнений и первой краевой задачи для двухмерного эллиптического уравнения.

В новом методе можно быс трее разобраться, если решить конкретную задачу. В качестве источников таких задач в пособии описаны задачи одномерной нес тационарной теплопроводности, двухмерной стационарной теплопроводности и задачи о колебаниях струн и стержней.

Для проведения вычислительного эксперимента, согласно алгоритму метода, выбранного для решения конкретной задачи, в пособии приведены постановки лабораторных работ, с реализацией этого метода в математически ориентированном пакете MathCAD.

Пособие предназначено для студентов вузов, изучающих специальные курсы современных вычислительных методов. Оно будет полезным для аспирантов и инженеров, применяющих численные методы к решению прикладных задач.

Программы по лабораторным работам Для выполнения лабораторных работ разработаны в прикладной системе MathCAD 2000 professional файлы ODU.mcd, Parab.mcd, Giperb.mcd, Ellipt.mcd.

В начале каждого файла расписано задание на лабораторную работу, а в конце приведены варианты заданий, номера которых студенты получают от преподавателя. Тексты программ подробно описывают все дейс твия и операции. Для удобства изучения сами тексты с разобранными тестирующими примерами приведены в приложениях А, Б, В, Г данного пособия.

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

1. Математическое моделирование физических задач 1.1. Вывод уравнений одномерной теплопроводности Пусть дано материальное тело, расположенное между точками x = a и x = b оси Ox, продольный размер которого значительно превосходит размеры поперечного сечения, например, тонкий стержень, длинный трубопровод и т. д.

В дальнейшем будем называть это тело стержнем. Будем считать площадь S (x ) поперечного сечения (перпендикулярного оси Ox ) настолько малой, что всем точкам одного сечения в момент времени t можно приписать одну и ту же температуру u ( x, t). Будем считать, что стержень теплоизолирован вдоль боковой поверхности, а внутри стержня нет источников или стоков (поглотителей) тепла.

Рассмотрим элемент с тержня между его сечениями с абсциссами x и x + dx. Найдем количес тво тепла, которое накапливается в элементе за время dt. Согласно закону Фурье интенсивность q( x, t ) теплового потока в сечении x определяется выражением:

где K (x) – коэффициент теплопроводнос ти (K ( x) 0). Тогда разность dQ между количеством тепла, вошедшим в элемент через сечение x и вышедшим через сечение x + dx за время dt, будет равна:

Используя формулу Тейлора первого порядка с остаточным членом в форме Пеано для функций K ( x + dx ), S ( x + dx ), u x ( x + dx, t ), имеем Напомним, что символом o(x) обозначается величина бесконечно малая более высокого порядка, чем x.

С другой стороны, за счет притока тепла температура в элементе изменяется, и количество тепла dQ, поглощаемое элементом за время dt, равно где C (x ) – теплоемкость; (x) – объемная плотность вещества стержня Откуда, на основании теоремы о среднем для определяемого интеграла, получаем равенство преобразуется к виду Приравнивая, на основании закона сохранения энергии, выражения (1.1), (1.2) и осуществляя предельный переход при dt 0, получаем одномерное уравнение теплопроводности в виде Предположим теперь, что внутри стержня происходит выделение или поглощение тепла (это имеет мес то, например, при прохождении по телу электрического тока или вследствие происходящих в нем химических реакций).

Тогда количество тепла, накопленное в элементе с тержня за время dt за счет внутренних источников, будет равно dQ o = SF (x, t )dxdt, где F ( x, t ) – плотнос ть тепловых источников внутри стержня. Уравнение теплопроводности с учетом внутренних источников тепла принимает вид dQ = dQ + dQ o или Предположим далее, что на боковой поверхности стержня происходит теплообмен с окружающей средой. Тогда тепловой поток, проходящий за время dt через боковую поверхность элемента, согласно закону Ньютона пропорционален разнос ти температур поверхности тела и окружающей среды и определяется выражением где T (x, t ) – температура внешней среды; * ( x ) – коэффициент теплообмена, зависящий от свойств материала стержня и внешней среды, режима взаимодействия (условий контакта) стержня с внешней средой, а также от теплопроводности с учетом внутренних источников тепла и теплообмена на боковой поверхности имеет вид или Заметим, что если тепло распространяется в жидкости, которая движется со скоростью V ( x, t ) параллельно оси x, то уравнение теплопроводности запишется следующим образом:

Если тело однородно, т. е. С,, K – постоянные, и площадь сечения S постоянна, то уравнение (1.4) можно записать в виде где a 2 = K / C – коэффициент температуропроводнос ти.

Аналогично уравнению (1.4) выводится уравнение, описывающее процесс распространения тепла в трехмерных телах или в развернутой форме Для однородных тел уравнение (1.5) удобно представить в виде Для двухмерных тепловых полей в плас тинах, тонких плитах уравнение (1.5) примет вид или для однородных плас тин Согласно (1.4) стационарное (установившееся во времени) распределение теплового поля в стержне постоянного поперечного сечения (S ( x) = const ) описывается уравнением В (1.7) введены обозначения Перечислим основные типы граничных условий (на примере левого конца стержня при x = a ).

а) Известная температура при x = a : y (a) = Ta.

б) Задана интенсивность теплового потока через торцевое сечение x = a :

K a y (a ) = q a, K a = K (a ). В частности, если стержень теплоизолирован при в) На конце x = a имеет место теплообмен с окружающей средой известной температуры Ta :

Здесь a – коэффициент теплообмена на конце x = a. Последнее условие (условие Ньютона) означает, что тепловой поток, передаваемый в единицу времени с единицы площади поверхнос ти в окружающую среду, пропорционален разности температур поверхности тела и окружающей среды.

Аналогичные краевые условия могут быть заданы и на правом конце стержня при x = b. Например, условие теплообмена при x = b имеет вид В таблице 1.1 приведены возможные варианты краевых условий для определения стационарного распределения температуры в стержне согласно уравнению (1.7).

Напомним еще раз используемые в таблице 1.1 обозначения:

K a = K (a), K b = K (b) – коэффициенты теплопроводности;

а, b – коэффициенты теплообмена на левом и правом концах стержня соответс твенно;

Tа, Tb – температуры, которые поддерживаются на концах с тержня при x = a и при x = b ;

qа, qb – интенсивности тепловых потоков при x = a и при x = b.

Очевидно, что все приведенные в таблице 1.1 варианты краевых условий можно записать в виде при соответствующем выборе значений коэффициентов ai, bi.

Например, для первого варианта условий из таблицы 1.1 имеем а для девятого – Таким образом, математическая задача одномерной стационарной теплопроводности формулируется следующим образом: требуется найти функцию y (x), удовлетворяющую на отрезке [a, b] обыкновенному линейному дифференциальному уравнению (1.7), а на концах отрезка – граничным условиям (1.8).

Варианты краевых условий для уравнения (1.5) 1.3. Постановка начально-краевой задачи одномерной В пункте 1.2 рассмотрена краевая задача для одномерного с тационарного уравнения теплопроводности (1.7), которая представляет собой краевую задачу для обыкновенного дифференциального уравнения второго порядка. В случае нестационарной теплопроводности к краевым (граничным) условиям (1.8) добавляется начальное условие в некоторый начальный момент времени t = t (обычно t = 0 ) и говорят, что задана начально-краевая задача для уравнения параболического типа (1.4).

Согласно (1.6) стационарное (установившееся во времени) распределение теплового поля в пластине описывается уравнением При решении краевых задач для уравнения эллиптического типа (1.10) наиболее часто используются три типа краевых условий.

а) Краевая задача с граничными условиями первого рода (первая краевая задача).

Требуется найти решение уравнения (1.10) в некоторой области D, принимающее на границе этой области заданные значения. Т. е. нужно найти (применительно к рассматриваемой задаче) стачионарное распределение температуры внутри области, если задана температура на границе этой области Здесь D – граница области D, g (x, y ) – известная функция.

б) Краевая задача с граничными условиями второго рода (вторая краевая задача).

Требуется найти решение уравнения (1.10) в некоторой области, на границе которой задана внешняя нормальная производная (т. е. на границе задана интенсивность теплового потока).

где q =. Здесь D – граница области D, q( x, y ) – интенсивность теплового потока. При этом, если q 0, то тепловой поток направлен наружу, а если q 0, то тепловой поток направлен внутрь облас ти. При q = 0 имеем условие теплоизоляции в) Краевая задача с граничными условиями третьего рода (третья краевая задача).

Требуется найти решение уравнения (1.9) в некоторой области, которое удовлетворяет на границе условию где =. Здесь D – граница области D, на которой задан теплообмен с окружающей средой, температура которой равна T ; – коэффициент теплообмена.

Если на различных частях границы D заданы условия различного рода, то такие условия и соответствующие им задачи называют смешанными.

1.5. Вывод уравнений поперечных колебаний струны Рассмотрим тонкую гибкую упругую нить (струну), которая в положении равновесия занимает отрезок [a, b] оси Ox и концы которой закреплены.

Полагая струну тонкой, пренебрегаем весом струны по сравнению с внутренними силами натяжения и внешней нагрузкой. Полагая струну гибкой, считаем, что внутренние усилия, возникающие в струне, направлены по касательной к мгновенному профилю в каждой точке, т. е. струна не сопротивляется изгибу. Предполагаем также, что внешние силы лежат в вертикальной плоскости, в которой совершают колебания точки струны.

Рассмотрим элемент струны между точками x и x + dx (рис. 1.1) и обозначим смещение точек струны через u ( x, t ), а длину элемента струны через откуда, предполагая смещение струны u ( x, t ) малыми настолько, что получаем ds dx, т. е. в пределах принятой точности удлинения участков струны в процессе колебаний не происходит. Следовательно, согласно закону Гука величина натяжения в каждой точке струны не меняется со временем и является функцией только x, т. е. T = T (x).

Рис. 1.1. Иллюстрация к выводу уравнения колебаний струны Запишем условия динамического равновесия элемента струны, на который действуют в плоскости Oxu силы натяжения T1 = T (x ), T2 = T (x + dx ), внешняя распределенная по длине дуги с линейной плотнос тью F ( x, t ) поперечная сила и сила инерции, направленная вдоль оси Ou.

Проектируя силы на ось Ox, получаем Так как, согласно тождествам тригонометрии и геометрического смысла производной, то, учитывая условие (1.15), из (1.16) получим T (x + dx ) = T ( x). Откуда, в силу произвольности выбора точек x и x + dx, следует, что величина натяжения не зависит и от x, т. е. является постоянной, T (x ) = T0 = const.

Проектируя теперь все силы на ось Ou, получаем где (x ) – линейная плотность струны.

Аналогично формулам (1.16) устанавливаем откуда, согласно условию (1.14), имеем Теперь, применяя для входящих в формулу (1.17) интегралов теорему о среднем, а для u ( x + dx, t ) – формулу Тейлора первого порядка с остатком в форме Пеано, получаем 1 и 2 принадлежат отрезку [x, x + dx ]. Почленно деля последнее где равенство на dx и осуществляя предельный переход при dx 0, получаем уравнение колебания струны следующего вида:

Если струна дополнительно по всей длине связана с вязкоупругим основанием, то для описания ее колебания можно получить уравнение где ( x, t ), ( x, t ) – коэффициенты жесткос ти и демпфирования основания;

Q ( x, t, ) – ядро релаксации, учитывающее изменение с течением времени физико-механических свойств материала основания (т. е. его старение).

Заметим, что при выводе уравнения (1.19) предполагалось, что реакция основания пропорциональна его деформации (модель Винклера).

В статистических задачах профиль струны u = u (x) определяется согласно (1.12), решением уравнения 1.6. Вывод уравнения продольных и крутильных колебаний стержня Для вязкоупругого тела при одномерном растяжении (сжатии) связь между деформацией (относительным удлинением) ( x, t ) и напряжением ( x, t ) представляется формулой где E – модуль упругости; R – ядро релаксации, учитывающее старение материала тела; – коэффициент внутреннего трения. Заметим, если R 0 и 0, то получаем закон Гука для упругого тела.

Рассмотрим элемент с тержня (рис. 1.2), заключенный между поперечным сечением с координатами x и x + dx.

Рис. 1.2 Иллюстрация к выводу уравнения продольных колебаний стержня В сечении « x » на элемент действует сила N ( x, t ) = (x, t )S (x ), где S (x ) – площадь сечения, в сечении « x + dx » – сила N ( x + dx, t ) = ( x + dx, t )S ( x + dx )).

Предполагая, что на стержень действует внешняя нагрузка, распределенная по длине стержня с объемной плоскостью F ( x, t ), аналогично выводу уравнения (1.18) получаем уравнение продольных колебаний стержня следующего вида:

где (x) – объемная плотнос ть материала с тержня; u ( x, t) – продольное смещение сечения стержня с координатой x в момент времени t от положения, которое занимало это сечение, когда стержень находился в ненапряженном состоянии.

Учитывая, что и подставляя (1.21) в (1.22), имеем Если боковая поверхность стержня скреплена с вязкоупругим основанием (модель Винкера), то приходим к следующему уравнению:

где (x, t ), ( x, t) – коэффициенты жесткос ти и демпфирования основания;

Q ( x, t, ) – ядро релаксации основания. Заметим, что форма записи уравнения (1.23) не изменится, если считать S и зависящими от времени t.

определяются, согласно (1.23), решением уравнения Для вязкоупругого стержня, находящегося в состоянии кручения (рис. 1.3), связь между напряжением, вызванным сдвигом образующей на угол, и этим углом может быть предс тавлена формулой где G – модуль сдвига; R – ядро релаксации стержня; – коэффициент внутреннего трения.

Заметим, если R 0, 0, то получаем известный закон сдвига для упругого тела.

Если обозначить через u ( x, t) угол поворота сечения с координатой x, то (см. рис. 1.3) из равенства rdu = dx, имеем Рис.1.3. Иллюстрация к выводу уравнения крутильных колебаний стержня.

соответс твующем координате x, определяется формулой Отсюда, используя выражения (1.25), (1.26), получаем где J 0 = r 2 dS – полярный момент инерции сечения.

Рассмотрим элемент стержня, заключенный между поперечными сечениями с координатами x и x + dx (рис. 1.3). В сечении « x » действует крутящий момент M ( x, t), в сечении « x + dx » – M ( x + dx, t ). Предполагая, что на с тержень действует крутящий момент внешних сил, распределенный по длине стержня с линейной плотностью F ( x, t ), из уравнения динамического равновесия получаем где – плотность с тержня; 1 и 2 – принадлежат [x, x + dx ]. Откуда аналогично уравнению (1.18) получаем уравнение крутильных колебаний стержня которое, с учетом (1.27), принимает вид Если боковая поверхность стержня скреплена с вязкоупругим основанием (модель Винклера), то для описания крутильных колебаний приходим к уравнению где,, Q – коэффициенты жес ткости, демпфирования и ядро релаксации основания.

Заметим, что форма записи уравнения (1.27) не изменится, если считать и J 0 функциями двух переменных x и t.

Статические углы поворота u (x) сечений стержня при кручении определяются, согласно (1.28), решением уравнения 1.7. Постановка статических краевых задач для струны и стержня В с татическом варианте профиль струны, продольные и угловые перемещения сечений стержня, согласно (1.20), (1.24) и (1.29), определяются решением уравнения K ( x) = T0, g (x ) = F ( x), если рассматривается задача (1.20);

K ( x) = J 0 (x )G( x), g ( x) = F ( x), если – задача (1.29).

Перечислим основные типы граничных условий при x = a для уравнений (1.20), (1.24), (1.29) в обозначениях уравнения (1.30).

а) y (a) = 0 ; это условие соответствует жесткому закреплению левого конца струны и стержня.

б) K (a ) y (a ) = qa ; это условие соответствует заданию на левом конце стержня продольной силы N (a ) = q a для задачи (1.24) и заданию крутящего момента M (a ) = qa в случае задачи (1.29). В частности, если левый конец свободен, то qa = 0.

в) K (a ) y (a ) = a y (a ) ; это условие соответс твует упругому закреплению левого конца стержня, когда qa = a y (a) ( qa или равно N (a ), или – M (a ) ), где a – соответс твующий задаче (1.24) или (1.29) коэффициент закрепления.

Аналогичные краевые условия могут быть заданы и на правом конце струны или стержня при x = b. Очевидно, что все возможные варианты краевых условий для уравнения (1.30) можно получить из условий (1.8) при соответс твующем выборе значений коэффициентов ai, bi.

Таким образом, рассматриваемые статические краевые задачи для струны и стержня математически формулируется так же, как и задача стационарной теплопроводности из раздела 1.2.

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

Уравнение продольных колебаний струны Уравнение продольных колебаний стержня Уравнение крутильных колебаний стержня Уравнения (1.31)–(1.33) являются уравнениями гиперболического типа.

Рассмотрим гармонические колебания упругих тел. В этом случае решение уравнений (1.31)–(1.33) и приложенную внешнюю нагрузку F ( x, t ) предс тавим в виде:

где (частота колебаний) и – постоянные. Тогда для u * (x ) = y (x ) получим уравнение (1.30), в котором F (x ) следует заменить на F* (x ), а (x) – на * (x ), где * ( x) = ( x) (x ) 2 соответствует уравнению (1.31), Приведем основные типы граничных условий при x = a.

а) u ( x, t ) = a (t ); это условие соответствует движению левого конца струны или стержня по закону a (t ).

стержня продольной силы N (a, t ) = q a (t ) для задачи (1.32) и заданию крутящего момента M (a, t ) = qa (t ) в случае задачи (1.33). В частнос ти, если левый конец свободен, то qa = 0.

закреплению левого сечения стержня, движущегося (вращающегося) по закону Предполагая функции a (t ), a (t), q a (t ) периодическими во времени, аналогично (1.34) положим где a, a, qa – постоянные. Тогда для u * (x ) = y (x ) будем иметь граничные условия следующего вида:

Условия на правом конце x = b задаются аналогично.

Замечание. Аналогичные краевые задачи получим в случае, когда где n, n, n, n, q n, q0 – постоянные; u 0 ( x ) – решение с тационарных краевых задач, описанных в (1.5); а u n (x ) – решение краевых задач, рассмотренных выше в этом параграфе.

Следует иметь в виду, что частота колебаний n являются в общем случае неизвестными величинами, определяемыми в процессе решения задачи.

2. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка методом Рассмотрим следующую краевую задачу: требуется на отрезке [a, b] найти решение Y (x ) дифференциального уравнения удовлетворяющее условиям где p( x), q( x), f (x ) – заданные функции, непрерывные на [a, b]; a0, a1, a 2, b0, b1, b2 – заданные действительные числа, причем a0 + a12 0, b0 + b12 0.

Напомним, что в отличие от имеющей всегда единс твенное решение задачи Коши для уравнения (2.1), краевая задача (2.1), (2.2) может иметь или одно решение, или бесконечно много решений, или, наконец, может совсем не иметь решений.

Везде далее будем предполагать существование единс твенного решения Y (x ) поставленной краевой задачи, что часто вытекает из физического смысла того явления или процесса, математическое моделирование которого привело к задаче (2.1), (2.2).

В методе Галеркина для нахождения приближенного решения рассматриваемой задачи строится функциональная последовательнос ть {y n (x)} из пробных решений y n (x ) следующим образом.

Задаемся на отрезке [a, b] некоторой системой дважды непрерывно дифференцируемых функций u 0 ( x), u1 ( x),..., un (x ) таких, что u0( x) удовлетворяет краевым условиям (2.2), а функции u1 (x ), u 2 ( x),..., un ( x), называемые пробными функциями, линейно независимы на [a, b] и удовлетворяют однородным краевым условиям Составляем функцию с неизвестными пока постоянными коэффициентами C1, C 2,..., C n. Подчеркнем, что в силу линейнос ти условий (2.2), функция (2.4) при любых значениях C1,..., C n удовлетворяет этим условиям. Подставляя функцию y n (x ) из (2.4) вместо y (x) в уравнение (2.1), получаем функцию которая называется невязкой. Как видно из (2.5), невязка линейно зависит от параметров C1, C 2,..., C n и является характеристикой уклонения функции (2.4) от точного решения Y (x ) задачи (2.1), (2.2). Во всяком случае, если при некоторых значениях параметров C1, C 2,..., C n невязка на [a, b ] тождественно равна нулю, то Y (x ) yn (x ) в силу единственности Y (x ).

Однако в общем случае невязка оказывается отличной от нуля. Поэтому подбираем значения параметров C1,..., C n так, чтобы невязка в каком-то смысле была бы наименьшей. В обобщенном методе Галеркина значения параметров C1,..., C n определяются из системы уравнений где а W1 ( x),...,Wn (x ) – заданные непрерывные и линейно независимые на [a, b] функции, часто называемые поверочными функциями. Заметим, что если в качестве поверочных функций взять пробные, то получится метод Галеркина в авторском варианте [1]. Заметим также, что если W1 ( x),...,Wn (x ) входят в полную сис тему функций, то при n равенства (2.6) свидетельс твуют об ортогональнос ти невязки всем элементам полной системы [3]. Значит, невязка сходится при n к нулю в среднем, и можно ожидать сходимости последовательности (2.4) к точному решению Y (x ) в среднем, т. е.

Записав условие (2.6) в развернутом виде, для определения значений параметров получаем неоднородную систему линейных алгебраических уравнений n-го порядка где Решив систему (2.8) и подставив определяемые этим решением значения параметров C1,..., C n в (2.4), заканчиваем построение пробного решения y n (x ).

Опишем теперь возможный алгоритм приближенного решения задачи (2.1), (2.2) методом Галеркина, предполагая, что y n (x ) сходится к Y (x ) при 1. Подготовленный шаг алгоритма. На этом шаге выбираем функцию u 0 ( x ), пробные функции u1 (x ),..., u n ( x ) и поверочные функции W1 ( x),...,Wn (x ).

Находим функцию R 0 ( x ) = L[u0 ] f (x ), т. е. невязку от подстановки u 0 ( x ) в уравнение (2.1). Если x [a, b ] : R0 ( x) = 0, то u 0 ( x ) = Y ( x), и вычисления заканчиваем. Если же R 0 ( x) 0, то переходим к следующему шагу алгоритма.

2. Первый шаг алгоритма. Строим y1 (x ) = u0 ( x) + C1u1 ( x), определив значение C1 из решения системы (2.8) при n = 1. Находим невязку Если x [a, b ] : R (C1, x ) = 0, то Y (x ) = y1 ( x), и задача решена, если же R (C, x ) 0, то находим или Если 11 1 или 12 2, где 1 и 2 заданные меры точнос ти приближенного решения, то полагаем Y (x ) y1 ( x) и вычисления заканчиваем, если же 11 или 12 2, то переходим к вычислениям на следующем шаге и т. д.

Таким образом, на m -м (m 1) шаге алгоритма строим функцию определив значения C1,..., C m из решения системы (2.8) при n = m, и определяем невязку заканчиваем.

Если R (C1,..., C m, x) 0, то находим – переходим к (m + 1) -му шагу.

2.2.Построение систем пробных и поверочных функций Известно, что степенные функции 1, x, x 2,..., x n,... линейно независимы на всей числовой прямой R и, следовательно, на любом ее отрезке [a, b] R.

Покажем, что на любом отрезке [a, b] линейно независима любая система многочленов последовательных с тепеней. Рассмотрим произвольную систему многочленов:

LLLLLLLLLLLLLLLLLL

и решим относительно неизвестных 0, 1,..., n определенное на R тождество Из условий тождественного равенства нулю многочлена n -й степени (равенство нулю коэффициентов при всех степенях x ) последовательно получаем

Таким образом, условие (2.10) выполняется тогда и только тогда, когда 0 = 1 =... = n = 0, т. е. система многочленов P0 ( x ),..., Pn ( x) и любая подсистема из них линейно независима на R и, следовательно, на любом [a, b] R.

Для построения u 0 ( x ) и линейно независимой на [a, b], системы пробных функций u1 (x ),..., un (x ), являющихся многочленами, можно применить метод неопределенных коэффициентов.

Например, предположим u 0 = A = P0 ( x), из условий (2.2) получаем систему линейных алгебраических уравнений относительно A В том случае, когда эта система совместна, коэффициент A определяется. Если система не совместна, то ищем аналогичным образом u 0 ( x ) в виде u 0 ( x ) = A + Bx = P1 (x ) и т. д., до тех пор, пока не будет найдена u 0 ( x) = Pr0 (x ), удовлетворяющая условиям (2.2).

Далее, используя условия (2.3), методом неопределенных коэффициентов определяем последовательно так же, как и u 0 ( x ),

LLLLLLLLLL

Пример 1. Построить u 0 ( x ) и систему из пяти пробных функций для задачи с краевыми условиями Решение. Пусть u0 (x) = A, тогда u = 0 и условия (2.11) дают несовместную систему из уравнений A = 1 и A = 4.

Пусть u 0 = A + Bx, тогда u 0 = B и условия (2.11) дают Определяем u1 (x ). Если u1 = A или u1 = A + Bx, то однородные условия, соответс твующие условиям (2.11), выполняются, если u1 0, что невозможно из-за требования линейной независимости пробных функций.

условий, соответствующих (2.11), получаем систему Решая ее методом Гаусса, имеем Видим, что система имеет множество решений Выбираем одно решение из G при =, тогда Аналогично, используя формулу u k = A0 + A1x +... + Ak +1 x k + 1, находим Пример 2. Построить u 0 ( x ) и систему из трех пробных функций для задачи с краевыми условиями Решение. Если u 0 ( x ), то условия (2.12) приводят к несовместной системе Предположим, что u 0 = A + Bx, тогда u = B и условия (2.12) дают тоже несовместную систему.

Полагаем u 0 = A + Bx + Cx 2, тогда u = B + 2Cx и условия (2.12) дают которая несовместна.

(2.12) имеем Решаем полученную систему методом Гаусса в матричной форме, чтобы найти все решения системы.

Прямой ход метода:

ABCD ABCD ADCB

Видим, что система совместна, ибо ранг матрицы системы (rg ) равен рангу расширенной матрицы и равен 2. Так как число неизвестных системы 4 больше rg = 2, то система неопределена, и все множество решений G 0 системы получаем обратным ходом метода Гаусса, придавая двум неизвестным C и B произвольные значения. Получаем Выбираем решение из G 0 при 1 = 2 = 0. Тогда Определяем теперь u1 (x ). Если u1 (x ) = A 0, то однородные условия, соответс твующие условиям (2.12), выполняются при A = 0, что недопустимо.

Пусть u1 (x ) = A + Bx, u1 (x ) = B, и из однородных условий, соответс твующих условиям (2.12), имеем Эта система неопределена, ее множество решений Выбираем одно ненулевое решение при = 1, тогда u1 (x ) =1 x.

однородные условия дают систему Решая ее методом Гаусса, находим множество решений Выбирая одно ненулевое решение ( C 0 ), при 1 = 2 = 1, получаем u 3 ( x) = B + 2Cx + 3Dx 2, и из однородных условий имеем систему которая противоречит условию D 0.

u 3 ( x) = B + 2Cx + 3Dx 2 + 4 Ex 3, и из однородных условий получаем систему Решая ее методом Гаусса, получаем множество решений Выбирая одно ненулевое решение ( E 0 ) при 1 = 1, 2 = 1, 3 = 1, имеем Применяя метод неопределенных коэффициентов, можно строить системы функций, используя другие системы линейно независимых на R функций, такие как Важным ис точником для пос троения ортогональных на [a, b] пробных функций является множество решений задачи, называемой задачей на собственные значения для дифференциального оператора L[ y ] = y.

Рассмотрим конкретный пример такой задачи.

Пример 3. Требуется найти действительные значения параметра, при которых существуют нетривиальные решения дифференциального уравнения удовлетворяющие однородным условиям Решение. Пус ть = 0, тогда общее решение уравнения (2.13) будет иметь вид y = C1 x + C 2. Пытаясь удовлетворить условиям (2.14), получаем Таким образом, = 0 не является собственным значением, так как ему соответс твует единс твенное тривиальное ( y 0 ) решение задачи (2.13), (2.14).

системе уравнений Следовательно, среди отрицательных действительных чисел собственных значений задачи (2.13), (2.14) нет.

Пусть теперь 0. Тогда y = C1 cos(x) + C 2 sin(x), =, и краевые условия (2.14) дают Видим, что существуют нетривиальные решения задачи (2.13), (2.14), если Таким образом, множество собственных значений определяется формулой n = (n ) 2, n = 1, 2,..., а множество собственных функций, соответс твующих собственному значению n, имеет базисную функцию y n ( x), y m ( x ) (n m), достаточно проверить, что При выборе систем поверочных функций полезно вспомнить и о других системах функций, ортогональных на некотором отрезке. Например, известно [3], что многочлены Лежандра, определяемые формулой ортогональны на [–1,1]. Так что, если в качестве поверочных функций Wk (x) решено взять, например, первые пять многочленов Лежандра, ортогональных на [a, b], то в первые пять выражений из (2.15):

следует подставить Методом Галеркина найти наиболее точное приближенное решение краевой задачи построенное при помощи системы из n пробных функций – многочленов и двух систем поверочных функций, одна из которых составлена из пробных функций, а вторая – из многочленов Лежандра. За меру точности выбрать (по указанию преподавателя) или Варианты заданий, определяемые различными наборами заданий параметров d 0, d1, d 2, a 0, a1, a 2, b0, b1, b 2, a, b задачи (2.16) приведены в таблице 2.1.

Лабораторная работа выполняется с использованием прикладной системы MathCAD, в которой реализуется алгоритм построения пробных решений y m (x) методом Галеркина.

Перед обращением к программе необходимо подготовить числовые и строчные данные, вводимые в процессе диалога с клавиатуры дисплея.

Числовые данные:

a, b – концы отрезка интегрирования;

n – максимальное число параметров C i в пробном решении.

Значение параметра n задает преподаватель.

Строчные данные:

Аналитические выражения для функций u 0 ( x ),..., u n ( x) и для функций W1 ( x),..., Wn ( x ) ;

Аналитические выражения набираются клавишами дисплея по определенным правилам, которые программа при необходимости может напомнить.

В результате расчета программа выводит на экран дисплея значения коэффициентов С i, таблице всех пробных решения и их невязок.

Анализируя данные этих таблиц, необходимо найти обоснованный ответ на поставленную в лабораторной работе задачу.

вар.

вар.

2.4. Выполнение работы в компьютерном классе 1. Прежде чем начать выполнение лабораторной работы на ЭВМ, внимательно ознакомьтесь с данной инструкцией.

2. При необходимости включите сами (или попросите лаборанта) питание компьютера. После того, как система загрузится, запускаем двойным щелчком левой кнопки мыши на рабочем столе программу Mathcad, если же ярлык отсутс твует, тогда открываем программу через кнопку «Пуск» (Программы Mathsoft Mathcad).

3. Узнайте у лаборанта расположение файла ODU.mcd и откройте его (File Open или, если программа русифицирована, Файл Открыть). При любой ошибке ввода программы нужно обратиться к лаборанту.

4. Прочитайте в начале файла задание на лабораторную работу и просмотрите пример выполнения работы, для которого исследование уже проведено. Программа файла ODU.mcd состоит из шести пунктов «1.

Постановка задачи», «2. Получение точного решения», «3. Получение приближенного решения методом Галеркина», «4. Получение приближенного решения вариационным методом Ритца», «5. Получение приближенного решения интегральным методом наименьших квадратов», «6. Выводы». Для выполнения лабораторной работы «Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка методом Галеркина» необходимо использовать пункты 1, 2 и 3. Цели и задачи каждого из пунктов описаны ниже.

5. Для набора функций нужно либо воспользоваться всплывающим меню инструментов «Calculator», либо ввести ее с клавиатуры, используя следующие символы арифметических действий и стандартных функций: сложение – ‘+’;

вычитание – ‘–‘; умножение – ‘*’; деление – ‘/’; возведение в степень – ‘^’;

квадратный корень – ‘\’; синус – sin(x); косинус – cos(x); экспонента – exp(x);

натуральный логарифм – ln(x). При вводе числовых данных, являющихся десятичными дробями, целую и дробную части нужно разделять точкой (например, 0.5, 1.5 и т. д.).

6. Порядок выполнения работы Вам укажет программа подсказками и заданиями, выделенными красным цветом.

2.5. Порядок выполнения лабораторной работы Рекомендуется следующий порядок выполнения лабораторной работы.

1. Изучить разделы 1.1–1.2, 2.1–2.4 и подготовить ответы на контрольные вопросы из раздела 2.6.

2. Пройти собеседование с преподавателем, получить допуск к выполнению работы в диалоге с ПЭВМ, номер варианта задания и значение параметра n.

3. В соответс твии с вариантом задания выполнить подготовительный шаг алгоритма метода Галеркина и подготовить, если u 0 ( x) не является точным решением задачи, все числовые и строчные исходные данные для расчетов на ПЭВМ.

4. Выполнить основную расчетную часть лабораторной работы в системе MathCAD. Следует переписать с экрана дисплея значения коэффициентов C i пробных решений, составить итоговые таблицы пробных решений и их невязок.

5. Оформить и защитить отчет по лабораторной работе, который должен содержать титульный лист, математическую постановку задачи, результаты выполнения подготовительных расчетов, основные результаты расчетов на ПЭВМ и итоговый отчет.

Методом Галеркина найти на [0, 1] приближенное решение краевой задачи 1. Запускаем программу MathCAD. Открываем файл ODU.mcd (текст программы приведен в приложении A). В пункте «Пос тановка задачи» вводим числовые параметры a0, a1, a 2, b0, b1, b2, a, b и функции p( x), q (x ), f (x ), входящие в задачу (2.17) Замечание. Для задачи (2.16) необходимо еще ввести числовые параметры d 0, d 1, d 2, входящие в функции p( x), q (x ).

2. В пункте «Получение точного решения в системе MathCAD» записываем дифференциальное уравнение (2.17) в виде нормальной системы дифференциальных уравнений второго порядка Далее с помощью функции bvalfit (см. п. 8.2) краевую задачу приводим к задаче Коши, получая начальные условия Решая полученную задачу Коши для нормальной системы второго порядка с помощью функции rkfixed (см. п. 8.2), находим решение дифференциального уравнения, разбив отрезок [0, 1] на N = 100 частей, т. е. с шагом интегрирования h = 0,01 (в дальнейшем будем называть его точным решением).

Выписываем полученные значения этого решения y (x) в промежуточных точках (табл. 2.2), чтобы в дальнейшем оценить дейс твительную точнос ть пробных решений.

y 0,846776 0,86522 0,886703 0,905491 0,914371 0, y 0,863987 0,779347 0,632877 0,402849 0, График точного решения имеет вид Заметим, что поставленная задача имеет единс твенное точное решение вида которое получено аналитическим методом, известным из теории линейных дифференциальных уравнений с постоянными коэффициентами.

Если найти значения этого решения y (x) в тех же промежуточных точках, то получим значения таблицы 2.2, т. е. компьютерное решение найдено с точностью 10 6. Если сравнить точное аналитическое решение с решением, полученным в системе MathCAD, видим, что локальная погрешность не превышает 1,77 10 8. Следовательно, решение найдено достаточно точно.

Замечание. Компьютерное решение будет тем точнее, чем больше число точек разбиения введено в функцию rkfixed. Например, при N = 1000, локальная погрешность не превышает 3,45 10 12.

3. В пункте «Получение приближенного решения методом Галеркина»

вводим порядок пробного решения n := 5. В качестве пробных функций u 0 (x), u1 ( x),..., u5 (x ) используем функции, построенные в примере 1 раздела 2.2, Тогда 3.1. Воспользовавшись указаниями из раздела 2.2, в качестве поверочных функций возьмем пробные u1 (x ),..., u5 ( x) В результате расчета по программе при n = 5 получим вектор коэффициентов следовательно, пробное решение имеет вид Основные результаты расчета при n 5 (т. е. подставляя последовательно n = 0,1, 2,3, 4,5 ) представлены в таблицах 2.3 и 2.4. В приложении A приведен пример при n = 5.

0,0 25,00 – 5,904624 2,837772 – 0,583105 0,090419 – 0, 0,1 24,58 – 4,179480 1,126703 – 0,049648 – 0,014480 0, 0,2 24,12 – 2,567052 – 0,123462 0,191588 – 0,031421 0, 0,3 23,62 – 1,067341 – 0,943584 0,225923 – 0,010841 – 0, 0,7 21,22 3,804335 – 0,540820 – 0,213391 – 0,012818 0, 0,8 20,52 4,740462 0,326387 – 0,147810 – 0,029362 – 0, Анализируя таблицы, видим, что наилучшее приближение к точному решению дает пробное решение y 5 (x ), для которого 3.2. Воспользовавшись указаниями из раздела 2.2, в качестве поверочных функций возьмем многочлены Лежандра В результате расчета по программе при n = 5 получим вектор коэффициентов С = (1,135995 2,510843 2, 638116 0,080296 1, 220556), следовательно, пробное решение имеет вид Основные результаты расчета при n 5 представлены в таблицах 2.5 и 2.6.

В приложении A приведен пример при n = 5.

0,0 25,00 – 7,542857 3,206936 – 0,578802 0,080526 – 0, 0,1 24,58 – 5,704000 1,374150 – 0,034754 – 0,019452 0, 0,2 24,12 – 3,981714 0,021919 0,210140 – 0,031797 0, 0,3 23,62 – 2,376000 – 0,881665 0,242606 – 0,007511 – 0, 0,4 23,08 – 0,886857 – 1,368509 0,143171 0,018995 – 0, 0,5 22,50 0,485714 – 1,470520 – 0,013838 0,028644 0, 0,6 21,88 1,741714 – 1,219607 – 0,160293 0,016448 0, 0,7 21,22 2,881143 – 0,647676 – 0,234266 – 0,009752 0, 0,8 20,52 3,904000 0,213364 – 0,180028 – 0,030533 – 0, 0,9 19,78 4,810286 1,331607 0,051950 – 0,016143 – 0, Анализируя решения при n = 0,1,...,5, видим, что наилучшее приближение к точному решению дает пробное решение y 5 (x ), для которого 4. Анализируя полученные погрешности, видим, что для задачи (2.17) многочлены Лежандра, как поверочные функции, дают лучшее приближение решения этой краевой задачи.

1. Найдите решение краевой задачи (2.17) аналитическим методом.

2. Каковы отличия краевой задачи от задачи Коши?

3. Каким условиям должны удовлетворять пробные функции в методе Галеркина?

4. Как находится функция, названная в методе Галеркина невязкой пробного решения?

5. Какими свойствами должны обладать поверочные функции в методе Галеркина?

6. Как в методе Галеркина строится система линейных алгебраических уравнений для определения коэффициентов пробного решения? Проверьте истиннос ть формул (2.8), (2.9).

7. В каком случае невязка пробного решения сходится к нулю в среднем 8. Опишите алгоритм приближенного решения краевой задачи для линейного дифференциального уравнения второго порядка методом Галеркина.

9. Приведите пример пос троения пробных функций методом неопределенных коэффициентов.

10.Напишите два многочлена Лежандра, ортогональные [2,4], и проверьте их ортогональнос ть.

11.Напишите уравнение и краевые условия задачи на собственные значения.

12.Напишите две собственные функции задачи (2.13), (2.14) и проверьте их ортогональнос ть.

13.Приведите физические интерпретации изучаемой краевой задачи.

3. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка Рассмотрим краевую задачу в следующей постановке.

Требуется на отрезке [a, b] найти решение Y (x ) дифференциального уравнения удовлетворяющее двум краевым (или граничным) условиям где K (x), K (x), (x ), g (x ) – заданные непрерывные на [a, b] функции ( K ( x) 0 ); a 0, a1, a 2, b0, b1, b2 – заданные действительные числа, причем Заметим, что краевая задача (2.1), (2.2) может быть сведена к задаче (3.1), (3.2) после умножения уравнения (2.1) на положительный множитель Идея вариационного метода состоит в замене краевой задачи (3.1), (3.2) равносильной задачей об отыскании дважды непрерывно дифференцируемой на [a, b] функции Y (x ), доставляющей экстремум следующему функционалу причем значения параметров a, b, q a, qb, Ta, Tb в этом функционале определяются в зависимости от значений a0, a1, a 2, b0, b1, b2 по таблице 3.1.

В методе Ритца для нахождения приближенного решения краевой задачи (3.1), (3.2) строится функциональная последовательнос ть {y n (x )} из пробных решений y n (x ) следующим образом.

Как и в методе Галеркина, задаемся на [a, b] функцией u 0 ( x ) и пробными функциями u1 (x ),..., u n ( x), такими, что u 0 ( x ) удовлетворяет условиям (3.2), а u1 (x ),..., u n ( x) удовлетворяют однородным условиям (2.3), и составляем функцию где С i ( i = 1, n ) – некоторые постоянные. Значения постоянных С i (i = 1, n ) подберем так, чтобы функция (3.6) доставляла экстремум функционалу (3.5).

Подс тавляя y ( x) = y n ( x) в (3.5), получаем квадратичную функцию переменных C 1,..., C n Необходимые условия экстремума функции (3.7), как известно из математического анализа, имеют вид:

Записав условия (3.8) в развернутом виде, для определения значений переменных C1, C 2,..., Cn получаем неоднородную систему линейных алгебраических уравнений n -го порядка где Решив систему (3.10) и подставив определяемые этим решением значения постоянных C1, C 2,..., Cn в (3.6), завершаем построение пробного решения Опишем теперь возможный алгоритм приближенного решения задачи (3.1), (3.2) методом Ритца, предполагая, что {y n (x)}0 сходится к Y (x ) при 1. Подготовительный шаг алгоритма. На этом шаге определяем значения параметров функционала (3.5) в соответствии с таблицей 3.1.

Выбираем функции u 0 ( x), u1 ( x),..., un (x ) и находим функцию т. е. невязку от подстановки u 0 ( x) в уравнение (3.1). Если x [a, b] : R 0 ( x) = 0, то u 0 ( x) = y 0 ( x) – искомое решение и вычисления заканчиваем. Если же R 0 ( x) 0, то переходим к следующему шагу алгоритма.

2. Первый шаг алгоритма. Строим функцию y1 ( x ) = u0 ( x ) + C1u1 ( x), определив значение C1 из решения системы (3.9) при n = 1.

Находим невязку Если x [a, b] : R (C1, x) = 0, то Y (x ) = y1 (x ), и задача решена.

Если R (C1, x ) 0 на [a, b], то находим или приближенного решения, то полагаем Y (x ) y1 (x ) и вычисления заканчиваем.

Если же 11 1 или 12 2, то переходим к вычислениям на следующем шаге.

Таким образом, на m -м шаге ( m 1 ) алгоритма сначала строим функцию определив значения C1, C 2,..., C m из решения системы (3.9) при n = m, а затем находим невязку переходим к (m + 1) -му шагу, и т. д.

Некоторые методы подбора пробных функций были приведены в разделе 2.2.

Подчеркнем здесь, что если пробные функции выбираются на множестве многочленов, то их всегда можно найти методом неопределенных коэффициентов, причем неоднозначно. Одним из возможных наборов пробных функций u i (x ) ( i 1 ) будут многочлены или где Напомним также, что пробные функции u i (x ), i 1 можно выбрать из собственных функций соответствующей задачи на собственные значения, представляемых в виде Приведем теперь еще несколько примеров построения пробных функций для некоторых вариантов граничных условий (3.2).

Пример 1. Построить u 0 ( x ) и систему пробных функций для задачи (3.1) с краевыми условиями y (0) = 2, y (1) = 3.

Решение. Пусть u 0 = Bx + Cx 2. Тогда из граничных условий находим Функции u k (x), k 1 будем искать в виде (3.13) sin( k x) = 0. Отсюда имеем k = k, 2 = ( k ) 2, k = 1,2,.... Тогда, обозначая C k = B k, находим u k ( x) = C k cos( kx). Таким образом пробное решение можно искать в виде Если же функции u k (x) ( k 1) взять в виде (3.11), то если же в виде (3.12), то Пример 2. Построить u 0 ( x ) и систему пробных функций для краевой задачи с условиями y (0) = 2, y (1) = 4.

Решение. Положим u 0 ( x ) = A + Bx. Удовлетворяя граничным условиям, Зададим Согласно однородным граничным условиям имеем Вводя обозначение C k = Ak, получаем Тогда Используя же пробные функции вида (3.11), (3.12), получаем или Пример 3. Построить u 0 ( x ) и систему пробных функций для задачи с граничными условиями 2 y (0) = y (0) 1, y (1) = 3.

Решение. Ищем искомые функции в виде u 0 ( x ) = A + Bx, Потребуем, чтобы функция u 0 ( x ) удовлетворяла неоднородным граничным условиям. Тогда 2B = A 1, A + B = 3, т. е. A = 7 / 3, B = 2 / 3, следовательно, Из соответс твующих однородных условий находим Откуда получаем 2 k Bk = Ak и трансцендентное уравнение для определения собственных значений.

Последнее уравнение имеет счетное множес тво дейс твительных корней 1, 2,..., что подтверждает рисунок 3.1.

Корни уравнения (3.14) определяются приближенными численными методами, например, метод хорд, метод Ньютона, методом итерации или методом половинного деления. В системе MathCAD корни уравнений отыскиваются с помощью стандартной функции root (см. п. 8.2).

Таким образом, обозначая C k = Ak, имеем тогда Рис. 3.1. Геометрическая иллюстрация корней уравнения (3.14) Если использовать пробные функции вида (3.11), (3.12), то получаем или Методом Ритца найти наиболее точное приближенное решение краевой задачи (2.16), построенное при помощи системы из n пробных функций – многочленов и системы из n пробных функций вида (3.13).

За меру точности выбрать (по указанию преподавателя) Варианты заданий приведены в таблице 2.1.

Лабораторная работа выполняется с использованием прикладной системы MathCAD, в которой реализуется алгоритм построения пробных решений y m (x) методом Ритца.

Перед обращением к программе необходимо подготовить числовые и строчные данные, вводимые в процессе диалога с клавиатуры дисплея.

Числовые данные:

a, b – концы отрезка интегрирования [a, b ] ;

n – максимальное число параметров C1,..., C n в пробном решении.

Значение параметра n задает преподаватель.

Строчные данные:

аналитические выражения для пробных функций u 0 ( x ),..., um ( x).

В результате расчета программа выводит на экран дисплея значения коэффициентов C1,..., C n и таблицы всех пробных решений y1 (x ),..., y n ( x) и их невязок. Анализируя данные этих таблиц, надо найти обоснованный ответ на поставленную задачу лабораторной работы.

3.4. Выполнение работы в компьютерном классе 1. Прежде чем начать выполнение лабораторной работы на ЭВМ, внимательно ознакомьтесь с данной инструкцией.

2. При необходимости включите сами (или попросите лаборанта) питание компьютера. После того, как система загрузится, запускаем двойным щелчком левой кнопки мыши на рабочем столе программу Mathcad, если же ярлык отсутс твует, тогда открываем программу через кнопку «Пуск» (Программы Mathsoft Mathcad).

3. Узнайте у лаборанта расположение файла ODU.mcd и откройте его (File Open или, если программа русифицирована, Файл Открыть). При любой ошибке ввода программы нужно обратиться к лаборанту.

4. Прочитайте в начале файла задание на лабораторную работу и просмотрите пример выполнения работы, для которого исследование уже проведено. Программа файла ODU.mcd состоит из шести пунктов «1.

Постановка задачи», «2. Получение точного решения», «3. Получение приближенного решения методом Галеркина», «4. Получение приближенного решения вариационным методом Ритца», «5. Получение приближенного решения интегральным методом наименьших квадратов», «6. Выводы». Для выполнения лабораторной работы «Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка вариационным методом Ритца» необходимо использовать пункты 1, 2 и 4. Цели и задачи каждого из пунктов описаны ниже.

5. Для набора функций нужно либо воспользоваться всплывающим меню инструментов «Calculator», либо ввести ее с клавиатуры, используя следующие символы арифметических действий и стандартных функций: сложение – ‘+’;

вычитание – ‘–‘; умножение – ‘*’; деление – ‘/’; возведение в степень – ‘^’;

квадратный корень – ‘\’; синус – sin(x); косинус – cos(x); экспонента – exp(x);

натуральный логарифм – ln(x). При вводе числовых данных, являющихся десятичными дробями, целую и дробную части нужно разделять точкой (например, 0.5, 1.5 и т. д.).

6. Порядок выполнения работы Вам укажет программа подсказками и заданиями, выделенными красным цветом.

3.5. Порядок выполнения лабораторной работы Рекомендуется следующий порядок выполнения лабораторной работы.

1. Изучить разделы 3.1–3.4 настоящей главы и подготовить ответы на контрольные вопросы из раздела 3.7.

2. Пройти собеседование с преподавателем, получить допуск к выполнению работы на ПЭВМ, номер варианта задания лабораторной работы и значение параметра n.

3. В соответствии с полученным задания выполнить подготовительный шаг алгоритма метода Ритца и подготовить, если u 0 ( x) не является точным решением задачи, числовые и строчные исходные данные для расчетов на ПЭВМ.

4. Выполнить основную расчетную часть лабораторной работы в диалоге с ПЭВМ. В процессе диалога следует переписать с экрана дисплея значения коэффициентов C i пробных решений, а в конце диалога – итоговые таблицы пробных решений и их невязок.

5. Оформить и защитить отчет по лабораторной работе, который должен содержать титульный лист, математическую постановку задачи, результаты выполнения подготовительных расчетов, основные результаты расчетов на ПЭВМ и обоснованный итоговый отчет.

Методом Ритца найти на [0, 1] приближенное значение краевой задачи (2.17).

Сводим задачу (2.17) к задаче (3.1), (3.2), определяя, согласно (3.4), Получаем задачу где K (x) = 3e 3 x ; (x) = 2e3 x; g(x) = 2x2 6x + 2 e3x. Напомним, что задача (3.16) имеет точное решение (2.18), значения которого представлены в таблице 2.2.

Для выполнения лабораторной работы используется тот же файл ODU.mcd, что и в разделе 2.5. В пункте «Постановка задачи» вводим числовые параметры и функции, входящие в задачу (2.17). В пункте «Получение точного решения в системе MathCAD» получаем таблицу 2.2 точного решения задачи.

В пункте «Получение приближенного решения вариационным методом Ритца» в качестве пробных функций используем те же функции, что и в разделе 2.5:

причем R (u 0 ) 0.

Определяем параметры функционала (3.5). Так как a0 = a1 = b0 = b1 = 1, то в соответс твии с таблицей 3.1 имеем В результате расчета по программе ODU.mcd при n = 5 получим вектор коэффициентов С = (1,140938 2,564241 2,466128 0,133009 1,130778), следовательно, пробное решение имеет вид Основные результаты расчета по программе ODU.mcd при n представлены в таблицах 3.2 и 3.3. В приложении A приведен пример при n = 5.

Анализируя таблицы, устанавливаем, что наилучшее приближение к точному решению Y (x ) дает решение y 5 (x ), для которого 0,0 6,0 1,041255 0,835765 0,847800 0,846705 0, 0,1 5,5 1,020600 0,860600 0,865221 0,865260 0, 1. Опишите алгоритм решения краевой задачи для линейного дифференциального уравнения второго порядка аналитическим методом.

2. Каким образом уравнение (2.1) свести к равносильному уравнению типа (3.1)?

3. В чем основная идея вариационного подхода к решению краевой задачи (3.1), (3.2)?

4. Проверьте правильность данных, представленных в таблице 3.1.

5. Какими свойствами должны обладать пробные функции в методе Ритца?

6. Как в методе Ритца находится невязка пробного решения?

7. Докажите, что ортогональная на [a, b] система функций, среди которых нет тождес твенно равной нулю, линейно независима.

8. Как в методе Ритца строится система алгебраических уравнений для определения коэффициентов пробного решения? Проверьте справедливость соотношений (3.9), (3.10).

9. Опишите алгоритм приближенного решения краевой задачи (3.1),(3.2) методом Ритца.

10.Приведите пример пробных функций для решения задачи методом Ритца.

11.Проверьте, что функции (3.13) являются собственными функциями задачи (2.13).

12.Опишите алгоритм метода хорд-касательных для приближенного вычисления корней уравнения f (x ) = 0.

13.Опишите алгоритм метода половинного деления для приближенного вычисления корней уравнения f (x ) = 0.

14.Опишите алгоритм метода итераций для приближенного вычисления корней уравнения f (x ) = 0.

4. Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка интегральным Снова рассмотрим краевую задачу: найти на отрезке [a, b ] решение Y (x ) дифференциального уравнения удовлетворяющее условиям где p( x), q( x), f (x ) – заданные функции, непрерывные на отрезке [a, b ] ; a0, a1, a2, b0, b1, b2 – заданные числа, причем a0 + a12 0, b02 + b12 0.

Заметим здесь, что краевая задача (3.1), (3.2) может быть сведена к задаче (4.1), (4.2). Для этого достаточно разделить обе части уравнения (3.1) на K (x) и ввести обозначение Для нахождения приближенного решения задачи (4.1), (4.2) интегральным последовательность {y n (x )}0 из пробных решений вида где u 0 ( x), u1 ( x),..., un (x ) – функции, удовлетворяющие таким же условиям и требованиям, что и аналогичные функции в методах Галеркина и Ритца.

Подс тавляя пробное решение (4.3) вместо y (x) в уравнение (4.1), получим невязку Напомним, что функция (4.4), линейно зависящие от параметров C1,..., C n, является характеристикой уклонения пробного решения (4.3) от точного решения задачи Y (x ). Поэтому подберем значения C1,..., C n так, чтобы они доставляли глобальный минимум следующей функции переменных C1,..., C n Заметим, что, так как (C1,..., C n ) из (4.5) неотрицательная квадратичная функция n переменных, то глобальный минимум ее существует и совпадает с локальным.

Необходимое условие локального минимума функции (4.5) дают откуда Записав условия (4.6) в развернутом виде, для определения значений переменных С1,..., C n получаем неоднородную систему уравнений n -го порядка где Решив систему (4.7) и подставив определяемые этим решением значения параметров С1,..., C n в (4.4), завершаем построение пробного решения y n (x).

Этапы возможного алгоритма приближенного решения задачи (4.1), (4.2) интегральным методом наименьших квадратов качес твенно полностью совпадают с этапами алгоритма решения задачи методом Галеркина. Имеется только одно количес твенное различие, связанное с тем, что параметры С1,..., C n пробного решения на первом и последующих этапах определяются решением системы (4.7), а не системы (2.8), как было в методе Галеркина.

Подчеркнем, что пробные функции можно подбирать так же, как и в методах Галеркина и Ритца.

Интегральным методом наименьших квадратов найти наиболее точное решение краевой задачи (2.16), построенное при помощи системы из n пробных функций многочленов. За меру точнос ти выбрать (по указанию преподавателя) или 1, или 2 из (3.15). Варианты заданий приведены в таблице 2.1.

Лабораторная работа выполняется с использованием прикладной системы MathCAD, в которой реализуется алгоритм построения пробных решений y m (x) интегральным методом наименьших квадратов.

Перед обращением к этой программе необходимо подготовить следующие числовые и строчные данные.

Числовые данные:

a, b – концы отрезка интегрирования;

n – максимальное число параметров С i в пробном решении.

Значение параметра n дает преподаватель.

Строчные данные:

аналитические выражения для функций u 0 ( x),..., u n ( x ).

4.3. Выполнение работы в компьютерном классе 1. Прежде чем начать выполнение лабораторной работы на ЭВМ, внимательно ознакомьтесь с данной инструкцией.

2. При необходимости включите сами (или попросите лаборанта) питание компьютера. После того, как система загрузится, запускаем двойным щелчком левой кнопки мыши на рабочем столе программу Mathcad, если же ярлык отсутс твует, тогда открываем программу через кнопку «Пуск» (Программы Mathsoft Mathcad).

3. Узнайте у лаборанта расположение файла ODU.mcd и откройте его (File Open или, если программа русифицирована, Файл Открыть). При любой ошибке ввода программы нужно обратиться к лаборанту.

4. Прочитайте в начале файла задание на лабораторную работу и просмотрите пример выполнения работы, для которого исследование уже проведено. Для выполнения лабораторной работы «Решение краевой задачи для линейного обыкновенного дифференциального уравнения второго порядка интегральным методом наименьших квадратов» необходимо использовать пункты 1, 2 и 5, 6. Цели и задачи каждого из пунктов описаны ниже.

5. Для набора функций нужно либо воспользоваться всплывающим меню инструментов «Calculator», либо ввести ее с клавиатуры, используя следующие символы арифметических действий и стандартных функций: сложение – ‘+’;

вычитание – ‘–‘; умножение – ‘*’; деление – ‘/’; возведение в степень – ‘^’;

квадратный корень – ‘\’; синус – sin(x); косинус – cos(x); экспонента – exp(x);

натуральный логарифм – ln(x). При вводе числовых данных, являющихся десятичными дробями, целую и дробную части нужно разделять точкой (например, 0.5, 1.5 и т. д.).

6. Порядок выполнения работы Вам укажет программа подсказками и заданиями, выделенными красным цветом.

4.4. Порядок выполнения лабораторной работы Рекомендуется такой порядок выполнения лабораторной работы.

1. Изучить разделы 4.1–4.3 настоящей главы, повторить разделы 2.2, 3.2 и подготовить ответы на контрольные вопросы из раздела 4.6.

2. Пройти собеседование с преподавателем; получить допуск к выполнению работы в диалоге ПЭВМ, номер варианта задания и значение параметра n .

3. Выполнить подготовительный шаг алгоритма метода наименьших квадратов и подготовить, если u 0( x) не является точным решение задачи, все числовые и строчные исходные данные для расчетов на ПЭВМ.

4. Выполнить расчетную часть работы на ПЭВМ, переписав с экрана дисплея значения коэффициентов C i пробных решений и итоговые таблицы пробных решений и их невязок.

5. Оформить и защитить отчет по лабораторной работе такого же содержания, что и отчеты по работам предыдущих глав.

Интегральным методом наименьших квадратов найти приближенное решение краевой задачи (2.17), используя пробные функции из раздела (2.5).

Для выполнения лабораторной работы используется тот же файл ODU.mcd, что и в разделах 2.5 и 3.5. В пункте «Пос тановка задачи» вводим числовые параметры и функции, входящие в задачу (2.17). В пункте «Получение точного решения в системе MathCAD» получаем таблицу 2.2 точного решения задачи.

В пункте «Получение приближенного решения интегральным методом наименьших квадратов» в качестве пробных функций используем те же функции, что и в разделе 2.5:

В результате расчета по программе при n = 5 получим вектор коэффициентов С = (1,137789 2,526360 2,595082 0, 032496 1,202066), следовательно, пробное решение имеет вид Основные результаты расчета по программе ODU.mcd при n представлены в таблицах 4.1 и 4.2. В приложении A приведен пример при n = 5.

Наилучшее приближенное к точному решению дает пробное решение y 5 (x ), для которого 0,3 23,62 – 1,340378 – 0,961014 0,195123 – 0,002561 – 0, 1. Как в методе наименьших квадратов строится система линейных уравнений для определения параметров пробного решения?

2. Получите самостоятельно развернутый вид условий (4.6), проверив тем самым справедливость формул (4.7), (4.8).

3. Опишите алгоритм приближенного решения краевой задачи (4.1), (4.2) интегральным методом наименьших квадратов.

4. Приведите примеры пробных функций для решения задачи (4.1), (4.2) интегральным методом наименьших квадратов.

5. Решение начально-краевой задачи для одномерного параболического уравнения методом Галеркина.

Рассмотрим следующую начально-краевую задачу. Требуется в плоской области найти решение U (x, t ) дифференциального уравнения удовлетворяющее двум краевым или граничным условиям и начальному условию где K ( x, t), K (x, t ), ( x, t ), g ( x, t ), a 2 (t ), b2 (t ) – заданные, непрерывные на D функции (K ( x, t ) 0); a0, a1, b0, b1 – заданные действительные числа, причем a0 + a12 0, b0 + b12 0 ; f (x ) – заданная функция, непрерывная на [a, b] вместе с f (x) и такая, что Напомним, что в такой форме может быть поставлена задача одномерной нестационарной теплопроводнос ти [1]. Например, типичная задача о нестационарной теплопередаче путем теплопроводности в однородном стержне единичной длины, концы которого поддерживаются при температурах T1 и T2, при начальном распределении температуры вдоль стержня по закону получается как частный случай сформулированной задачи В методе Галеркина для нахождения приближенного решения задачи (5.1)– (5.4) строится функциональная последовательнос ть {un (x, t )} из пробных решений u n ( x, t ) следующим образом.

Задаемся в области D некоторой системой дважды дифференцируемых функций u 0 ( x, t), u1 ( x),..., un (x ) таких, что u 0 ( x, t) удовлетворяет краевым условиям (5.2), а пробные функции u i (x) (i 1) являются линейно независимы на [a, b] и удовлетворяют однородным краевым условиям Составляем функцию с неизвестными пока функциями v1 (t ), v2 (t ),..., vn (t ), зависящими только от аргумента t.

Подчеркнем, что в силу линейности условий (5.2) и (5.6), функция (5.7) удовлетворяет условиям (5.2) при любых функциях v1 (t ),...,v n (t ). Значит, следует так определить vi (t ) (i 1) и количество (n) этих функций, чтобы u n ( x, t ) из (5.7) удовлетворяла уравнению (5.1) и начальному условию (5.3) с заданной точностью.

Подставляя u n ( x, t ) вместо u ( x, t) в уравнение (5.1), получаем невязку или Подставляя u 0 ( x,0), в полученную из (5.7) при t = 0, в (5.3), находим невязку Невязки R1 и R2 являются характеристиками уклонения функции (5.7) от точного решения U (x,t ) задачи (5.1)–(5.4). Во всяком случае, если при некотором наборе функций v1 (t ),...,v n (t ) R1 0 и R2 0, то функция u n ( x, t ) из (5.7) – точное решение U (x,t ).

В общем случае эти невязки оказываются отличными от нуля. Поэтому накладываем дополнительные условия на функции vk (t ) и их начальные значения vk (0) так, чтобы невязки в каком-то смысле были бы наименьшими.

В обобщенном методе Галеркина эти условия определяются системой уравнений:

где w1 (x ),..., wn ( x) – заданные линейно независимые на [a, b] поверочные функции; а Напомним здесь, что если поверочные функции w1 (x ),..., wn ( x) входят в полную на [a, b] систему функций, то можно ожидать сходимости последовательности {un (x,t )} в среднем к точному решению U (x,t ) [1].

Запишем условия (5.10) в развернутом виде или или где Если ввести в рассмотрение матрицы то система (5.12) в матричном виде запишется так Покажем, что матрица A всегда невырожденная, т. е. det A 0.

Рассмотрим однородную линейную алгебраическую систему уравнений относительно неизвестных 1, 2,..., n Если det A = 0, то система (5.17) имеет множество ненулевых решений.

Пусть одним из таких решений является совокупность 1, 2,..., n, где, например, m 0. Подставляя это решение в уравнение системы (5.17), суммируя все получившиеся при этом равенства и используя свойства скалярного произведения, получаем Так как функции wk (x ) линейно независимы, то w1 +... + wn 0. Значит,/ должно выполнятся тождество 1u1 +... + nu n 0, m 0. Но это невозможно из-за линейной независимости функций u1,..., u n. Значит, ненулевых решений у системы (5.17) нет, а для этого необходимо и достаточно, чтобы det A 0.

Таким образом, матрица A невырожденная и, следовательно, имеет обратную матрицу A1.

Теперь из (5.16) получаем Таким образом, функции v j (t ) должны удовлетворять нормальной системе линейных обыкновенных дифференциальных уравнений n -го порядка.

Заметим, что если функции K ( x, t), ( x,t ) зависят только от x, то система (5.18) – система с постоянными коэффициентами. Заметим так же, что если в качестве поверочных функций выбраны пробные, которые ортогональны, то матрицы A и A1 являются диагональными матрицами.

Запишем теперь в развернутом виде условия (5.11). Получаем или или где akj определяются формулами (5.13), а Если ввести матрицу D = (d k )n,1, то из (5.19) получаем Таким образом, для нахождения функций Vk (t ), k = 1, n, определяющих пробное решение (5.7), получаем задачу Коши для нормальной системы (5.18) линейных обыкновенных дифференциальных уравнений n -го порядка с начальными условиями (5.20). Решив указанную задачу Коши и подставив определяемые этим решением функции v k (t) в (5.7), заканчиваем построение пробного решения u n ( x,t ).

Опишем возможный алгоритм построения приближенного решения задачи (5.1)–(5.3) методом Галеркина, предполагая, что последовательность {un (x,t )} сходится равномерно к точному решению U (x,t ).

1. Подготовительный шаг алгоритма. На этом шаге выбираем функцию u 0 ( x, t) и находим невязку R10 ( x,t ) = L[u0 ] g ( x, t ) от подстановки функции u 0 ( x, t ) в уравнение (5.1). Находим невязку R20 ( x ) = u0 ( x,0) f ( x) для условия (5.3). Определяем Если 10 1 и 20 2, где 1 и 2 заданные меры точности приближенного решения, то полагаем U (x,t ) u0 (x,t ). В противном случае переходим к следующему шагу алгоритма, предварительно выбрав u j (x ) и поверочные wk (x ) функции.

2. Первый шаг алгоритма. Определив функцию v1 (t ) из решения задачи Коши (5.18), (5.20) при n = 1, строим функции u1 (x,t ) = u0 + v1(t )u1( x). Находим по формулам (5.8), (5.9) невязки R11 (v1 (t), x,t ), R21 (v1(0), x ) и определяем max R11 (v1, x, t ) = 11 и max R21(v1(0), x) = 21. Если 11 1 и 21 2, то полагаем U (x,t ) ~ u1( x, t), и вычисления заканчиваем. В противном случае переходим к вычислениям на втором шаге алгоритма и т. д.

Таким образом, на m-ом (m 1) шаге алгоритма строим функцию определив предварительно функции v1 (t ),...,v m (t ) из решения задачи Коши (5.18), (5.20) при n = m. Находим по формулам (5.8), (5.9) невязки а затем вычисляем max R1m (v1 (t ),...,vm (t ), x, t ) = 1m и max R2m (v1(0),...,vm (0), x) = 2m.

Если 1m 1 и 2m 2, то полагаем U (x,t ) u m ( x,t ), в противном случае переходим к (m + 1) -му шагу алгоритма.

Пробные и поверочные функции можно выбирать так же или такими же методами, как описано в предыдущих главах.

Поэтому обсудим здесь только возможность построения функции u 0 ( x, t) в виде многочлена относительно x с коэффициентами, зависящими от t, и рассмотрим несколько примеров, иллюстрирующих эту возможность.

Например, положив u 0 ( x, t) = A(t ), из условий (5.2) получаем систему функциональных уравнений и если a0b2 b0a2, то система совместна и A(t ) = 2. Если же a0b2 b0 a 2, то система несовместна, и ищем u 0 ( x, t) в виде Для определения A(t ) и B(t ) из условий (5.2) получаем систему функциональных уравнений которую можно исследовать, используя теорему Кронекера-Капелли, как линейную неоднородную алгебраическую систему относительно неизвестных функций A(t ) и B(t ).

Если 1 = a0 (b0b + b1 ) b0 (a 0a + a1 ) 0, то система совместна и определена при этом и функция u 0 ( x, t) = P ( x,t ) определяется однозначно. Если 1 = 0, то система несовместна, и ищем u 0 (x, t ) в виде Для определения A(t ) и B(t ) из условий (5.2) получаем систему неопределена, причем B(t ) можно придавать произвольные значения. Если 2 = 0, то система несовместна, и ищем u 0 ( x, t) в виде Условия (5.2) приводят к системе b0 A(t ) + (b0b + b1 ) B(t ) + (b0b 2 + 2b1b) C(t ) + (b0b3 + 3b1b 2 ) D(t ) = b2 (t).

Покажем, что это система всегда совместна и, следовательно, неопределена.

Для этого надо доказать, что для любых значений параметров a,b, a0,b0, a1,b1 не могут выполняться все условия несовместности, отмеченные выше, одновременно Введем обозначения x1 = a0b0, x2 = a0b1, x3 = b0a1 и заметим, что, в силу ограничений на параметры (a0 + a12 0, b02 + b12 0) и последнего из выписанных условий несовместности, переменные x1, x 2 и x 3 одновременно в ноль обратиться не могут. Тогда первые три условия несовместности можно записать в виде линейной однородной системы относительно x1, x 2 и x которая должна иметь ненулевое решение. Для этого необходимо и достаточно, чтобы определитель третьего порядка Последнее невозможно, так как 3 = (b a )4 и a b.

Таким образом, при любых значениях параметров a, b, a0, a1,b0,b1 всегда найдется хотя бы одна функция вида u 0 ( x,t ) = P3 (x,t ), удовлетворяющая условиям (5.2).

Пример 1. Построить u 0 ( x, t) для задачи с краевыми условиями Решение. Пусть u 0 ( x, t) = A(t ), тогда условия (5.21) дают т. е. – несовместную систему. Если теперь u 0 ( x, t) = A(t ) + B(t ) x, то условия (5.21) приводят к системе Следовательно, в качестве функции u 0 ( x, t) можно взять функцию Пример 2. Построить функцию u 0 ( x, t) для задачи с краевыми условиями Решение. Пусть u 0 ( x, t) = A(t ), тогда условия (5.22) дают т. е. – несовместную систему. Если теперь u 0 ( x, t) = A(t ) + B(t ) x, то условия (5.22) приводят также к несовместной системе Полагая u 0 ( x, t) = A(t ) + B(t ) x + С(t) x 2, снова получаем несовместную систему (5.23). Ищем поэтому u 0 ( x, t) в виде u 0 ( x, t) = A(t ) + B(t ) x + С(t) x 2 + D(t )x 3. Из условий (5.22) имеем Получили совместную систему. Одним из решений ее будет, например, следующая совокупность функций A(t) =et, B(t) = 0, C(t) = 0, D(t) = 0,5e2t + 0,25et.

Таким образом, u 0 ( x, t) = e t + (0,25e t 0,5e 2t )x 3.

Пример 3. Построить u 0 ( x, t) и систему из пяти пробных функций для задачи с краевыми условиями Решение. Если u 0 ( x) = A, то получаем из (5.24) несовместимую систему Если u 0 ( x) = A + Bx, то условия задачи (5.24) дают Таким образом, В качестве пробных функций можно взять Рассматривается начально-краевая задача. Требуется в плоской области найти решение u ( x, t) дифференциального уравнения удовлетворяющее условиям где c1,c 2, c3, c4 – некоторые заданные постоянные величины.

Заметим, что эта задача получается как частный случай задачи (5.1)–(5.3) b1 = 1, b2 = c3.

Варианты заданий, определяемые различными наборами значений постоянных c1,c 2, c3, c4 задачи (5.25)–(5.27) и параметра T, приведенные в таблице 5.1.

Лабораторная работа выполняется с использованием прикладной системы MathCAD, которая реализует алгоритм построения пробных решений u m ( x, t ) задачи (5.25)–(5.27) методом Галеркина.

Перед обращением к программе необходимо подготовить числовые и строчные данные, вводимые в процессе диалога с ЭВМ с клавиатуры дисплея.

Числовые данные:

l – правый конец отрезка изменения переменной x ;

c1 – числовой параметр уравнения (5.25);

T – значение параметра T задачи.



Pages:   || 2 | 3 |
Похожие работы:

«Федеральное агентство по образованию Сыктывкарский лесной институт – филиал государственного образовательного учреждения высшего профессионального образования Санкт-Петербургская государственная лесотехническая академия имени С. М. Кирова КАФЕДРА ФИЗИКИ ФИЗИКА САМОСТОЯТЕЛЬНАЯ РАБОТА СТУДЕНТОВ Методические указания для подготовки дипломированного специалиста по специальностям 280201 Охрана окружающей среды и рациональное использование природных ресурсов, 230201 Информационные системы и...»

«Л.С. Атанасян, В.Т. Базылев Геометрия в двух частях Допущено Министерством образования и науки РФ   в качестве учебного пособия   для студентов физико-математических факультетов   педагогических вузов Часть 2 КНОРУС • МОСКВА • 2013 УДК 514.1(075.8) ББК 22.151.1я73 А92 Рецензенты: Л.Е. Евтушик, д-р физ.-мат. наук, В.И. Близникас, проф. Атанасян Л.С. А92 Геометрия : в 2 ч. — Ч. 2 : учебное пособие / Л.С. Атанасян, В.Т. Базылев. — 2-е изд., стер. — М. : КНОРУС, 2013. — 424 с....»

«ПРИОРИТЕТНЫЙ НАЦИОНАЛЬНЫЙ ПРОЕКТ ОБРАЗОВАНИЕ РОССИЙСКИЙ УНИВЕРСИТЕТ ДРУЖБЫ НАРОДОВ Г.А. ЛУКИЧЕВ, В.М. ФИЛИППОВ СИСТЕМЫ ФИНАНСИРОВАНИЯ ВЫСШЕГО ОБРАЗОВАНИЯ В ЗАРУБЕЖНЫХ СТРАНАХ Учебное пособие Москва 2008 Инновационная образовательная программа Российского университета дружбы народов Создание комплекса инновационных образовательных программ и формирование инновационной образовательной среды, позволяющих эффективно реализовывать государственные интересы РФ через систему экспорта образовательных...»

«Владимирский государственный университет ЛАБОРАТОРНЫЕ РАБОТЫ ПО ДИСЦИПЛИНЕ ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ Методические указания в двух частях Часть 1 Владимир 2004 Министерство образования Российской Федерации Владимирский государственный университет Кафедра технологии переработки пластмасс ЛАБОРАТОРНЫЕ РАБОТЫ ПО ДИСЦИПЛИНЕ ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ Методические указания в двух частях Часть 1 Составитель Н.А. КОЗЛОВ Владимир УДК 678.64 (076.5) Рецензент Кандидат химических наук, доцент...»

«ПРИОРИТЕТНЫЙ НАЦИОНАЛЬНЫЙ ПРОЕКТ ОБРАЗОВАНИЕ РОССИЙСКИЙ УНИВЕРСИТЕТ ДРУЖБЫ НАРОДОВ Л.Е. РОССОВСКИЙ КАЧЕСТВЕННАЯ ТЕОРИЯ ДИФФЕРЕНЦИАЛЬНЫХ И ФУНКЦИОНАЛЬНОДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ Учебное пособие Москва 2008 Инновационная образовательная программа Российского университета дружбы народов Создание комплекса инновационных образовательных программ и формирование инновационной образовательной среды, позволяющих эффективно реализовывать государственные интересы РФ через систему экспорта образовательных...»

«Федеральное агентство по образованию АМУРСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ГОУВПО АмГУ УТВЕРЖДАЮ И.о. зав.кафедрой ТиЭФ Е.А. Ванина _2007г. ИСТОРИЯ И МЕТОДОЛОГИЯ ФИЗИКИ УЧЕБНО-МЕТОДИЧЕСКИЙ КОМПЛЕКС ПО ДИСЦИПЛИНЕ для специальности 010701 – Физика Составитель: Е.А. Ванина Благовещенск 2007 г. Печатается по решению редакционно-издательского совета инженерно-физического факультета Амурского государственного университета Е.А. Ванина Учебно-методический комплекс по дисциплине История и методология...»

«СЫКТЫВКАРСКИЙ ЛЕСНОЙ ИНСТИТУТ КАФЕДРА ВЫСШЕЙ МАТЕМАТИКИ ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА САМОСТОЯТЕЛЬНАЯ РАБОТА СТУДЕНТОВ Методические указания для подготовки дипломированных специалистов по направлению 654700 Информационные системы специальности 230201 Информационные системы и технологии СЫКТЫВКАР 2007 ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ СЫКТЫВКАРСКИЙ ЛЕСНОЙ ИНСТИТУТ – ФИЛИАЛ ГОСУДАРСТВЕННОГО ОБРАЗОВАТЕЛЬНОГО УЧРЕЖДЕНИЯ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ САНКТ-ПЕТЕРБУРГСКАЯ ГОСУДАРСТВЕННАЯ...»

«П ПРАКТИКУМ В ДЛЯ ВУЗОВ ПРАКТИКУМ ПО БИОФИЗИКЕ Учебное пособие для студентов высших учебных заведений Издание второе, исправленное и дополненное Москва 2004 ББК 28.071я73 П69 А в т о р ы: В.Ф. Антонов, А.М. Черныш, В.И. Пасечник, С.А. Вознесенский, Е.К. Козлова Практикум по биофизике: Учеб. пособие для студ. высш. П69 учеб. заведений. — М.: Гуманит. изд. центр ВЛАДОС, 2001. — 352 с. ISBN 5 691 00698 3. Пособие является составной частью учебного комплекта Био физика и служит практическим...»

«ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ РФ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ МИФИ Ю.Н. Громов Пособие по физике Колебания и волны В помощь учащимся 10 – 11 классов Москва 2009 УДК 534.1(075) ББК 22.32я7 Г 87 Громов Ю.Н. УЧЕБНОЕ ПОСОБИЕ ПО ФИЗИКЕ КОЛЕБАНИЯ И ВОЛНЫ. В помощь учащимся 10 – 11 классов. – М.: МИФИ, 2009. – 48 с. Дано систематизированное изложение основного содержания школьного курса физики по разделу Колебания и волны в соответствии с требованиями образовательного...»

«Кафедра физики полупроводников А.В. Бурмистров Проектирование баз данных в MS Access МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ПРАКТИЧЕСКИМ ЗАНЯТИЯМ ПО ДИСЦИПЛИНЕ СИСТЕМЫ УПРАВЛЕНИЯ БАЗАМИ ДАННЫХ Саратов 2012 Содержание 1. Лабораторная работа № 1 3 2. Лабораторная работа № 2 7 3. Лабораторная работа № 3 12 4. Лабораторная работа № 4 14 5. Лабораторная работа № 5 17 6. Лабораторная работа № 6 24 7. Литература 27 Лабораторная работа № 1 1. Построить базу данных, основанную на двух таблицах. Пусть одна таблица...»

«БЕЛОРУССКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ХИМИЧЕСКИЙ ФАКУЛЬТЕТ Кафедра физической химии А. В. Блохин ТЕОРИЯ ЭКСПЕРИМЕНТА Курс лекций В двух частях Часть 1 МИНСК 2002 Автор-составитель Блохин А.В., кандидат химических наук. Рецензенты: кандидат химических наук Н.Н. Горошко; Л.М. Володкович. Утверждено на заседании Ученого совета химического факультета 29 марта 2002 г., протокол № 5. 2 ПРЕДИСЛОВИЕ Учебное пособие представляет собой лекции по курсу Теория эксперимента для студентов IV курса...»

«9435 УДК 519.711; 378.4 ОПЫТ ПРЕПОДАВАНИЯ ТЕОРИИ УПРАВЛЕНИЯ СТУДЕНТАМ ФИЗИЧЕСКОГО ФАКУЛЬТЕТА А.Ю. Ощепков Пермский государственный национальный исследовательский университет Россия, 614990, Пермь, Данщина ул., 19 E-mail: aos57@mail.ru Ключевые слова: система автоматического управления, преподавание теории управления, физические исследования, применение теории управления в физике, Аннотация: В докладе излагается опыт преподавания теории автоматического управления студентам физического факультета...»

«ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ Факультет радиофизики и электроники Кафедра интеллектуальных систем КУРС ЛЕКЦИЙ по специальному курсу Теория принятия решений и распознавания образов Учебное пособие для студентов факультета радиофизики и электроники Минск 2005 1 УДК 681.31:621.38 ББК 32.841я43+32.85я43 ISBN 5-06-0004597 Рецензенты доктор технических наук В. А. Зайка кандидат технических наук, доцент А. А. Белый Рекомендовано Ученым советом факультета радиофизики и электроники 2003 г., протокол №_...»

«ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ Сыктывкарский лесной институт – филиал государственного образовательного учреждения высшего профессионального образования Санкт-Петербургская государственная лесотехническая академия имени С. М. Кирова КАФЕДРА ХИМИИ АНАЛИТИЧЕСКАЯ ХИМИЯ И ФИЗИКО-ХИМИЧЕСКИЕ МЕТОДЫ АНАЛИЗА Раздел Аналитическая химия Методические указания и контрольные задания для студентов специальности 240406 Технология химической переработки древесины заочной формы обучения Самостоятельное...»

«Белорусский государственный университет Химический факультет Кафедра физической химии Л.А.Мечковский Л.М.Володкович Развернутая программа дисциплины “Физическая химия” с контрольными вопросами и заданиями Учебно-методическое пособие для студентов химического факультета специальности Н 03.01.00—химия Минск 2004 1 УДК. ББК. Рецензенты Кандидат химических наук доцент Г.С. Петров Кандидат химических наук доцент А.Ф. Полуян Мечковский Л.А., Володкович Л.М. Развернутая программа дисциплины...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ САМАРСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ Межфакультетская кафедра истории отечества МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ИЗУЧЕНИЮ КУРСА “ОТЕЧЕСТВЕННАЯ ИСТОРИЯ” Издательство “Самарский университет” 2003 Печатается по решению Редакционно-издательского совета Самарского государственного университета Методические указания содержат программу, планы семинарских занятий, тематику контрольных работ, список литературы и рекомендации по работе над материалами курса....»

«ИНСТИТУТ ПРИКЛАДНОЙ МАТЕМАТИКИ ИМ. М.В.КЕЛДЫША РАН МОСКОВСКИЙ ФИЗИКО-ТЕХНИЧЕСКИЙ ИНСТИТУТ Методическое пособие к курсу МЕТОДЫ НАВИГАЦИИ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ (второй семестр) Профессор А.К.Платонов Аспирант Д.С. Иванов Москва 2013 г. Пособие разработано в процессе чтения лекций на кафедре МФТИ Прикладная математика по специализации Управление динамическими системами, направленных на подготовку студентов-магистров. Цель курса – освоение студентами фундаментальных знаний в области...»

«Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования Владимирский государственный университет Кафедра технологии переработки пластмасс ЛАБОРАТОРНЫЕ РАБОТЫ ПО ДИСЦИПЛИНЕ ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ Методические указания В двух частях Часть 2 Составитель Н.А. КОЗЛОВ Владимир 2006 1 УДК 678.64 (076.5) ББК 32.81 Л12 Рецензент Кандидат химических наук, доцент Владимирского государственного университета М.В. Ольшевский Печатается по...»

«Ho IL М А К А Р К И Н И. М. ШАРАНОВ Н. Ф. Д Ю Д Я Е В в ; Ф. Б А Й Н Е В ИНФОРМАТИКА УЧЕБНОЕ ПОСОБИЕ Рекомендовано Министерством общего и профессионального образования Российской Федерации в качестве учебного пособия для студентов высших учебных заведений, обучающихся по экономическим и гуманитарным специальностям САРАНСК ИЗДАТЕЛЬСТВО МОРДОВСКОГО УНИВЕРСИТЕТА 1998 УДК 338.242 ББК У.ф М151 Рецензенты: кафедра информационно-вычислительных систем Саранского кооперативного института Московского...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ МИФИ Л.Н. ДЕМИНА МЕТОДЫ И СРЕДСТВА ИЗМЕРЕНИЙ, ИСПЫТАНИЙ И КОНТРОЛЯ Рекомендовано УМО Ядерные физика и технологии в качестве учебного пособия для студентов высших учебных заведений Москва 2010 УДК 006.91(075) ББК 30.10я7 Д 30 Демина Л.Н. Методы и средства измерений, испытаний и контроля: Учебное пособие. – М.: НИЯУ МИФИ, 2010. – 292 с. В учебном пособии изложены основные понятия, методы и...»






 
© 2013 www.diss.seluk.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Методички, учебные программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.