CC BY
28
Выводы
Как показывает предварительный этап исследования, при выборе почвообрабатывающих агрегатов предпочтение может быть отдано агрегатам большей ширины захвата как имеющим меньшее тяговое сопротивление. При этом необходимо учесть тот факт, что удельный расход топлива минимален для агрегатов меньшей ширины. Возникает противоречивая ситуация, разрешение которой требует проведения дополнительных исследований.
Список литературы
1. Протоколы испытаний / ФГУ «Сибирская государственная зональная машиноиспытательная станция». - Сосновское, 2002-2010 гг.
2. Крейтер, С.В. Основы конструирования и агрегатирования : учеб. пособие / С.В. Крейтер, А.Р. Нестеров, В.В. Данилевский. - М. : Изд-во стандартов, 1983. - 224 с.
3. Карпенко, А.Н. Сельскохозяйственные машины / А.Н. Карпенко, В.М. Халанский. - М., Колос, 1983. -
495 с.
SUMMARY
G. V. Redreev, S.N. Shhetinina
On the question of technical characteristics of the study soil cultivating units
The article posed the problem of choosing from among soil cultivating machines on the market of agricultural machinery. A review of soil cultivating units tested in the machine-zonal station. Proposed solution to the problem of choosing soil cultivating units.
Key words: soil cultivating machines, technical specifications, research, problem selection.
УДК 621.318
М.В. Семеняк, В.К. Федоров
ВЫБОР ИЗМЕРИТЕЛЬНЫХ УСТРОЙСТВ ЕСТЕСТВЕННОЙ ОСВЕЩЁННОСТИ, ИХ МАТЕМАТИЧЕСКИХ ОПИСАНИЙ И МЕТОДИК ИССЛЕДОВАНИЯ
В данной статье на основе обработки и анализа экспериментальных дискретных значений E(n) исследовали динамические процессы естественной освещённости. Кроме того, отражены математические модели процесса изменения естественной освещённости E(t).
Ключевые слова: естественное освещение, динамический процесс, случайный сигнал, математическая статистика, случайная составляющая, сглаживающий фильтр, математическое ожидание.
Введение
В статье описаны подходы для решения задач, связанных с уменьшением влияния случайных факторов на результаты измерений.
Объекты и методы исследования Построение моделей принятия решений и исследование динамического процесса естественной освещённости E(t) проводилось на основе обработки и анализа экспериментальных дискретных значений E(n), полученных с помощью информационно-измерительного устройства в следующем составе: фоторезистивный датчик, модуль текущих интегральных телеизмерений (ТИТ), контроллер, персональный компьютер (ПК). Схема измерительного комплекса приведена на рис. 1.
© Семеняк М.В., Федоров В.К., 2012
Рис. 1. Схема измерительного комплекса естественной освещённости: АЦП - аналогово-цифровой преобразователь; ПУ - пункт управления
В качестве датчика был использован фоторезистор ФР764. Подключение производилось таким образом, чтобы ток фоторезисторного датчика был пропорционален уровню естественного освещения. Модуль ТИТ обеспечивает измерение унифицированных сигналов постоянного тока, поступающих на измерительные каналы ТИТ, и передачу обработанных значений измерений по межмодульному интерфейсу CAN. В качестве контроллера применялся модуль МП-03, предназначенный для использования в качестве центрального процессорного устройства для обеспечения в автоматическом режиме сбора, обработки принятой информации, её хранения и передачи результатов измерений в ПУ. Автоматизированное рабочее место, или пункт управления (ПУ), представляет собой персональный компьютер или сервер с установленным программным обеспечением верхнего уровня для связи с контроллером, сбора и сохранения измеренных значений в базе данных. Контроллер и модуль ТИТ являются законченными устройствами, обеспечивающими работу телемеханической автоматизированной системы диспетчерского управления. Таким образом, данные об измеренных значениях тока, пропорционального естественной освещённости, поступают на персональный компьютер (АРМ), где собираются в базу данных SQL. Дальнейшая статистическая обработка производилась используя полученные базы данных.
Результаты исследований График изменения естественной освещённости, построенный по измеренным значениям Е(п) с шагом дискретизации Т, равным пяти минутам, приведён на рис. 2.
у.е
О, DC!
0,034 0.029
_! 1 ПШ
---U— ллг Щfc
1 --------— -UP ^"nfliL —
- -J г.. 4Ti -
Jftff m /
1, мм.......
м Г ~f nl—-[ --
1 1 lit ■ ----
i ■ ч-
9:10 10:22 11:34 12:4fi 13:S3 1S10 16;22 17:34 10:^5
21:10
Рис. 2. График освещённости полного сигнала
Анализ построенной зависимости Е(^ позволяет выбрать в качестве математической модели процесса изменения естественной освещённости Е(^ аддитивную зависимость вида
e ^) = e0(t) + Е($) . (1)
В выражении (1) под E0(t) понимается трендовое изменение освещённости, e ^)- её случайная составляющая. Аналогичную модель можно записать и для дискретных значений
e^) = eo(n) + е(п) . (2)
Статистическая обработка случайных процессов сводится к получению трендовой составляющей E0 ^) по результатам наблюдаемых значений e ^), а затем вычислению случайной e ^) по формулам:
e(t) = e ^) - eo(t); (3)
e(t) = e (п) - eo(n). (4)
Задача восстановления или определения динамической функции тренда является оцениванием процессов, под которым в математической статистике принято называть обработку данных измерений (наблюдений) с целью уменьшения влияния случайных факторов и ошибок. В настоящее время для решения задач оценивания используются три подхода, а именно:
• применение методов оценивания математической статистики;
• приближение функций;
• фильтрация процессов.
В первой группе при выборе вероятностных критериев оптимизации приходят к методам максимальной апостериорной вероятности или максимального правдоподобия, но в вычислительной практике чаще используется критерий в виде функционалов невязок, или ошибок оценивания, которые приводят к методам наименьших квадратов (МНК), расстояний (МНР), произведений (МНП) и т.д. При обработке случайных процессов в режиме реального времени рекомендуется применять рекуррентные алгоритмы, в частности наиболее популярен рекуррентный метод наименьших квадратов (РМНК).
Во втором подходе оцениваются степени приближения исследуемого динамического
процесса к некоторой заданной или выбранной функции. Естественная освещённость e^), как это видно из графика (рис. 2), является периодической, а простейшая колебательная функция - это синусоида. Проверку возможности такого описания тренда E0 ^) можно осуществить по соответствующей эмпирической гистограмме теоретического закона распределения - арксинуса.
Для подхода, связанного с фильтрацией, разработано достаточно большое число алгоритмов оценивания как непрерывных, так и дискретных процессов, которые различаются типами моделей, применяемых для описания процессов и условий наблюдения, а также вероятностных характеристик случайных сигналов. Если оцениваемый процесс описывается линейным уравнением и наблюдения также линейны, то применяется достаточно известный линейный непрерывный или дискретный фильтр Калмана - Бьюси (ФКБ). Алгоритм состоит из дифференциальных, или разностных, уравнений для вектора восстанавливаемого процесса и ковариационной матрицы ошибок оценивания. В вычислительном отношении такой дискретный ФКБ совпадает с РМНК. Кроме описания процессов во временной области переменных состояния, в теории фильтрации используются и комплексные преобразования, то есть спектры сигналов, позволяющие проектировать и реализовывать достаточно простые аналоговые и цифровые фильтры.
Простейшей операцией, уменьшающей влияние ошибок, является усреднение или, с точки зрения математической статистики, вычисление математического ожидания. Примем следующие обозначения. Допустим, что исследуется непрерывный сигнал u(t) и известны
его значения и(п) , полученные из и(V), с шагом дискретизации Т, а результатом обработки или вычислений являются дискретные величины х(п). Оценку математического ожидания или среднее значение вычисляют по формуле
1 Ы
х(п) = и (n), (5)
П п=1
где N - объём обрабатываемой информации или количество измерений. Выражение (5) можно записать и в рекуррентном виде
х(п) = х(п-1)+-Ч(п); х(0) = о,п = ^ (б)
Операция усреднения, или простейшая процедура сглаживания, является частным случаем более общей цифровой фильтрации, описываемой разностным уравнением
Р ч
х(п) + ^ а(к) • х(п - к) = Ъ(0) • и(п) + ^ Ъ(к) • и(п - к) . (7)
к=1 к=1
Если коэффициенты а(к) при к = 1, р и Ь(к) при к = 1, ч - постоянные, то модель (7) и соответствующий цифровой фильтр линейны и можно записать передаточную функцию через коэффициенты в виде
7) Ъ(°) + 1 Ъ(к) В(2) _ к=1
-к
2
Н (2) = ^^ =-Р^- , (8)
А( 2) 1 + ]Га(к) • 2-к
к=1
или через импульсную характеристику
ад
н (2) = 2 И(к) • 2-к . (9)
к=0
В частном случае при а(к) = 0 и к = 1, р формулы (8) и (9) определяют выражения передаточной функции
Н (2) = В(2) = ]Г Ъ(к) • 2-к = ¿Л(к) • 2-к . (10)
к=0 к=0
Импульсная характеристика будет конечной и такие цифровые устройства называют КИХ-фильтрами, для них справедливо, что Ь(к) = Ь(к), к = 0, ч. В общем случае в формулу (9) входит бесконечное число значений Цк) и передаточная функция (8) описывает БИХ-фильтр.
Основной частотной характеристикой является амплитудная, для которой выполняется
Н(I)|2 = Н(2)Н(2-1) = 2ПТ . (11)
1 1 I 12
Частотные характеристики дискретных процессов и систем периодичны частоте дискретизации /0 = 1/ Т, поэтому, как правило, они рассматриваются либо в основной полосе
I е
/0 . 1
либо в полуполосе / е
0;^ 2
. Также вводится и нормированная частота
2 2
/ при Т = 1, для неё / е [0;0,5].
При проектировании сглаживающих фильтров вводится идеальная амплитудная час-
тотная характеристика
Нп (I)
, значения которой в полосе пропускания равны единице, а
остальные - нулю. Для М значений нормированной частоты / е [0;0,5] вводится квадратичный функционал
м ~ ~ if L(a, b) = £ Hn f ) - H f, a, b) ] (12)
i=0
и решается задача минимизации
(a, b) = arg min L(a, b) . (13)
Как правило, одним из численных методов оптимизации функции L(a, b) находят век-
I |Т I \T
торы a = a(1),a(2),...,a(p)\ и b = b(1),b(2),...,b(q)| .
Можно при проектировании использовать средства автоматизации, так как в состав пакета цифровой обработки Siqnal Processing Toolbox программной системы MATLAB входит процедура yulewalk, осуществляющая расчёт коэффициентов БИХ-фильтра по заданной амплитудной частотной характеристике [1].
Линейные разностные уравнения вида (7) достаточно широко используются и в математической статистике под названием моделей авторегрессии скользящего среднего, или АРСС (p, q), и авторегрессии АР(р).
Простейшими являются модели СС (q) и их коэффициенты b(k), определяющие веса, с
которыми учитываются значения входного сигнала u(n - к), к = 1, q в формировании выхода x(n). Если все коэффициенты одинаковы, то алгоритм СС аналогичен формуле (8) вычисления среднего значения, или оценке математического ожидания. Более естественным является выбор коэффициентов b(k) таким образом, чтобы ближайшие к x(n) значения u(n - к) учитывались с большим весом. В качестве такой убывающей функции может быть выбрана экспонента.
В работе предлагается использовать известный алгоритм экспоненциального сглаживания
x(n) = у x(n -1) + (1 -у) • u (n), (14)
который, как и рекуррентная процедура (8), является простейшим БИХ-фильтром первого порядка, или моделью АР(1).
При обработке экспериментальной информации коэффициент у выбирается из диапазона
0 <у< 1, (15)
но, если сравнивать (14) с формулой (5), то численное значение у следует ожидать более близким к единице.
Передаточная функция для уравнения (14) имеет вид
1 - Y
H о( z) = --^. (16)
1 - у • z
Дискретным моделям экспоненциального сглаживания (14) и (16) соответствует непрерывный аналог в виде простейшего фильтра нижних частот (ФНЧ) с передаточной и весовой функциями:
t
1 1 -Г W0(s) = --- и w(t) = — • е .
ТФ •s +1 ТФ
Характеристика w(t) изменяется по экспоненциальному закону, отсюда и название алгоритма. Можно показать, что параметры передаточных функций H0 (z) и W0 (s) связаны выражением
у = —Тф—. (17)
Тф + Т V '
Введём в рассмотрение относительные постоянную времени фильтра Тф = —ф и часто-
Т
Л Т
ту сопряжения Еф = -Е—, тогда формулу (17) можно переписать в виде
Т
7 =
/ф +1
(18)
для которой график зависимости параметра у от /ф для основной полуполосы нормированной частоты /ф е [0;0,5] приведён на рис. 3.
Рис. 3. Зависимость параметра у алгоритма экспоненциального сглаживания
В соответствие с рис. 3 диапазон (15) для коэффициента у изменится и будет определяться неравенством
0,67 <у< 1. (19)
На этапе предварительных исследований динамических процессов можно выбрать середину основной полуполосы /ф = 0,25 и соответствующее ему значение параметра экспоненциального сглаживания у = 0,8.
В работе предлагается измеренные (наблюдаемые) дискретные значения Е(п) процесса естественной освещённости Е(^) с целью определения тренда (постоянной составляющей) Е0(п) обрабатывать алгоритмом экспоненциального сглаживания (15), тогда
Е0(п) = у Е,(п -1) + (1 -у) - Е(п),
(20)
а для вычисления случайной составляющей использовать формулу
Е(п) = Е (п) - Е0(п) . (21)
Алгоритму цифровой обработки соответствует фильтр нижних частот Н0 (г) и W0 (*) , а для выражения (20) можно записать
~ 1 - 2-1 ~ Тф - S
Н(г) = 1 - Н 0 (г) = ^--1 и Щ*) = 1 - Wo(s) = у ф
1 -у - г-
Тф - * +1
то есть уравнению вычисления случайной составляющей Е(п) соответствует простейший фильтр верхних частот, поэтому алгоритмы (17) и (18) можно использовать для цифровой обработки динамического процесса естественной освещённости.
1
Вывод
В работе предлагается задачу исследования динамических процессов естественной освещённости E (t), E0 (t) и E (t) проводить в два этапа. На первом определить моментные характеристики (математическое ожидание, дисперсию, ассиметрию и эксцесс) и эмпирические гистограммы с целью проверки соответствия их теоретическим законам распределения, в первую очередь для тренда E0 (t) . На втором этапе оценить частотные характеристики и провести спектральный анализ для подтверждения достоверности предлагаемых в работе алгоритмов определения тренда E0 (t) и случайной составляющей E (t) процесса естественной освещённости.
Список литературы
1. Оптенгейм, А. Цифровая обработка сигналов / А. Оптенгейм. - М. : Техносфера, 2007. - 856 с.
SUMMARY
M. V. Semenjak, V.K. Fyodorov
The choice of measuring devices of natural illumination, its mathematical description
and the methods of exploration
In given article, on the basis of processing and the analysis of experimental discrete values E (n) investigated dynamic processes of natural light exposure. Besides, here mathematical models of process of change of natural light exposure E (t) are reflected. In article approaches for the decision of the problems connected with reduction of influence of random factors on results of measurements are described.
Key words: daylight lightning, dynamic process, random signal, mathematical statistics, random component, ripple filter, mathematical expectation.
