Использование дисперсионного анализа в методических исследованиях
Министерство образования и науки Российской Федерации
ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
__________________ ____________________
наименование кафедры
Допускаю к защите
Руководитель_________________
_ __
И.О. Фамилия
______________________________
наименование темы
Курсовая работа по дисциплине
_ Математическое моделирование эксперимента________
наименование дисциплины
Выполнил студент группы _________
Курсовая работа защищена с оценкой ___________________________
Иркутск 2013 г.
Министерство образования и науки Российской Федерации
ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
ЗАДАНИЕ
НА КУРСОВОЕ ПРОЕКТИРОВАНИЕ (КУРСОВУЮ РАБОТУ)
По курсу ____Математического моделирования
эксперимента_____________ ______________________________
Студенту __________________________
(фамилия, инициалы)
Тема проекта (работы) Использование_дисперсионного_
Исходные данные_______________
1)Заполнить таблицу;______________________
2)Найти_вспомогательные_суммы;
3)Проверить_однородность_
4)Рассчитать_относительные_
5)Сделать_вывод_об_
Рекомендуемая литература ______________________________
2)Гмурман В.Е. Теория вероятностей и математическая статистика. – М.: Высшая школа, 2003.-523с.
3)Никаноров А.В. Математическое моделирование эксперимента, учебное пособие_2008. -108 с,
______________________________
Графическая часть на ______________ листах.
Дата выдачи задания
“______” __________________________2013
Дата представления проекта руководителю “______” ___________2013 г.
Руководитель курсового проектирования (курсовой работы) _____
Содержание
- Введение 4
- Основная часть 5
2.1)Однофакторный
2.2)Двухфакторный
2.3)Многофакторный
- Задача для курсовой работы 18
3.1)Решение задачи для курсовой работы 20
- Заключение 22
- Список используемых источников 23
Введение
Проникновение математических методов в самые разнообразные, подчас неожиданные сферы человеческой деятельности дает возможность пользоваться новыми, как правило, весьма плодотворными средствами исследования. Рост математической культуры специалистов в соответствующих областях приводит к тому, что изучение общих теоретических положений и методов вычислений уже не встречает серьезных трудностей. Вместе с тем на практике оказывается, что одних лишь математических познаний далеко не достаточно для решения той или иной прикладной задачи – необходимо еще получить навыки в переводе исходной формулировки задачи на математический язык. В этом и состоит проблема овладения искусством математического моделирования.
В рамках же нашей работы стоит цель овладеть искусством дисперсионного анализа. В чем же оно заключается?
Оценивая разброс результатов исследований вокруг среднего значения, определяем дисперсию воспроизводимости методики исследования (анализа), характеризующую сумму погрешностей, вносимых на различных этапах исследования (анализа).
Экспериментатору, когда он проводит исследования с целью повышения точности исследования или анализа, необходимо знать, на каких этапах вносится наибольший вклад в суммарную дисперсию. Например, при рентгеноспектральном анализе образца погрешности могут быть внесены при отборе анализируемой навески (погрешности неоднородности материала), при изготовлении излучателя, измерении интенсивности аналитической линии, расчете содержания.
В зависимости от выбранного способа пробоподготовки второй этап можно разбить на ряд подэтапов. При получении излучателей в виде стеклянных дисков можно выделить погрешности, вносимые при взвешивании и перемешивании пробы и флюса, сплавлении, отливке излучателя.
Для повышения воспроизводимости
методики исследования прежде всего
следует совершенствовать тот этап,
на котором вносится наибольшая погрешность.
Основная часть
Дисперсионный анализ (от латинского Dispersio – рассеивание) –статистический метод, позволяющий анализировать влияние различных факторов на исследуемую переменную. Метод был разработан биологом Р. Фишером в 1925 году и применялся первоначально для оценки экспериментов в растениеводстве. В дальнейшем выяснилась общенаучная значимость дисперсионного анализа для экспериментов в психологии, педагогике, медицине и др.
Целью дисперсионного анализа является проверка значимости различия между средними с помощью сравнения дисперсий. Дисперсию измеряемого признака разлагают на независимые слагаемые, каждое из которых характеризует влияние того или иного фактора или их взаимодействия. Последующее сравнение таких слагаемых позволяет оценить значимость каждого изучаемого фактора, а также их комбинации.
При истинности нулевой гипотезы (о равенстве средних в нескольких группах наблюдений выбранных из генеральной совокупности), оценка дисперсии, связанной с внутригрупповой изменчивостью, должна быть близкой к оценке межгрупповой дисперсии.
При проведении исследования рынка часто встает вопрос о сопоставимости результатов. Например, проводя опросы по поводу потребления какого-либо товара в различных регионах страны, необходимо сделать выводы, на сколько данные опроса отличаются или не отличаются друг от друга. Сопоставлять отдельные показатели не имеет смысла и поэтому процедура сравнения и последующей оценки производится по некоторым усредненным значениям и отклонениям от этой усредненной оценки. Изучается вариация признака. За меру вариации может быть принята дисперсия. Дисперсия σ2– мера вариации, определяемая как средняя из отклонений признака, возведенных в квадрат.
На практике часто возникают задачи более общего характера – задачи проверки существенности различий средних выборочных нескольких совокупностей. Например, требуется оценить влияние различного сырья на качество производимой продукции, решить задачу о влиянии количества удобрений на урожайность с/х продукции.
Иногда дисперсионный анализ применяется, чтобы установить однородность нескольких совокупностей (дисперсии этих совокупностей одинаковы по предположению; если дисперсионный анализ покажет, что и математические ожидания одинаковы, то в этом смысле совокупности однородны). Однородные же совокупности можно объединить в одну и тем самым получить о ней более полную информацию, следовательно, и более надежные выводы.
При проведении дисперсионного анализа
должны выполняться следующие статистические
допущения: независимо от уровня фактора
величины отклика имеют нормальный (Гауссовский)
закон распределения и одинаковую дисперсию.
Такое равенство дисперсий называется
гомогенностью. Таким образом, изменение
способа обработки сказывается лишь на
положении случайной величины отклика,
которое характеризуется средним значением
или медианой. Поэтому все наблюдения
отклика принадлежат сдвиговому
семейству нормальных распределений.
Говорят, что техника дисперсионного анализа
является "робастной". Этот термин,
используемый статистиками, означает,
что данные допущения могут быть в некоторой
степени нарушены, но несмотря на это,
технику можно использовать.
При неизвестном законе распределения
величин отклика используют непараметрические
(чаще всего ранговые) методы анализа.
В основе дисперсионного анализа лежит
разделение дисперсии на части или компоненты.
Вариацию, обусловленную влиянием фактора,
положенного в основу группировки, характеризует
межгрупповая дисперсия σ2. Она является
мерой вариации частных средних по группам
вокруг общей средней
и определяется по формуле:
,
где k - число групп;
nj - число единиц в j-ой группе;
- частная средняя по j-ой группе;
- общая средняя по совокупности единиц.
Вариацию, обусловленную влиянием прочих
факторов, характеризует в каждой группе
внутригрупповая дисперсия σj2.
.
Между общей дисперсией σ02,
внутригрупповой дисперсией σ2
и межгрупповой дисперсией
существует соотношение:
σ02 =
+ σ2.
Внутригрупповая дисперсия объясняет
влияние неучтенных при группировке факторов,
а межгрупповая дисперсия объясняет влияние
факторов группировки на среднее значение
по группе.
Однофакторный дисперсионный
анализ
Однофакторная дисперсионная модель
имеет вид:
xij = μ + Fj
+ εij,
где хij – значение исследуемой переменой,
полученной на i-м уровне фактора (i=1,2,...,т)
c j-м порядковым номером (j=1,2,...,n);
Fi – эффект, обусловленный влиянием
i-го уровня фактора;
εij – случайная компонента,
или возмущение, вызванное влиянием неконтролируемых
факторов, т.е. вариацией переменой внутри
отдельного уровня.
Основные предпосылки дисперсионного
анализа:
1)математическое ожидание возмущения εij равно нулю для любых i, т.е.
M(εij) = 0; (2)
2)возмущения εij взаимно независимы;
3)дисперсия переменной xij (или возмущения εij) постоянна для
любых i, j, т.е. D(εij) = σ2; (3)
4)переменная xij (или возмущение εij)
имеет нормальный закон
распределения N(0;σ2).
Влияние уровней фактора может
быть как фиксированным или систематическим
(модель I), так и случайным (модель II).
Пусть, например, необходимо выяснить,
имеются ли существенные различия между
партиями изделий по некоторому показателю
качества, т.е. проверить влияние на качество
одного фактора - партии изделий. Если
включить в исследование все партии сырья,
то влияние уровня такого фактора систематическое
(модель I), а полученные выводы применимы
только к тем отдельным партиям, которые
привлекались при исследовании. Если же
включить только отобранную случайно
часть партий, то влияние фактора случайное
(модель II). В многофакторных комплексах
возможна смешанная модель III, в которой
одни факторы имеют случайные уровни,
а другие – фиксированные.
Пусть имеется m партий изделий. Из каждой
партии отобрано соответственно n1,
n2, …, nm изделий (для простоты
полагается, что n1=n2=...=nm=n).
Значения показателя качества этих изделий
представлены в матрице наблюдений:
x11 x12 … x1n
x21 x22 … x2n
………………… = (xij), (i = 1,2, …,
m; j = 1,2, …, n).
xm1 xm2 … xmn
Необходимо проверить
Если полагать, что элементы строк матрицы
наблюдений – это численные значения
случайных величин Х1,Х2,...,Хm,
выражающих качество изделий и имеющих
нормальный закон распределения с математическими
ожиданиями соответственно a1,а2,...,аm
и одинаковыми дисперсиями σ2, то
данная задача сводится к проверке нулевой
гипотезы Н0: a1=a2 =...= аm,
осуществляемой в дисперсионном анализе.
Усреднение по какому-либо индексу обозначено
звездочкой (или точкой) вместо индекса,
тогда средний показатель качества изделий
i-й партии, или групповая средняя для i-го
уровня фактора, примет вид:
где
i* – среднее значение по столбцам;
ij – элемент матрицы наблюдений;
n – объем выборки.
А общая средняя:
Сумма квадратов отклонений наблюдений
хij от общей средней
** выглядит так:
2=
2+
2+
+2
2. (6)
или
Q = Q1 + Q2 + Q3.
Последнее слагаемое равно нулю
=0. (7)
так как сумма отклонений значений
переменной от ее средней равна нулю,
т.е.
2=0.
Первое слагаемое можно записать в виде:
В результате получается тождество: Q
= Q1 + Q2,
где
- общая, или полная, сумма квадратов отклонений;
- сумма квадратов отклонений групповых
средних от общей средней, или межгрупповая
(факторная) сумма квадратов отклонений;
- сумма квадратов отклонений наблюдений
от групповых средних, или внутригрупповая
(остаточная) сумма квадратов отклонений.
В разложении (8) заключена основная идея
дисперсионного анализа. Применительно
к рассматриваемой задаче равенство (8)
показывает, что общая вариация показателя
качества, измеренная суммой Q, складывается
из двух компонент – Q1
и Q2, характеризующих изменчивость
этого показателя между партиями (Q1)
и изменчивость внутри партий (Q2),
характеризующих одинаковую для всех
партий вариацию под воздействием неучтенных
факторов.
В дисперсионном анализе анализируются
не сами суммы квадратов отклонений, а
так называемые средние квадраты, являющиеся
несмещенными оценками соответствующих
дисперсий, которые получаются делением
сумм квадратов отклонений на соответствующее
число степеней свободы.
Число степеней свободы определяется
как общее число наблюдений минус число
связывающих их уравнений. Поэтому для
среднего квадрата S12, являющегося
несмещенной оценкой межгрупповой дисперсии,
число степеней свободы k1=m-1, так
как при его расчете используются m групповых
средних, связанных между собой одним
уравнением (5). А для среднего квадрата
S22, являющегося несмещенной
оценкой внутригрупповой дисперсии, число
степеней свободы k2=mn-m, т.к. при ее
расчете используются все mn наблюдений,
связанных между собой m уравнениями (4).
Таким образом:
= Q1/(m-1),
= Q2/(mn-m).
Если найти математические ожидания средних
квадратов
и
, подставить в их формулы выражение xij
(1) через параметры модели, то получится:
т.к. с учетом свойств математического
ожидания
а
(10)
Для модели I с фиксированными уровнями
фактора Fi(i=1,2,...,m) – величины неслучайные,
поэтому
M(S
) =
2 /(m-1) +σ2.
Гипотеза H0 примет вид Fi = F*(i
= 1,2,...,m), т.е. влияние всех уровней фактора
одно и то же. В случае справедливости
этой гипотезы
M(S
)= M(S
)= σ2.
Для случайной модели II слагаемое
Fi в выражении (1) – величина
случайная. Обозначая ее дисперсией
получим из (9)
и, как и в модели I
M(
)= σ2.
В таблице №1 представлен общий
вид вычисления значений, с помощью дисперсионного
анализа.
Таблица №1 – Базовая таблица дисперсионного
анализа
Компоненты дисперсии |
Сумма квадратов |
Число степеней свободы |
Средний квадрат |
Математическое ожидание среднего квадрата |
Межгрупповая |
|
m-1 |
= Q1/(m-1) |
|
Внутригрупповая |
|
mn-m |
= Q2/(mn-m) |
M(S
)= σ2 |
Общая |
|
mn-1 |
Гипотеза H0 примет вид σF2
=0. В случае справедливости этой гипотезы
M(
)= M(
)= σ2.
В случае однофакторного комплекса как
для модели I, так и модели II средние квадраты
S2 и S2, являются несмещенными
и независимыми оценками одной и той же
дисперсии σ2.
Следовательно, проверка нулевой гипотезы
H0 свелась к проверке существенности
различия несмещенных выборочных оценок
S2
и S
дисперсии σ2.
Гипотеза H0 отвергается, если фактически
вычисленное значение статистики
F = S
/S
больше критического Fα:K1:K2, определенного
на уровне значимости α при числе степеней
свободы k1=m-1 и k2=mn-m, и принимается,
если F < Fα:K1:K2
.
F- распределение Фишера (для x > 0) имеет
следующую функцию плотности (для
= 1, 2, ...;
= 1, 2, ...):
где
- степени свободы;
Г - гамма-функция.
Применительно к данной задаче опровержение
гипотезы H0 означает наличие существенных
различий в качестве изделий различных
партий на рассматриваемом уровне значимости.
Для вычисления сумм квадратов Q1,
Q2, Q часто бывает удобно использовать
следующие формулы:
т.е. сами средние, вообще говоря, находить
не обязательно.
Таким образом, процедура однофакторного
дисперсионного анализа состоит в проверке
гипотезы H0 о том, что имеется одна
группа однородных экспериментальных
данных против альтернативы о том, что
таких групп больше, чем одна. Под однородностью
понимается одинаковость средних значений
и дисперсий в любом подмножестве данных.
При этом дисперсии могут быть как известны,
так и неизвестны заранее. Если имеются
основания полагать, что известная или
неизвестная дисперсия измерений одинакова
по всей совокупности данных, то задача
однофакторного дисперсионного анализа
сводится к исследованию значимости различия
средних в группах данных.
Двухфакторный дисперсионный анализ
Принимается аддитивная и независимая модель действия факторов:
, причем , . (15)
Величины aj и bi называются вкладами факторов. Последние два условия всегда можно выполнить масштабированием величин aj и bi за счет изменения величины c.
Для каждого наблюдения из рассматриваемой совокупности справедливо уравнение:
xij = c + aj + bi + e ij, i =1, ..., n; j =1, ..., k. (16)
Обычно наблюдения представляют структурной таблицей статистического комплекса. Приведем простейший двухфакторный комплекс, в которой каждому сочетанию (Aj, Bj) уровней (градаций) факторов, т.е. одной клетке таблицы, соответствует одно наблюдение (в таблице сочетание символов “( )^” обозначает статистическую оценку групповых средних): Разложение результатов измерения при двухфакторном дисперсионном анализе представлены в таблице №2.
Таблица №2:
Фактор B |
Фактор A |
Средние по строкам |
B1 B2 ... Bn |
x11 x12 ... x1k x21 x22 ... x2k ... xn1 xn2 ... xnk |
x1· =(c+b1)^ x2· =(c+b2)^ ... xn· =(c+bn)^ |
Средние по столбцам |
x· 1= x· 2= x· k= |
x· · =c^ |
Основное тождество
Оценки c, bi, aj могут быть получены с помощью метода наименьших квадратов (МНК) минимизацией суммы (17)
Основываясь на МНК-оценках
, (18)
введем следующие обозначения:
- для сумм квадратов отклонений под влиянием k уровней фактора А и n уровней фактора В:
, ; (19)
- для остаточной суммы квадратов:
; (20)
- для полной суммы квадратов наблюдений относительно общего среднего :
. (21)
Тогда справедливо следующее
, (22)
т.е. полная сумма квадратов отклонений является суммой квадратов вкладов по факторам и квадратов случайных отклонений (или остатков ). Другими словами, полное рассеяние есть сумма вариации под влиянием факторов и случайной составляющей.
Проверка нулевых гипотез
По имеющимся наблюдениям
HA: a1 = a2 = . . . = ak = 0
HB: b1 = b2 = . . . = bn = 0 .
Основой процедуры проверки гипотезы является сравнение двух статистически независимых оценок дисперсии s 2.
Одна из них, s 2* оценивает дисперсию вне зависимости от того, верна или нет гипотеза HA (или HВ) и основана на сумме квадратов случайных отклонений:
. (23)
Другая, s 2** оценивает дисперсию, если HA (или HВ) верна. Для гипотезы HA эта дисперсия основана на сумме квадратов разностей “между столбцами”, т.е. по уровням фактора A:
. (24)
Если гипотеза HA верна, то отношение
(25)
имеет F-распределение Фишера с (k -1) и r степенями свободы. Если
FA ³ F1-a , (26)
где F1-a – квантиль этого распределения порядка 1- a , a – выбранный уровень значимости, то гипотеза HA отклоняется.
Вместо (5.54) можно использовать эквивалентную процедуру: гипотеза HA отклоняется, если
P{ F ³ FA } £ a, (27)
где P{ F ³ FA } – вероятность при справедливости HA получить значение, большее, чем FA; F – случайная величина, имеющая распределение Фишера.
Для проверки гипотезы HВ используют сумму квадратов разностей "между строками", то есть по уровням фактора B: . (28)
Аналогичным образом, если отношение велико, то гипотеза HB отклоняется.
Многофакторный дисперсионный
анализ
Следует сразу же отметить, что
принципиальной разницы между многофакторным
и однофакторным дисперсионным анализом
нет. Многофакторный анализ не меняет
общую логику дисперсионного анализа,
а лишь несколько усложняет ее, поскольку,
кроме учета влияния на зависимую переменную
каждого из факторов по отдельности, следует
оценивать и их совместное действие. Таким
образом, то новое, что вносит в анализ
данных многофакторный дисперсионный
анализ, касается в основном возможности
оценить меж-факторное взаимодействие.
Тем не менее, по-прежнему остается возможность
оценивать влияние каждого фактора в отдельности.
В этом смысле процедура многофакторного
дисперсионного анализа (в варианте ее
компьютерного использования) несомненно
более экономична, поскольку всего за
один запуск решает сразу две задачи: оценивается
влияние каждого из факторов и их взаимодействие.
Общая схема двухфакторного эксперимента,
данные которого обрабатываются дисперсионным
анализом имеет вид:
Зависимая переменная xi |
|
Взаимодействие факторов A и B |
Прочие не учитываемые (случайные) факторы |
Фактор B: |
Фактор А: |
|
Рисунок 1.1 – Схема двухфакторного эксперимента
Данные, подвергаемые многофакторному
дисперсионному анализу, часто обозначают
в соответствии с количеством факторов
и их уровней.
Предположив, что в
рассматриваемой задаче о качестве различных
m партий изделия изготавливались на разных
t станках и требуется выяснить, имеются
ли существенные различия в качестве изделий
по каждому фактору:
А - партия изделий;
B - станок.
В результате получается переход к задаче
двухфакторного дисперсионного анализа.
Все данные представлены в таблице №3,
в которой по строкам - уровни Ai фактора
А, по столбцам — уровни Bj фактора
В, а в соответствующих ячейках, таблицы
находятся значения показателя качества
изделий xijk (i=1,2,...,m; j=1,2,...,l; k=1,2,...,n).
Таблица №3 – Показатели качества изделий
B1 |
B2 |
… |
Bj |
… |
Bl | |
A1 |
x11l,…,x11k |
x12l,…,x12k |
… |
x1jl,…,x1jk |
… |
x1ll,…,x1lk |
|
A2 |
x21l,…,x21k |
x22l,…,x22k |
… |
x2jl,…,x2jk |
… |
x2ll,…,x2lk |
|
… |
… |
… |
… |
… |
… |
… |
Ai |
xi1l,…,xi1k |
xi2l,…,xi2k |
… |
xijl,…,xijk |
… |
xjll,…,xjlk |
|
… |
… |
… |
… |
… |
… |
… |
Am |
xm1l,…,xm1k |
xm2l,…,xm2k |
… |
xmjl,…,xmjk |
… |
xmll,…,xmlk |
Двухфакторная дисперсионная модель
имеет вид:
где xijk - значение наблюдения в ячейке
ij с номером k;
μ - общая средняя;
Fi - эффект, обусловленный влиянием
i-го уровня фактора А;
Gj - эффект, обусловленный влиянием
j-го уровня фактора В;
Iij - эффект, обусловленный взаимодействием
двух факторов, т.е. отклонение от средней
по наблюдениям в ячейке ij от суммы первых
трех слагаемых в модели (29);
εijk - возмущение, обусловленное
вариацией переменной внутри отдельной
ячейки.
Предполагается, что εijk имеет нормальный
закон распределения N(0; с2), а все
математические ожидания F*,
G*, Ii*, I*j равны нулю.
Групповые средние находятся по
формулам:
- в ячейке:
,
по строке:
по столбцу:
общая средняя:
В таблице №4 представлен общий вид вычисления
значений, с помощью дисперсионного анализа.
Таблица №4 – Базовая таблица дисперсионного
анализа
Компоненты дисперсии |
Сумма квадратов |
Число степеней свободы |
Средние квадраты |
Межгрупповая (фактор А) |
|
m-1 |
|
Межгрупповая (фактор B) |
|
l-1 |
|
Взаимодействие |
|
(m-1)(l-1) |
|
Остаточная |
|
mln - ml |
|
Общая |
|
mln - 1 |
Проверка нулевых гипотез HA, HB, HAB об отсутствии влияния на рассматриваемую переменную факторов А, B и их взаимодействия AB осуществляется сравнением отношений
Если n=1, т.е. при одном наблюдении в ячейке, то не все нулевые гипотезы могут быть проверены так как выпадает компонента Q3 из общей суммы квадратов отклонений, а с ней и средний квадрат
С точки зрения техники вычислений для нахождения сумм квадратов Q1, Q2, Q3, Q4, Q целесообразнее использовать формулы:
Q3 = Q – Q1 – Q2 – Q4.
Отклонение от основных предпосылок дисперсионного анализа — нормальности распределения исследуемой переменной и равенства дисперсий в ячейках (если оно не чрезмерное) — не сказывается существенно на результатах дисперсионного анализа при равном числе наблюдений в ячейках, но может быть очень чувствительно при неравном их числе. Кроме того, при неравном числе наблюдений в ячейках резко возрастает сложность аппарата дисперсионного анализа. Поэтому рекомендуется планировать схему с равным числом наблюдений в ячейках, а если встречаются недостающие данные, то возмещать их средними значениями других наблюдений в ячейках. При этом, однако, искусственно введенные недостающие данные не следует учитывать при подсчете числа степеней свободы.
Задача для курсовой работы

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