1. ИСПОЛЬЗОВАНИЕ ЧИСЛЕННЫХ МОДЕЛЕЙ В ЗАДАЧАХ НАНОТЕХНОЛОГИЙ
1.1. О ПОСТРОЕНИИ АНАЛИТИЧЕСКОГО ЗАДАНИЯ
W W
k-ЗНАЧНОЙ ПОРОГОВОЙ ФУНКЦИИ
Бурделёв Александр Владимирович, ст. преподаватель кафедры Математического моделирования и анализа данных факультета Прикладной математики и информатики. Белорусский государственный университет. E-mail: [email protected]
Никонов Владимир Глебович, доктор технических наук, доцент, член президиума РАЕН
Аннотация: Задача: В работе [2] предложен ряд подходов к решению задачи нахождения коэффициентов линейной формы пороговой булевой функции. Эти подходы предполагают использование характеристического вектора булевой функции в качестве первого приближения коэффициентов линейной формы и окончательное их уточнение с помощью нескольких итеративных алгоритмов. В данной работе рассмотрен вопрос построения порогового представления k-значной пороговой функции.
Модель: Для решения вопроса построения порогового представления k-значной пороговой функции предлагается несколько трактовок близости двух k-значных функций: мультипликативные, разностные и квадратичные коэффициенты, коэффициенты роста и коэффициенты возрастания. Рассматривается возможность аппроксимации коэффициентов линейной формы данными коэффициентами и возможность дальнейшей коррекции.
Выводы: На основании примеров сделано заключение о том, что для первого приближения коэффициентов линейной формы предпочтение стоит отдать использованию коэффициентов возрастания. При этом аналогично булевому случаю подтверждается предположение о необходимости введения итеративной процедуры. Предложен итеративный алгоритм нахождения коэффициентов линейной формы k-значной пороговой функции на основе коэффициентов возрастания
Ключевые слова: пороговая k-значная функция, нахождение линейной формы пороговой функции, пороговая логика
ABOUT CONSTRUCTION OF ANALYTICAL DEFINITION OF k-VALUE THRESHOLD FUNCTION
Burdeliov Alexander V., Senior Lecturer at the Department of Mathematical Modeling and Data Analysis at the Faculty of Applied Mathematics and Computer Science, Belarusian State University. E-mail: aburd2011@mail. ru
Nikonov Vladimir G., Doctor of Sciences, Docent, member of presidium of Russian Academy of Natural Sciences.
Abstract: Task: In [2] proposed a few methods for finding coefficients of linear form of Boolean threshold function. These methods are founding on using coefficients of characteristic vector, as first approximation of coefficients of linear form, and then a few algorithms for correction of this approximation. In this paper submitted for consideration the question of finding coefficients of linear form of k-value threshold function.
Model: In this paper submitted a few interpretations of closeness of two k-value functions by definition of multiplication, differential and quadratic coefficients, also expansion coefficients and increase coefficients. Considered potential of these coefficients to approximate the coefficients of linear form and possibility of further correction.
Findings: In this paper made the conclusion that expansion coefficients and increase coefficients are better for approximation the coefficients of linear form. Submitted algorithm for finding coefficients of linear form of k-value threshold function funding on increase coefficients.
Index terms: k-value threshold function, finding coefficients of linear form of threshold function, threshold logic.
1. ВВЕДЕНИЕ
Одно из актуальных направлений развития современных вычислительных технологий связано с использованием бионических принципов моделирования нейронных сетей живых организмов, что привело к построению так называемых нейрокомпьютеров. Базовыми элементами нейрокомпьютера являются формальные нейроны, функционирование которых описывается пороговыми функциями, - булевыми, или как показали последующие исследования к-значными [6]. В представленной статье рассматриваются некоторые подходы к решению задачи построения порогового задания к-значной функции.
Переход от операций в булевой области к преобразованиям в к-значной логике соответствует современным тенденциям развития систем связи, а именно увеличению быстродействия и возрастанию объёмов перерабатываемой информации. Поэтому изучение свойств и характеристик различных классов к-значных функций представляется актуальным направлением прикладной дискретной математики. Особый интерес приобретает изучение пороговых к-значных функций, обладающих простотой технической реализации как в существующей, так и в перспективной элементной базе, в частности, в оптической вычислительной среде. В области синтеза пороговые соотношения дают возможность компактно реализовать в арифметической вычислительной среде достаточно сложные в традиционном понимании к-значные функции большого числа переменных.
Определение 1. Булева функция f(x1,■..,xп), для которой существует линейная форма
¿(Х], ■■■, Х~П) = а]Х] + ^2X2 + "* + ЯпХп
с вещественными коэффициентами и вещественный порог Т такая, что
/<Х.....хп) =
1,если ^ о , л:, > Г
! = 1
И
^ < Т
называется пороговой булевой функцией. Вектор (а1,.,ап) называется вектором структуры.
При инвертировании переменных х1,...,хп или их перестановке свойство принадлежности функции классу пороговых сохранится. Такое преобразование для булевых функций называется преобразованием однотипности.
Определение 2. Функция к-значной логики fk(x1,■..,xп), для которой существует линейная форма
¿(Х],..., Хп) = а]Х] + а2Х2 + "* + апХп
с вещественными коэффициентами и набор вещественных порогов Ь0, Ь1,..., Ьк такие, что для всех I е 0, к — 1 выполняется импликация
называется пороговой к-значной функцией.
При переходе в к-значную область понятие однотипности сохраняется с заменой булевой операции инвертирования на к-значное преобразование с помощью стрелки Лукашевича: при преобразовании переменных Х1, . , Хп с помощью стрелки Лукашевича или их перестановке свойство принадлежности функции классу пороговых сохранится.
2. ПРОБЛЕМА ХАРАКТЕРИЗАЦИИ ПОРОГОВОЙ ФУНКЦИИ В БУЛЕВОМ СЛУЧАЕ
К числу важных задач пороговой логики является распознавание принадлежности функции классу пороговых функций и, в случае положительного решения, - нахождение коэффициентов линейной формы и порога. Известно, что, несмотря на простоту формулировки, обе эти задачи оказалась достаточно сложными.
В работе [2] предложен ряд подходов к решению задачи нахождения коэффициентов линейной формы и порога в булевом случае. Данные подходы предполагают нахождение характеристического вектора булевой функции f(x), представленной в базисе {-1,1}.
0, если
Определение 3. Пусть функция fix), представлена в базисе {-1,1}. Упорядоченный набор весов с0, с1,..., сп, где
= ^ xjix), i — 1,п
XEVf
с0= ^ f{x),
XEV{l
называется характеристическим вектором с функции fix).
Определение 4. Будем говорить, что линейная форма Lixi,..., x~n) — x1 + + Q^x^ + * * * + Unxn дает чистое разделение областей значений функции fk(£1,... ,£п), если для любого i — 0, к — 2 выполняется строгое неравенство
max {Qjfi + ■
• + вп«п} < min {ajfj +
fk('i.....«n)=i+l
•+ an£nl
В случае если построенная с помощью коэффициентов линейная форма дает чистое разделение областей значений функции /к(£1,..., £п), границы Ь0,Ь1,...,Ьк можно определить, например, следующим способом:
£>, - min {а1£1 + ■■■ + апЕ,J, i = 0,к - 1:
fЧ<t
Ъь =
та х
{fl + ■■■ + fl„f1 J + 1.
В работе [2] коэффициенты характеристического вектора используются в качестве модели первого приближения для коэффициентов а1,...,ап искомого порогового задания
а1х1 + а2х2 + —+ апхп > Ь булевой функции /(х1,.,хп) как а^ « с. Но коэффициенты характеристического вектора не всегда точно аппроксимируют коэффициенты искомого порогового задания. Покажем это на двух примерах.
Для удобства изложения рассмотрим коэффициенты, характеризующие число знако-перемен по направлению соответствующей переменной х^, которые эквивалентны коэффициентам характеристического вектора, но в отличие от них не требуют перехода в
базис {-1,1}. Число знакоперемен для пороговой функции в точности равно числу ребер, на концах которых функция принимает разные значения.
Определение 5. Числом знакоперемен по переменной х^, (7 — 1,п) функции f(х1,... ,хп) называется
(*1.....*!-
I
/(.V1,..., Л'г-i, 1 A'i + i ■••*„)
fix i,
Легко показать, что коэффициенты с1, ...,сп характеристического вектора функции в базисе {-1,1} и коэффициенты а1,.,ап, характеризующие число знакоперемен в базисе
{0,1} эквивалентны и имеют соотношение
*
Cj — —— 1,п.
1 2
Рассмотрим несколько примеров использования числа знакоперемен для нахождения коэффициентов разделяющей плоскости.
Пример 1. Рассмотрим функцию
заданную геометрически на Рис. 1.
f> Xi
- /(*!, x2l = О • - f{xt, хг, х2) = 1
Рисунок 1. Функция f(х1,х2,х3).
Для данной функции * — 1, *2 — 1, — 3.
Подставив а1,а2,аъ в качестве коэффициентов разделяющей плоскости, получаем линейную форму Ь(х1,х2,х3) — х1 + х2 + 3х3, для которой выполняется неравенство
/(ж
max, {.Vj + ,v2 + Зл'з} < , min , {.Vj + v2 + 3.v3} = 4.
1.хг-Хз)~0 f(x1.xz,x2)=l
Таким образом, функция f(x1,x2,xз) реализуется пороговой функцией и в качестве коэффициентов а1,а2,аз искомого порогового задания можно использовать числа знакоперемен а1, *2, *з.
Следующий пример показывает, что число знакоперемен не всегда точно аппроксимирует коэффициенты искомого порогового задания.
Пример 2. Рассмотрим функцию
f(x1,x2,x3lx4) =
ест и 3.Vj — 3.i'2 + 5х3 — 3.v4 > 2
заданную геометрически на Рис. 2.
Рисунок 2. Функция f(x1,x2,xз,x4). Для данной функции
*1 = 2, *2 = —2, *з = 6, = —2.
Подставив числа знакоперемен в качестве коэффициентов разделяющей плоскости, получаем линейную форму ¿(Х1,Х2,Хз,Х.) = 2Х1 — 2Х2 + 6Хз — 2Х4, которая не разделяет точки (1,0,0,0) и (0,1,1,1):
¿(1,0,0,0) = ¿(0,1,1,1) = 2. Из Примера 2 видно, что число знакопере-мен можно использовать в качестве модели первого приближения для коэффициентов а1, а2,..., ап, но для получения разделяющей
плоскости необходимо провести некоторую процедуру корректирования.
Кроме характеристического вектора и числа знакоперемен для характеризации пороговой булевой функции используются коэффициенты Чоу. Коэффициенты Чоу вычисляются по тем же формулам, что и коэффициенты характеристического вектора (Определение 3), только для функции, представленной в базисе {0,1} и имеют с ними следующие соотношения:
с0 = 20 — 2
( mi\
i
1,п,
где 0, I = 0,# - коэффициенты Чоу. Коэффициенты характеристического вектора и коэффициенты Чоу в одинаковой степени содержат информацию о принадлежности булевой функции классу пороговых функций. Чоу принадлежат две соответствующие теоремы.
Первая теорема Чоу. Если коэффициенты Чоу у двух функций равны, то либо они обе принадлежат классу пороговых функций, либо обе ему не принадлежат.
Вторая теорема Чоу. Если коэффициенты Чоу у двух функций равны и одна из них принадлежит классу пороговых функций, то эти функции равны.
Коэффициенты Чоу, как и коэффициенты характеристического вектора и число знако-перемен, могут быть использованы для первого приближения в поиске коэффициентов а1, а2,..., ап искомого порогового задания а1х1 + а2х2 + —+ апхп > Ь булевой функции f(x1,... ,хп). В работе [2] изложено несколько итеративных алгоритмов вычисления коэффициентов а1,а2,... ,ап с помощью коэффициентов Ь1,.,Ьп характеристического вектора функции: метод минимизации функционала [2, глава 3] и итеративный метод синтеза порогового элемента [2, глава 4]. Таким образом, задача нахождения параметров пороговой булевой функции является объективно более сложной математической задачей, чем задача определение принадлежности функции классу по-
роговых функций, и сводится к итеративной процедуре.
3. НАХОЖДЕНИЕ АНАЛИТИЧЕСКОГО
ПРЕДСТАВЛЕНИЯ к-ЗНАЧНЫХ ПОРОГОВЫХ ФУНКЦИЙ
Рассмотрим вопрос построения порогового представления к-значной пороговой функции. Число знакоперемен в булевом случае можно трактовать как меру близости функций fk(x1,... ,хп) и хь1 е 1,п. Переход в к-значную область ставит первый важный вопрос о возможных вариантах трактовки близости или отличия двух к-значных функций, а в случае нахождения аналитического представления к-значных функций - её близости к функциям х^. Можно предложить несколько таких трактовок.
Определение 6. Мультипликативным коэффициентом переменной х^ функции fk(x1, ...,хп) называется
5Хб = Х1~ f (х1,..-,хп).
(х1,.,хп)
Определение 7. Разностным коэффициентом переменной х^ функции fk(x1,... ,хп) называется
9х6 = ^ к— ^О^.^хЛ
(х1,-,хп)
Определение 8. Квадратичным коэффициентом переменной х^ функции fk(x1,... ,хп) называется
8х1= ^ {х^— ^(Х1,...,Хп)) .
(х1,...,хп)
Определение 9. Коэффициентом роста по переменной х^ функции fk(x1,... ,хп) называется
Определение 10. Коэффициентом возрастания по переменной х^ функции fk(x1, ...,хп) называется
1С—2 К—1
X XX .....+ .....*»>-
> Х- . Д'т-'Л'-.- 1—0 1 + 1
Изучение этих коэффициентов позволяют заключить, что все они характеризуют меру близости функции fk к х^ и поэтому они введены и будут рассматриваться далее.
Для некоторых пороговых к-значных функций коэффициенты роста и возрастания приводят к прямому нахождению коэффициентов а1,..., ап. Рассмотрим пример.
Пример 3. Рассмотрим функцию fз(x1,x2,xз): Жз хЖз хЖз ^ Ъз, задающуюся линейной формой ¿(х1,х2,хз) = 3х1 + 2х2 + хз и системой порогов следующим образом:
- р{х±,х2,х-{] - О • - р(х1,х2,х-{) = 1 □ ~/3(х1,х2,х3) — 2
Рисунок 3. Функция fз(x1,x2,xз). Вычислим для данной функции введённые выше коэффициенты:
1. Мультипликативные коэффициенты равны: 5х1=41,5х2 = 35,5х3 = 31.
2. Разностные коэффициенты равны:
9х8 = 8, 9х8 = 16, 9х3 = 20.
3. Квадратичные коэффициенты равны: 5х1 = 8, ¿х2 = 20, 8х3 = 28.
4. Коэффициенты роста равны: АХ1 — 14, АХ2 — 8, — 4.
5. Коэффициенты возрастания равны:
ЛХ1 — 28> ЛХ2 — 16, ЛХз — 8
Из представленных 5 наборов коэффициентов лишь два набора дают разделяющую плоскость - это коэффициенты роста и коэффициенты возрастания.
Действительно, подставив коэффициенты роста АХ1 — 14, АХ2 — 8, АХз —4 в линейную форму 1У(х1,х2) — 14х1 + 8х2 + 4х3 получаем чистое разделение областей значений функции f3(х1 , х2, х, ). Вычислив границы Ь0 — 0, Ь1 — 20, Ь2 — 34, Ь3 — 53, приведенным выше способом, получим реализацию функции f3(х1,х2,х3):
/Ч*1,*2.*з) = 0 « 0 < ¿Хх1,х2,х3) < 20 /3Сч,х2,х3) =1 ~ 20 < Г(х1гх2,х3) < 34 /3(х1,х21х3) = 2 => 34 < ЬХх1,х2,х3) < 53.
Подставив коэффициенты возрастания ЛХ1 — 28, ЛХ2 — 16, ЛХз — 8, в линейную форму Ь^^(х1,х2) — 28х1 + 16х2 + 8х3 получаем чистое разделение областей значений функции f3(х1 , х2, хз ). Вычислив границы Ь0 — 0, Ь1 — 40, Ь2 — 68, Ь3 — 103, приведенным выше способом, получаем реализацию функции f3(х1,х2,х3):
г3(хих2,х3) = 0 => 0 < ¿"(*1,*2.*з) < 40
/3(*1,л:2,.*3) = 1 <=» 40 < ¡Т&ъхг.х3) < 68
^(х1,х2,х3) = 2 => 68 < Г(Х1,Х2,Х3) < 103
В то же время мультипликативные, разностные и квадратичные коэффициенты не дают построение разделяющей плоскости. В дальнейшем будем рассматривать только коэффициенты роста и возрастания.
4. ИСПОЛЬЗОВАНИЕ КОЭФФИЦИЕНТОВ РОСТА И ВОЗРАСТАНИЯ ДЛЯ ПОСТРОЕНИЯ ПОРОГОВОГО ПРЕДСТАВЛЕНИЯ ^ЗНАЧНОЙ ПОРОГОВОЙ ФУНКЦИИ
Использование коэффициентов роста не всегда приводит к непосредственному нахождению коэффициентов а1,..., ап.
Пример 4. Рассмотрим функцию ^.(х1,х2):Ж4 хЖ. ^ Ж., задающуюся линейной формой Ь(х1,х2) — 2х1 + 5х2 и системой порогов следующим образом:
3 п-2 »
1
-р-
-р
-f*(:xllx2lx3) = 0
• -f*{x1,x2,x3) = 1 □ -f*(xL,x2,x3) = 2
♦ -f*(Xi,x2,x3) = 3
Рисунок 4. Функция f4(х1,х2). Вычислим для данной функции коэффициенты роста и возрастания: 1. Коэффициенты роста равны: Ах1 — 5, АХ2 — 10.
Х2
2. Коэффициенты возрастания равны:
Лх1 — 15, ЛХ2 — 33.
Полученная с помощью коэффициентов роста АХ1 — 5, АХ2 — 10 линейная форма Ь^(х1,х2) — 5х1 + 10х2 не дает чистое разделение областей значений функции
f4(x1, х2у.
шах {5л-! + 10х2} = Г (3,2) = 35 =
Подставив коэффициенты возрастания ЛХ1 — 15 и ЛХ1 — 33 в линейную форму Ъ^(х1,х2) — 15х1 + 33х2 получаем чистое разделение областей значений функции f3(х1,х2). Вычислив границы
Ь0 — 0, Ь1 — 45, Ь2 — 78, Ь3 — 114, Ь. — 145, приведенным выше способом, получаем реализацию функции f3(х1,х2):
Таким образом, линейная форма, построенная с помощью коэффициентов роста, не всегда дает чистое разделение областей значений функции. На основании примера можно заключить, что при нахождении линейной формы предпочтение стоит отдать использованию коэффициентов возрастания, которые дают более точное приближение коэффициентов линейной формы.
Следующий пример показывает, что и коэффициенты возрастания не всегда дают чистое разделение областей значений функции. Такое замечание следует признать ожидаемым т.к. в более простом - булевом - случае подобная задача решается итеративно.
Пример 5. Рассмотрим функцию f.(x1,x2):Ж.xЖ. ^ Ж., задающуюся линейной формой ¿(х1, х2) = 5х1 + х2 и системой порогов следующим образом:
/4(х1,х2) = О 0 < их1,х2) < 15,5 /4(Х1,Х2) = 1 о 15,5 < их1,х2) < 16,5 /4(*1,*2) = 2 <=> 16,5 < 1(х1,х2) < 17,5 ГЧх1,х2) = 3 = 17,5 < Ь(х1,х2)-
-л ч>
О 1
Рисунок 5. Функция f.(x1,x2).
Вычислим для данной функции коэффициенты роста и возрастания:
1. Коэффициенты роста равны:
Д*8 = 6, Д*2 = 3.
2. Коэффициенты возрастания равны: DXl = 18, DXz = 10.
Полученная с помощью коэффициентов ро-
ста
Axi = 6,
ДХ2 = 3
-/4(:х1,х2,х3) = 0
• -f*iXi,x 2lx3) = 1 □ -f*(_xL,x2,x3) = 2
♦ -f*iXi,x 2,.y3) = 3
линейная форма Lv(*i>*2) = 6% + 3x2 не дает чистое разделение областей значений функции
/4(х1; х2):
min {бл-i + 3.v2} = Г(ЗД) = 21 =
/J0riji2)=l
Полученная с помощью коэффициентов возрастания DXl = 18 и DXl = 10 линейная форма Lvv(xx,х2) = 18% + 10х2 также не дает чистое разделение областей значений функции /4(хх, х2):
min {18.fi + 10х2} = Г(ЗД) = 64 < 66 = ¿42,3) =
Л2')=1
= max, {IS.Vi + 10.v2}.
fibclje2)=0
Таким образом, среди предложенных вариантов трактовки близости или отличия двух k-значных функций для характеризации пороговых функций можно выделить коэффициенты возрастания. При этом аналогично булевому случаю подтверждается предположение о необходимости введения итеративной процедуры для нахождения коэффициентов линейной формы.
5. ИТЕРАТИВНАЯ ПРОЦЕДУРА ДЛЯ НАХОЖДЕНИЯ КОЭФФИЦИЕНТОВ ЛИНЕЙНОЙ ФОРМЫ
В Примере 5 можно увидеть необходимость введения итеративной процедуры для нахождения коэффициентов линейной формы. На идейном уровне такая процедура может быть изложена следующим образом.
Предположим, что некоторая линейная форма ...,хп), построенная с помощью коэффициентов роста, не даёт чистое разделение областей значений функции /fc(*i>..., хп), то есть существует i G 0, к — 2 такое, что выполняется неравенство
/Ч
max {f^fi н-----ь atl£fl} >
> min ia1sl + "' + an€n).
fk(ev....£n)=i+l
Рассмотрим точки
А = (ц±.....гО 6 М: = {.V: /(л) = £} и
В = Оч,...,!',,) € М,+1 = {*:/(*) = I + 1},
такие, что линейная форма 1(х1,.,хп) принимает на них максимальное значение из множества Мг и минимальное значение из множества М[+1 соответственно. Заметим, что таких точек может быть несколько. Для них выполняется неравенство Ь(А) > Ь(В). Уравнение Ь(х1,.,хп) — Ь, Ь £ Е задает гиперплоскость, рассекающую п-мерный куб.
Для новой линейной формы L\x1,.,xn) подсчитываются её значения на множествах
Рисунок. 6.
Алгоритм корректировки должен наклонить множество гиперплоскоскостей, задающихся линейной формой Ь(х1г... ,хп) так, чтобы при смещении гиперплоскость отсекала точки А и В в другом порядке, чем до корректировки.
Отклонение «в нужную сторону» можно обеспечить прибавлением к линейной форме Ь(х1г... ,хп) двух линейных форм, соответствующих двум рассматриваемым точкам А и В. Линейные формы образованы координатами точек и берутся с разными знаками: в первую линейную форму подставляются координаты точки В, во вторую - координаты точки А с отрицательными знаками.
М^ = [x:f(x) = i], iE0,k-l, и определяется, даёт ли она чистое разделение областей значений функции fk(x1,... ,xn). В случае, если нет, то процедура повторяется. Если же после корректировки линейная форма даёт чистое разделение областей значений функции fk{x1,.,xn), то алгоритм заканчивает работу: коэффициенты линейной формы найдены.
Возможно, что в случае зацикливания алгоритма необходимо будет прибавлять корректирующие линейные формы с коэффициентами меньшими единицы:
Пример 6. Рассмотрим коррекцию линейной формы
L(x1,x2) = l8x1 + l0x2, полученной с помощью коэффициентов возрастания в Примере 5. Чистому разделению мешают точки (3,1) и (2,3). Прибавим соответствующие линейные формы:
Полученная линейная форма дает чистое разделение областей значений функции f.{x1,x2) и исправляет нарушение разделения в корректируемых точках:
min
f*
{19.1
г)=1\3,1)
42,3)
= max {19.VJ + 8х2 }.
f*(x
лч)=0
В заключение необходимо отметить, что с теоретической точки зрения задача нахождения порогового представления к-значной функции f(х1, ...,хп) сводится к решению системы неравенств, вообще говоря, двухсторонних для каждого значения I вида:
^ а-^х-^ + а2х2 + """ + апхп < где х1,х2,.,хп - известные координаты векторов, на которых функция принимает соответствующие значения, а параметры а1,а2,... ,ап,Ь0,Ь1,... ,Ьк - неизвестные. Для решения системы такого вида применим полиномиальный алгоритм Хачияна [7], что
обеспечивает принципиальное решение рассматриваемой задачи. Однако высокая вычислительная сложность и большое число необходимых операций алгоритма Хачияна оставляет актуальным поиск прямых конструктивных путей решения данной задачи, в частности тех, которые изложены в данной статье.
Список литературы:
1. Бутаков Е.А. Методы синтеза релейных устройств из пороговых элементов. Москва, Энергия. 1970.
2. Дертоузос М. Пороговая логика. Москва, Мир. 1967.
3. Зуев А.Ю. Пороговые функции и пороговые представления булевых функций. // «Математические вопросы кибернетики», выпуск 5, 1994 г.
4. Никонов В.Г. Пороговые представления булевых функций // «Обозрение прикладной и промышленной математики», 1994, Т. 1, вып. 3.
5. Никонов В.Г., Никонов Н.В. Особенности пороговых представлений к-значных функций // «Труды по дискретной математике», 2008, Т. 11.
6. Вальцев В.Б., Григорьев В.Р., Никонов В.Г. Некоторые структурные принципы организации высших функций мозга // «Нейрокомпьютер как основа мыслящих ЭВМ. - М.: Наука, 1993.
7. Хачиян Л.Г. Полиномиальные алгоритмы в линейном программировании // «Журнал вычислительной математики и математической физики», 1980, Т 20.
8. Золотых Н.Ю. Расшифровка пороговой функции, заданной расширенным оракулом// Вестник Нижегородского университета им. Н.И. Лобачевского, 2012, №3(1).