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

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

Поскольку внешние условия, в которых работает комбайн, существенно влияют на установочные параметры, то к информации о внешних факторах предъявляются высокие требования в части достоверности и адекватности. Информация поступает в виде оценок экспертов, а также в результате обработки данных датчиков и анализа изображений, что говорит о ее нечетком характере. Взаимосвязи регулируемых параметров комбайна с внешними условиями еще более очевидно можно отнести к нечетким в силу неконкретного описания внешних условий, а также неоднозначности и наличия сложных корреляционных зависимостей. Для описания функционирования подобных систем традиционные методы математического моделирования, такие как построение регрессионных моделей и другие статистические методы, опирающиеся на экспериментальные данные и цифровую обработку [2–6], практически непригодны, поскольку применимы только в рассматриваемом диапазоне величин, а оперативное перестроение модели, необходимое в полевых условиях, невозможно без анализа новых экспериментальных данных. Помимо этого, наличие множества неоднозначных связей между внешними условиями, регулируемыми параметрами и показателями качества уборки приводит к трудно оптимизируемым, громоздким математическим конструкциям в рамках множественного регрессионного анализа и делает само применение этого подхода неэффективным.

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

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

Обзор литературы

Для описания процессов принятия решений по управлению технологическими процессами в сложных системах используется математический аппарат нечеткой логики. В последние годы он достаточно широко применяется для решения задач агропромышленного комплекса, среди которых идентификация сорняков, оценка урожайности, оценка качества почвы, повышение продуктивности культур и эффективности использования удобрений и др. [7–10]. Обосновано применение данного подхода к оценке значений регулируемых параметров комбайна [11–13] и к проблеме оперативной корректировки технологических настроек в случае обнаружения нарушений качества уборки [14].

Задача предварительной настройки относится к классу задач принятия решений в нечетких условиях, поскольку характеристики внешних факторов имеют нечеткие размытые границы, неоднозначно влияние различных комбинаций параметров настройки на показатели качества уборки, а устанавливаемые априорно взаимосвязи могут оказаться неадекватными для реальной ситуации. Все эти факты послужили основанием для выбора способа решения данной задачи – создания интеллектуальной системы, использующей для получения решения нечеткий логический вывод, основанный на информации экспертного, оценочного и эвристического характера [11]. Источниками такой информации являются квалифицированные специалисты в данной области, накопленные знания о взаимосвязях показателей качества уборки, внешние условия и параметры настройки комбайна. Такой подход был использован для вывода важнейших регулируемых параметров комбайна: скорости движения, частоты вращения молотильного барабана, частоты вращения вентилятора очистки [11; 15].

Наряду с перечисленными параметрами для качественной уборки урожая существенное значение имеет правильный выбор параметров жатвенной части комбайна, таких как вынос мотовила по вертикали и горизонтали и другие. Предварительная настройка регулируемых параметров комбайна, соответствующих условиям уборки, охватывает весь комплекс параметров, и настоящее исследование посвящено разработке метода выбора оптимальных значений параметров жатвенной части. В проведенных ранее исследованиях представлено лингвистическое описание внешних факторов, рассмотрены урожайность, влажность, засоренность и соломистость зерна [11; 14; 15]. На основе алгоритма были установлены оптимальные модели лингвистических переменных, соответствующих перечисленным факторам [16]. Для настройки параметров жатвенной части являются актуальными и другие параметры, такие как высота, полеглость хлебостоя и другие. Их лингвистическое описание представлено в другой работе [17].

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

Материалы и методы

Предлагаемый подход к формированию метода предварительной настройки регулируемых параметров жатвенной части базируется на математическом аппарате нечеткой логики и в общем виде состоит из этапов фаззификации, композиции и дефаззификации [11; 15]. Далее остановимся подробно на каждом из них.

На этапе фаззификации формулируется лингвистическое описание условий задачи [18], определяется множество { X } входных переменных задачи в виде совокупности значений факторов внешней среды, существенно влияющих на величину выходной переменной Y – регулируемого параметра машины. В нашем случае при наличии нескольких регулируемых параметров выводы их значений рассматриваются как параллельные процессы.

Модели входных и выходных признаков X , Y представляются в виде семантических пространств и соответствующих им функций принадлежности (ФП) [18].

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

При построении ФП используются нормальные нечеткие множества с верхней границей, равной 1 ( sup ц_А ( x ) = 1 ), причем как унимодаль-

X € E ные (например, треугольной или экспоненциальной формы), так и имеющие область толерантности (например, трапециевидные). Среди процедур построения ФП выберем представление с помощью стандартных функций треугольной и трапециевидной форм, оценка параметров которых не вызывает затруднений у экспертов и отвечает специфике неопределенности.

Согласованность экспертной информации оценивается следующими характеристиками [16; 18]:

– матрицей K парной согласованности моделей экспертов;

– матрицей D нечеткости;

– аддитивным k и мультипликативным k показателями.

Первые две характеристики определят парную согласованность моделей, а последние – согласованность моделей в целом.

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

Том 32, № 4. 2022

экспертной информации реальным условиям. Согласованность экспертной информации является одним из критериев ее качества. Низкая согласованность может возникать из˗за недостаточной квалификации кого˗то из экспертов. Для установления этого факта анализируют матрицы нечеткости D и парной согласованности K моделей и выполняют корректировку согласно алгоритму [16]. Общая рассогласованность требует более углубленного анализа модели, положенной в основу формального описания.

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

i : A ^ B ;( F ) , где i – имя нечеткой продукции; A ⇒ B – ядро продукции, в котором А – антецендент, включающий элементарные предложения, соединенные логическими связками «и», а В – консеквент, выражающий заключение; ⇒ – оператор логического следования; F – коэффициент достоверности продукции². В словесной форме ядро продукции имеет вид «ЕСЛИ А , ТО В ». Для проверки адекватности сформулированных правил строят поверхности, отражающие зависимости выходной ЛП от различных комбинаций входных ЛП, взятых попарно.

В общем случае развернутая форма нечеткого логического вывода для системы продукций имеет вид3

M b -=\ ⁽ M a • ( ^x ) a M r ( ^x, y )). (1) x e X

На этапе дефаззификации для конкретной комбинации внешних факторов выводятся методом центра тяжести значения выходных параметров4. Вывод реализуется в пакете Fuzzy Logic Toolbox (MatLab).

Результаты исследования

На основе систематизации описанных выше этапов и их приложения к исследованию предметной области «предварительная настройка жатвенной части комбайна» предложен метод решения задачи выбора оптимальных начальных значений параметров рабочих органов агрегата. Иерархия составляющих метода включает несколько уровней и приведена на рисунке 1. Далее рассмотрим подробно отдельные элементы данной структурной схемы.

В результате использования метода концептуализации выполняется определение входных и выходных факторов, их отношений и типов данных. Рассматриваются внешние факторы, оказывающие влияние на регулируемые параметры жатвенной части. Среди них наиболее существенными являются урожайность, засоренность, влажность, высота, полеглость, густота, спутанность хлебостоя. При определении лингвистической шкалы ЛП, соответствующих перечисленным факторам, необходимо учитывать тип убираемой культуры. Может потребоваться и дифференциация такого фактора, как урожайность. Так, для пшеницы рассматривают урожайность примерно 60, 50, 40 ц/га и т. д.

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

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

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

Приведем лингвистическое описание нескольких внешних факторов и регулируемых параметров для одной из культур – пшеницы (табл. 1).

Выбор границ термов проиллюстрирован на рисунке 3 для ЛП «положение мотовила по высоте».

Указано количество термов ЛП, выбранных на основе критериев согласованности при обработке информации четырех экспертов [16]. Характеристики согласованности вычислены с помощью программного комплекса [19] и представлены в таблице 2 для нескольких ЛП.

На рисунке 4 приведено несколько графиков ФП рассматриваемых ЛП.

На основе частных ФП построена обобщенная ФП для всех ЛП, которые используются при разработке базы продукционных правил.

2-

Т а б л и ц а 1

T a b l e 1

Описание лингвистических переменных Tuples of the linguistic variables

Кортежи ЛП / Linguistic variable tuples

Условное обозначение / Unit designation

<Урожайность 40 ц/га {Менее 40, Примерно 40, Более 40}, 34–46> /

<Влажность хлебостоя,% {Сухой; Нормальный, Влажный}, 0–20> /

<Соломистость, % {Малая, Нормальная}, 40–70> /

<Засоренность хлебостоя, % {Низкая, Большая}, 0–40> /

<Высота хлебостоя, см {Низкорослый, Ниже среднего, Средний, Высокий}, 30–120> /

<Полеглость хлебостоя, % {Прямостоящий, Низкая, Средняя, Высокая}, 0–40> / < Lodging of corn crops, % {Right, Low, Medium, High}, 0–40>

<Спутанность хлебостоя, % {Нормальный, Спутанный}, 0–40> / < Confusion of corn crops, % {Normal, Confused}, 0–40>

<Густота хлебостоя, стеб./м2 {Разреженный, Ниже среднего, Средний, Густой}, 100–1 000> / < Density of corn crops, stems/m2 {Fragile, Below average, Medium, Dense}, 100–1,000>

<Положение мотовила по высоте, см {Очень низкое, Низкое, Среднее, Высокое}, -10–110> / < Reel height position, cm {Very low, Low, Medium, High}, -10–110>

<Положение мотовила по горизонтали, см {Незначительный, Ниже среднего, Средний, Дальний}, 20–70> / < Reel horizontal position, cm {Significant, Below average, Medium, Far}, 20–70>

УР ={УРМ40, УРП40, УРБ40}

ВС ={СХ, НОРХ, ВЛХ}

СОЛ ={МСОЛ, НСОЛ}

ЗХ ={НЗХ, БЗХ}

ВысХ = {НХ, НСрХ, СрХ, ВХ}

ПХ = {ППХ, НПХ, СПХ, ВПХ}

СХ = {НСХ, ССХ}

ГХ = {ГРХ, НСрХ, СГХ, ГГХ}

МПВ = {ОН, Н, С, В}

МПГ = {Н, НС, С, Д}

Положение мотовила по высоте / Height position of reel: -10–35 очень низкое / -10–35 very low 10–60 низкое / 10–60 low 40–95 среднее / 40–95 average 70–110 высокое / 70–110 high

-10 0 10 20 30 40 50 60  70 80 90 100 110 см / cm

Р и с. 3. Границы термов ЛП «положение мотовила по высоте»

F i g. 3. Limits of the terms of the linguistic variable “Reel vertical position”

Т а б л и ц а 2

T a b l e 2

Матрица K парной согласованности и показатели к и k общей согласованности Matrix K of pair consistency and Indices k and k of overall consistency

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

База знаний сформулирована для каждой культуры и регулируемого параметра, коэффициент достоверности правил равен 1. Например, для пшеницы с урожайностью примерно 40 ц/га Technologies, machinery and equipment

и выходной ЛП «положение мотовила по высоте» база знаний содержит 128 правил вида:

ЕСЛИ <высота хлебостоя «низкорослый», и полеглость хлебостоя «прямостоящий», и спутанность хлебостоя «нормальный», и густота хлебостоя «разреженный», ТО положение мотовила по высоте «очень низкое»>.

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

ЛП / Linguistic Variable		Матрица K / Matrix K			k	k
«Соломистость хлебостоя» /	1,000	0,667	0,667	0,964	0,699	0,677
“Strawness of corn crops”	0,667	1,000	0,946	0,667
	0,667	0,946	1,000	0,667
	0,964	0,667	0,667	1,000
«Высота хлебостоя» / “Height of	1,000	0,945	0,901	0,812	0,812	0,812
corn crops”	0,945	1,000	0,951	0,862
	0,901	0,951	1,000	0,911
	0,812	0,862	0,911	1,000
«Полеглость хлебостоя» /	1,000	0,868	0,856	0,929	0,820	0,812
“Lodging of corn crops”	0,868	1,000	0,912	0,839
	0,856	0,912	1,000	0,927
	0,929	0,839	0,927	1,000
«Густота хлебостоя» / “Density	1,000	0,958	0,896	0,872	0,787	0,779
of corn crops”	0,958	1,000	0,857	0,832
	0,896	0,857	1,000	0,867
	0,872	0,832	0,867	1,000
«Положение мотовила по	1,000	0,868	0,868	0,874	0,850	0,840
высоте» / “Reel height positionˮ	0,868	1,000	1,000	0,747
	0,868	1,000	1,000	0,747
	0,874	0,747	0,747	1,000
«Положение мотовила по	1,000	0,860	0,860	0,892	0,880	0,870
горизонтали» / “Reel horizontal	0,860	1,000	0,960	0,810
position”	0,860	0,960	1,000	0,920
	0,892	0,810	0,920	1,000

Ц ⁽ х )

1,0

0,5

0,0

30   35   40   45   50   55   60   65   70   75   80   85   90   95  100 105 110 115 120

x

Низкорослый / Low-growth Ниже среднего / Below average Средний / Average Высокий / High

a)

b)

Ц ( х )   1,0

0,5

0,0

100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1 000

x

Разреженный / Sparse Ниже среднего / Below average Средний / Average Густой / Dense ^с)

d)

e)

Р и с. 4. Графики функции принадлежности термов ЛП: a) «Высота хлебостоя»; b) «Полеглость хлебостоя»; c) «Густота хлебостоя»; d) «Положение мотовила по высоте»;

• e) «Положение мотовила по горизонтали»

F i g. 4. Graphs of membership functions of the terms of the linguistic variables: a) “Height of corn crops”; b) “Lodging of corn crops”; c) “Density of corn crops”; d) “Reel vertical position”; e) “Reel horizontal position”

Положение мотовила

УР ВлХ   СОЛ ЗХ   ВысХ   ПХ    СХ    ГХ

по высоте /

Reel height position

Р и с. 5. Система взаимосвязей входных ЛП и выходной ЛП «положение мотовила по высоте».

Условные обозначения, используемые на рисунке, расшифрованы в таблице 1

F i g. 5. The system of interrelationships of input linguistic variables and output linguistic variable “Reel vertical position”. The symbols used in the figure are shown in Table 1

Обсуждение и заключение

Разработан метод выбора начальных параметров жатвенной части комбайна. В рамках предложенного подхода подробно рассмотрены элементы формально˗логической схемы вывода решений. Даны лингвистические описания значимых внешних факторов и параметров настройки, выбраны оптимальные модели на основе показателей согласованности, создана база нечетких экспертных знаний. Полученные результаты, совместно с аналогичными разработками по настройке параметров

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

Поступила 25.08.2022; одобрена после рецензирования 14.11.2022; принята к публикации 21.11.2022

Об авторах:

Все авторы прочитали и одобрили окончательный вариант рукописи.

Submitted 25.08.2022; approved after reviewing 14.11.2022; accepted for publication 21.11.2022

All authors have read and approved the final manuscript.