Методы математического моделирования и обработки данных. Практикум
.pdfМИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ
РОССИЙСКОЙ ФЕДЕРАЦИИ
Федеральное государственное бюджетное образовательное учреждение высшего образования
«Воронежский государственный технический университет»
П. С. Куприенко, А. В. Кочегаров, Е. П. Вялова, Н. В. Ильина
МЕТОДЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ И ОБРАБОТКИ ДАННЫХ
Практикум
Воронеж 2025
1
УДК 519.8:519.2:004.42(075.8) ББК 87.256.631.0я7
М54
Рецензенты:
кафедра информационных систем и технологий Воронежского института высоких технологий;
С. В. Дорохин, д-р техн. наук, доцент, декан автомобильного факультета Воронежского государственного лесотехнического университета
им. Г. Ф. Морозова
Методы математического моделирования и обработки данных:
практикум [Электронный ресурс]/ П. С. Куприенко, А. В. Кочегаров, Е. П. Вялова, Н. В. Ильина. - Электрон. текстовые и граф. данные М54 (1,5 Мб). – Воронеж: ФГБОУ ВО «Воронежский государственный технический университет», 2025. – 1 электрон. опт. диск (DVD-ROM):
цв. – Сиcтем. требования: ПК 500 и выше; 256 Мб ОЗУ; Windows 10; SVGA c разрешением экрана 1024х768; Adobe Acrobat; CD-ROM
дисковод; мышь. – Загл. с экрана.
ISBN 978-5-7731-1237-2
В практикуме рассмотрены способы решения практических задач по курсу «Методы математического моделирования и обработки данных» с помощью пакета математических вычислений Scilab и Matlab. Показаны методы постановки и решения задач математического моделирования и статистического анализа, приведены варианты для индивидуального решения задач.
Практикум предназначен для студентов направления 20.04.01 «Техносферная безопасность» всех форм обучения, а также для широкого круга заинтересованных читателей.
Ил. 9. Табл. 10. Библиогр.: 10 назв.
УДК 519.8:519.2:004.42(075.8) ББК 87.256.631.0я7
Издается по решению редакционно-издательского совета Воронежского государственного технического университета
ISBN 978-5-7731-1237-2 © Куприенко П. С., Вялова Е. П.,
Ильина Н. В., Кочегаров А. В., 2025
©ФГБОУ ВО «Воронежский государственный технический университет», 2025
2
ОГЛАВЛЕНИЕ |
|
Введение......................................................................................................... |
4 |
Практическая работа № 1. Определение характеристик статистического |
|
распределения…………………………………………………………………. |
5 |
Практическая работа № 2. Статистические оценки параметров |
|
распределения…………………………………………………………………. |
13 |
Практическая работа № 3. Построение линейной регрессии по |
|
несгруппированным данным………………………………………………… |
19 |
Практическая работа № 4. Анализ распределения данных: проверка |
|
соответствия нормальному закону…………………………………………… |
27 |
Практическая работа № 5. Обработка экспериментальных данных методом |
|
наименьших квадратов ………………………………………………………… |
34 |
Практическая работа № 6. Полный факторный эксперимент в |
|
регрессионном анализе………………………………………………………………. |
41 |
Практическая работа № 7. Исследование алгоритма нечеткой |
|
кластеризации…………………………………………………………………. |
56 |
Заключение........................................................................................................ |
61 |
Библиографический список............................................................................ |
62 |
3
ВВЕДЕНИЕ
Современный мир немыслим без анализа данных и математического моделирования, которые стали неотъемлемой частью научных исследований, инженерных разработок и принятия решений в различных сферах деятельности. Умение работать с данными, интерпретировать их и строить математические модели является ключевым навыком для специалистов в области естественных наук, техники, информационных технологий и техносферной безопасности.
Данный практикум «Методы математического моделирования и обработки данных. Практикум» предназначен для студентов, магистров, аспирантов и всех, кто интересуется применением математических методов в практических задачах. В практикуме рассматриваются основные подходы к обработке данных, методам математической статистики и построению моделей с использованием среды Scilab – мощного инструмента для научных вычислений.
Особое внимание уделено практической стороне: каждая тема сопровождается примерами и заданиями, которые помогут закрепить теоретические знания и развить навыки работы с реальными данными. Практикум также включает рекомендации по выбору методов анализа, интерпретации результатов и визуализации данных.
Цель практикума – помочь читателю не только освоить инструменты Scilab, но и научиться применять их для решения актуальных задач, что делает его полезным как для учебного процесса, так и для профессиональной деятельности.
4
ПРАКТИЧЕСКАЯ РАБОТА №1 ОПРЕДЕЛЕНИЕ ХАРАКТЕРИСТИК СТАТИСТИЧЕСКОГО
РАСПРЕДЕЛЕНИЯ
Цель работы
1.Приобрести практические навыки расчета основных числовых характеристик выборки:
o выборочного среднего o моды и медианы
o выборочного среднего квадратического отклонения
2.Освоить методы визуализации статистических данных:
o построение гистограмм распределения
o создание графиков эмпирической функции распределения
Исходные данные
Статистические данные о количестве лесных пожаров в различных областях за отчетный период, сгруппированные по интервалам значений.
Количество |
Менее 20 |
20–40 |
40–60 |
60–80 |
80–100 |
Более 100 |
пожаров |
|
|
|
|
|
|
Число областей |
3 |
6 |
17 |
25 |
10 |
4 |
Задачи исследования
1.Рассчитать основные статистические характеристики распределения:
o среднее арифметическое количество пожаров
o показатели центра распределения (моду и медиану)
o меру разброса данных (среднее квадратическое отклонение)
2.Выполнить графическое представление данных:
o построить гистограмму частот распределения
o создать график эмпирической функции распределения 3. Провести анализ полученных результатов:
o определить характер распределения данных
o сделать выводы о вариабельности количества пожаров по областям
Методика выполнения
Работа выполняется с использованием пакета прикладных программ для статистического анализа (например, Scilab, Excel). Все расчеты сопровождаются заполнением расчетной таблицы с промежуточными и итоговыми результатами.
Порядок выполнения работы
Все задания мы будем выполнять в системе Scilab – это программная среда для математических вычислений, предназначенная для решения инженерных и научных задач. Scilab также предоставляет широкие возможности для создания и редактирования различных типов графиков и поверхностей.
5
Система Scilab включает в себя множество встроенных команд, операторов и функций. При этом пользователь может самостоятельно создавать новые команды или функции и использовать их вместе со встроенными инструментами.
Scilab – это мощная программная среда для выполнения инженерных и научных расчетов, которая сочетает в себе простоту освоения с широкими вычислительными возможностями. Она активно применяется для анализа данных, математического моделирования, а также визуализации результатов в виде двумерных и трёхмерных графиков, диаграмм и поверхностей.
Ключевые особенности системы 1. Богатый набор инструментов:
o Встроенные математические функции, операторы и библиотеки для решения задач линейной алгебры, статистики, оптимизации и других областей. o Возможность создания пользовательских функций и интеграции их со
стандартными командами.
2. Интуитивный интерфейс:
o После запуска Scilab открывается главное окно с меню, панелью инструментов и рабочей областью.
o Готовность системы к работе отображается знаком приглашения --> в
командной строке.
3.Процесс выполнения команд:
o Ввод команд осуществляется с клавиатуры непосредственно в командной строке.
o Для выполнения команды и отображения результата достаточно нажать Enter (пример на рис. 1.1).
Практические преимущества
Гибкость настройки под конкретные задачи за счет комбинации встроенных и пользовательских функций.
Удобная визуализация данных, включая построение графиков в реальном времени.
Поддержка скриптов для автоматизации сложных вычислений.
Среда Scilab идеально подходит как для учебных проектов, так и для профессиональных инженерных решений, делая процесс работы с данными наглядным
После запуска Scilab на экране отображается главное окно программы. Оно состоит из меню, панели инструментов и рабочей области. О том, что система готова к выполнению команд, свидетельствует наличие знака приглашения (-->). Рабочую область с этим знаком называют командной строкой. Ввод команд в Scilab осуществляется с клавиатуры. Нажатие клавиши Enter приводит к выполнению команды и выводу результата на экран
(рис. 1.1).
6
Рис. 1.1. Рабочая область Scilab
Записи, расположенные за знаком // являются комментариями и в вычислениях не участвуют, приводятся для пояснения работы.
В ходе работы необходимо провести необходимые вычисления и заполнить расчётную таблицу (табл. 1.1)
Таблица 1.1
Расчётная таблица
Интервалы |
|
ai |
xi |
ni |
xi ni |
|
wi |
fi |
sni |
swi |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
0–20 |
0 |
20 |
|
3 |
|
|
|
|
|
|
20–40 |
20 |
40 |
|
6 |
|
|
|
|
|
|
40–60 |
40 |
60 |
|
17 |
|
|
|
|
|
|
60–80 |
60 |
80 |
|
25 |
|
|
|
|
|
|
80–100 |
80 |
100 |
|
10 |
|
|
|
|
|
|
100–120 |
100 |
120 |
|
4 |
|
|
|
|
|
|
Итого
//Листинг программы
// Ввод данных: частоты ni, левые границы интервалов ai–1, правые границы ai. -->n=[3 6 17 25 10 4];
-->left_a=[0 20 40 60 80]; -->right_a=[20 40 60 80 100 120]; -->x=(left_a+right_a)./2;
-->xn=x.*n // результаты вычислений внести в столбец 6
7
-->x2n=x.^2.*n // результаты вычислений внести в столбец 7 -->sum1=sum[x2n] // результаты вычислений внести в столбец 7, в
последнюю ячейку
-->sum2=sum[xn] // результаты вычислений внести в столбец 6, в последнюю ячейку
//Вычислим характеристики выборки: -->h=right_a-left_a; //вычисление ширины интервала -->xsr=sum2/sum[n];//значения выборочного среднего -->x2sr=sum1/sum[n];//квадрат выборочного среднего -->D=x2sr-xsr^2; // выборочная дисперсия
-->sigma =sqrt(D); // среднеквадратичное отклонение //Раcсчитаем относительные частоты wi
-->w=n./sum(n) // результаты вычислений внести в столбец 8 //Раcсчитаем плотности относительной частоты fi -->f=w./h // результаты вычислений внести в столбец 9
//вычислим накопленные частоты sni и накопленные относительные частот
swi,
-->sn=[0 f.*h] // результаты вычислений внести в столбец 10 -->sw=[0 w.*h] // результаты вычислений внести в столбец 11 // вычислим значения моды и медианы.
-->[max_value, index_max]=max(n) //находим интервал с наибольшей частотой
-->index=ans; -->Mo=left_a[index]+((n[index]-n[index-1])/(n[index]-n[index-1]+n[index]-
n[index+1]))*h[index] //вычисление моды -->k=sum[n]/2;
-->[max_value, index_max]=max(n) //находим интервал с наибольшей частотой
-->indexMe=ans -->Me=left_a[index_Me]+(k-sn[indexMe-1]/n[indexMe])*(right_a[indexMe]-
left_a[index_Me])
//Строим гистограмму относительных частот вариационного ряда. -->hf=w;
-->histogram(x, hf)
Результаты построения представлены на рис. 1.2
8
Рис. 1.2. Гистограмма
//Строим эмпирическую функцию распределения
-->for nx<x -->F= nx/sum(n); -->plot(x, F)
Результаты построения представлены на рис. 1.3
Рис 1.3. График функции по эмпирическим данным
По итогам выполнения работы
1.Будет создан файл с результатами Lab_1.scl.
2.Все вычисления необходимо заносить в расчётную таблицу.
Требования к отчёту
Титульный лист работы;
Таблицы с исходными данными и результатами расчётов;
Гистограмма;
График эмпирической функции распределения.
9
Варианты индивидуального задания
В ходе обследования 1000 пожароопасных объектов региона была сформирована случайная выборка из 100 объектов. На основе полученных данных выявлено распределение нарушений техники безопасности:
Вариант 1
|
Время работы |
|
|
|
|
|
|
|
|
|
|
30–40 |
|
|
|
40–50 |
|
50–60 |
|
|
|
60–70 |
|
||||||||||||
|
технологической линии, час |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
Количество инцидентов |
|
|
|
|
|
|
8 |
|
|
|
|
12 |
|
|
|
|
22 |
|
|
|
|
10 |
|
|
||||||||||
|
Вариант 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время работы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
технологической |
|
|
|
20–30 |
30–40 |
|
40–50 |
|
|
50–60 |
|
|
60–70 |
70–80 |
|
|
|||||||||||||||||||
линии, час |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Количество |
|
|
|
|
8 |
|
19 |
|
|
|
28 |
|
|
|
|
|
32 |
|
|
|
42 |
|
|
|
21 |
|
|
||||||||
инцидентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Вариант 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время работы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
технологической |
50 |
|
50–100 |
100–150 |
|
150–200 |
200–250 |
|
250–300 |
|
|
||||||||||||||||||||||||
линии, час |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Количество |
10 |
|
|
20 |
|
|
|
30 |
|
|
40 |
|
|
|
|
60 |
|
|
|
|
|
|
80 |
|
|
||||||||||
инцидентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Вариант 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время работы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30 и |
|
||
технологической |
|
|
До 5 |
|
|
5–10 |
10–15 |
15–20 |
|
|
20–25 |
|
|
25–30 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
более |
|
|||||||||||||||||||||||||
линии, час |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Количество |
|
|
1 |
|
|
1 |
|
|
2 |
|
|
4 |
|
|
|
|
6 |
|
|
|
6 |
|
|
|
8 |
|
|
||||||||
инцидентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Вариант 5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время работы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Более |
|
||
технологической |
|
|
|
До 2 |
|
2–5 |
|
5–7 |
|
|
7–9 |
|
|
9–11 |
|
|
11–13 |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
13 |
|
|
|||||||||||||||||||||
линии, час |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Количество |
|
|
12 |
|
|
|
|
13 |
|
14 |
|
10 |
|
|
|
5 |
|
|
|
3 |
|
|
|
1 |
|
|
|||||||||
инцидентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Вариант 6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время работы |
|
|
До |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Более |
|
|||
технологической |
|
|
|
|
20–50 |
|
50–80 |
|
80–110 |
|
110–140 |
|
|
140–170 |
|
|
|||||||||||||||||||
|
|
20 |
|
|
|
|
|
|
|
|
170 |
|
|
||||||||||||||||||||||
линии, час |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Количество |
|
|
25 |
|
21 |
|
|
18 |
|
|
14 |
|
|
|
|
10 |
|
|
|
3 |
|
|
|
|
7 |
|
|
||||||||
инцидентов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
10
