Аппроксимация МНК

Автор работы: Пользователь скрыл имя, 14 Октября 2009 в 09:59, Не определен

Описание работы

Аппроксимация функции методом наименьших квадратов является простой и легко реализуемой как в среде Pascal, так и в MathCAD. МНК «сглаживает» функцию, выбирая промежуточные значения, что является выгодным решением.

Файлы: 1 файл

Оглавление.doc

— 934.50 Кб (Скачать файл)

Оглавление 

 

Введение

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

     Как известно, между величинами может  существовать точная (функциональная) связь, когда одному значению аргумента  соответствует одно определенное значение.

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

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

     В данной курсовой работе программа на языке Pascal разработана в оболочке PascalABC версия 1.0 beta.  Решение задачи в среде MathCad производили в Mathcad версия 14.0.0.163.

 

Теоретическая часть

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

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

  1. Разработать алгоритм нахождения коэффициентов трёх аппроксимирующих полиномов (многочленов) вида для табулированной функции y=f(x):
x 0,15 -0,01 0,46 1,71 3,94 7,4 12,41 19,47 29,21 42,49
f(x) 7,26 7,95 8,77 9,72 10,78 12,45 14,74 18,18 23,36 30,99
x 60,43 84,42 116,18 157,8 211,78 281,04 369 479,59 617,3 787,21
f(x) 41,97 57,38 78,45 106,65 143,64 191,28 251,71 327,27 420,58 534,52
x 995,04 1247,17 1550,71 1913,5 2344,18 2852,21 3447,91 4142,52 4948,19 5878,09
f(x) 672,24 837,2 1033,14 1264,13 1534,54 1849,11 2212,91 2631,36 3110,26 3665,8
x 995,04 1247,17 1550,71 1913,5 2344,18 2852,21 3447,91 4142,52 4948,19 5878,09
f(x) 672,24 837,2 1033,14 1264,13 1534,54 1849,11 2212,91 2631,36 3110,26 3665,8
x 6946,36 8768,24 9560,03 11139,19 12924,34 14935,29 17193,14 19720,24 22540,29 25678,32
f(x) 4274,55 4973,5 5760,63 6641,98 7627,6 8725,62 9945,22 11296,05 12788,24 14432,44

для степени полиномов n=2, 4, 5.

  1. Построить блок-схему алгоритма.
  2. Создать программу на языке Pascal, реализующую разработанный алгоритм.
  3. Рассчитать среднеквадратичные отклонения для каждого из трех случаев по формуле: 
  4. Построить графики 3-х полученных приближающих функций в одной системе координат. На графике должны содержаться и исходные точки i,yi).
  5. Решить задачу средствами MathCAD.

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

Построение  эмпирических формул  методом наименьших квадратов

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

     При аналитическом исследовании взаимосвязи  между двумя величинами x и y производят ряд наблюдений и в результате получается таблица значений:

x ¼ ¼
y ¼ ¼

     Эта таблица обычно получается как итог каких-либо экспериментов, в которых  (независимая величина) задается экспериментатором, а получается в результате опыта. Поэтому эти значения будем называть эмпирическими или опытными значениями.

     Между величинами x и y существует функциональная зависимость, но ее аналитический вид  обычно неизвестен, поэтому возникает практически важная задача - найти эмпирическую формулу

                                                   (1)

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

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

     Если  в эмпирическую формулу (1) подставить исходные , то получим теоретические значения , где .

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

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

                                                 (2)

будет минимальной.

     Поясним геометрический смысл метода наименьших квадратов.

     Каждая  пара чисел  из исходной таблицы определяет точку на плоскости . Используя формулу (1) при различных значениях коэффициентов можно построить ряд кривых, которые являются графиками функции (1). Задача состоит в определении коэффициентов таким образом, чтобы сумма квадратов расстояний по вертикали от точек до графика функции (1) была наименьшей.

     Построение  эмпирической формулы состоит из двух этапов: выяснение общего вида этой формулы и определение ее наилучших параметров.

     Если  неизвестен характер зависимости между данными величинами x и y, то вид эмпирической зависимости является произвольным. Предпочтение отдается простым формулам, обладающим хорошей точностью. Удачный выбор эмпирической формулы в значительной мере зависит от знаний исследователя в предметной области, используя которые он может указать класс функций из теоретических соображений. Большое значение имеет изображение полученных данных в декартовых или в специальных системах координат (полулогарифмической, логарифмической и т.д.). По положению точек можно примерно угадать общий вид зависимости путем установления сходства между построенным графиком и образцами известных кривых.

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

     Для того, чтобы найти набор коэффициентов  , которые доставляют минимум функции S , определяемой формулой (2), используем необходимое условие экстремума функции нескольких переменных - равенство нулю частных производных.  В результате получим нормальную систему для определения коэффициентов :

                                           (3)

     Таким образом, нахождение коэффициентов  сводится к решению системы (3).

     Эта система упрощается, если эмпирическая формула (1) линейна относительно параметров , тогда система (3) - будет линейной.

     Конкретный  вид системы (3) зависит от того, из какого класса эмпирических формул мы ищем зависимость (1). В случае линейной зависимости система (3) примет вид:

                           (4)

     Эта линейная система может быть решена любым известным методом (методом  Гаусса, простых итераций, формулами  Крамера).

     В случае квадратичной зависимости  система (.3) примет вид:

                                            (5)

 

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

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

     При этом матрица СЛАУ приводится  треугольному виду, где ниже главной диагонали  располагаются только нули.

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

     Имеем Ax=b  (6), где A=[aij] - матрица размерности n*n, det A 0, b=(b1, b2, …, bn)T.

     В предположении, что a11 0, первое уравнение системы (6):

делим на коэффициент a11, в результате получаем уравнение

.

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

Информация о работе Аппроксимация МНК