Referat.me

Название: Метод наименьших квадратов в случае интегральной и дискретной нормы Гаусса

Вид работы: реферат

Рубрика: Математика

Размер файла: 212.25 Kb

Скачать файл: referat.me-215108.docx

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

Метод наименьших квадратов в случае интегральной и дискретной нормы Гаусса

Метод наименьших квадратов в случае интегральной и дискретной нормы Гаусса

1. Постановка задачи

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

Например: может быть, в частности, и непрерывной функцией на , а соответствующая - алгебраическим или тригонометрическим многочленом, который «достаточно хорошо» приближает функцию .

Например: всякую функцию из можно представить приближённо соответствующим многочленом степени с помощью формулы Тейлора:

(1)

т.е.

; (2)

где , - многочлен степени , приближающий функцию , - остаточный член. Ясно, что

(3)

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

Известно также, что можно приблизить с помощью тригонометрического многочлена – отрезка ряда Фурье.

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

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

Если - многочлен степени , то рассматриваемый процесс приближения называется параболическим интерполированием или процессом построения интерполяционного многочлена (частным примером является многочлен Лагранжа, т.е. );

б) функцию можно выбрать так, чтобы норма - отклонения невязки – достигала минимального значения, причём норма может быть определена по-разному, и разным нормам соответствуют различные степени приближения.

В функциональном пространстве Гильберта , норме невязки имеет вид (интегральная норма Гаусса):

(4)

часто, в качестве нормы рассматривают Чебышевскую норму (Т – первая буква фамилии Чебышева на немецком языке):

(5)


При использовании нормы (5) говорят о равномерном приближении функции , функцией .

Подробная теория Т-приближений была развита в работах немецкого математика Л. Коллатца.

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

(6)

Ясно, что метод наименьших квадратов (6) – является дискретным аналогом функции Гаусса (4).

Принципиальную возможность приближения любой непрерывной функции многочленом даёт теорема Вейерштрасса: Если , тогда , - многочлен, что имеет место неравенство:

(7)

2. Метод наименьших квадратов в случае приближения функции

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

Как известно, для линейной независимости системы функций необходимо и достаточно, чтобы определитель Грама этой системы был отличен от нуля, т.е.

(8)

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

(9)

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

2.1 Квадратичное приближение таблично заданной функции по дискретной норме Гаусса

Рассмотрим задачу приближения функции в случае использования невязки в форме (6). Т.е. используем дискретную норму Гаусса:

(10)

где неизвестная функция аппроксимируется функцией из (9). Для известны лишь значения в различных точках , т.е. , где . Таким образом, для определения имеем задачу: найти точку минимума - невязки функции Гаусса - для таблично заданной функции , если

, (где ). (11)

Очевидно, что условия минимума дискретной функции невязки Гаусса - имеют вид:

, (12)

Эти условия для (11) преобразуются к виду:

, (13)


Раскрывая систему (13) получаем систему уравнений для определения коэффициентов разложения в виде:

(14)

Нетрудно увидеть, что вводя скалярные произведения в соответствующем функциональном пространстве в виде:

(15)

систему (14) можно переписать в нормальном виде Гаусса:

(16)

Ясно, что эта система имеет единственное решение, т.к. определитель системы (16) совпадает с определителем

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

Найдя из системы (16) и подставляя в (9) мы получаем функцию:


(17)

которая является приближением к функции в смысле минимума квадратичного отклонения Гаусса (10) по норме индуцированной скалярным произведением (15), действительно:

(18)

а дискретная норма Гаусса невязки имеет вид:

(19)

2.2 Интегральное приближение функции заданной аналитически

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

Рассмотрим теперь случай, когда аналитически заданную, на интервале , функцию - надо аппроксимировать обобщённым многочленом:

(20)


так, чтобы минимизировалась интегральная норма невязки Гаусса :

(21)

иначе говоря, нам нужно минимизировать интеграл

(22)

Для решения этой задачи подставим (20) в (22), тогда функционал (22) превратится в функцию многих переменных, т.е. . Условия же минимума функции многих переменных имеют вид:

, (23)

Эти условия приобретают вид:

(24)

т.е.

(25)


Определитель этой системы представляет собой определитель Грама для функций , в , поэтому система (25) имеет единственное решение . Подставляя эти значения в разложение (20) имеем приближение для . Характер приближения оценивается соответствующей нормой невязки .

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

2.3 Числовые примеры на применение метода наименьших квадратов Гаусса для приближения функций заданных таблично или аналитически

а) Рассмотрим пример в случае табличного задания функции :

Пример 1: пусть функция задана таблично:

0.5 1.0 1.5 2.0 2.5 3.0
0.31 0.82 1.29 1.85 2.51 3.02

с помощью метода наименьших квадратов аппроксимировать эту функцию в классе линейных функций. Т.е. допускаем, что . Для нахождения коэффициентов , составляем невязку по дискретной норме Гаусса:

(26)

Необходимые условия минимума для имеют вид:


(27)

Из (27) – получаем нормальные уравнения Гаусса:

(28)

Решение имеет вид:

(29)

т.е.

(30)

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

Пример 2: Функцию , заданную на интервале аппроксимировать линейной функцией , определив параметры и по методу Гаусса (используем интегральную норму невязки Гаусса).

Решение: интегральная норма невязки для данной функции имеет вид:

(31)


Необходимые условия минимума для - имеют вид:

(32)

т.е.

(33)

(34)

Из уравнений (33) и (34) находим

(35)

аппроксимирующий многочлен имеет вид:

(36)

или

(37)

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

Литература

1). К. Ректорис. Вариационные методы в математической физике и механике. Мир, М.,1995

2). С.Г. Михлин. Численная реализация вариационных методов, М., Наука, 1996

3). Л.А. Кальницкий, Д.А. Добротин, В.Ф. Жевердеев. Специальный курс высшей математики для втузов. М., ”Высшая математика”, 1996

4). Т. Шуп. Решение инженерных задач на ЭВМ. Мир, М., 1982

5). Л. Коллатц. Функциональный анализ и вычислительная математика. Мир, М., 1999

6). Р. Варга. Функциональный анализ и теория аппроксимации в численном анализе. Мир, М., 1994

7). Л. Коллатц, Ю. Альбрехт. Задачи по прикладной математике. Мир, М.,1998.

Похожие работы

  • Задачи и примеры их решения по теории вероятности

    Вариант 3. 1. Решите уравнение Решение По определению Тогда и уравнение принимает вид откуда получаем Так как m может быть только натуральным числом, то значение

  • Метод Гаусса

    Методические рекомендации по выполнению заданий методом гауса. Примеры выполнения заданий.

  • Аппроксимация функций

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

  • Численное интегрирование функций

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

  • Передаточная функция дискретной системы

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

  • Методы решения систем линейных уравнений

    1. Решение системы линейных уравнений методом Гаусса Задачи аппроксимации функции, а также множество других задач прикладной математики м вычислительной физики сводятся к задачам о решении систем линейных уравнений. Самым универсальным методом решения системы линейных уравнений является метод последовательного исключения неизвестных, называемый методом Гаусса.

  • Аппроксимация функции методом наименьших квадратов

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

  • Визуализация численных методов

    Общая характеристика и особенности двух методов решения обычных дифференциальных уравнений – Эйлера первого порядка точности и Рунге-Кутта четвёртого порядка точности. Листинг программы для решения обычного дифференциального уравнения в Visual Basic.

  • Линеаризация без метода наименьших квадратов

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

  • Использование пакетов прикладных программ для решения задач линейного программирования

    ГОУ ВПО Уфимский государственный авиационный технический университет Кафедра вычислительной математики и кибернетики КУРСОВАЯ РАБОТА по дисциплине