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

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

МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ

Учреждение образования

«Гомельский государственный технический университет имени П.О. Сухого»

Факультет автоматизированных и информационных систем

Кафедра «Информационные технологии»

направление специальности 1-40 01 02-01 «Информационные системы и технологии в проектировании и производстве»

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

к курсовой работе

по дисциплине «Основы алгоритмизации и программирования»

на тему: «Определение параметров регрессионной зависимости полной фактической себестоимости топлива и электроэнергии, расходуемой на технологические цели»

Исполнитель: студент гр. ИТ-12

Дедков Р.А.

Гомель 2015

Содержание

Введение

. Теоретические сведения

.1 Функции ввода-вывода строк и символов языка С+

.2 Функции ввода информации

. Постановка задачи и исходные данные

3. Графическая схема алгоритма

3.1 Вычисление среднего значения, дисперсии, среднеквадратических отклонений х и у

3.2 Вычисление коэффицента парной кораляции

3.3 Вычисление регресии для первой и второй функции

3.4 Вычисление остаточных дисперсий

. Расчет параметров регрессионных зависимостей с помощью надстройки Пакет анализа MSExcel

4.1 Расчет параметров регрессионных зависимостей

.2 График с исходными данными и регрессионными зависимостями

. Описание функционирования программного комплекса

Заключение

Список использованных источников

Приложение А

Приложение Б

Введение

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

Цели регрессионного анализа:

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

предсказание значения зависимой переменной с помощью незаисимой;

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

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

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

В ходе регрессионного анализа решаются две основные задачи:

построение уравнения регрессии, т.е. нахождение вида зависимости между результатным показателем и независимыми факторами x1, x2, ..., xn;

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

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

факторная дисперсия, характеризующаяся вариацию результативного признака, объясняемую только признаком фактором;

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

общая дисперсия, складывающаяся за счет влияния всех факторов;

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

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

Регрессия может быть представлена в виде суммы не случайной и случайной составляющих.

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

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

Применяется регрессионный анализ главным образом для планирования, а также для разработки нормативной базы.

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

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

В процессе выполнения работы решаются следующие задачи:

разработка алгоритма вычисления статистических характеристик по заданным формулам;

запись алгоритма в виде блок-схемы и программы на языке Си;

подготовка данных для отладки программы средствами табличного процессора MSExcel;

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

1. 
Теоретические сведения

1.1    Функции ввода-вывода строк и символов языка С

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

Декларации функций ввода/вывода, как уже упоминалось, приведены в заголовочном файле stdio.h.

Для вывода информации на экран монитора (дисплей) в языке Си чаще всего используются функции: printf() и puts().

Формат функции форматного вывода на экран:( управляющая строка , список объектов вывода);

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

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

Количество и порядок следования форматов должен совпадать с количеством и порядком следования выводимых на экран объектов.

Функция printf выполняет вывод данных в соответствии с указанными форматами, поэтому формат может использоваться и для преобразования типов выводимых объектов.

Для чисел long добавляется символ l, например, %ld - длинное целое, %lf - число вещественное с удвоенной точностью - double.

Если нужно напечатать сам символ %, то его нужно указать 2 раза:("Только %d%% предприятий не работало. \n",5);

Получим: Только 5% предприятий не работало.

Управляют выводом специальные последовательности символов: \n - новая строка; \t - горизонтальная табуляция; \b - шаг назад; \r - возврат каретки; \v - вертикальная табуляция; \\ - обратная косая; \' - апостроф; \" - кавычки; \0 - нулевой символ (пусто).

Пример:

#define PI 3.14159

. . .number = 5;bat = 255;cost = 11000;

. . .(" %d студентов съели %f бутербродов. \n", number, but);(" Значение числа pi равно %f. \n", pi);

printf(" Стоимость этой вещи %d %s. \n", cost, "Руб.");

. . .

В модификаторах формата функции printf после символа % можно указывать число, задающее минимальную ширину поля вывода, например, %5d - для целых, %4.2f - для вещественных - две цифры после запятой для поля шириной 4 символа. Если указанных позиций для вывода целой части числа не хватает, то происходит автоматическое расширение.

Если после «%» указан знак «минус», то выводимое значение будет печататься с левой позиции поля вывода, заданной ширины, например: % - 10d.

Использование функции printf для преобразования данных:

1) printf("%d", 336.65);        получим:   336;

) printf("%o", 336);       получим:   520, т.е. 5*8**2+2*8+0*1 = 336;

) printf("%x", 336);      получим:   150 (шестнадцатеричное).

Можно использовать функцию printf для нахождения кода ASCII некоторого символа:

printf (" %c - %d\n", 'a', 'a');

получим десятичный код ASCII символа а:    a - 65 .       

Функция puts(ID строки); выводит на экран дисплея строку символов, автоматически добавляя к ней символ перехода на начало новой строки (\n).

Аналогом такой функции будет: printf(“%s \n”, ID строки);

Функция putchar() выдает на экран дисплея один символ без добавления символа ‘\n’.

1.2 Функции ввода информации

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

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

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

Если нужно ввести значение строковой переменной, то использовать символ & не нужно, т.к. строка - это массив символов, а ID массива является адресом его первого элемента. Например:

course;grant;name[20];(" Укажите курс, стипендию, имя \n ");("%d %lf %s", &course, &grant, name);

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

Функция scanf() использует практически тот же набор модификаторов форматов, что и printf(); отличия от функции вывода следующие: отсутствует формат %g, форматы %e,%f - эквивалентны. Для ввода коротких целых чисел введен модификатор формата %h.

Внимание. Функцией scanf() по формату %s строка вводится только до первого пробела.

Для ввода фраз, состоящих из слов, разделенных пробелами, используется функция gets (ID строковой переменной);

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

Если же использовать ее в правой части операции присваивания, например:

char c;

. . .

c = getch();

то символьная переменная с получит значение кода нажатой клавиши.

С началом работы любой программы автоматически открываются стандартные потоки для ввода (stdin) и вывода данных (stdout), которые по умолчанию связаны с клавиатурой и экраном монитора соответственно.

Ввод данных функциями gets(), getch() выполняется с использованием потока stdin. Если указанная функция не выполняет своих действий (проскакивает), перед использованием необходимо очистить поток (буфер) ввода с помощью функции fflush(stdin).

2. Постановка задачи и исходные данные

Исходными данными задачи являются:

1)      набор значений фактора  и набор значений результата , указанных в таблице 1.1 , где n - количество значений;

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

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

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

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

. средние значения этих наборов

, ;

. дисперсии

, ;

3. среднеквадратические отклонения:

, ;

. коэффициент парной корреляции:

;

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

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

6. вычислить значения результата по регрессионным зависимостям:

, ;

. вычислить остаточные дисперсии:

, ;

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

, ;

9.      в заключение алгоритм должен сравнить полученные значения  и  - лучшей регрессионной зависимостью является та, для которой коэффициент Фишера имеет большее значение;

В таблице 1.1 представлены факторы и результаты

Таблица 1.1 - Таблица факторов и результатов

Формула для расчета первой регрессионной зависимости:

Формула для расчета второй регрессионной зависимости:

Формула для расчета по первой регрессии:

Формула для расчетапо первой регрессии:

Где , n - количество факторов       

Формула для расчета по второй регрессии:

Формула для расчетапо второй регрессии:

Где , n - количество факторов

3. Графическая схема алгоритма решения задачи

.1 Вычисление среднего значения, дисперсии, среднеквадратических отклонений х и у (блок AB)

Для вычисления среднего значения, дисперсии, среднеквадратических отклонений фактических массивов x и y нужно два раза обратиться к алгоритму с разными параметрами. Схема алгоритма блока АВ представлена на рисунке 3.2 :

3.2 Вычисление коэффицента парной кораляции (блок СD):

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

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

Рисунок 3.1 - Обобщенная схема алгоритма решения задачи

Рисунок 3.2 - Вычисление среднего значения, дисперсии, среднеквадратических отклонений х и у (детализация блока AB)

Рисунок 3.3 - Вычисление коэффицента парной кораляции ( детализация блока СD )

3.3 Вычисление регресии для первой и второй функции (блок EF)

Для вычисления результатов по регрессионным зависимостям, нужно вычислить коэффициенты a₀ и a₁ для двух регрессионных зависимостей по следующему алгоритму представленном на рисунке 3.4 :

Рисунок 3.4 - Вычисления коэффициентов a₀ и a₁ для регрессионных зависимостей (детализация блока EF)

В качестве параметра fi в алгоритме на рисунке(3.4) будут передаваться функции Func1 и Func2. Схема функций представлена на рисунке 3.5:

Рисунок 3.5 - Функции Func1 и Func2

3.4 Вычисление остаточных дисперсий (блок GI)

Для вычисления остаточных дисперсий нам нужно два раза обратиться к алгоритму с разными параметрами. Схема алгоритма блока GI будет иметь вид (рисунок 3.6):

Рисунок 3.6 - Вычисление остаточных дисперсий (детализация блока GI)

4. Расчет параметров регрессионных зависимостей с помощью надстройки Пакет анализа MS Excel

.1 Расчет параметров регрессионных зависимостей

В таблице 4.1.1 отображены исходные данные: факторы и результаты. Таблицы 4.1.2 и 4.1.3 содержат значения X и Y для первой и второй регрессий. В таблице 4.1.4 находятся значения a0 и a1 для первой и второй регрессий. Таблицы 4.1.5 и 4.1.6 содержат значения регрессионная статистика и дисперсионного анализа по первой регрессии, 4.1.7 и 4.1.8 - по второй регрессии соответственно.

Таблица 4.1.1 - Факторы и результаты

Таблица 4.1.2 - Значения функции от X для первой регрессии, значения Y для первой регрессии

Таблица 4.1.9 - Значения функции от X для первой регрессии, значения Y для первой регрессии (записанные в формульном виде)

Таблица 4.1.3 - Значения функции от X для второй регрессии, значения Y для второй регрессии

Таблица 4.1.10 - Значения функции от X для второй регрессии, значения Y для второй регрессии(записанные в формульном виде)

Таблица 4.1.4 - Значения a0 и a1 для первой и второй регрессий

Таблица 4.1.5 - Регрессионная статистика по первой регрессии

Таблица 4.1.6 - Дисперсионный анализ по первой регрессии

Таблица 4.1.7 - Регрессионная статистика по второй регрессии

Таблица 4.1.8 - Дисперсионный анализ по второй регрессии

.2 график с исходными данными и регрессионными зависимостями

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

График 4.2.1 - График с исходными данными и регрессионными зависимостями

5. Описание функционирования программного комплекса

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

В случае отсутствия файла с введенным именем будет выдана ошибка (рисунок 5.3). Если такой файл существует, то вам нужно выбрать вы хотите вывести на экран результат выполнения программы или записать результаты в отдельный файл (рисунок 5.4).Далее выводится таблица входных данных из трех столбцов: названия райпо, факторы X, результаты Y(рисунок 5.5). Далее будет выведены средние значение фактора и результата, дисперсия по X иY, среднеквадратичное отклонение Y, коэффициент парной корреляции, коэффициенты a0 и a1 для первой и второй регрессии (рисунок 5.6).

В случае, если коэффициент парной корреляции меньше чем 0.5, будет выведено сообщение о слабой связи. Если же коэффициент парной корреляции больше 0.5, то будет выведена таблица результатов из пяти столбцов: названия райпо, факторы X, результаты Y, значения Y по первой регрессии, значения Y по второй регрессии (рисунок 5.7). После этого будут выведены результаты: остаточные дисперсии по Y, коэффициенты Фишера для первой и второй регрессий, и вывод о том, какая регрессия более полно отражает результаты (рисунок 5.8).

Рисунок 5.1 - Информация о работе

Рисунок 5.2 - Задачи

Рисунок 5.3 - Ошибка открытия файла

Рисунок 5.4 - Меню выбора действия с файлом

Рисунок 5.5 - Таблица факторов и результатов

Рисунок 5.6 - Вывод средних значений факторов и результатов, дисперсии и среднеквадратического отклонения по X и Y, коэффициента парной корреляции, коэффициенты a0 и a1 для первой и второй регрессии

Рисунок 5.7 - Таблица выходных результатов работы программы (значений Y по первой и второй регрессиям)

Рисунок 5.8 - Вывод результатов (остаточные дисперсии, коэффициенты Фишера, вывод о большем соответствии первой или второй регрессии исходным результатам)

дисперсия корреляция регрессия

Заключение

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

Кроме того были высчитаны:

1)      Средние значения заданных наборов фактора и результата;

2)      Дисперсии;

)        Среднеквадратические отклонения;

)        Коэффициент парной корреляции;

)        Характеристики регрессионных зависимостей;

)        Значения результатов по регрессионным зависимостям;

)        Остаточные дисперсии;

)        Коэффициенты Фишера.

Список использованных источников

1. Регрессионный анализ 

. Язык программирования Си

3. Мелещенко, А. Основы программирования на языке /А.А. Мелещенко. - Мн.: БГУИР, 2004. - 232 с.

. Бусько, В Основы алгоритмизации и программирования: Конспект лекций для студ. всех спец. и форм обуч. БГУИР / В.Л.Бусько, А.Г.Корбит, Т.М. Кривоносова. - Мн.: БГУИР, 2004. - 103 с.

5. Программирование ввода-вывода данных и линейных вычислительных алгоритмов на языке Си: практическое пособие к выполнению лабораторных и контрольных работ по дисциплине “Вычислительная техника и программирование” для студентов технических специальностей дневной и заочной форм обучения (м/УК №3089) / О.А. Кравченко, А.М. Мартыненко. - Гомель: ГГТУ им. П.О. Сухого, 2005. - 33 с.

Приложение А

Код программы на языке C файла kursovaya.cpp:

#include"bibliot.h"//подключаем головную библиотеку

/*подключаем русский язык*/bufRus[256];

char *Rus(const char * text)

{(text, bufRus);bufRus;

}()//главная функция

{("color F0");// Цвет экрана и шрифта (F - белый фон, 0 - чёрный шрифт)("\n");(Rus("\t МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ"));(Rus("\t УЧРЕЖДЕНИЕ ОБРАЗОВАНИЯ ГОМЕЛЬСКИЙ ГОСУДАРСТВЕННО ТЕХНИЧЕСКИЙ"));(Rus("\t УНИВЕРСИТЕТ ИМ. П.О.СУХОГО"));

puts("\n\n\n\n\n\n\n");(Rus("\t Факультет автоматизированных и информационных систем"));(Rus("\t Кафедра ''Информационные технологии''"));(Rus("\t направление специальности 1-40 01 02-01 Информационные системы и"));(Rus("\t технологии в проектировании и производстве\n\n\n\n\n"));(Rus("\t КУРСОВАЯ РАБОТА"));(Rus("\t по дисциплине ''Основы алгоритмизации и программирования''"));(Rus("\tна тему:''Определение параметров регрессионной зависимости валовых"));(Rus("\t доходов от торговли от амортизации основных средств торговли''"));

puts("\n\n\n\n\n\n\n\n");(Rus("\t Исполнитель: студент гр. ИТ-12"));

puts(Rus("\t  Ясонов В.А."));(Rus("\t Руководитель: Лозовская Е.В."));("\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n");

puts(Rus("\t Гомель 2012")); ();("cls");(Rus("Задание:"));(Rus("Целью курсовой работы является получение студентами навыков по алгоритмизации и программированию.\nВ процессе выполнения работы студент решает следующие задачи:\n1)разработка алгоритма вычисления статистических характеристик по заданным формулам;\n2)запись алгоритма в виде блок-схемы и программы на языке Си;\n3)подготовка данных для отладки программы средствами табличного процессора MS Excel;\n4)изучение теоретических вопросов использования системы программирования.\n\nНа основе указанных исходных данных требуется написать программу на языке C, которая для заданных наборов фактора и результата рассчитывает следующие статистические характеристики:\n1)средние значения этих наборов\n2)дисперсии\n3)среднеквадратические отклонения\n4)коэффициент парной корреляции\n5)регрессионные зависимости\n6)значения результата по регрессионным зависимостям\n7)остаточные дисперсии\n8)коэффициенты Фишера\n"));

getch();

system("cls");nazv[50], s[100], vib;// Название райпо, вспомогательная строка, переменная для выбора - на экран выводить результаты или в файлx[100], y[100];// x - прочие производственне расходы, y - полная фактическая себестоимостьSx, Sy, Srx, Sry;// Sx, Sy - суммы элементов массивов x и y, Srx, Sry - Средние значения элементов массивов x и yi, n;// i - номер райпо, параметр цикла, n - количество райпоSomx, Somy, r;// Среднеквадратичные отклонения, коэффициент парной корреляцииregr1[100], regr2[100], a0, a1;//значения первой и второй регрессий, коэффициенты a0, a1OstDisp1, OstDisp2, Fish1, Fish2;//Значения остаточных дисперсий, коэффициенты Фишераfile[100];*f;// Описание файловой переменной*fout;//Файловая переменная для записи результатов(Rus("Введите название файла"));(file);=fopen(file,"r");// Открытие файла в режиме "r" - только чтение("cls");(f==NULL)// Проверяем, открылся ли файл, если значение NULL - ошибка открытия.

{("color 4F");// Цвет экрана и шрифта (4 - красный фон, F - белый шрифт)(Rus("\n\n\n\n\n\n\n\n\n\n\t\t\tОшибка открытия файла!"));// Вывод сообщения об ошибке

}// Иначе...

{(Rus("1)Для вывода на экран нажмите клавишу Enter\n2)Для записи данных в файл нажмит любую другую клавишу"));("%c", &vib);("cls");//Очистка экрана

if(vib=='\n')// Нажата клавиша Enter?=fopen("con","a");// con - консоль (вывод на экран)// Иначе...

{=fopen("fail.dat","w");//Открытие файла в режиме "w" - запись данных, все предыдущие данные стирются( Rus("Файл записан,для выхода нажмите любую клавишу...\n"));//Вывод сообщения о том что информация записана в файл

}

/*Выводим шапку таблицы*/(fout,"\t\xDA\xc4\xc4\xc4\xc2\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc2\xc4\xc4\xc4\xc4\xc4\xc4\xc2\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xbf\n");=0;// Счетчик количества строк в файле

/*Считываем данные из файла,одновременно выводя их в таблицу*/(!feof(f))// Пока не конец файла...

{(s, 80, f);// Считываем строку, не более 80 символов. (Считает до ближайшего перехода на следующую строку - символа '\n').(nazv, s, 21);//Из строки s копируем первые 21 символ в строку navz[19]='\0';// Последнему символу присваиваем значение '\0' - конец строки.[i]=atof(&s[20]);//Начиная с 21 символа считываем значение х (все цифры до пробела)[i]=atof(&s[25]);//Начиная с 25 символа считываем значение у (все цифры до конц строки)(fout,"\t\xb3%3d\xb3\%-21s\xb3%6.1f\xb3%7.1f\xb3\n", i+1, Rus(nazv), x[i], y[i]); //Выводим полученные значения.++;// Увеличиваем счетчик элементов на единицу

}=i;// Запоминаем количество элементов массива

/*Выводим низ таблицы*/(fout,"\t\xc0\xc4\xc4\xc4\xc1\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xc1\xc4\xc4\xc4\xc4\xc4\xc4\xc1\xc4\xc4\xc4\xc4\xc4\xc4\xc4\xd9\n");

/*Вычисляем статистические величины - сумму, среднее значение, среднеквадратичное отклонение*/

statistic(n, x, &Sx, &Srx, &Somx);(n, y, &Sy, &Sry, &Somy);

/* Вывод сумм, средних значений и среднеквадратичных отклонений*/(fout,Rus("Результаты вычислений в массиве x\nСреднее значение: %-15.f\nСреднеквадратичное отклонение: %.3f\n\n"), Sx, Srx, Somx);(fout,Rus("Результаты вычислений в массиве y\nСреднее значение: %-15.f\nСреднеквадратичное отклонение: %.3f\n"), Sy, Sry, Somy);=top(n, x, Srx, y, Sry)/((n-1)*Somx*Somy);//Вычисление коэффициента парной корреляции(fout,Rus("\nКоэффициент парной корреляции %.3f\n\n"), r);//Вывод коэффициента парной корреляции(fabs(r)<=0.5)//Если коэффициент парной коррелдяции по модулю меньше 0.5

{("color FC");// Цвет экрана и шрифта (F - белый фон, C - красный шрифт)(fout, Rus("Связь между x и y слабая!"));//Вывод сообщения о слабой связи между x и у.

}//Иначе

{(n, Sy, fi1, x, &a0, &a1, y);//Вычисление коэффициентов a0 и a1 для первой регрессии

/*Вычисление значений первой регрессии*/(i=0; i<n; i++)

regr1[i]=a0+a1*fi1(x[i]);

fprintf(fout,Rus("Коэффициенты а0 и а1 для первой регрессии:\n a0=%.3f\n a1=%.3f\n\n"),a0,a1);(n, Sy, fi2, x, &a0, &a1, y);//Вычисление коэффициентов a0 и a1 для второй регрессии

/*Вычисление значений второй регрессии*/(i=0; i<n; i++)

regr2[i]=a0+a1*fi2(x[i]);

fprintf(fout,Rus("Коэффициенты а0 и а1 для второй регрессии:\n a0=%.3f\n a1=%.3f\n\n"),a0,a1);

/*Вывод значений первой и второй регрессии*/=fopen(file,"r");(fout,"\xDA\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC2\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC2\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC2\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC2\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xBF");//Шапка таблицы

fprintf(fout,"\n\xB3");(fout, Rus("Название "));(fout,"\xB3 X \xB3 Y \xB3 reg1 \xB3 reg2 \xB3\n");(fout,"\xC3\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xB4");

i=0;(!feof(f))

{(s, 80, f);// Считываем строку, не более 80 символов. (Считает до ближайшего перехода на следующую строку - символа '\n').(nazv, s, 21);//Из строки s копируем первые 21 символ в строку navz

nazv[19]='\0';(fout,"\n\xB3%-20s\xB3%-8.3f\xB3%-8.3f\xB3%-10.3f\xB3%-10.3f\xB3\n",Rus(nazv),x[i],y[i],regr1[i],regr2[i]);(fout,"\xC3\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC5\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xB4");++;

}(fout,"\r\xC0\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC1\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC1\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC1\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC1\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xC4\xD9\n");//Низ таблицы

n=i;

/*Вычисление остаточных дисперсий первой и второй регрессии*/

OstDisp1=disp (n, y, regr1);=disp (n, y, regr2);

/*Вывод остаточных дисперсий*/(fout, Rus("\n\nОстаточная дисперсия первой регрессии: %-10.1f\nОстаточная дисперсия второй регрессии: %-10.1f\n"), OstDisp1, OstDisp2);

/*Вычисление коэффициентов Фишера*/

Fish1=Somy*Somy/OstDisp1;

Fish2=Somy*Somy/OstDisp2;

fprintf(fout,Rus("\nКоэффициент Фишера для первой регрессии: %.3f\nКоэффициент Фишера для второй регрессии: %.3f\n\n"), Fish1, Fish2);//Вывод коэффициентов Фишера(Fish1>Fish2)//Если коэффициент Фишера для первой регрессии больше, чем для второй...(fout,Rus("Первая регрессия \"лучше\""));//Выводим сообщение//Иначе...(fout,Rus("Вторая регрессия \"лучше\"\n"));//Выводим сообщение

}(f);//Закрытие файла с исходными данными(fout);//Закрытие файла с результатами

}

getch();

return(0);

}

Приложение Б

Код библиотеки bibliot.h:

#ifndef _bibliot_H_

#define _bibliot_H_s_

#include <stdio.h>

#include <conio.h>

#include <windows.h>

#include <string.h>

#include <math.h>float (*func)(float);

/*Вычисление среднего х и у, дисперсий x и у, среднеквадратичное отклонение х и у*/

void statistic(int n, float a[], float *S, float *Sr, float *Som)

{i;D;

*S=0;(i=0; i<n; i++)

*S=*S+a[i];

*Sr=*S/n;=0;(i=0; i<n; i++)=D+(a[i]-*Sr)*(a[i]-*Sr);=D/(n-1);

*Som=sqrt(D);

}

/*Вычисление коэффициента парной корреляции*/

float top (int n, float a[], float Sra, float b[], float Srb)

{i;S;=0;(i=0;i<n;i++)=S+(a[i]-Sra)*(b[i]-Srb); (S);

}

/*Вычисление результатов по регрессионным зависимостям*/

void regr(int n, float Sy, func fi, float a[], float *a0, float *a1, float b[])

{i;S1, S2, S3, d;=0;(i=0; i<n; i++)=S1+fi(a[i])*fi(a[i]);=0;(i=0; i<n; i++)=S2+fi(a[i]);=0;(i=0; i<n;i++)=S3+b[i]*fi(a[i]);=n*S1-S2*S2;

*a0=(Sy*S1-S2*S3)/d;

*a1=(n*S3-S2*Sy)/d;

}

/*Функция x=sqrt(x)*(log(sqrt(x)));*/

float fi1 (float x)

{=sqrt(x)*(log(sqrt(x)));(x);

}

/*Функция x=(x*x)/exp(1/x)*/

float fi2 (float x)

{=(x*x)/exp(1/x);(x);

}

/*Вычисление остаточных дисперсий*/

float disp (int n, float a[], float b[])

{i;S;D;=0;(i=0; i<n; i++)=S+(a[i]-b[i])*(a[i]-b[i]);=S/(n-1);(D);

}

#endif