Решение линейных уравнений. Контрольная работа. Математика.
⚡ 👉🏻👉🏻👉🏻 ИНФОРМАЦИЯ ДОСТУПНА ЗДЕСЬ ЖМИТЕ 👈🏻👈🏻👈🏻
Вы можете узнать стоимость помощи в написании студенческой работы.
Помощь в написании работы, которую точно примут!
Похожие работы на - Решение линейных уравнений
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Скачать Скачать документ
Информация о работе Информация о работе
Нужна качественная работа без плагиата?
Не нашел материал для своей работы?
Поможем написать качественную работу Без плагиата!
1.1 Основные виды линейных интегральных уравнений
Под интегральным уравнением понимается уравнение, содержащее
неизвестную функцию у (х) под знаком определенного интеграла. В дальнейшем мы
ограничимся рассмотрением линейных интегральных уравнений, в которые
неизвестная функция входит лишь в первой степени (линейно).
Приведем некоторые наиболее часто встречающиеся типы
интегральных уравнений. Уравнение вида
где (ядро) и f(х) - известные функции, называется интегральным уравнением
Фредгольма первого рода. Уравнение вида
где - числовой параметр, носит название
интегрального уравнения Фредгольма второго рода.
Параметр вводится по следующим соображениям: при
данном значении интегральное уравнение (2) не всегда
имеет решения. Варьируя параметр , можно добиться того, чтобы решение уравнения (2) существовало.
Параметр можно также ввести в левую часть
уравнения Фредгольма первого рода (1).
Если в (2) f(х) 0, то получается однородное уравнение
допускается нулевое (тривиальное) решение у≡0. Те
значения параметра λ, при которых однородное
интегральное уравнение (3) имеет нетривиальные решения, называются собственными
значениями (собственными числами) ядра К ( x , s ) , или соответствующего
уравнения (2), а отвечающие им нулевые решения - собственными функциями.
Основной результат теории следующий (теорема Фредгольма): 1) если λ не есть собственное значение ядра К ( x , s ) , то соответствующее
неоднородное интегральное уравнение Фредгольма (2) с регуляторным ядром К ( x , s ) и непрерывным свободным
членом f ( x ) имеет единственное непрерывное решение у(х)
(а≤х≤ b ) , 2), если же λ есть собственное значение, то уравнение (2) или не имеет решений,
или же допускает бесчисленное множество их.
В приложениях важную роль играют уравнения Фредгольма второго
рода с симметричным ядром К ( x , s ) , т.е. таким, что
Симметрическое ядро обладает следующими свойствами:
· для всякого симметрического ядра
существует по меньшей мере одно собственное значение;
· все собственные значения симметрического
ядра действительны;
· собственные функции (х) и (х) симметрического ядра, соответствующие различным собственным
значениям и ( ) ,
ортогональны между собой на основном промежутке (а, b ), т.е.
Пример 1. Пусть простой замкнутый кусочно-гладкий контур
ограничивает конечную область G. Тогда функция u (х, у), дающая решение
соответствующей задачи Дирихле т.е. такая, что
(f, (t) - известная функция), может быть представлена в виде
и функция удовлетворяет интегральному уравнению
В приложениях встречаются также интегральные уравнения вида
которые носят названия интегральных уравнений Вольтерра
соответственно первого и второго рода. Вводя функцию
уравнения Вольтерра (4) и (5) можно записать в виде
соответствующих уравнений Фредгольма с ядром . Таким образом, теория уравнений Вольтерра сводится к теории
уравнений Фредгольма; однако в некоторых случаях уравнения Вольтерра полезно
изучать независимо.
Примером уравнения Вольтерра первого рода является обобщенное
уравнение Абеля
где - известная непрерывно дифференцируемая
функция. Решение уравнения (6) дается формулой
в чем можно убедиться непосредственно.
Заметим, что если ядро и - непрерывно дифференцируемые функции,
причем 0 при , то уравнение Вольтерра первого рода (4)
сводится к уравнению Вольтерра второго рода (5). Действительно, дифференцируя
уравнение (4) по x, будем иметь
К линейным интегральным уравнениям может быть приведено
большое количество задач математической физики.
В основном мы это будут интегральное уравнение Фредгольма
второго рода и частично уравнение Вольтерра второго рода (для краткости в
дальнейшем мы их будем именовать просто уравнениями Фредгольма и Вольтерра).
Основными проблемами здесь являются следующие:
1. нахождение приближенного или точного решения неоднородного
интегрального уравнения при заданном значении параметра ;
. нахождение собственных значений и соответствующих
собственных функций однородного интегрального уравнения.
1.2 Метод последовательных приближений
Будем искать решение в форме степенного ряда
Подставляя выражение (2) в интегральное уравнение (1) и
приравнивая коэффициенты при одинаковых степенях параметра , будем иметь
Из формул (3) по индукции получаем Поэтому сходимость ряда (2) будет обеспечена, если
мы получим приближенное решение интегрального уравнения (1) с
погрешностью
Формула (2) дает аналитическое относительно решение уравнения Фредгольма (1) в
окрестности точки = 0. Из формул (3) вытекает, что решение
(2) можно
Коэффициенты , так называемые итерированные ядра, могут
быть последовательно найдены по формулам
Функция называется резольвентой уравнения (1) и
при малых определяется степенным рядом (6). Пользуясь
аналитическим продолжением, резольвенту можно продолжить на всю комплексную плоскость параметра , за исключением собственных значений (особые точки), которые являются
полюсами резольвенты. Тогда формула (5) дает решение интегрального
уравнения (1) при любом
Рассмотрим теперь соответствующее уравнение Вольтерра
Следовательно, ряд (8) сходится при любом и дает единственное решение уравнения
(7). Погрешность приближенного решения
1.3 Решение интегрального уравнения методом
конечных сумм
Метод основывается на приближенном вычислении определенного
интеграла с помощью некоторой квадратурной формулы
где - абсциссы точек отрезка - числовые коэффициенты, не зависящие от
выбора функции , и - остаточный член (ошибка) формулы (1).
Например, в случае равноотстоящих точек , где , будем иметь:
3. для общей формулы
Симпсона при :
Пусть теперь дано интегральное уравнение Фредгольма второго
рода
На основании формулы (1) будем иметь
где - соответствующие ошибки. Отбрасывая в
системе (3) величины для приближенных значений решения в узлах получим линейную алгебраическую систему
то система (4’) имеет единственное решение которое можно найти методом Гаусса или
другими методами, разработанными для решения систем алгебраических линейных
уравнений.
Найдя , для решения получаем из уравнения (2) приближенное
аналитическое выражение
Различные между собой корни алгебраического уравнения представляют собой, вообще говоря, приближения собственных
значений ядра . Если - соответствующие ненулевые решения однородной системы
Метод конечных сумм может быть применен также к интегральному
уравнению Фредгольма первого рода
В этом случае приближенные значения решения в узлах будут определяться из системы
Особенно просто применение метода конечных сумм для решения
интегрального уравнения Вольтерра второго рода
которое можно рассматривать как уравнение Фредгольма второго рода.
Здесь = 0 при j>i, и, следовательно, соответствующая
система (4) имеет вид
Получилась линейная система с треугольной матрицей. Если
Будем искать приближенное решение уравнения (1) в виде
функции определенного вида
со свободными параметрами (неопределенные
коэффициенты). Подставляя выражение (2) в уравнение (1), получим невязку
Если у является точным решением, то, очевидно, невязка =0. Поэтому стараются подобрать параметры так, чтобы n невязка была в определенном смысле возможно
малой. Минимизировать невязку можно различными способами. Обычно для простоты выкладок берут
функцию , линейно зависящую от параметров . Найдя параметры , получают приближенное решение (2).
Отметим одно обстоятельство. Если невязка получилась малой, то она близка к невязке
=0, даваемой точным решением y. Но если два оператора и R[y] имеют
близкие значения, то отсюда, вообще говоря, не следует, что функции и у близки между собой в обычном смысле.
Поэтому возникает математическая задача: по известной невязке оценить погрешность приближенного решения . Но это трудная проблема, связанная с
глубокими теоремами функционального анализа, и мы ее оставим без рассмотрения.
Если справедливо равенство (4), то данным методом решение у можно
найти с любой степенью точности, взяв достаточно большое число параметров .
Перейдем теперь к изложению одного из конкретных методов
построения приближенного решения .
где - известные функции (координатные функции)
и - неопределенные коэффициенты, причем
функции линейно независимы. Заметим, что, в
частности, можно полагать . Подставляя выражение (5) в левую часть
уравнения (1), получим невязку
Согласно методу коллокации требуем, чтобы невязка обращалась в нуль в заданной системе
точек из отрезка [а, b] (точки коллокации),
т.е. полагаем, что
Отсюда на основании формулы (6) для определения коэффициентов получаем алгебраическую линейную систему
уравнений
Если определитель этой системы , то из системы можно однозначно определить величины и, следовательно, найти приближенное
решение по формуле (5).
Приравнивая нулю определитель , получим уравнение =0, из которого, вообще говоря, можно найти приближенные значения первых собственных чисел ядра .
Если положить то вместо нашей системы будем иметь однородную систему
Найдя ненулевые решения ( ) системы (8), получим для ядра приближенные собственные функции
отвечающие его собственному значению .
где - известные функции и -неопределенные коэффициенты, причем линейно независимы.
Подставляя (2) в левую часть уравнения (1), получим невязку
Заметим, что если , то и, следовательно, . Так как то матрица системы (8) симметрическая.
Вместо интегрального метода наименьших квадратов можно воспользоваться точечным
способом наименьших квадратов.
Метод наименьших квадратов применяется также для приближенного
нахождения собственных значений и собственных функций ядра , аналогично тому как это делается для
метода коллокации. А именно, полагая f(x) 0 и , откуда , определяем приближенные значения
собственных чисел из алгебраического уравнения.
Аналогично предыдущему, будем искать приближенное решение
уравнения (1) в виде конечной суммы
интегральный уравнение линейный коллокация
где -некоторые известные линейно независимые
функции (координатные функции) и - неопределенные коэффициенты. Подставляя выражение (2) в левую
часть уравнения (1), получим невязку
Согласно методу моментов коэффициенты определяются из условия ортогональности
невязки ко всем координатным функциям . Это дает систему уравнений
Если определитель системы (4) отличен от нуля, то из этой системы можно однозначно определить
коэффициенты . Тогда формула (2) даст приближенное
решение интегрального уравнения (1). Из уравнения =0 приближенно находятся собственные
значения ядра К (х, s).
Найдя ненулевые решения однородной линейной системы
легко построить приближенные собственные функции , отвечающие данному собственному значению . Заметим, что метод моментов по идее
совпадает с методом Галеркина.
Пример 1. Найти первые два собственных значения интегрального
уравнения
Ортогонализируя невязку R[Y], будем иметь систему
Приравнивая пулю определитель системы (6), получим уравнение
дня определения собственных значений:
полученные из решения соответствующей краевой задачи:
Таким образом, погрешность равна примерно ,
В заключение заметим, что рассмотренные методы минимизации
невязки применимы также к решению нелинейных интегральных уравнений. Имеются
также другие методы решения интегральных уравнений, например метод Монте-Карло.
Решить интегральное уравнение методом конечных сумм и методом
моментов. В методе моментов использовать функции , k = 0, 1, 2,…, n:
где , а = 0, в = 1, f(x) = 1/8 (x + 8), ядро К = 1 + .
1) задаемся числом разбиений n = 6, тогда
2) выбираем равноотстоящие узлы
x 1 = 0, x 2 = 0,2, x 3
= 0,4 x 4 =0,6, x 5 =0,8, x 6 =1, t 1 = 0,
t 2 = 0,2, t 3 = 0,4, t 4 =0,6, t 5 =0,8,
t 6 =1.
3) Рассчитываем значения К ij = K(x i , t j ), i, j, =
) Рассчитываем значения f i = f(x i ), i = 1, n
) Решаем систему относительно i = 1, 6
) Приближенное решение находим по формуле
) Также находим приближенное решение для числа разбиений n = 4 и n = 6
10) Находим погрешность по формуле
11) Находим погрешность по формуле
) Сравнивая и , делаем вывод: если < , то задача решена правильно.
В результате выполнения курсовой работы мною был изучен
теоретический материал по решению линейных интегральных уравнений. Были
подробно рассмотрены такие численные методы решения линейных интегральных
уравнений как метод последовательных приближений, метод конечных сумм, метод
вырожденных ядер, метод коллокации, метод наименьших квадратов и метод моментов.
1.
Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение
функций, дифференциальные и интегральные уравнения. М.: Наука, 1967
.
Калиткин Н.Н. Численные методы. М: Наука, 1978.
.
Сборник задач по математике для втузов. В 4-х частях. Под ред. Ефимова А.В.,
Поспелова А.С. (2001-2003)
.
Турчак Л.И. Основы численных методов. М: Физматлит, 2002.
.
Кунцман Ж. Численные методы. М.: Наука, 1979
Похожие работы на - Решение линейных уравнений Контрольная работа. Математика.
Вегетативная Нервная Система Ликвор Неврология Реферат
Создание Собственного Дела Реферат
Дипломная работа по теме Товароведная характеристика, экспертиза качества обуви и организационно-экономические условия её реализации
Реферат: Wasting Time Essay Research Paper Wasting Death
Контрольная Работа По Истории 1945 1964
Курсовая работа: Сравнительная оценка качества растительных масел отечетвенного и импортного производства. Скачать бесплатно и без регистрации
Реферат: День взятия Бастилии
Реферат: Закон РФ о таможенном тарифе. Скачать бесплатно и без регистрации
Умственное Развитие Младшего Школьника Курсовая Работа
Товар И Его Свойства Реферат По Экономике
Реферат Технология Консервирования Мяса Птицы
Курсовая работа по теме Титул главы государства, способы передачи власти, права верховной власти с XVIII века по 1917 год
Реформа генриха ii плантагенета
Реферат: Особо охраняемые территории Урала
Реферат На Тему А.С. Пушкин Как Прозаик, Драматург, Историк
Дипломная работа по теме Вопросы воспитания средствами театральной педагогики
Курсовая По Транспортной Безопасности
Реферат Архитектура
Курсовая работа по теме Социально-демографический состав потребителей, типологии потребителей в СО сервисе
Дипломная работа по теме Программирование на языке С#
Курсовая работа: Логистическая оптимизация деятельности предприятия на примере ООО Совлит
Похожие работы на - Разработка укрупненного технологического процесса
Похожие работы на - Пути продвижения бренда компании на материалах ОАО 'Винап'