Некоторые дополнительные вычислительные методы

Министерство образования и науки РФ

ГОУ ВПО “УГТУ-УПИ”

Курсовая работа

по “Вычислительной математике”

на тему: “Некоторые дополнительные вычислительные методы”

Семестр № 3

Преподаватель

Кочнев В.П.

Студент гр. № р-23021д

Логиновских М.А.

Номер зачетной книжки 17309013

Екатеринбург

2004

____________________________________________________________________________

Домашнее задание по ________________________________ № ________________

№ записи в книге регистрации __________________ дата регистрации

___________200_г.

Преподаватель _________________________________________

Студент _________________________________________ группа № ________________

Деканат ФДО _______________

СОДЕРЖАНИЕ

1. Решение систем линейных уравнений …………………………………………………… 3

а) Схема Халецкого ……………………………………………………………………....... 3

б) Метод Зейделя и условия сходимости ………………………………………………… 5

2. Методы решения нелинейных уравнений ……………………………………………….. 6

а) Метод хорд ………………………………………………………………………………. 7

б) Метод Ньютона (метод касательных) …………………………………………………. 8

в) Метод итерации ………………………………………………………………………… 9

3. Интерполирование и экстраполирование ……………………………………………….. 11

а) Интерполирование с помощью многочленов ………………………………………… 11

б) Интерполяционный многочлен Лагранжа ……………………………………………. 12

в) Интерполяционные многочлены Стирлинга и Бесселя ……………………………… 13

г) Тригонометрическое интерполирование …………………………………….………... 15

д) Интерполяция сплайнами ……………………………………………………..………... 15

4. Численное дифференцирование и интегрирование ……………………………….…….. 16

а) Постановка задачи численного интегрирования ……………………………………... 16

б) Составные квадратурные формулы ………………………………………………….… 17

5. Приближенное вычисление обыкновенных дифференциальных уравнений …………..

а) Метод Рунге-Кутта ……………………………………………………………………… 18

б) Экстраполяционные методы Адамса ………………………………………………….. 20

в) Метод Милна ……………………………………………………………………………. 20

г) Краевые задачи для обыкновенных дифференциальных уравнений ………………...

6. Приближенные методы решения дифференциальных уравнений с частными

производными ………………………………………………………………………………… 21

а) Классификация дифференциальных уравнений второго порядка …………………… 22

б) Постановка краевых задач ……………………………………………………………… 23

в) Метод конечных разностей (метод сеток) …………………………………………….. 24

г) Разностные схемы для решения уравнения теплопроводности ……………………… 25

д) Разностные схемы для решения уравнения колебания струны ……………………… 26

7. Список литературы ………………………………………………………………………… 27

1. Решение систем линейных уравнений

Системы линейных уравнений (СЛУ) имеют в вычислениях очень большое

значение, так как к ним может быть приведено приближенное решение широкого

круга задач. Так, основными источниками возникновения СЛУ являются теория

электрических цепей, уравнения балансов и сохранения в механике, гидравлике

и т.д. Существует несколько способов решения таких систем, которые в

основном делятся на два типа: 1) точные методы, представляющие собой

конечные алгоритмы для вычисления корней системы, 2) итерационные методы,

позволяющие получать корни системы с заданной точностью путем сходящихся

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

приближенными из-за неизбежных округлений. Для итерационных процессов также

добавляется погрешность метода.

Пример системы линейных уравнений: [pic]

Или в матричном виде: [pic],

где [pic]матрица коэффициентов системы;

[pic] - вектор неизвестных; [pic] - вектор свободных членов.

Схема Халецкого

Запишем систему линейных уравнений в матричном виде: [pic],

где A=[aij] – квадратная матрица порядка n и

[pic], [pic] - векторы-столбцы.

Представим матрицу A в виде произведения нижней треугольной матрицы B=[bij]

и верхней треугольной матрицы C=[cij] с единичной диагональю [pic], где

[pic] и [pic].

Тогда элементы bij и cij определяются по формулам

[pic] и [pic]

Отсюда искомый вектор x может быть вычислен из уравнений [pic] и [pic].

Так как матрицы B и C – треугольные, то системы легко решаются:

[pic] и [pic]

Из этих двух формул видно, что числа yi выгодно вычислять вместе с

коэффициентами cij. Этот метод получил название схемы Халецкого. В схеме

применяется обычный контроль с помощью сумм. Если матрица A –

симметрическая aij=aji, то [pic]

Пример. Решить систему [pic]

Решение.

В первый раздел таблицы впишем матрицу коэффициентов системы, ее свободные

члены и контрольные суммы. Далее так как [pic] [pic], то первый столбец из

раздела 1 переносится в первый столбец раздела II. Чтобы получить первую

строку раздела II, делим все элементы первой строки раздела I на

элемент[pic], в нашем случае на 3.

Имеем: [pic]; [pic]; [pic]; [pic]; [pic].

Переходим к заполнению второго столбца раздела II, начиная со второй

строки. Пользуясь формулами, определяем [pic]: [pic]; [pic]; [pic].

Далее определяя по формулам, заполняем вторую сетку для раздела II:

[pic][pic]

[pic]

Затем переходим к третьему столбцу, вычисляя его элементы [pic] и [pic] по

формулам и т.д., пока не будет заполнена вся таблица раздела II. Таким

образом, заполнение раздела II происходит способом “елочки”: столбец -

строка, столбец - строка и т.д.

В разделе Ш, пользуясь формулами, определяем [pic] и [pic].

Текущий контроль осуществляется с помощью столбца S, над которым

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

|0 |1,200|0,000|0,000|

| |0 |0 |0 |

|1 |1,200|1 |0,948|

| |0 |,0600|0 |

|2 |0,999|1,005|0,999|

| |2 |4 |1 |

|3 |0,999|1.000|1,000|

| |6 |1 |1 |

|4 |1 |1,000|1,000|

| |,0000|0 |0 |

|5 |1 |1,000|1,000|

| |,0000|0 |0 |

Точные значения корней: [pic].

2. Методы решения нелинейных уравнений

Как известно, далеко не всякое уравнение f(x)=0 можно решить точно, т.е. не

всегда можно найти число [pic] такое что f([pic])?0. В первую очередь это

относится к трансцендентным уравнениям. Кроме того, даже для алгебраических

уравнений степени выше четвертой не существуют формулы, выражающей их

решения через коэффициенты уравнения при помощи арифметических операций и

извлечение корней. Для уравнений третьей и четвертой степени формулы для

отыскания корней существуют, но они настолько сложны, что практически не

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

уравнения f(x)=0. Для этого существует множество методов некоторые, из

которых мы рассмотрим.

Метод хорд

Пусть дано уравнение f(x)=0, где функция f(x) определена и непрерывна на

интервале

[a, b] и f(a)f(b)0. Разделим

отрезок [a, b] в отношении - f(a):f(b). Это даст нам приближенное значение

корня x1 = a + h1, где

[pic].

Далее этот прием применяем к одному из отрезков [a, x1] или [x1, b], на

концах которого функция f(x) имеет противоположные знаки. Аналогично

находим второе приближение x2 и т.д. Геометрически этот способ эквивалентен

замене кривой y = f(x) хордой, проходящей через точки А[a, f(a)] и B[b,

f(b)].

f(b)

f(a)

? x3 x2 x1 b=x0 a=x0 x1 x2

f(b)

Действительно, уравнение хорды АВ имеет вид [pic]

При х = х1 и y = 0, получим [pic]

Полагая, что на отрезке [a, b] вторая производная f''(x) сохраняет

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

Из рис. 1 видно, что конец а неподвижен и последовательные приближения:

x0=b;

[pic]

образуют ограниченную монотонно убывающую последовательность, причем

a0.

Пример. Вычислить методом Ньютона отрицательный корень уравнения

[pic] с пятью верными знаками.

Решение. Полагая в левой части уравнение [pic] получим [pic].

Следовательно, искомый корень [pic] находится в интервале [pic]. Сузим

найденный интервал. Так как [pic] то [pic]. В этом последнем интервале

[pic] и [pic].Так как [pic] и [pic], то можем принять за начальное

приближение [pic]. Последовательные приближения [pic] вычисляем по

следующей схеме:

|[pi|[pic] |[pic] |[pic] |[pic] |

|c] | | | | |

|0 |-11 |3453 |-5183 |0,7 |

|1 |-10,3 |134,3 |-4234 |0,03 |

|2 |-10,27 |37,8 |-4196 |0,009 |

|3 |-10,261|0,2 |- |- |

Останавливаясь на [pic], проверяем знак значения [pic]. Так как [pic], то,

[pic] и любое из этих чисел дает искомое приближение.

Метод итерации

Заменим уравнение f(x)=0 эквивалентным x=?(x). Выберем некоторое начальное

приближение x0 и вычислим дальнейшие приближения по формулам x1= ?(x0), x2=

?(x1), …, xn= ?(xn-1). Если последовательность xn имеет предел, то

итерационный процесс

xn= ?(xn-1) (n=1, 2, …) называется сходящимся. Пусть функция ?(x)

непрерывна. Переходя к пределу в равенстве xn= ?(xn-1), получим

[pic] Следовательно, [pic] является корнем уравнения x=?(x) и может быть

вычислен по формуле xn= ?(xn-1) (n=1, 2, …) с любой точностью. Для

данного метода существуют две теоремы:

Теорема 1. Пусть корень [pic] уравнения x=?(x), а также его

последовательные приближения x0, xn= ?(xn-1) (n=1, 2, …) содержатся в

интервале [a, b] и [pic] на [a, b]. Тогда справедливы утверждения:

1. итерационный процесс xn= ?(xn-1) сходится к корню уравнения [pic];

2. [pic];

3. [pic].

Следствие 1. Если требуется найти корень с точностью ?, то кончаем

итерационный процесс тогда, когда [pic]0 на

(a, b), то последовательные приближения xn= ?(xn-1) (n=1, 2, …) сходятся

к корню монотонно; если ?'(x)1, то итерационный процесс расходится. Найдем функцию

[pic], обратную функции ?(x). Так как [pic], то итерационный процесс [pic]

будет сходится. [pic], результаты вычислений приведены в таблице:

| n| |[pic] |[pic] |[pic] |

| |zn | | | |

| 0 |0,5 | | 0,503 | |

|1 |0,50|-0,688 |0,504 |------ |

|2 |3 |-0,686 |------ |0,0015 |

| |0,50|------ | |0,0005 |

| |3 | | | |

3. Интерполирование и экстраполирование

Задача интерполирования состоит в том, чтобы по значениям функции f(x) в

нескольких точках отрезка восстановить ее значения в остальных точках

данного отрезка. Разумеется, такая постановка задачи допускает сколь угодно

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

когда известны результаты измерений yk = f(xk) некоторой физической

величины f(x) в точках xk, k = 0, 1,…, n и требуется определить ее

значение в других точках. Интерполирование используется также при

необходимости сгущения таблиц, когда вычисление значений f(x) по точным

формулам трудоемко. Иногда возникает необходимость приближенной замены

(аппроксимации) данной функции (обычно заданной таблицей) другими

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

(экспериментальных) зависимостей, результаты обычно представлены в

табличном или графическом виде. Задача заключается в аналитическом

представлении искомой функциональной зависимости, то есть в подборе

формулы, корректно описывающей экспериментальные данные.

Интерполирование с помощью многочленов

Пусть функциональная зависимость задана таблицей y0 = f(x0); …, y1= f(x1);

…, yn = f(xn). Обычно задача интерполирования формулируется так: найти

многочлен P(x) = Pn(x) степени не выше n, значения которого в точках xi

(i = 0, 1 2,…, n) совпадают со значениями данной функции, то есть

P(xi) = yi. Геометрически это означает, что нужно найти алгебраическую

кривую вида [pic] проходящую через заданную систему точек Мi(xi, yi) (см.

рис. 4). Многочлен Р(х) называется интерполяционным многочленом. Точки xi

(i = 0, 1, 2,…, n) называются узлами интерполяции.

Для любой непрерывной функции f(x) сформулированная задача имеет

единственное решение. Действительно, для отыскания коэффициентов а0, а1, а2

,…, аn получаем систему линейных уравнений [pic] определитель которой

отличен от нуля, если среди точек xi (i = 0, 1, 2,…, n) нет совпадающих.

Решение системы можно записать различным образом. Однако наиболее

употребительна запись интерполяционного многочлена в форме Лагранжа или в

форме Ньютона.

Интерполяционный многочлен Лагранжа

Пусть на отрезке [a,b] некоторая функция f(x) задана лишь в некоторых

точках [pic], т.е. известны ее значения [pic], которые, собирают в таблицу:

|x |x0 |x1 |...|xn |

|f(x|y0 |y1 |...|yn |

|) | | | | |

Кроме того, пусть задана некоторая точка [pic]. Построим по таблице

следующий многочлен: [pic].

Этот многочлен называется многочленом Лагранжа.

Его основные свойства:

1) это - многочлен степени [pic];

2) [pic], т.е. многочлен Лагранжа имеет в точках [pic] те же

значения, что и функция [pic];

3) если фиксировать любое число [pic] то окажется выполненным

неравенство [pic]

где [pic] на участке [pic], т.е. число [pic] ограничивает производную

[pic]го порядка функции [pic].

Сказанное означает, что если функция [pic] задана своей таблицей и

требуется найти значение [pic] где-то в промежуточной точке c, то можно по

таблице построить многочлен Лагранжа и его значение в этой точке принять за

значение функции. Отыскание промежуточного значения функции называется

интерполяцией; когда это делается с помощью многочлена Лагранжа, то говорят

об интерполяционном многочлене Лагранжа или об интерполяции по Лагранжу.

Пример. Построить интерполяционный многочлен Лагранжа для функции заданной

таблицей

|x |1 |2 |3 |5 |

|y |1 |5 |14 |81 |

И найти значение функции при x=4.

Решение. Используя формулу Лагранжа найдем: [pic]

После некоторых преобразований получим [pic] Тогда f(4)?L3(4)=36,5.

Интерполяционные многочлены Стирлинга и Бесселя

Взяв среднее арифметическое первой и второй интерполяционных формул Гаусса

[pic] [pic]

[pic] и [pic] [pic]

[pic], получим формулу Стирлинга [pic]

где [pic].

Легко видеть, что [pic] при [pic].

Кроме формулы Стирлинга, часто употребляется формула Бесселя. Для вывода

этой формулы воспользуемся второй интерполяционной формулой Гаусса [pic]

[pic]

[pic].

Возьмем [pic] равностоящих узлов интерполирования [pic] с шагом [pic], и

пусть [pic] — заданные значения функции [pic].

Если выбрать за начальные значения [pic] и [pic], то, используя узлы [pic],

будем иметь:

[pic][pic].

Примем теперь за начальные значения [pic] и [pic] и используем узлы [pic].

Тогда [pic], причем соответственно индексы всех разностей в правой части

предыдущей формулы возрастут на единицу. Заменив в правой части этой

формулы [pic] на [pic] и увеличив индексы всех разностей на 1, получим

вспомогательную интерполяционную формулу:

[pic][pic].

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

интерполяционную формулу Бесселя

[pic]

где [pic].

Интерполяционная формула Бесселя, как следует из способа получения ее,

представляет собой полином, совпадающий с данной функцией [pic] в [pic]

точках [pic].

Тригонометрическое интерполирование

Пусть функция f(х) представлена на некотором отрезке [0, 2?] таблицей

значений f(хi) в

равноотстоящих узлах хi =2?(i-1)/(2N+1), i =1, 2, ..., 2N+1. Тогда

тригонометрическим интерполирующим многочленом назовем многочлен степени m

вида:

[pic].

Задача тригонометрической интерполяции состоит в построении

тригонометрического полинома, который бы наиболее полно удовлетворял

условиям Рm (хi)= f(хi ) для любого i=1, 2, ..., 2 N+1.

Можно показать, что решением этой задачи является полином именно того вида,

коэффициенты которого вычисляют по следующим формулам:

[pic];

[pic].

Интерполяция сплайнами

Пусть отрезок [a, b] разбит на n равных частей [xi, xi+1], где xi=a+ih,

i=0, ..., n, xn=b,

h=(b-a)/n.

Сплайном называется функция, которая вместе с несколькими производными

непрерывна на всем заданном отрезке [a, b], а на каждом частичном отрезке

[xi, xi+1] в отдельности является некоторым алгебраическим многочленом.

Максимальная по всем частичным отрезкам степень многочленов называется

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

непрерывной на [a,b] производной - дефектом сплайна.

На практике широкое распространение получили сплайны третьей степени,

имеющие на [a, b] непрерывную, по крайней мере, первую производную. Эти

сплайны называются кубическими и обозначаются S3(x).

Пусть на отрезке [a, b] в узлах сетки ( заданы значения некоторой функции

fi =f(xi), i=0, ..., n.

Интерполяционным кубическим сплайном S3(x) называется сплайн

S3(x)=аi0 +аi1(x - xi)+аi2(x - xi)2 +аi3(x - xi)3, x([xi, xi+1],

удовлетворяющий условиям

S3(xi)=f(xi), i=0, ..., n.

Данный сплайн на каждом из отрезков [xi, xi+1], i=0, ..., n-1 определяется

четырьмя коэффициентами, и поэтому для его построения на всем промежутке

[a, b] необходимо определить 4n коэффициентов. Для их однозначного

определения необходимо задать 4n уравнений.

Условие S3(xi)=f(xi), i=0, ..., n дает 2n уравнений, при этом функция

S3(xi), удовлетворяющая этим условиям, будет непрерывна во всех внутренних

узлах.

Условие непрерывности производных сплайна [pic], r=1,2 во всех внутренних

узлах xi, i=1, ..., n-1 сетки ( дает 2(n-1) равенств.

Вместе получается 4N-2 уравнений.

Два дополнительных условия обычно задаются в виде ограничений на значение

производных сплайна на концах промежутка [a, b] и называются краевыми

условиями.

Наиболее употребительны следующие типы краевых условий:

а) S'3(а)=f'(а), S'(b)=f'(b);

б) S"3(а)=f"(а), S"(b)=f"(b);

в) [pic];

г) S'''3(x p+0)=S'''3(x p-0), р =1, n-1.

4. Численное дифференцирование и интегрирование

Если функция f(x) заданна аналитически ее первообразная F(x) является

элементарной функцией, то [pic] вычисляется по формуле Ньютона-Лейбница:

[pic] В тех случаях, когда функция f(x) задана аналитически, но ее

первообразная не является элементарной функцией или отыскать ее сложно, а

также в случае, когда функция f(x) задана графически или таблично, для

вычисления [pic] применяются приближенные методы.

Постановка задачи численного интегрирования

Задача численного интегрирования функции заключается в вычислении

определенного интеграла на основании ряда значений подынтегральной функции.

Численное вычисление однократного интеграла называется механической

квадратурой. Обычный прием механической квадратуры состоит в том, что

данную функцию f(x) на рассматриваемом отрезке [a, b] заменяют

интерполирующей или аппроксимирующей функцией ?(x) простого вида, а затем

приближенно полагают: [pic] Функция ?(x) должна быть такова, чтобы интеграл

[pic] вычислялся непосредственно. Если функция f(x) заданна аналитически,

то ставится вопрос об оценке погрешности. Пусть для функции y=f(x) известны

в n+1 точках x0, x1, x2, …, xn отрезка [a, b] соответствующие значения

f(xi)=yi (i=0, 1, 2, …, n). Требуется приближенно найти [pic] По

заданным значениям yi построим полином Лагранжа [pic], где

Пn+1(x)=(x-x0)(x-x1)…(x-xn), причем Ln(xi)=yi (i=0, 1, 2, …, n). Заменяя

функцию f(x) полиномом Ln(x), получим равенство [pic] где Rn[f] – ошибка

квадратурной формулы. Отсюда получаем приближенную квадратурную формулу

[pic] где [pic] (i=0, 1, 2, …, n). Для вычисления Ai заметим, что

1) коэффициенты Ai при данном расположении узлов не зависят от выбора

функции f(x);

2) для полинома степени n полученная формула – точная, так как тогда

Ln(x)=f(x); следовательно, формула [pic] - точная при y=xk (k=0, 1, 2,

…, n), т.е. Rn[xk]=0 при k=0, 1, …, n. Полагая y=xk (k=0, 1, 2, …, n),

получим линейную систему из n+1 уравнений [pic] - где [pic] (k=0, 1, …,

n), из которой можно определить коэффициенты A0, A1, …, An.

Составные квадратурные формулы

Приведем ряд простейших квадратурных формул, используемых в практике

численного интегрирования функции f(x) на некотором интервале [a, b],

разбитого на n равных отрезков точками a0=a, a1=a+h, a2=a+2h, …, an=a+nh+b,

где n=0,1, …, k и [pic] Положим f(xn)=yn=f(a+nh).

Формула прямоугольников: [pic]

Погрешность формулы определяется выражением

[pic] где [pic]

Формула трапеций: [pic]

Погрешность формулы определяется выражением

[pic] где [pic]

Формула Симпсона: [pic] где [pic]

Погрешность формулы определяется выражением

[pic] где [pic]

Если длина интервала [a, b] велика для применения простейших квадратурных

формул, то поступают следующим образом:

1) интервал [a, b] разбивают точками xi, [pic] на n интервалов по

некоторому правилу;

2) на каждом частичном интервале [xi, xi+1] применяют простейшую

квадратурную формулу, находят приближенное значение интеграла [pic]

3) из полученных выражений Qi составляют (отсюда и название составная

формула) квадратурную формулу для всего интервала [a, b];

4) абсолютную погрешность R составной формулы находят суммированием

погрешностей Ri на каждом частичном интервале.

5. Приближенное вычисление обыкновенных дифференциальных уравнений

Обыкновенным дифференциальным уравнением называется равенство [pic], в

котором [pic] - независимая переменная, изменяющаяся в некотором отрезке

[pic], а [pic] - неизвестная функция от [pic], которую и надо найти.

Различают два типа обыкновенных дифференциальных уравнений - уравнения без

начальных условий и уравнения с начальными условиями. Уравнения без

начальных условий - это как раз то, что было только что определено. А

уравнение с начальными условиями - это записанное выше уравнение

относительно функции [pic], но в котором требуется найти лишь такую функцию

[pic], которая удовлетворяет при некотором [pic] следующим условиям:

[pic], т.е. в точке [pic] функция [pic]и ее первые [pic] производных

принимают наперед заданные значения. В этой ситуации число [pic]называется

порядком уравнения.

Метод Рунге-Кутта

Изложим идею метода на примере: [pic]

Интегрируя это уравнение в пределах от x до x + h (0 < h 0 – коэффициент

теплопроводности; ?0, ?1, ?2 – заданные на Г функции, причем ?2 есть

произведение коэффициента теплопроводности на температуру внешней среды,

соприкасающейся с телом.

Таким образом, краевая задача заключается в том, чтобы найти классическое

решение уравнения Пуассона или Лапласа, удовлетворяющее одному из граничных

условий.

3. Смешанная краевая задача. Рассмотрим задачу распространения тепла в

тонком

стержне единичной длины. Поместим один из концов в точку x=0, а другой – в

точку x=1. Распределение температуры в таком стержне в течение некоторого

интервала времени 0

[pic], а для единственности решения в этом случае необходимо еще задать

температурный режим на концах стержня. Это можно сделать с помощью

граничных условий, аналогичных тем, которые были сформулированы для

уравнений Пуассона и Лапласа.

Граничное условие I рода (на конце стержня x=0 заданна температура): [pic].

Граничное условие II рода (на конце стержня x=0 задан тепловой поток):

[pic].

Граничное условие III рода: [pic].

Для другого конца стержня x=1 правые части граничных условий заменяются

соответственно на ?0(t), ?1(t), ?2(t). Заметим, что начальное и граничное

условия должны удовлетворять так называемым условиям сопряжения, т.е. при

условии I рода u0(0)=?0(0), при условии II рода u0x(0)=?1(0), при условии

III рода -u0x(0)+?u0(0)=?2(0). Аналогичные условия сопряжения должны

выполнятся и на другом конце стержня x=1.

Сформулируем одну из возможных краевых задач. Найти классическое решение

уравнения [pic], удовлетворяющее начальному условию [pic] и следующим

граничным условиям [pic]. Эта задача обычно называется первой краевой

задачей для уравнения теплопроводности. Соответственно краевые задачи с

граничными условиями II роди или III называются второй и третьей краевой

задачей для уравнения теплопроводности.

Метод конечных разностей (метод сеток)

Численные методы, основанные на разностной аппроксимации производных

называется разностным методом, методом конечных разностей или методом

сеток.

Пусть заданно линейное дифференциальное уравнение, записанное в

символическом виде: [pic]. Здесь u – искомое решение уравнения; L –

некоторый дифференциальный оператор, сокращенно обозначающий

соответствующую дифференциальную операцию; f – правая часть уравнения

(заданная функция).

Для единственного решения данного уравнения к нему необходимо присоединить

краевые условия: [pic].

Разностный метод решения этих двух задач можно представить в виде двух

этапов:

1. построение разностной схемы, аппроксимирующей данную непрерывную

задачу;

2. получение решения разностной задачи и оценка погрешности этого

решения.

Для построения разностной схемы первым шагом является замена области [pic]

непрерывного изменения аргументов областью дискретного их изменения –

сеточной областью [pic], т.е.множеством точек (xn, ym), называемых узлами

сетки. Для квадрата [pic]сеточную область можно построить следующим

образом. Проведем прямые [pic]. Множество точек пересечения этих прямых и

составит сеточную область, а сами точки образуют узлы сетки. Всякая функция

[pic], определенная на ссеке [pic], называется сеточной функцией и

обозначается [pic].

Второй шаг в построении разностной схемы состоит в аппроксимации

дифференциального выражения Lu некоторым разностным выражением, а функцию

непрерывного аргумента f – сеточной функцией, т.е. в построение некоторого

разностного аналога для данного уравнения, при данных краевых условиях.

Такая аппроксимация приводит к системе алгебраических уравнений

относительно значений некоторой сеточной функции [pic]. Эту систему можно

записать в следующем виде:

[pic]

Где Lh и ?h – разностные операторы, аппроксимирующие соответственно L и l;

vh – искомая сеточная функция, аппроксимирующая решение u; fh, ?h –

заданные сеточные функции, аппроксимирующие f и ?.

Совокупность разносных уравнений, аппроксимирующих исходную задачу – есть

разностная схема. Рассмотрим их подробнее на примерах уравнения

теплопроводности и колебания струны.

Разностные схемы для решения уравнения теплопроводности (параболический

тип)

Рассмотрим первую краевую задачу для уравнения теплопроводности в

прямоугольнике [pic]. Требуется найти непрерывное в [pic] решение задачи:

[pic]

В области [pic] введем прямоугольную равномерную сетку [pic]{xn, tk} с

шагом h=1/N по координате x и с шагом ?=T/M по координате t:

[pic].

Производные левой части уравнения [pic] аппроксимируем следующим

разностными выражениями: [pic]

В соответствии с данной аппроксимацией построим два разностных аналога

уравнения [pic] с неизвестной сеточной функцией vh?: [pic]

Здесь [pic] - значение некоторой сеточной функции fh?, соответствующей

правой части уравнения [pic]. Для первой разностной схемы [pic], а для

второй - [pic] .

Начальное и граничное условия для первой краевой задачи аппроксимируются

точно: [pic]

Для второй и третьей краевых задач граничные условия аппроксимируются на

основе разностных выражений.

Полагая r=?/h2 получим [pic] - для первой разностной схемы, [pic] - - для

второй разностной схемы.

Анализ показывает, что погрешность аппроксимации схем есть [pic].

Разностные схемы для решения уравнения колебания струны (гиперболический

тип)

Рассмотрим первую краевую задачу для уравнения колебания струны в

прямоугольнике [pic]. Требуется найти непрерывное в [pic] решение задачи:

[pic]

Применение метода конечных разностей к решению задачи по существу мало чем

отличается от его применения к уравнению теплопроводности. Область [pic]

покрывается сеткой [pic]. Отличие заключается в приближении второй

производной по переменной t:

[pic].

Разностная аппроксимация принимает вид

[pic].

Начальные условия аппроксимируются следующим образом: [pic].

Граничные условия аппроксимируются точно так же, как и для уравнения

теплопроводности: [pic].

Значение [pic] является фиктивным неизвестным, которое можно определить по

формуле: [pic], где ?=?/h.

Анализ показывает, что погрешность аппроксимации схем есть [pic].

Список литературы

1. Демидович Б.П., Марон И.А. Основы вычислительной математики. Наука,

1970.

2. Минкова Р.М., Вайсбурд Р.А. Методы вычислительной математики. УПИ,

1981.

3. Боглаев Ю.П. Вычислительная математика и программирование. Высшая

школа, 1990.

4. Кацман Ю.Я. Прикладная математика. Численные методы. ТПУ, 2000.

-----------------------

0 x

0 f(a) x

Рис. 1. Рис. 2.

Рис. 3.

0 f(a) x2 x1 x

a ?

Рис. 4.

y = P(x)

y = f(x)

x0 x1 xn x

b=x0 x1’

МЕНЮ

Некоторые дополнительные вычислительные методы

Некоторые дополнительные вычислительные методы

ИНТЕРЕСНОЕ