Обработка случайных выборок

Размещено на http://www.allbest.ru

Размещено на http://www.allbest.ru

Содержание

Введение

1. Обработка одномерной случайной выборки

1.1 Нахождение точечных оценок для не сгруппированной выборки

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

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

1.4 Расчёт критерия Пирсона

1.5 Расчёт критерия Колмогорова

2. Обработка двумерной случайной выборки

2.1 Построение поля рассеивания, гипотеза о виде корреляционной зависимости

2.2 Построение корреляционной таблицы

2.3 Расчёт коэффициентов уравнения прямой регрессии

2.4 Нахождение выборочного коэффициента корреляции

2.5 Расчёт коэффициентов уравнения криволинейной регрессии

2.6 Нахождение корреляционного отношения

2.7 Расчёт критерия Фишера

Заключение

Литература

Приложения

Введение

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

Например, для изучения физических явлений производят наблюдения или опыты. Их результаты обычно регистрируют в виде значений некоторых наблюдаемых величин. При повторении опытов мы обнаруживаем разброс их результатов. Например, повторяя измерения одной и той же величины одним и тем же прибором при сохранении определенных условий (температура, влажность и т.п.), мы получаем результаты, которые хоть немного, но все же отличаются друг от друга. Даже многократные измерения не дают возможности точно предсказать результат следующего измерения. В этом смысле говорят, что результат измерения есть величина случайная. Еще более наглядным примером случайной величины может служить номер выигрышного билета в лотерее. Можно привести много других примеров случайных величин. Все же и в мире случайностей обнаруживаются определенные закономерности. Математический аппарат для изучения таких закономерностей и дает теория вероятностей. Таким образом, теория вероятностей занимается математическим анализом случайных событий и связанных с ними случайных величин.

1 Обработка одномерной случайной выборки

1.1 Нахождение точечных оценок для не сгруппированной выборки

Математимческое ожидамние -- среднее значение случайной величины, распределение вероятностей случайной величины. Для несгруппированной выборки находится по формуле:

где - значение из выборки, - величина выборки.

Дисперсия - мера разброса данной случайной величины, то есть её отклонения от математического ожидания. Выборочная дисперсия равна:

Среднеквадратичное отклонение (СКО) - показатель рассеивания значений случайной величины относительно её математического ожидания.

Ассиметрия - величина, характеризующая асимметрию распределения данной случайной величины. Коэффициент асимметрии A статистического распределения определяется по формуле:

Коэффициент эксцесса - мера остроты пика распределения случайной величины. Эксцесс E статистического распределения определяется по формуле:

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

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

Количество интервалов, на которые необходимо разбить полученную выборку определяется по формуле: k = (1 + 3.32 * log₁₀ (n));

Данные о выборке Х:

Выборка: 100

Величина: 0,475

Минимум = 0,0191

Максимум = 3,82Количество интервалов: 8

Интервалы: 1) - 12; 2) - 24; 3) - 28; 4) - 23; 5) - 6; 6) - 4; 7) - 2; 8) - 0;

О выборке У:

Выборка: 100

Величина: 0,406

Минимум = -0,076

Максимум = 3,17

Количество интервалов: 8

Интервалы: 1) - 2; 2) - 2; 3) - 5; 4) - 13; 5) - 32; 6) - 28; 7) - 9; 8) - 8;

Математическое ожидание для сгруппированной выборки:

где где , - границы интервала, - вероятность попадания в этот интервал.

Дисперсия:

Среднеквадратическое отклонение:

Ассиметрия:

Эксцесс:

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

ДМ= |М_нс - М_с|, ДD= |D_нс - D_с|, Д ? = | ? _нс - ? _с|, ДA= |A_нс - A_с|, ДE= |E_нс - E_с|.

Таблица со значениями расчётов представлены на рисунке 1.1. и рисунке 1.2.

Рисунок 1.1 Таблица расчета точечных оценок для выборки Х



Рисунок 1.2 Таблица расчета точечных оценок для выборки У

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

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

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

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

Для проверки гистограммы, использовались следующие законы распределения.

Гауссовский нормальный:

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

Гаусовское стандартное распределение:

Экспоненциальный:

Сдвинутый экспоненциальный:

Вейбула:

Релея:

Равномерный:

Нормированный равномерный:

Проанализируем результаты, построенных гистограмм.

Выборка Х.

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

Рисунок 1.3 Гауссовский нормальный закон распределения для выборки Х

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

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



Рисунок 1.3 Гауссовский стандартный закон распределения для выборки Х

Из рисунка (рисунок 1.3. Гауссовский стандартный закон распределения) видно, что данный закон распределения не подходит.

И последнее подходящее распределение - это сдвинутое экспоненциальное.

Рисунок 1.3 Сдвинутый экспоненциальный закон распределения для выборки Х

График (рисунок 1.3. Сдвинутый экспоненциальный закон распределения для выборки Х) не совсем точно накладывается на гистограмму.

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

Выборка У.

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

Рассмотрим гауссовский нормальный закон.

Рисунок 1.3 Гауссовский нормальный закон распределения для выборки У

График распределения на рисунке 1.3. Гауссовский нормальный закон распределения для выборки У почти совпадает с гистограммой.

Рисунок 1.3 Гауссовский стандартный закон распределения для выборки У

Как видим из рисунка 1.3. Гауссовский стандартный закон распределения для выборки У гауссовский стандартный не подходит, как и в предыдущей выборке.

Рассмотрим сдвинутый экспоненциальный закон.

Рисунок 1.3 Сдвинутый экспоненциальный закон распределения для выборки У

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

1.4 Расчёт критерия Пирсона

Критерий согласия К. Пирсона , который можно представить как сумму отношений квадратов расхождений между f' и f к теоретическим частотам:

Вычисленное значение критерия _расч необходимо сравнить с табличным (критическим) значением _табл. Табличное значение определяется по специальной таблице, оно зависит от принятой вероятности Р и числа степеней свободы k (при этом k = m - 3, где m - число групп в ряду распределения для нормального распределения). При расчете критерия согласия Пирсона должно соблюдаться следующее условие: достаточно большим должно быть число наблюдений (n 50), при этом если в некоторых интервалах теоретические частоты < 5, то интервалы объединяют для условия > 5.

Если _{расч табл}, то расхождения между эмпирическими и теоретическими частотами распределения могут быть случайными и предположение о близости эмпирического распределения к нормальному не может быть отвергнуто.

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

К_Ром,

где m - число групп; k = (m - 3) - число степеней свободы при исчислении частот нормального распределения.

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

1.5 Расчёт критерия Колмогорова

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

где D - максимальное значение разности между накопленными эмпирическими и теоретическими частотами; - сумма эмпирических частот.

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

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

Расчёт критериев Пирсона и Колмогорова выборки Х:

Критерий Пирсона: 6,24

Критерий Колмагорова: 0,0894

Расчёт критериев Пирсона и Колмогорова выборки У:

Критерий Пирсона: 183

Критерий Колмагорова: 0,344

Вся информация о выборках также приведена на рисунке 1.5.Обработка однмерной выборки Х в приложении Б и рисунке Рисунок 1.5. Обработка однмерной выборки У в приложении В.

2. Обработка двумерной случайной выборки

2.1 Построение поля рассеивания, гипотеза о виде корреляционной зависимости

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

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

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

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

Если р = 0, то значения, x_i, y_i, полученные из двумерной нормальной совокупности, располагаются на графике в координатах х, у в пределах области, ограниченной окружностью (рисунок 2.1. Виды зависимости а). В этом случае между случайными величинами Х и Y отсутствует корреляция и они называются некоррелированными. Для двумерного нормального распределения некоррелированность означает одновременно и независимость случайных величин Х и Y.

Рисунок 2.1 Виды зависимости

Если р = 1 или р = -1, то между случайными величинами Х и Y существует линейная функциональная зависимость (Y = c + dX). В этом случае говорят о полной корреляции. При р = 1 значения x_i, y_i определяют точки, лежащие на прямой линии, имеющей положительный наклон (с увеличением x_i значения y_i также увеличиваются), при р = -1 прямая имеет отрицательный наклон (рисунок 2.1. Виды зависимости б).

В промежуточных случаях (-1 < p < 1) точки, соответствующие значениям xi_, y_i, попадают в область, ограниченную некоторым эллипсом (рисунок 2.1. Виды зависимости в, рисунок 2.1. Виды зависимости г), причем при p > 0 имеет место положительная корреляция (с увеличением x_i значения y_i имеют тенденцию к возрастанию), при p < 0 корреляция отрицательная. Чем ближе р к ±1, тем уже эллипс и тем теснее экспериментальные значения группируются около прямой линии.

Здесь же следует обратить внимание на то, что линия, вдоль которой группируются точки, может быть не только прямой, а иметь любую другую форму: парабола, гипербола и т. д. В этих случаях мы рассматривали бы так называемую, нелинейную (или криволинейную) корреляцию (рисунок 2.1. Виды зависимости д). Таким образом, визуальный анализ корреляционного поля помогает выявить не только наличия статистической зависимости (линейную или нелинейную) между исследуемыми признаками, но и ее тесноту и форму. Это имеет существенное значение для следующего шага в анализе ѕ выбора и вычисления соответствующего коэффициента корреляции.

Модель корреляционного поля, полученного для данной двумерной выборки представлено на рисунке 2.1.Корреляционное поле.

2.1 Корреляционное поле

Т.к. коэффициент корреляции для данной двумерной выборки равен +0,15, т.е. корреляционная зависимость очень слабая. Чему соответствует вид поля - широкий эллипс, почти круг, имеющий слабую тенденцию возрастания.

2.2 Построение корреляционной таблицы

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

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

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

Корреляционная таблица

y/x |0,019-0,49|0,49-0,97|0,97- 1,4| 1,4- 1,9| 1,9- 2,4| 2,4- 2,9| 2,9- 3,3| 3,3- 3,8|Ny

-0,076-0,33| 0| 0| 1| 1| 5| 3| 2| 0| 12| 0,33-0,74 | 0| 1| 1| 6| 7| 5| 2| 2| 24| 0,74- 1,1 | 2| 0| 1| 3| 10| 7| 1| 4| 28| 1,1- 1,5 | 0| 1| 2| 3| 6| 8| 2| 1| 23| 1,5- 2 | 0| 0| 0| 0| 1| 3| 1| 0| 5| 2- 2,4 | 0| 0| 0| 0| 3| 1| 0| 0| 4| 2,4- 2,8 | 0| 0| 0| 0| 0| 0| 1| 1| 2| 2,8- 3,2 | 0| 0| 0| 0| 0| 0| 0| 0| 0| Nx | 2| 2| 5| 13| 32| 27| 9| 8|

2.3 Расчёт коэффициентов уравнения прямой регрессии

Выбрав вид функции регрессии, т.е. вид рассматриваемой модели зависимости Y от Х (или Х от У), например, линейную модель y_x=a+bx, необходимо определить конкретные значения коэффициентов модели.

При различных значениях а и b можно построить бесконечное число зависимостей вида y_x=a+bx т.е на координатной плоскости имеется бесконечное количество прямых, нам же необходима такая зависимость, которая соответствует наблюдаемым значениям наилучшим образом. Таким образом, задача сводится к подбору наилучших коэффициентов.

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

Обозначим: Y_i - значение, вычисленное по уравнению Y_i=a+bx_i. y_i - измеренное значение, е_i=y_i-Y_i - разность между измеренными и вычисленными по уравнению значениям, е_i=y_i-a-bx_i.

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

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

Если разделить обе части нормальных уравнений на n, то получим:

Учитывая, что

Получим , отсюда, подставляя значение a в первое уравнение, получим:

При этом b называют коэффициентом регрессии; a называют свободным членом уравнения регрессии.

2.4 Нахождение выборочного коэффициента корреляции

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

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

В качестве измерителей степени тесноты парных связей между количественными переменными используются коэффициент корреляции (или то же самое "коэффициент корреляции Пирсона") и корреляционное отношение.

Пусть при проведении некоторого опыта наблюдаются две случайные величины X и Y, причем одно и то же значение x встречается n_x раз, y - n_y раз, одна и та же пара чисел (x,y) наблюдается n_xy раз. Все данные записываются в виде таблицы, которую называют корреляционной.

Выборочная ковариация k(X,Y) величин X и Y определяется формулой

где а x^*, y^*- выборочные средние величин X и Y

Выборочный коэффициент корреляции находится по формуле:

где - выборочные средние квадратические отклонения величин X и Y.

Выборочный коэффициент корреляции r(X,Y) показывает тесноту линейной связи между X и Y: чем ближе r(X,Y) к единице, тем сильнее линейная связь между X и Y.

2.5 Расчёт коэффициентов уравнения криволинейной регрессии

Это уравнение вида Y = b₀ +b₁X₁ + b₂X₂;

1) Найтин еизвестные b₀, b₁,b₂ можно, решим систему трехлинейных уравнений с тремя неизвестными b₀,b₁,b₂:

Для решения системы можете воспользоваться решение системы методом Крамера

2) Или использовав формулы^

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

Здесь z'_jj - j-тый диагональный элемент матрицы Z-1 =(XTX)-1.

При этом:



где m - количество объясняющих переменных модели.

2.6 Нахождение корреляционного отношения

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

Эмпирическое корреляционное отношение определяется по формуле:

, где

- межгрупповая дисперсия;

- общая дисперсия.

2.7 Расчёт критерия Фишера

Критерием Фишера (F-критерием, ц*-критерием) -- называют любой статистический критерий, тестовая статистика которого при выполнении нулевой гипотезы имеет распределение Фишера (F-распределение).

Статистика теста так или иначе сводится к отношению выборочных дисперсий (сумм квадратов, деленных на "степени свободы"). Чтобы статистика имела распределение Фишера необходимо, чтобы числитель и знаменатель были независимыми случайными величинами и соответствующие суммы квадратов имели распределение Хи-квадрат. Для этого требуется, чтобы данные имели нормальное распределение. Кроме того, предполагается, что дисперсия случайных величин, квадраты которых суммируются, одинакова.

Тест проводится путем сравнения значения статистики с критическим значением соответствующего распределения Фишера при заданном уровне значимости. Известно, что если , то ,. Кроме того, квантили распределения Фишера обладают свойством F_1-б=1/F_б. Поэтому обычно на практике в числителе участвует потенциально большая величина, в знаменателе - меньшая и сравнение осуществляется с "правой" квантилью распределения. Тем не менее тест может быть и двусторонним и односторонним. В первом случае при уровне значимости используется квантиль F_б/2, а при одностороннем тесте F_б.

Более удобный способ проверки гипотез - с помощью p-значения p(F) - вероятностью того, что случайная величина с данным распределением Фишера превысит данное значение статистики. Если p(F) (для двустороннего теста - 2p(F) меньше уровня значимости б, то нулевая гипотеза отвергается, в противном случае принимается.

Статистика теста для двумерной выборки будет вычисляться по формуле:

Если статистика больше критического, то дисперсии не одинаковы, в противном случае дисперсии выборок одинаковы.

Все значения, рассчитанные в данном разделе:

Линейное уравнение прямой регрессии: y=1,94-0,00267*x

Коэффициент корреляции: 0,151

Корреляционное отношение: 2,5

Критерий Фишера: 1,41

Также вся информация о двумерной выборке приведена на рисунке 2.7. Обработка двумерной выборки.

Заключение

В данной курсовой работе была смоделирована и решена задача по обработке одномерной и двумерной случайной выборки при помощи вычислительной техники в приложении Microsoft Visual Studio на языке программирования C#. Результаты машинних экспериментов были наглядно проиллюстрированы графиками различных распределений для одномерной выборки и корреляционным полем для двумерной.

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

Литература

случайная выборка корреляционный

1. Вентцель Е.С. Теория вероятностей. - М.: Наука, 1969. - 576с.

2. Гмурман В.Е. Теория вероятностей и математическая статистика. - М.: Высш.школа, 1992. - 368с.

3. Гмурман В.Е. Руководство к решению задач по теории вероятностей и математической статистики. - М.: Высш.школа, 1979. - 400с.

4. Гурский Е.И. Теория вероятностей с элементами математической статистики. - М.: Высш.школа, 1991. - 328с.

Приложение А

Программный код

Form1.

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Linq;

using System.Text;

using System.Windows.Forms;

namespace kurs

{

public partial class Form1: Form

{

public Form1()

{

InitializeComponent();

Text = "= Теория Вероятности =";

MenuStrip Menu = new MenuStrip();

ToolStripMenuItem p1 = new ToolStripMenuItem("Одномерная");

Menu.Items.Add(p1);

ToolStripMenuItem x = new ToolStripMenuItem("Выборка x");

p1.DropDownItems.Add(x);

ToolStripMenuItem y = new ToolStripMenuItem("Выборка y");

p1.DropDownItems.Add(y);

ToolStripMenuItem p2 = new ToolStripMenuItem("Двумерная");

Menu.Items.Add(p2);

ToolStripMenuItem p3 = new ToolStripMenuItem("Выход:)");

Menu.Items.Add(p3);

x.Click += X;

y.Click += Y;

p2.Click += D;

p3.Click += Close;

Controls.Add(Menu);

MainMenuStrip = Menu;

}

private void Form1_Load(object sender, EventArgs e)

{

}

void X(object who, EventArgs e)

{

Form2 f = new Form2("X.txt");

f.Show();

}

void Y(object who, EventArgs e)

{

Form2 f = new Form2("Y.txt");

f.Show();

}

void D(object who, EventArgs e)

{

Form6 f = new Form6();

f.Show();

}

void Close(object who, EventArgs e)

{

DialogResult result = MessageBox.Show("Вы хотите выйти?", "Выход:)", MessageBoxButtons.YesNo);

if (result == DialogResult.Yes) Application.Exit();

}

public int s { get; set; }

}

}

Form2.

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Linq;

using System.Text;

using System.Windows.Forms;

using System.IO;

namespace kurs

{

public partial class Form2: Form

{

int kol;

int k;

double h;

double xmax,xmin;

string name;

double[] pogr = new double[5];

double[] Et = new double[5];

double[] Pr = new double[5];

double[] prak;

double[] y;

double[] Em;

double[] zn;

double[] Points;

public void proverka()

{

if (xmin < 0 || xmin > 0) linkLabel3.Enabled = linkLabel5.Enabled = linkLabel6.Enabled = linkLabel7.Enabled = linkLabel8.Enabled=false;

if (xmin > 0 && xmax<1||xmin==0 || xmax==1) linkLabel1.Enabled = linkLabel2.Enabled = linkLabel4.Enabled = linkLabel8.Enabled = false;

}

void file(string s) { name = s; }

public Form2(string s)

{

InitializeComponent();

Text = "Одномерная выборка";

file(s);

}

private void Form2_Load(object sender, EventArgs e)

{

string[] lines = File.ReadAllLines(name);

int i;

double[] X = new double[lines.Length];

kol = X.Length;

label1.Text = "Выборка: " + X.Length;

for (i = 0; i < kol; i++)

{

X[i] = Convert.ToDouble(lines[i]);

}

k = Convert.ToInt32(1 + 3.32 * Math.Log10(Convert.ToDouble(kol)));

prak = new double[k];

Em = new double[k];

y = new double[k];

zn = new double[k];

Points = new double[k+1];

double t = xmin;

double K = k;

unsort(ref X);

sort(ref X);

for (i = 0; i < Points.Length; i++)

{

Points[i] = t;

t += h;

}

for (i = 1, zn[0] = xmin; i < k; i++)

zn[i] = zn[i - 1] + h;

for (i = 0; i < 5; i++)

pogr[i] = Math.Abs(Et[i] - Pr[i]);

K = 0;

for (i = 0; i < k; i++)

{

y[i] = prak[i] / h;

Em[i] = K;

K += prak[i];

}

chart1.Series[0].Points.DataBindXY(zn, y);

chart2.Series[0].Points.DataBindXY(zn, Em);

proverka();

}

private void button1_Click(object sender, EventArgs e)

{

Form5 f = new Form5();

f.Show();

f.Table(Et, Pr, pogr);

}

public delegate double F(double x, double m, double sko);

public double square (double a, double b, double m, double sko, F f)

{

double hh = b - a, S1, S2, x;

S2 = f(a, m, sko) * hh;

do

{

hh = hh / 2;

S1 = S2;

S2 = 0;

x = a;

while (x < b - hh / 2)

{

S2 += f(x, m, sko);

x += hh;

}

S2 *= hh;

}

while (Math.Abs(S2 - S1) > 0.0001);

return S2;

}

public void krit(double m, double sko, F f)

{

double pr = 0, kl = 0;

double [] Y = new double[k];

double [] teor = new double [k];

double[] R = new double[k];

double P = 0;

int i;

for (i = 0; i < k; i++)

{

Y[i] = square(Points[i], Points[i + 1], m, sko, f);

pr += Math.Pow((prak[i] - Y[i]), 2) / Y[i];

teor[i] = P;

R[i] = Math.Abs(teor[i] - Em[i]);

P += Y[i];

Y[i] = Y[i] / h;

}

pr *= kol;

Array.Sort(R);

kl = R[R.Length - 1] / Math.Sqrt(Em[Em.Length - 1]);

label6.Text = "Критерий Пирсона: " + pr;

label7.Text = "Критерий Колмагорова: " + kl;

chart1.Series[1].Points.DataBindXY(zn, Y);

chart2.Series[1].Points.DataBindXY(zn, teor);

}

public void clear()

{

for (int i=1; i <= 8; i++)

chart1.Series[i].Points.Clear();

}

public double f1(double x, double m, double sko)

{

return (1 / (sko * (Math.Sqrt(2 * Math.PI)))) * Math.Exp(-((Math.Pow(x - m, 2) / (2 * Math.Pow(sko / 1, 2)))));

}

public double f3(double x, double m, double sko)

{

return Math.Sqrt(1 / Math.Pow(sko, 2)) * Math.Exp(-Math.Sqrt(1 / Math.Pow(sko, 2)) * x);

}

public double f4(double x, double l, double c)

{

return (l / 2) * Math.Exp(-l * Math.Abs(x - c));

}

public double f5(double x, double l, double c)

{

return c * l * Math.Pow(x, l - 1) * Math.Exp(-c * Math.Pow(x, l));

}

public double f6(double x, double m, double sko)

{

return (x / Math.Pow(sko, 2)) * (Math.Exp(-Math.Pow(x, 2) / (2 * Math.Pow(sko, 2))));

}

public double f7(double x, double m, double sko)

{

return 1 / (xmax - xmin);;

}

public double f8(double x, double m, double sko)

{

return xmax / (1 - xmin);

}

private void linkLabel1_LinkClicked(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(Et[0], Et[2], f1);

}

private void linkLabel2_LinkClicked(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(0, 1, f1);

}

private void linkLabel3_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(Et[0], Et[2], f3);

}

private void linkLabel4_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(1, 1, f4);

}

private void linkLabel5_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(1, 1, f5);

}

private void linkLabel6_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(Et[0], Et[2], f6);

}

private void linkLabel7_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(Et[0], Et[2], f7);

}

private void linkLabel8_LinkClicked_1(object sender, LinkLabelLinkClickedEventArgs e)

{

krit(Et[0], Et[2], f8);

}

void unsort(ref double [] X)

{

int i;

double M = 0, D = 0, SKO = 0, A = 0, E = 0;

for (i = 0; i < kol; i++)

M += X[i];

M = M / (double)kol;

for (i = 0; i < kol; i++)

{

D += Math.Pow((X[i] - M), 2);

A += Math.Pow((X[i] - M), 3);

E += Math.Pow((X[i] - M), 4);

}

D = D / (kol - 1);

SKO = Math.Sqrt(D);

A = A / kol / Math.Pow(SKO, 3);

E = E / kol / Math.Pow(SKO, 4) - 3;

Et[0] = M; Et[1] = D; Et[2] = SKO; Et[3] = A; Et[4] = E;

}

void sort(ref double [] X)

{

int i,j;

double M_=0,D_=0,SKO=0,t,H,A_=0,E_=0;

for(i=0;i<kol-1;i++)

for(j=0;j<kol-1;j++)

if(X[j+1]>X[j])

{

t=X[j];

X[j]=X[j+1];

X[j+1]=t;

}

double [] KOL=new double [k];

double [] M=new double [k];

double [] D=new double [k];

double [] A=new double [k];

double [] E=new double [k];

for(i=0,xmin=X[0],xmax=X[0];i<kol;i++)

{

if(xmin>X[i])

xmin=X[i];

if(xmax<X[i])

xmax=X[i];

}

h=(xmax-xmin)/(double)k;

label2.Text = "Величина: " + h;

label3.Text = "Диапазон: xmin = " + xmin + " xmax = " + xmax;

label4.Text = "Количество интервалов: " + k;

for(i=0,H=xmin;i<k;i++,H+=h)

for(j=0;j<kol;j++)

if(X[j]>=H&&X[j]<=H+h)

KOL[i]++;

label5.Text = "Интервалы: ";

for (i = 0; i < k; i++)

label5.Text += (i + 1) + ") - " + KOL[i] + "; ";

for(i=0,H=xmin;i<k;i++,H+=h)

M_+=M[i]=(H+H+h)/2*KOL[i]/(double)kol;

for (i = 0, H = xmin; i < k; i++, H += h)

{

D_+=D[i]=Math.Pow(((H+H+h)/2-M_),2)*KOL[i]/(double)kol;

A_+=A[i]=Math.Pow(((H+H+h)/2-M_),3)*KOL[i]/(double)kol;

E_+=E[i]=Math.Pow(((H+H+h)/2-M_),4)*KOL[i]/(double)kol;

SKO=Math.Sqrt(D[i]);

prak[i]=KOL[i]/(double)kol;

}

Pr[0]=M_;

Pr[1]=D_;

Pr[2]=SKO=Math.Sqrt(D_);

Pr[3]=A_=A_/Math.Pow(SKO,3);

Pr[4]=E_=E_/Math.Pow(SKO,4)-3;

}

}

}

Form4.

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Linq;

using System.Text;

using System.Windows.Forms;

using System.IO;

using System.Windows.Forms.DataVisualization.Charting;

namespace kurs

{

public partial class Form4: Form

{

public Form4()

{

InitializeComponent();

Text = "Корреляционное поле";

}

private void Form4_Load(object sender, EventArgs e)

{

string [] vib = File.ReadAllLines("X.txt");

string[] vib2 = File.ReadAllLines("Y.txt");

int kol = vib.Length,i;

double[] X = new double[kol];

double[] Y = new double[kol];

for (i = 0; i < kol; i++)

{

X[i] = Convert.ToDouble(vib[i]);

Y[i] = Convert.ToDouble(vib2[i]);

}

for (i = 0; i < kol; i++)

{

chart1.Series[0].Points.AddXY(X[i],Y[i]);

}

}

private void chart1_Click(object sender, EventArgs e)

{

}

}

}

Form5.

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Linq;

using System.Text;

using System.Windows.Forms;

namespace kurs

{

public partial class Form5: Form

{

public Form5()

{

InitializeComponent();

Text = "Таблица точечных оценок";

}

private void Form5_Load(object sender, EventArgs e)

{

}

public void Table(double[] Et, double[] Pr, double[] pogr)

{

dataGridView1.AllowUserToAddRows = false;

DataGridViewRow row = new DataGridViewRow();

DataGridViewCell Cell0 = new DataGridViewTextBoxCell(); Cell0.Value = "Негруппированная";

DataGridViewCell Cell1 = new DataGridViewTextBoxCell(); Cell1.Value = Convert.ToString(Et[0]);

DataGridViewCell Cell2 = new DataGridViewTextBoxCell(); Cell2.Value = Convert.ToString(Et[1]);

DataGridViewCell Cell3 = new DataGridViewTextBoxCell(); Cell3.Value = Convert.ToString(Et[2]);

DataGridViewCell Cell4 = new DataGridViewTextBoxCell(); Cell4.Value = Convert.ToString(Et[3]);

DataGridViewCell Cell5 = new DataGridViewTextBoxCell(); Cell5.Value = Convert.ToString(Et[4]);

row.Cells.AddRange(Cell0, Cell1, Cell2, Cell3, Cell4, Cell5);

dataGridView1.Rows.Add(row);

DataGridViewRow row1 = new DataGridViewRow();

DataGridViewCell Cell0a = new DataGridViewTextBoxCell(); Cell0a.Value = "Cгруппированная";

DataGridViewCell Cell1a = new DataGridViewTextBoxCell(); Cell1a.Value = Convert.ToString(Pr[0]);

DataGridViewCell Cell2a = new DataGridViewTextBoxCell(); Cell2a.Value = Convert.ToString(Pr[1]);

DataGridViewCell Cell3a = new DataGridViewTextBoxCell(); Cell3a.Value = Convert.ToString(Pr[2]);

DataGridViewCell Cell4a = new DataGridViewTextBoxCell(); Cell4a.Value = Convert.ToString(Pr[3]);

DataGridViewCell Cell5a = new DataGridViewTextBoxCell(); Cell5a.Value = Convert.ToString(Pr[4]);

row1.Cells.AddRange(Cell0a, Cell1a, Cell2a, Cell3a, Cell4a, Cell5a);

dataGridView1.Rows.Add(row1);

DataGridViewRow row2 = new DataGridViewRow();

DataGridViewCell Cell0b = new DataGridViewTextBoxCell(); Cell0b.Value = "Погрешность";

DataGridViewCell Cell1b = new DataGridViewTextBoxCell(); Cell1b.Value = Convert.ToString(pogr[0]);

DataGridViewCell Cell2b = new DataGridViewTextBoxCell(); Cell2b.Value = Convert.ToString(pogr[1]);

DataGridViewCell Cell3b = new DataGridViewTextBoxCell(); Cell3b.Value = Convert.ToString(pogr[2]);

DataGridViewCell Cell4b = new DataGridViewTextBoxCell(); Cell4b.Value = Convert.ToString(pogr[3]);

DataGridViewCell Cell5b = new DataGridViewTextBoxCell(); Cell5b.Value = Convert.ToString(pogr[4]);

row2.Cells.AddRange(Cell0b, Cell1b, Cell2b, Cell3b, Cell4b, Cell5b);

dataGridView1.Rows.Add(row2);

}

}

}

Form6.

using System;

using System.Collections.Generic;

using System.ComponentModel;

using System.Data;

using System.Drawing;

using System.Linq;

using System.Text;

using System.Windows.Forms;

using System.IO;

namespace kurs

{

public partial class Form6: Form

{

int k1;

int k2;

double h1;

double h2;

public void table(double[] X, double[] Y, int kol1, int kol2, out double xmin2, out double xmax2)

{

int i, j, l, sum=0;

double xmin1, xmax1, h, hh;

for (i = 0, xmin1 = X[0], xmax1 = X[0]; i < kol1; i++)

{

if (xmin1 > X[i])

xmin1 = X[i];

if (xmax1 < X[i])

xmax1 = X[i];

}

for (i = 0, xmin2 = Y[0], xmax2 = Y[0]; i < kol2; i++)

{

if (xmin2 > Y[i])

xmin2 = Y[i];

if (xmax2 < Y[i])

xmax2 = Y[i];

}

k1 = Convert.ToInt32(1 + 3.32 * Math.Log10(Convert.ToDouble(kol1)));

k2 = Convert.ToInt32(1 + 3.32 * Math.Log10(Convert.ToDouble(kol2)));

h1 = (xmax1 - xmin1) / k1;

h2 = (xmax2 - xmin2) / k2;

int[,] Mas = new int[k1, k2];

int[] Sum = new int[k1];

for (i = 0; i < k1; i++)

{

for (j = 0; j < k2; j++)

Mas[i, j] = 0;

Sum[i] = 0;

}

for (l = 0; l < kol1; l++)

for (i = 0, h = xmin1; i < k1; i++, h += h1)

if (X[l] <= (h + h1) && X[l] >= h-0.0001)

for (j = 0, hh = xmin2; j < k2; j++, hh += h2)

if (Y[l] <= (hh + h2) && Y[l] >= hh)

Mas[i, j]++;

h=xmin1;

DataGridViewTextBoxColumn Column0 =new DataGridViewTextBoxColumn();

Column0.HeaderText = "y/x";

DataGridViewTextBoxColumn Column1 =new DataGridViewTextBoxColumn();

Column1.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column2 =new DataGridViewTextBoxColumn();

Column2.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column3 =new DataGridViewTextBoxColumn();

Column3.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column4 =new DataGridViewTextBoxColumn();

Column4.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column5 =new DataGridViewTextBoxColumn();

Column5.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column6 =new DataGridViewTextBoxColumn();

Column6.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column7 =new DataGridViewTextBoxColumn();

Column7.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column8 = new DataGridViewTextBoxColumn();

Column8.HeaderText = Convert.ToString(h) + "-" + Convert.ToString(h + h1); h += h1;

DataGridViewTextBoxColumn Column9 =new DataGridViewTextBoxColumn();

Column9.HeaderText = "Ny";

dataGridView1.Columns.Add(Column0);

dataGridView1.Columns.Add(Column1);

dataGridView1.Columns.Add(Column2);

dataGridView1.Columns.Add(Column3);

dataGridView1.Columns.Add(Column4);

dataGridView1.Columns.Add(Column5);

dataGridView1.Columns.Add(Column6);

dataGridView1.Columns.Add(Column7);

dataGridView1.Columns.Add(Column8);

dataGridView1.Columns.Add(Column9);

for (i = 0,h=xmin2; i <k2; i++, sum = 0,h+=h2)

{

DataGridViewCell Cell0 = new DataGridViewTextBoxCell();

DataGridViewCell Cell1 = new DataGridViewTextBoxCell();

DataGridViewCell Cell2 = new DataGridViewTextBoxCell();

DataGridViewCell Cell3 = new DataGridViewTextBoxCell();

DataGridViewCell Cell4 = new DataGridViewTextBoxCell();

DataGridViewCell Cell5 = new DataGridViewTextBoxCell();

DataGridViewCell Cell6 = new DataGridViewTextBoxCell();

DataGridViewCell Cell7 = new DataGridViewTextBoxCell();

DataGridViewCell Cell8 = new DataGridViewTextBoxCell();

DataGridViewCell Cell9 = new DataGridViewTextBoxCell();

DataGridViewRow row = new DataGridViewRow();

Cell0.Value = Convert.ToString(h) + "-" + Convert.ToString(h+h2);

Cell1.Value = Convert.ToString(Mas[i, 0]); sum += Mas[i, 0]; Sum[0] += Mas[i, 0];

Cell2.Value = Convert.ToString(Mas[i, 1]); sum += Mas[i, 1]; Sum[1] += Mas[i, 1];

Cell3.Value = Convert.ToString(Mas[i, 2]); sum += Mas[i, 2]; Sum[2] += Mas[i, 2];

Cell4.Value = Convert.ToString(Mas[i, 3]); sum += Mas[i, 3]; Sum[3] += Mas[i, 3];

Cell5.Value = Convert.ToString(Mas[i, 4]); sum += Mas[i, 4]; Sum[4] += Mas[i, 4];

Cell6.Value = Convert.ToString(Mas[i, 5]); sum += Mas[i, 5]; Sum[5] += Mas[i, 5];

Cell7.Value = Convert.ToString(Mas[i, 6]); sum += Mas[i, 6]; Sum[6] += Mas[i, 6];

Cell8.Value = Convert.ToString(Mas[i, 7]); sum += Mas[i, 7]; Sum[7] += Mas[i, 7];

Cell9.Value = Convert.ToString(sum);

row.Cells.AddRange(Cell0, Cell1, Cell2, Cell3, Cell4, Cell5, Cell6, Cell7, Cell8,Cell9);

this.dataGridView1.Rows.Add(row);

}

DataGridViewCell Cell00 = new DataGridViewTextBoxCell();

DataGridViewCell Cell01 = new DataGridViewTextBoxCell();

DataGridViewCell Cell02 = new DataGridViewTextBoxCell();

DataGridViewCell Cell03 = new DataGridViewTextBoxCell();

DataGridViewCell Cell04 = new DataGridViewTextBoxCell();

DataGridViewCell Cell05 = new DataGridViewTextBoxCell();

DataGridViewCell Cell06 = new DataGridViewTextBoxCell();

DataGridViewCell Cell07 = new DataGridViewTextBoxCell();

DataGridViewCell Cell08 = new DataGridViewTextBoxCell();

DataGridViewCell Cell09 = new DataGridViewTextBoxCell();

DataGridViewRow row_ = new DataGridViewRow();

Cell00.Value = "Ny";

Cell01.Value = Convert.ToString(Sum[0]);

Cell02.Value = Convert.ToString(Sum[1]);

Cell03.Value = Convert.ToString(Sum[2]);

Cell04.Value = Convert.ToString(Sum[3]);

Cell05.Value = Convert.ToString(Sum[4]);

Cell06.Value = Convert.ToString(Sum[5]);

Cell07.Value = Convert.ToString(Sum[6]);

Cell08.Value = Convert.ToString(Sum[7]);

Cell09.Value = " ";

row_.Cells.AddRange(Cell00, Cell01, Cell02, Cell03, Cell04, Cell05, Cell06, Cell07, Cell08, Cell09);

this.dataGridView1.Rows.Add(row_);

}

public void ur_pr(double[] X, double[] Y, int kol1,out double a,out double b)

{

int i;

double Xi = 0, Yi = 0, XiYi = 0, Xi2 = 0;

for (i = 0; i < kol1; i++)

{

Xi += X[i];

Xi2 += X[i] * X[i];

Yi = Y[i];

XiYi += X[i] * Y[i];

}

b = (XiYi - Xi * Yi) / (Xi2 - Xi * Xi);

a = (Xi2 * Yi - Xi * XiYi) / (Xi2 - Xi * Xi);

}

public void koef_kor(double[] X, double[] Y, int kol1, int kol2, out double Dx, out double Dy, out double Y_, out double r)

{

int i;

Y_ = 0; Dx = 0; Dy = 0;

double X_ = 0, k = 0, SKOx, SKOy;

for (i = 0; i < kol1; i++)

{

X_ += X[i];

Y_ = Y_ + Y[i];

}

X_ /= kol1;

Y_ = Y_ / kol2;

for (i = 0; i < kol1; i++)

{

k += (X[i] - X_) * (Y[i] - Y_);

Dx = Dx + Math.Pow((X[i] - X_), 2);

Dy = Dy + Math.Pow((Y[i] - Y_), 2);

}

k /= kol1;

Dx = Dx / (kol1 - 1);

Dy = Dy / (kol1 - 1);

SKOx = Math.Sqrt(Dx);

SKOy = Math.Sqrt(Dy);

r = k / (SKOx * SKOy);

}

public void ur_kt(double[] X, double[] Y)

{

int i;

double o,a,b,c;

double sx=0,sy=0,sxy=0,ss=0,s2=0,s3=0,s4=0;

for(i=0;i<k1;i++)

{

sx+=X[i];

sy+=Y[i];

sxy+=X[i]*Y[i];

ss+=Math.Pow(X[i],2)*Y[i];

s2+=Math.Pow(X[i],2);

s3+=Math.Pow(X[i],3);

s4+=Math.Pow(X[i],4);

}

o=k1*s2*s4+sx*s3*s2+s2*sx*s3-s2*s2*s2-sx*sx*s4-k1*s3*s3;

a=sy*s2*s4+sx*s3*ss+s2*sxy*s3-s2*s2*ss-sxy*sx*s4-s3*s3*sy;

b=k1*sxy*s4+sy*s3*s2+sx*ss*s2-s2*sxy*s2-sx*sy*s4-s3*ss*k1;

c=k1*s2*ss+sx*sxy*s2+sx*s3*sy-sy*s2*s2-sx*sx*ss-s3*sxy*k1;

a=a/o;

b=b/o;

c=c/o;

label3.Text="y="+a+b+"x+"+c+"x^2";

}

public void kor_otn(double[] X, double[] Y, int kol2, double xmin, double xmax, double Dy, double Y_,out double N)

{

int i, j;

double h, Dy_gr = 0;

double[] kol_gr = new double[k2];

double[] sr_gr = new double[k2];

for (i = 0; i < k2; i++)

{

kol_gr[i] = 0;

sr_gr[i] = 0;

}

for (i = 0; i < kol2; i++)

for (j = 0, h = xmin; j < k2; j++, h += h2)

if (Y[i] <= (h + h2) && Y[i] >= h)

{

kol_gr[j]++;

sr_gr[j] += Y[i];

}

for (i = 0; i < k2; i++)

sr_gr[i] /= kol_gr[i];

for (i = 0; i < k2; i++)

Dy_gr += Math.Pow((sr_gr[i] - Y_), 2) * kol_gr[i] / kol2;

N = Math.Sqrt(Dy_gr / Dy);

}

public void fisher(double Dx, double Dy,out double F)

{

if (Dx > Dy)

F = Dx / Dy;

else

F = Dy / Dx;

}

public Form6()

{

InitializeComponent();

}

private void Form6_Load(object sender, EventArgs e)

{

int i;

double Dx, Dy, xmin, xmax, Y_, a, b, r, N, F;

string[] lines1 = File.ReadAllLines("X.txt");

double[] X = new double[lines1.Length];

for (i = 0; i < X.Length; i++)

{

X[i] = Convert.ToDouble(lines1[i]);

}

string[] lines2 = File.ReadAllLines("Y.txt");

double[] Y = new double[lines2.Length];

for (i = 0; i < Y.Length; i++)

{

Y[i] = Convert.ToDouble(lines2[i]);

}

int kol1 = X.Length;

int kol2 = Y.Length;

table(X,Y,kol1,kol2,out xmin,out xmax);

ur_pr(X, Y, kol1, out a, out b);

label1.Text = "Линейное уранение прямой регрессии: y=" + a + b + "*x";

koef_kor(X, Y, kol1, kol2,out Dx,out Dy,out Y_,out r);

label2.Text = "Коэффициент корреляции: " + r;

ur_kt(X,Y);

kor_otn(X, Y, kol2, xmin, xmax, Dy, Dx,out N);

label4.Text = "Корреляционное отношение: " + N;

fisher(Dx, Dy,out F);

label5.Text = "Критерий Фишера: " + F;

}

private void button2_Click(object sender, EventArgs e)

{

Form4 f = new Form4();

f.Show();

}

}

}

Размещено на Allbest.ru