Нахождение полиноминальной аппроксимации методом наименьших квадратов

Московский Авиационный Институт

(государственный технический университет)

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

«теории вероятностей и математической статистике»

на тему:

Нахождение полиноминальной аппроксимации методом наименьших квадратов

Вариант №2

Выполнила: Студентка  группы 05-202

Андреева Виктория

Принял:     Преподаватель кафедры 804

Молчанов Игорь Иванович

Москва

2010 г.


ЗАДАНИЕ (вариант № 2) :  Даны результаты измерений случайного процесса в равноотстоящие моменты времени (реализация временного ряда).

13,393 13,207 13,477 11,911 14,311 14,979 14,437 14,957 13,044 12,142

12,000 11,496 12,927 11,849 11,612 10,401 8,755 8,185 9,681 9,644

9,073 8,535 9,062 7,602 9,164 6,913 7,749 5,543 5,901 5,901

6,760 4,593 6,131 3,651 3,796 3,663 3,068 3,008 2,809 0,333

1,730 -0,072 0,479 -3,180 -2,962 -5,849 -6,153 -7,911 -10,134 -11,662

Измерения производятся с шагом по аргументу 0,08

ЦЕЛЬ РАБОТЫ:  Требуется найти полиноминальную аппроксимацию этого процесса методом наименьших квадратов.


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

Математическая статистика – наука о математических методах, позволяющих по статистическим данным, например по реализацией случайной величины (СВ), построить теоретико-вероятностную модель исследуемого явления. Задачи математической статистики являются, в некотором смысле, обратными к задачам теории вероятностей. Центральным понятием математической статистики является выборка.

Выборка

Однородной выборкой (выборкой) объема n при  называется случайный вектор , компоненты которого , называемые элементами выборки, являются независимыми СВ с одной и той же функцией распределения F(x). Будем говорить, что выборка соответствует функции распределения F(x).

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

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

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

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

Выборочное пространство может быть всем -мерным евклидовым пространством  или его частью, если СВ X непрерывна, а также может состоять из конечного или счетного числа точек из , если СВ X дискретна.

На практике при исследовании конкретного эксперимента распределения  СВ  редко бывают известны полностью. Часто априори (до опыта) можно лишь утверждать, что распределение  случайного вектора  принадлежит некоторому классу (семейству) .

Пара (S,F) называется статистической моделью описания серии опытов, порождающих выборку .

Если распределение  из класса F определены с точностью до некоторого векторного параметра , то такая статистическая модель называется параметрической и обозначается .

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

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

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


Кривая регрессии.

регрессия вероятность статистический опыт

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

Точечная оценка.

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

Оценка  параметра  называется несмещенной, если ее МО равно , т. е.  для любого .

Оценка  параметра  называется состоятельной, если она сходится по вероятности к , т. е.  при  для любого .

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

Очевидно, что если оценка сильно состоятельная, то она является также состоятельной.

Доверительный интервал.

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

Пусть для параметра  получена из опыта несмещенная оценка . Назначим некоторую достаточно большую вероятность  (например,  или 0,99) такую, что событие с вероятностью  можно считать практически достоверным, и найдем такое значение , для которого

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

рис 1

Вероятность  принято называть доверительной вероятностью, а интервал - доверительным интервалом. Границы интервала :  и  называются доверительными границами.

Интервальные оценки.

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

Интервал  со случайными концами, «накрывающий» с вероятностью , , неизвестный параметр , т. е.

,

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

Число  называется доверительной вероятностью или уровнем доверия.

Уровень значимости.

Уровнем значимости статистического критерия называется вероятность ошибки 1-го рода . Вероятность ошибки 1-го рода  может быть вычислена, если известно распределение .

Ошибки 1 и 2-го рода.

Ошибкой 1-го рода называется событие, состоящее в том, что гипотеза  отвергается, когда она верна.

Ошибкой 2-го рода называется событие, состоящее в том, что принимается гипотеза , когда верна гипотеза .

Проверка статистических гипотез.

Статистической гипотезой H или просто гипотезой называется любое предположение относительно параметров или закона распределения СВ , проверяемое по выборке .

Проверяемая гипотеза называется основной (или нулевой) и обозначается . Гипотеза, конкурирующая с , называется альтернативной и обозначается .

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

Статистическим критерием (критерием согласия, критерием значимости или решающим правилом) проверки гипотезы  называется правило, в соответствии с которым по реализации  статистики  гипотеза  принимается или отвергается.

Критической областью статистического критерия называется область реализаций  статистики , при которых гипотеза  отвергается.

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


Практическая часть.

Этап 1 (Вычисление оценок ,  неизвестных коэффициентов регрессии , ):

;

;

 - оценка полезного сигнала (кривая регрессии);

 - ошибка;

Формулируем все ошибки:

.

Находим наименьшую ошибку. Для этого продифференцируем уравнение по a и по , приравняем к 0, получив систему:

 - система нормальных уравнений.

Решаем систему методом Крамера:


Расчетная схема для оценок  по методу наименьших квадратов.

НомерYXy^2X*Yx^2δ^2=(y-at-b)^2
113,393-2179,37245-26,786484,52154547
213,207-1,92174,42485-25,35743,686477,33345969
313,477-1,84181,62953-24,79773,385663,01706699
411,911-1,76141,87192-20,96343,097679,54344995
514,311-1,68204,80472-24,04252,822435,20165139
614,979-1,6224,37044-23,96642,5621,89755599
714,437-1,52208,42697-21,94422,310421,49126227
814,957-1,44223,71185-21,53812,073612,46267515
913,044-1,36170,14594-17,73981,849623,59662672
10112,142-1,2812575,828-143,5421,63848991,966406
1112-1,2144-14,41,4422,3767305
1211,496-1,12132,15802-12,87551,254421,61124277
1312,927-1,04167,10733-13,44411,08166,928339322
1411,849-0,96140,3988-11,3750,92169,762864991
1511,612-0,88134,83854-10,21860,77447,705858746
1610,401-0,8108,1808-8,32080,6411,56902772
178,755-0,7276,650025-6,30360,518419,90687251
188,185-0,6466,994225-5,23840,409619,76777254
199,681-0,5693,721761-5,421360,31365,590769531
209,644-0,4893,006736-4,629120,23043,29736681
219,073-0,482,319329-3,62920,163,244500827
228,535-0,3272,846225-2,73120,10243,075233208
239,062-0,2482,119844-2,174880,05760,410905071
247,602-0,1657,790404-1,216320,02562,296447057
259,164-0,0883,978896-0,733120,00640,399692323
266,913047,789569001,067444888
279,7490,0895,0430010,779920,00645,704661232
285,5430,1630,7248490,886880,0256Определение дуальных и двойных чисел


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


Расчет показателей надежности и законов их распределения


Різницевий метод розв'язання крайових задач для звичайних диференціальних рівнянь


Фрактальна розмірність


Актуально: