ЗАДАЧА 1

advertisement
ЗАДАЧА 1
Приведение к каноническому виду линейных уравнений в частных производных 2-го
порядка в случае двух независимых переменных. Нахождение общего решения уравнения.
Решение задачи Коши для уравнений гиперболического и параболического типов.
Решение задачи приведения к каноническому виду линейных уравнений в частных
производных 2-го порядка в случае двух независимых переменных излагается в учебных
пособиях по уравнениям математической физики, например [1, 2, 3]. Приведенные к
каноническому виду уравнения гиперболического и параболического типов в простейших
случаях решаются. В предлагаемых задачах они решаются интегрированием с
понижением порядка уравнений. Получаются общие решения, содержащие произвольные
функции от канонических переменных. После подстановки выражений канонических
переменных через первоначальные x и y получаем аналитические формулы общих
решений заданные линейных уравнений, в которые входят произвольные функции от
промежуточных аргументов.
При решении задачи Коши в общее решение и его производную по аргументу x или
y подставляются заданные начальные условия. Получается система двух уравнений на две
неизвестные функции. Одно из уравнения дифференциальное. Из этой системы
определяются две функции (на промежуточном этапе - произвольные функции). Заменяя в
общем решении производные функции найденными конкретными функциями, получаем
решение задачи Коши.
Пример 1. Найти решение уравнения
удовлетворяющее заданным начальным условиям
По коэффициентам при вторых производных определяем тип уравнения
Уравнение (1) гиперболического типа. Составляем характеристическое уравнение
Решаем это квадратное уравнение
Находим общие интегралы каждого из полученных дифференциальных уравнений с
разделяющимися переменными
Левые части найденных общих интегралов берут в качестве новых переменных
Новые переменные и называют каноническими, или характеристическими.
Разрешимость уравнений (3) относительно x и у доказана для уравнений
гиперболического типа. Для преобразования заданного уравнения (1) к каноническому
виду вычисляем производные первого и второго порядков по переменным x и у, выражая
их через производимо по переменным
и , рассматриваемым как промежуточные
переменные, с применением правил дифференцирования сложных функций. Получаем
следующие выражения для производных первого порядка:
и, соответственно, для производных второго порядка
Подставляем полученные выражения первых и вторых производных по x и y в уравнение
(1) и группируем слагаемые с общими множителями
и т.д.
Получаем уравнение
Заметим, что коэффициенты при
должны обращаться в нуль. Их можно
и
было бы и не подсчитывать. Но при решении задачи возможны ошибки. Рекомендуется
подсчитывать также и коэффициенты при
. Необращение в нуль
хотя бы одного из этих коэффициентов показывает наличие ошибок. Следует просмотреть
ход решения задачи, отыскать и исправить ошибки.
Разделив все члены уравнения на коэффициент при
вид уравнения гиперболического типа (1)
Это уравнение легко решается. Обозначив
для функции U
получаем канонический
= U, получаем уравнение первого порядка
Интегрируем обе части равенства по переменному
произвольная функция другого переменного η)
(при этом войдет слагаемое -
Потенцируем равенство и заменяем U
Обе части полученного равенства интегрируем по переменному η
Интеграл от произвольной функции φ1(η) есть также произвольная функция
аргумента. Слагаемое
есть произвольная функция другого аргумента
образом, общее решение канонического уравнение
того же
. Таким
Подставляя выражения (3) канонических переменных, получим общее решение исходного
уравнения (1)
Остается выделить из всех решений, представляемых формулой (4),. то, которое
удовлетворяет начальным условиям (2). Дифференцируем функцию Z по переменному y
и подставляем в полученное равенство и равенство (4) y = 0. Получаем два условия
Это система двух дифференциальных уравнений на две неизвестные функции. Система
специфического вида. Первое уравнение содержит только функции и не содержит их
производных. Второе уравнение содержит и Функции и их производные. Дифференцируем первое уравнение
Умножив обе части полученного уравнения на 2 и сложив с соответственными частями
второго уравнения (5), исключим функцию Ψ и получим дифференциальное уравнение
для функции φ
Находим функцию
Функцию
интегрированием
находим теперь из первого уравнения (5)
Окончательно функции
и
получим, заменяя (X +1) через t,
Решение задачи Коши получим подстановкой в формулу (4) найденных функций, заменяя
при этом аргумент t у функции φ выражением (x + 2*y + cos y), а у функции
выражением (x - 2*y + cos y)
После раскрытия скобок и возможных упрощений получаем следующее выражение
функции, являющейся решением задачи Коши:
Правильность решения легко проверяется. Вычислив все производные первого и второго
порядков и подставив в уравнение (1), получаем тождество. Подставляя в функцию Z и ее
частную производную по y, значение y = 0, убеждаемся в выполнении начальных условий.
Решение правильное.
Пример 2. Найти решение уравнения
удовлетворяющее начальным условиям
Мы имеем уравнение параболического типа
Составляем характеристическое уравнение
Имеем только одно уравнение первого порядка (первое степени)
и, следовательно, один общий интеграл
Одно каноническое переменное определено
Второе каноническое переменное для уравнения параболического типа может быть
произвольной функцией от x и y. Но вместе функции
должны
разрешаться относительно x и y. Следует выбирать функцию
наиболее
простой, чтобы разрешимость функций
относительно x и y была
наиболее легкой. В данном примере наиболее простой будет функция
= y.
Производные первого порядка от функции Z по x и у выражаются через
производные по
равенствами
производные второго порядка – равенствами
Подставляя их в уравнение (6), получаем
Окончательный канонический вид уравнения подучается после замены в коэффициентах
уравнения x и y их выражениями через
. Каноническое уравнение
легко решается. Подстановкой
относительно функции U
= U оно приводится к уравнению первого порядка
Интегрируя обе части равенства по переменному
, имеем
Потенцируя последнее равенство и заменяя U черев
еще раз интегрируем по переменному η и находим общее решение
Заменив переменные
и их выражениями через x и y, получаем общее решение
первоначального уравнения (6)
Для выделения из него решения, удовлетворяющего заданным начальным условием (7),
дифференцируем функцию (8) по y
и подставляем в функцию (8) и ее производную, значение y -2. Учитывая условия (7),
получаем два соотношения
Дифференцируя первое равенство по X, умножаем обе части полученного равенства на 8 и
вычитаем соответственно из второго равенства. Получаем
Затем легко находим вторую функцию
В найденных функциях делаем подстановку x + 4 = t,
Решение задачи Коши получаем, подставляя в формулу (8) найденные функции
с заменой в них аргумента t суммой
и
Правильность решения проверяется. Найденная функция удовлетворяет начальным
условиям (7) и обращает уравнение (в) в тождество.
Домашнее задание. Задача 1
Найти решения линейных уравнений второго порядка, удовлетворяющие заданным
начальным условиям (задача Коши).
ЗАДАЧА 2
Решение методом Фурье смешанной задачи для линейного уравнения 2-го порядка
гиперболического типа в случае двух независимых переменных.
Рассматривается типовая задача. Можем считать, что это задача о плоских
колебаниях ограниченной однородной струны в среде без сопротивления. Концы струны
или закреплены или могут перемешаться в плоскости колебаний по перпендикулярам к
линии невозмущенного положения струны по определенным законам, зависящим от
времени. На струну действует непрерывно распределенная внешняя сила. Решение задачи
о движении точек струны сводится к решению линейного уравнения 2-го порядка
гиперболического типа от двух независимых переменных x и t вида
с заданными начальными условиями
и заданными граничными условиями
Функция U(x; t) есть смещение точки струны с абсциссой x в момент времени t в
направлении, перпендикулярном оси x, отрезком 0 <= x <= S которой изображается
ограниченная струна в невозмущенном состоянии. При решении задачи мы рассмотрим
только форму решения, не доводя решения до цифровых результатов. Поэтому некоторые
параметры задачи (длина струны S, плотность ρ, начальное натяжение Т0, а значит, и
) остаются выраженными буквенно, не задаются численно. Смещения точек
струны предполагаются достаточно малыми, соответственно малыми предполагаются и
функции f(x),
в начальных и граничных условиях, и свободный член g(x; t)
в уравнении. Малость этих функций выразим с помощью буквенных коэффициентов –
множителей α, β, μ, ν, λ. Эти коэффициенты имеют соответствующую размерность и в
конкретных задачах с числовыми данными малые числовые значения. Кроме того,
функции в начальных и граничных условиях должны удовлетворять условиям:
ЗАДАЧА 3
Численное приближенное решение первой краевой задачи для уравнения распространения
тепла в ограниченном стержне.
В уравнении теплопроводности для теплоизолированного по боковой поверхности
ограниченного стержня длины S
всегда можно сделать преобразование независимых переменных
после которого уравнение принимает вид
с коэффициентом a2 = 1. Длина стержня также равна 1 (0 <= y <= 1). Свободный член
уравнения изменяется
Соответственно видоизменяются функции в начальном и граничных условиях. Найдя
решение задачи
в новых переменных, легко можно перевести его в решение,
выраженное через старые переменные
В дальнейшем рассматриваем только уравнения теплопроводности с a2 = 1 и S = 1
Начальное распределение температур в стержне
На концах стержня поддерживается заданная температура
1. Окончательная формула решения уравнения (16), удовлетворяющего начальному
условию (17) и граничным условиям (18), записывается в виде ряда
где коэффициенты и функции вычисляются по формулам
bn = 2*
Для получения приближенного численного решения задачи разбиваем отрезок [0;1] на m
малых равных частей длины (шага) h. В точках разбиения xк = k:h (k = 0,1,2,.., m)
вычисляем “нулевую” строку таблицы решений из начального условия (17). Задавая
числовые значения tj = σ * j (j = 1,2,...) с малым промежутком σ между ними (шаг по t),
вычисляем приближенно последующие строки таблицы решений беря в решении по
формуле (19) кроме первых двух слагаемых еще несколько первых слагаемых ряда,
например четыре первых, отличных от нуля, члена.
Для примера зададим функции в правой части уравнения (16), в начальном условии
(17) и граничных условиях (18) в виде
По формулам (20) для заданных функций вычисляем
Решение задачи при заданных функциях (21) по формуле (18) представляется в виде
Числовые значения этого решения с шагом h = 0,1 по аргументу x вычисляем для t = 0 из
начального условия и для задаваемого t = 0,005 = 3*1/600 по формуле
ограничиваясь первыми четырьмя членами ряда в решении задачи. Два первых из
отбрасываемых членов (в скобках) в расчет не принимаются. Они выписаны для того,
чтобы примерно представлять точность приближения. Точность приближения 0,02.
Коэффициенты при синусах подсчитаны на формулы решения задачи. Например,
и т.д.
Таблица числовых значений решения задачи имеет вид:
X
0
0,1
0,2
0,3
0,4
U(x; 0)
0
0,099 0,192
0,273
0,338
U(x; 005) 0.3
0,1839 0,16573 0,24516 0,3378
0,5
0.6
0,7
0,8
0,9
1
0,375
0,384 0,357 0,288
0,171
0
0,3744
0,3583 0,3235 0,2681 0,1604 0
При вычислении значений тригонометрического многочлена с шагом h = 0,1
рекомендуется выписать значения sin18˚, sin36˚, sin54˚, sin72˚ из таблиц, а также
десятичные логарифмы синусов и коэффициентов. Произведения коэффициентов на
значения синусов указанных аргументов следует записывать и сохранять запись. В
процессе вычислений эти произведения встречаются повторно несколько раз со знаком
плюс или минус как результат применения формул приведения для синуса.
2. В разностном методе или методе сеток для уравнения теплопроводности (16)
область изменения независимых переменных 0 <=x <= 1, 0 <= t <= T разбивают на малые
прямоугольники сеткой прямых x=xi = h*i (i = 0,1,2,..., m) и прямых t = tj = τ*j (j = 0,1, 2,...,
p). Число h(mh = 1) называют шагом разбиения по аргументу x, а число τ (p*τ = T) шагом разбиения по аргументу t. Точки пересечения линий разбиения называют узлами
сетки. Для получения численного решения задачи нужно найти числовые значения
решения u(x;t) в узлах сетки
Явная схема метода сеток получается при замене в уравнении теплопроводности (16)
частных производных разностными отношениями
Уравнение (16) заменяется разностным уравнением
Числовые значения Ui,0(i=0,1,2,...,m) определяются из начального условия (17)
Значения решения U0,j и Um,j в крайних узлах сетки для каждого j = 1,2,..., p определяются
из граничных условий (18) задачи
Числа τg i,j = τ*g(h*i; τ*j) пpи выбранных h и τ рассчитываются для заданного свободного
члена g(x; t) уравнения (16) к представляются отдельной вспомогательной таблицей. По
формуле (22) можно последовательно вычислять числовые значения решения Ui,j+1 для i =
1,2,..., m-1 через известные уже, вычисленные на предыдущем шаге, значения решения
Ui+1,j, Ui,j, Ui-1,j и величины τ*gi,j. Геометрически схема изображена в виде
Вообще можно брать различные соотношения шагов h и τ в формуле (22), но устойчивое
решение методом сеток для явной схемы имеет место при условии
. Ошибка
аппроксимации, образующаяся в результате замены уравнения в частных производных
(16) разностным уравнением (22), также зависит от отношения
Наименьшая
ошибка имеет порядок h4 при
употребительна
. Получающаяся при этом формула наиболее
Для примера с функциями (21) берем h = 0,1 и, следовательно,
вспомогательную таблицу чисел
. Вычисляем
Ограничимся тремя строками j = 0,1,2. При атом для заданного примера первая строка
состоит из нулей. Значения i следует брать от 1 до m-1 (9 в нашем примере), так как g0,j и
gm,j в решении задачи не участвуют. Вспомогательная таблица чисел
имеет для
рассматриваемого примера следующий вид:
i
1
2
3
4
5
6
7
8
9
τ*gi,0 0
0
0
0
0
0
0
0
0
τ*gi,1 10-6*9 10-6*16 10-6*21 10-6*24 10-6*25 10-6*24 10-6*21 10-6*16 10-6*9
τ*gi,2 l0-6*18 10-6*32 10-6*42 10-6*48 10-6*50 10-6*48 10-6*42 10-6*32 10-6*18
Числовые значения решения задачи для примера с функциями (21 представлены
следующей таблицей:
xi
0 0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
ui,0 0 0,098 0,192 0,273 0,336 0,375 0,384 0,357 0,288 0,171 0
ui,1 0,1 0,098 0,190 0,270 0,332 0,370 0,378 0,350 0,280 0,162 0
ui,2 0,2 0,114 0,188 0,267 0,328 0,365 0,372 0,343 0,272 0,155 0
ui,3 0,3 0,140 0,189 0,264 0,324 0,360 0,366 0,336 0,264 0,148 0
В первой строке таблицы записаны значения аргумента x, а во второй - значения Ui,0.
вычисленные из начального условия. Вычисления по формуле (22') начинаются с третьей
строки. Например,
и т.д. Вычислены три строки. Последняя соответствует значениям
Ui,3 = U(0,1*i; 0,005). Сравнивая с соответствующими приближенными значениями U
(0,1*i; 0,005), вычисленными для решения методом Фурье, видим отклонения порядка
0,02 и более (для U1,3 отклонение 0,04), Если же учесть значения первых двух
отбрасываемых членов (в скобке) при решении по методу Фурье, то получим отклонения
значительно меньшие. Следовательно, значения решения по методу сеток более точны,
чем по методу Фурье с четырьмя первыми членами ряда. Заметим, что в рассматриваемом
примере на значения решения методом сеток очень мало влияли слагаемые τ*gi,j. Это
объясняется спецификой свободного члена g(x; t) в данном примере и малостью
рассмотренного отрезка времени
3. Неявная схема метода сеток для уравнения теплопроводности (16) получается при
замене частных производных разностными отношениями
Уравнение в частных производных (16) заменяется аппроксимирующим разностным
уравнением
Обозначая для краткости
, представим разностное уравнение в виде
Начальная строка таблицы числовых значений решения задачи определяется из
начального условия
Из граничных условий определяются значения
для каждого значения j = 1,2,…, p. Зная значения решения в предыдущей строке таблицы с
номером j - 1, находят значения решения для последующей строки таблицы с номером j из
системы линейных алгебраических уравнений (25) для
. Это система специфического “трехдиагонального” вида и
решается методом прогонки. Система может быть преобразована к виду
где вспомогательные величины ai,j и bi,j вычисляются последовательно по формулам
Вычисление этих величин называют "прямым ходом". С помощью найденных величин и
второго граничного условия задачи последовательно вычисляются значения решения для
строки с номером j
Эти вычисления называют "обратным ходом" прогонки. Получают числовые значения
решения задачи для строки с номером j. Для следующей строки процесс повторяется.
Отметим, что числа ai,j фактически не зависят от номера j. Эти числа вычисляются один
раз для всех строк. Числа bi,j для каждой строки вообще получаются различные. Строку
чисел ai,j и каждую строку чисел bi,j при расчетах удобно включать в таблицу числовых
значений решения Ui,j, причем строки чисел bi,j последовательно чередовать со строками
чисел Ui,j. В эту же таблицу можно вписать и строки чисел gi,j или чисел h2*gi,j, которые
входят слагаемыми в числа bi,j. Для примера с функциями (21) методом прогонки
вычислим одну строку чисел Ui,1, полагая h = 0,1 и
вычислений составляют следующую таблицу;
xi
0 0,1
0,2
0,8
0,4
0,5
Ui,0 0 0,099
0,182
0,273
0,386
0,875
h2gi,1 0 10-6 162 10-6 288 10-6 378 10-6 432 10-6 450
ai,1 - 0,25
0,2667 0,2678 0,2679 0,2679
bi,1 - 0,4982 0,5088 0,6321 0,8551 0,9788
Ui,1 0,3 0,1465 0,2077 0,2700 0,3258 0,3607
Продолжение таблицы
0,7
0,8
0,9
0,357
0,288
0,171
10-6 378 10-6 288 10-6 162
0,2679 0,2679 0,2070
0,9906 0,8417 0,5677
0,3368 0,2663 0,1521
. Результаты
0,6
0,384
10-6 432
0,267»
1,0309
0,3665
1
0
0
0
В первой строке таблицы записаны абсциссы узлов сетки, во второй строке - числовые
значения решения при t = 0, определяемые из начального условия задачи. В третьей
.
строке записаны значения h2*g(x; t) узлов сетки xi = h*i при h = 0,1 и
В четвертой и пятой строках записаны последовательно вычисленные значения числовых
коэффициентов (27) "прямого хода". Например, (ω = h2 / τ = 2)
В последней строке записаны числовые значения решения Ui,1 (t1 - τ = 0,005),
последовательно получаемые по формулам (28) "обратного хода". Например,
Сравнивая последнюю строку полученной таблицы с последней строкой таблицы метода
сеток явной схемы, видим отклонения, не превосходящие 0,01 (для U2,1 несколько
большие). Это соответствует порядку ошибки аппроксимации метода сеток неявной
схемы. Эта ошибка пропорциональна
Домашнее задание. Задача 3
Для приведенного уравнения теплопроводности (16) с начальными условиями (17) и
граничными условиями (18) найти численное приближенное решение с шагом h = 0,1 по
аргументу x для значения t = 0,005. Решение получить тремя способами.
1) Методом Фурье найти решение первой краевой задачи и подсчитать
приближенное решение, взяв 4 первых отличных от нуля члена ряда.
2) Методом сеток явной схемы подсчитать таблицу приближенных числовых
значений решения для
3) Методом сеток неявной схемы с применением метода прогонки подсчитать
таблицу числовых значений решения для t0 = 0 и t1 = 0,005 и значения вспомогательных
величин h2*gi,1, ai,1, bi,1. Функции
для каждого варианта задачи
представлены следующей таблицей:
№
g(x; t)
f(x)
φ(t)
Ψ(t)
по пор.
1
2
3
4
5
2
2
1
30x(1-x)
x (1-x)
30t
0
2
2
2
30x (1-X)
x(1-x)
0
60t
2
2
3
60x (1-x)
x(1-x)
90t
0
4
36x(1-x)2
x(1-x)
0
90t
2
5
36x (1-x)t
x
0
1
6
72x(1-x)t
1-x
1
0
7
72x2(1-x)2t
2x
0
2
2
8
36x(1-x) t
2(1-x)
2
0
2
9
30x
х(1-х )
60t
0
2
10
30x
x(1-х)
0
60t
2
2
11
60(1-x)
x (1-х)
30t
0
12
36x2t
x(1-х2)
0
30t
2
13
36(1-x) t
3х
0
4
2
14
72xt
(1-x)
1
0
2
15
72x t
2х
0
2
16
60x(1+x)
2х(1-x2)
60t
0
2
17
30x(2+x)
2x (1-x)
0
60t
18
90x(x-2)
3х(1-x)
30t
0
19
60x(x-3)
Зх(1-x2)
0
90t
2
20
90x(x-4)
Зх(1-x)
45t
0
21
30x(1-x)(2+x)t Зх
0
3
22
60x(1-x)(2-x)t 2х
0
2
23
60x(1-x)(3+x)t 2(1-x)
2
1
24
30x(1-x)(3-x)t З(1-x)
3
0
2
2
25
90x (1-x)t
(1-x)
1
0
2
26
60x
х (1-x2)
30t
0
27
30(1+x)
x(1-x2)
0
45t
2
28
60(1-x )
x(1-x)(2+x) 45t
0
29
30
30(1-x)(2+x)
60(1+x)(2-x)
x2(1-x)
x(1-x)(2-x)
0
90t
60t
0
Литература
1. Н.С. Кошляков, Э.Б. Глинер, М.М. Смирнов. Уравнения в частных производных
математической физики. М., "Высшая школа", 1970.
2. А.Н. Тихонов, А.А, Самарский. Уравнения математической физики. М., "Наука",
1966.
3. И.Г. Араманович, В.И. Левин. Уравнения математической физики. М., "Наука",
1964.
4. Б.П. Демидович, ИД. Марон, Э.З. Шувалова. Численные методы анализа, М.,
"Физматгиз", 1963.
Б. Н.В. Копченова, И.А. Марон. Вычислительная математика в примерах и
задачах, М., "Наука", 1972.
6. Б.И, Сегал, К.Л., Семендяев. Пятизначные математические таблицы. М.,
Физматгиз, 1962.
Download