ЛЕКЦИЯ № 12

ТИПОВЫЕСПОСОБЫ ПРИБЛИЖЕНИЯ ПРИ ПОМОЩИ МНК НЕЛИНЕЙНЫХ ФУНКЦИЙ В ЛИНЕЙНЫЕ

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

Аппроксимируемая Линейная Замена

функция функция

Вообще полиномы выше 6-ой степени при помощи МНК никогда не строят, т.к. появляются серьёзные ошибки округлений и раскачивания. На практике ограничиваются квадратической зависимостью.

МНК для системы линейно- независимых функций.

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

. Под линейно-независимой функцией понимаем такую систему, в которой ни одна из функций не может быть представлена в виде линейной комбинации остальных функций.

Задана f(x) таблично на [a;b] по системе узлов x_j,y_j=f(x_j)

Рассмотрим приближение f(x) при помощи обобщенного многочлена:

(12.1)

Необходимо найти неизвестные коэффициенты из (12.1)

(12.2)

Критерий (12.2) представляет собой квадратичную функцию относительно параметров b_i.

Запишем

(12.3)

Получим

(12.4)

Система (12.4) представляет собой СЛАУ относительно параметров b_i и может быть решена одним из известных методов.

Рассмотрим один из частных случаев этой системы, когда функции

являются ортогональными.

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

(12.5)

Линейно-независимая система функций

является ортогональной если

Для системы ортогональных функций решение системы (12.4) получается элементарно.

(12.6)

Коэффициенты (12.6) называются коэффициентами Фурье, а многочлен (12.1) называется обобщенным многочленом Фурье.

Тригонометрические ряды и полиномы Фурье в использовании МНК

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

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

f(x+KP)=f(x)

P-наименьший положительный период.

Пусть g(x) имеет P, тогда f(x)=g(Px/2π) будет иметь период 2π.

Пусть f(x) –функция, имеющая период 2π, тогда она может быть представлена рядом:

(12.7)

(12.8)

(12.8) тригонометрический ряд Фурье.

(12.9)

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

Пусть значения таблично-заданной функции известны в точках

Тригонометрические полиномы используются для тригонометрических процессов.

ЛЕКЦИЯ №13

ЧИСЛЕННЫЕ ДИФФИРИНЦИРОВАНИЕ И ИНТЕГРИРОВАНИЕ

Численное дифференцирование

1.1 Расчет производных аналитически заданной функции

1.2 Нахождение производных таблично заданной функции

Численное интегрирование

2.1 Формулы прямоугольников

2.2 Формулы Ньютона - Котеса

Формулы Симпсона и Ньютона

2.3 Формулы Чебышева и Гаусса

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

§1.1 Для аналитически заданных функций рассмотрим следующие способы численного дифференцирования:

1.) предел отношения приращений;

2.) при помощи центрированных разностях;

Предел отношения приращений

Строим последовательности {h_k}так, чтобы h_k→0 вычисляем предел последовательности {D_k}, где

D_k=

k=1,2..n

Вычисления проводят до некоторого n, при котором выполняется условие:

|D_n₊₁-D_n|≥|D_n-D_n_-1|

Шаг выбираем сами (обосновать).

Центрированные разности

Пусть наша функция трижды непрерывно дифференцируема на [a;b]:

c³[a;b] x-h, h, x

[a;b]

тогда

Эта приближенная формула имеет 2-ой порядок точности.

Ошибка: 0(h²). Разложим функцию в ряд Тейлора:

Вычтем из первого равенства второе:

f(x+h)-f(x-h)

ЛЕКЦИЯ №14

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

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

1) Формулы прямоугольника

Пусть y=f(x) непрерывна на [a,b]. Требуется вычислить

Разобьем отрезок интегрирования на n равных частей, точками x_i, i=0,n

x_i=a-i*h

шагразбиения

На отрезке [x_i_-1;x_i] возьмем произвольную точку x_i. Из определения наш интеграл равен

(14.1)

общая форма прямоугольника.

Геометрическая интерпритация формулы прямоугольника.

Площадь ограничена графиком функции y=f(x) на отрезке осью абцис и заменяется площадью прямоугольника с высотой равной f(x_i). Отметим частный случай формулы прямоугольника.

1) Пусть x_i это x_i. Из формулы 14.1 видно

Точность 0(h) порядка один. Формула левых прямоугольников

2) x_i=x_i.

формула правых прямоугольников 0(h).

3) x_i=1/2 (x_i_-1+ x_i). Получилась формула средних прямоугольников.

. Точность порядка два 0(h²).

2)Формула Ньютона-Котеса

-общая формула

Н-коэффициент Котеса.

x_k=a+b*hk=0,n

Пусть n=1 значит

Пусть n=2 значит

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

y_i=f(x_i), i=0,n. Точность порядка два 0(h²). Эта формула точна для многочленов первой степени.

Геометрическая интерпритация

На элементарном отрезке [x_i_-1;x_i] площадь под кривой полагают равной площади трапеции с основаниями y_i_-1и y_i и высотой h. h= x_i-x_i_-1

Формула Симсона

Пусть n=2m. Число разбиения отрезков четное. Тогда

точна для многочлена в третьей степени.

Численные методы (стр. 1 из 2)

ЛЕКЦИЯ №13

Геометрическая интерпритация

Формула Симсона