\documentclass{cmi}
\usepackage{upgreek}
\usepackage{multirow}
\usepackage{booktabs}

% Данные поля заполняются редакцией журнала
%\makeatletter
%\newcommand{\@volume}{}
%\newcommand{\@issue}{}
%\newcommand{\@num}{}
%\newcommand{\@pages}{}

% Подключение bib файлов
\addbblresources{biblio_ru.bib}{biblio_en.bib}

\begin{document}

% ЗАПОЛНЯЕТСЯ АВТОРАМИ: Релевантные коды УДК статьи
% (см. http://teacode.com/online/udc/)
\classify{УДК 004.81}

%\renewcommand*{\thefootnote}{\fnsymbol{footnote}} %to make a number in end of title delete this line

\title{ДИАГНОСТИКА НЕИСПРАВНОСТЕЙ ПОДШИПНИКОВ КАЧЕНИЯ \\
С ИСПОЛЬЗОВАНИЕМ ПИКОВ СПЕКТРА \\
И НЕЙРОННЫХ СЕТЕЙ
%	\footnote{Если статья рекомендована к публикации программным
%		комитетом научной конференции, это указывается в сноске к названию статьи. Благодарность за финансовую %поддержку в подготовке статьи необходимо поместить после заключения.}%
}

%\authors{}

%\organizations{}

%\newcommand{\emailsinsert}{}

% ЗАПОЛНЯЕТСЯ АВТОРАМИ: Дата поступления статьи в редакцию
\newcommand{\rdate}{20.04.2022}
%\emails{\emailsinsert}
\redactiondate{\rdate}

\renewcommand*{\thefootnote}{\arabic{footnote}}
\setcounter{footnote}{0}

% Оформление аннотации
\begin{abstract}%
Важнейшими составляющими деталей машин являются подшипники качения, контроль за состоянием которых необходим, так как возможные дефекты в их конструкции могут привести к неправильной работе или общему выходу машин из строя. Современные решения по диагностике неисправностей подшипников обычно используют сложные процессы  извлечения признаков, например, построение их изображений спектра Гильберта и дальнейшую мощную нейронную сеть для их классификации. В этой статье мы предлагаем простой, но, тем не менее, эффективный алгоритм решения данной задачи. Для выделения признаков из сигнала мы делим спектр сигнала на равные подинтервалы и находим максимум амплитуды и соответствующее значение частоты в каждом из них. В статье, на основе метода t-SNE, показано, что выделенные таким образом признаки, несмотря на свой небольшой размер, хорошо представляют разного типа сигналы. На втором этапе выделенные признаки поступают на вход простой нейронной сети классификатора. Предложенный метод обладает простотой в вычислительном отношении, как на этапе выделения признаков, так и на этапе обучения нейронной сети. Несмотря на это, метод дает 100\% точность для всех типов сигналов на коротких данных из набора данных IMS.


\keywords{Диагностика неисправностей; подшипник качения; спектр сигнала; нейронные сети; t-SNE}
\end{abstract}

%Образец цитирования статьи
\forcitation{Ибряева~О.Л., Мохаммад~М.Н. Диагностика неисправностей подшипников качения с использованием пиков спектра и нейронных сетей}

%Верхние колонтитулы
\markboth{О.Л.~Ибряева, М.Н.~Мохаммад}{Название статьи с многоточием в случае нехватки горизонтального\ldots}

% Заголовки разделов формируются при помощи команд  \section{}, \subsection{}, \subsubsection{}
\section*{Введение}
\label{sec-intro}
Подшипники качения играют ключевую роль практически в любой системе с механизмом вращения, такой как корабли, станы холодной прокатки и многие другие~\bblcite{Yu2017}. Ввиду чрезвычайной важности обнаружения неисправностей подшипников в области машиностроения к настоящему времени существует множество подходов к их диагностике. В погоне за большей точностью современные подходы часто используют вычислительно затратные методы, например, получают признаки с помощью построения спектральных изображений, на которых затем обучают мощную сверточную нейронную сеть~\bblcite{Xu2021,Han2021}.

Целью настоящей статьи явилась разработка вычислительно простого, но эффективного метода диагностики подшипников качения, который основан на применении быстрого преобразования Фурье (БПФ) к вибрационным сигналам подшипника. Извлеченные из спектра сигнала признаки используются далее в качестве входных данных нейронной сети для определения типа неисправности. Репрезентативность выделенных признаков показана с помощью метода t-Distributed Stochastic Neighbor (t-SNE) ~\bblcite{Hamad2019}. Предложенный метод протестирован на эталонном наборе данных IMS ~\bblcite{mnmWeb001} для диагностики подшипников и показал высокую эффективность.
Дополнительные преимущества нашего нового метода заключаются в простой процедуре выделения признаков, их небольшом количестве и простой архитектуры нейронной сети.

Основные ключевые моменты данной статьи:
\begin{itemize}
	\item[--]Разработан новый подход к выделению признаков из вибрационных сигналов подшипника.
	
	\item[--]Отличительной особенностью нового метода является простота процесса извлечения признаков.
	
	\item[--]Новый метод использует небольшое количество признаков и простую нейронную сеть по сравнению с другими современными подходами.
		
	\item[--]Эксперименты, проведенные на наборе данных IMS, доказывают высокую эффективность метода.
	
\end{itemize}

Структура данной работы такова. \S~\ref{sec-relatedWorks} содержит обзор аналогичных работ. В \S~\ref{sec-datset} описана экспериментальная установка и набор данных IMS. \S~\ref{sec-prop} содержит описание процесса извлечения признаков из вибрационных сигналов подшипника и структуры нейронной сети. Экспериментальные исследования представлены в \S~\ref{sec-res}. Последний раздел содержит заключение и выводы исследования.



\section{Обзор аналогов}
\label{sec-relatedWorks}

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

Обработка сигналов на основе вибрации является одним из наиболее часто используемых методов на первом этапе, обычно с использованием анализа во временной области, частотной области или частотно-временной области. Анализ в частотной области обычно производится с помощью быстрого преобразования Фурье (БПФ), а также методов на основе БПФ, методов спектрального анализа ~\bblcite{wang2015enhanced} и т. д. Анализ во временной области рассчитывает статистические параметры, такие как среднеквадратичное значение виброускорения, эксцесс, ударные импульсы ~\bblcite{zhen2008bearing} и т. д. Частотно-временной анализ используется для расширения возможностей анализа в частотной области на нестационарные вибрационные сигналы. и включает в себя такие методы, как краткосрочное преобразование Фурье ~\bblcite{wang1993early}, вейвлет-преобразование ~\bblcite{luo2003line}, разложение по эмпирическим модам преобразования Гильберта---Хуанга ~\bblcite{lei2013review} и т. д.

На втором этапе извлеченные признаков используются в качестве входных данных для методов машинного обучения с целью распознавания ошибок. Среди этих методов есть метод k-NN, метод опорных векторов и методы на основе дерева решений, но, конечно, преобладают нейронные сети, особенно сверточные и рекуррентные нейронные сети ~\bblcite{yuan2018intelligent,liang2019intelligent,Shao2017}. Например, Лю в ~\bblcite{Liu2018} предложил систему диагностики неисправностей подшипников качения с рекуррентной нейронной сетью (RNN) и автоматическим кодировщиком. В работе ~\bblcite{Xu2020} представлен метод, использующий разложение по вариационным модам на английском языке – Variational Mode Decomposition (VMD)  и глубокие сверточные нейронные сети для выполнения классификации неисправностей подшипников качения в ветряных турбинах.

В целом, в диагностике неисправностей подшипников наблюдается тенденция к более сложным и глубоким моделям машинного обучения. Но глубокое обучение обнажает проблемы большого размера выборки и длительного времени обучения, а в реальных условиях работы количество доступных данных с помеченными ошибками относительно невелико, поэтому модель сложной архитектуры с хорошей обобщающей способностью и высокой точностью трудно обучить. Некоторые исследователи уже говорят о необходимости других, более простых и прикладных методов ~\bblcite{Chen2020,Zhang2017a,Li2020}.


\section{Описание набора данных}
\label{sec-datset}

Для проверки работоспособности предлагаемого метода в данной статье используются экспериментальные данные IMS (Intelligent Maintenance Systems) ~\bblcite{mnmWeb001}, представленные Центром интеллектуальных сервисных систем Университета Цинциннати. Экспериментальная установка представлена на рис.~\ref{figOne}. На валу были установлены четыре двухрядных подшипника Rexnord ZA-2115 с двумя акселерометрами PCB 353B33 на каждом. Скорость вращения вала поддерживалась постоянной на уровне 2000 оборотов в минуту. Он приводился в движение двигателем переменного тока, который был соединен с валом фрикционными ремнями. К подшипникам и валу с помощью пружинного механизма добавлялась нагрузка в 6000 фунтов. Частота дискретизации была установлена равной 20480 Гц. После вращения более чем на 100 миллионов оборотов в системе проявились такие неисправности подшипников, как дефект внутреннего кольца, внешнего кольца и шарика.

\begin{figure}[h]%
	\centering%
	\includegraphics*[width=6.43in, height=2.36in,    keepaspectratio=false]{Figure 1.png}%
%	\includegraphics*[width=6.43in, height=2.36in,    keepaspectratio=false]{Figure 1.jpg}%
	\caption{Тестовая платформа IMS}%
	\label{figOne}%
\end{figure}%

IMS содержит наборы вибрационных сигналов подшипника с появившимися дефектами. Каждый набор содержит все записи от состояния исправного подшипника до его поломки.
Мы рассмотрели сигналы каждого типа продолжительностью 20 секунд, как это было сделано в работе [19]. Для сигналов, отвечающих нормальному состоянию подшипника, мы взяли первые 20 секунд из набора данных для подшипника 1, когда он (как и все подшипники) был еще исправен. Последние 20 с записи для того же подшипника уже соответстуют его дефекту внутреннего кольца. Для сигнала с дефектом шарика были выбраны последние 20 секунд записи для подшипника 4. Наконец, запись с 2004-04-17 21:52:55 до 2004-04-18 01:02:55 с подшипника 3 содержит данные с дефектом внешнего кольца.

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

Для расширения данных сигналы вибрации были сегментированы на перекрывающиеся кадры по 4260 отсчетов, что эквивалентно (4260)/20480=0.208 секунды, с процентом перекрытия 50\%, т. е. 2130 отсчетов, как показано на рис.~\ref{fig2}. В табл.~\ref{tabOne} представлены классы сигналов набора данных IMS и количество образцов в каждом.

\begin{figure}[h]%
	\centering%
	\includegraphics*[width=3.78in, height=3.16in, keepaspectratio=false]{Figure 2.png}%
	\caption{Сегментация сигнала}%
	\label{fig2}%
\end{figure}%

\begin{table}[h]%
	\centering
	\ttabbox[\FBwidth]{	\caption{Классы набора данных IMS}\label{tabOne}}
	{\begin{tabular}
{|p{0.7in}|p{0.8in}|p{1.4in}|p{1.3in}|} \hline
\centering \textbf{Класс} & \centering \textbf{Тип дефекта} & \centering \textbf{Число образцов в обучающей выборке} &  \multicolumn{1}{p{1.3in}|}{ \centering \textbf{Число образцов в тестовой выборке}} \\ \hline
 \centering 1 & нет & \centering  134 & \multicolumn{1}{p{1.3in}|}{ \centering 56} \\ \hline
 \centering 2 & дефект внутреннего кольца & \centering  134 & \multicolumn{1}{p{1.3in}|}{ \centering 56} \\ \hline
 \centering 3 & дефект внешнего кольца &  \centering  134 &  \multicolumn{1}{p{1.3in}|}{ \centering 56} \\ \hline
 \centering 4 & дефект шарика & \centering  134 &    \multicolumn{1}{p{1.3in}|}{ \centering 56} \\ \hline
\end{tabular}}
\end{table}

\section{Предлагаемый метод}
\label{sec-prop}

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

\subsection{Извлечение признаков}

На рис.~\ref{fig3} приведены примеры спектров сигналов разного типа. Можно видеть, что сигналы разного типа имеют существенно разные спектральные формы. Преобладающие спектральные компоненты каждого типа сигнала расположены в разных частотных диапазонах.


\begin{figure}[h]%
	\centering%
	\includegraphics*[width=6.53in, height=4.5in, keepaspectratio=false]{Figure 3.png}%
	\caption{Спектр различных сигналов; а) нормальный; б) дефект внутреннего кольца; в) дефект внешнего кольца; г) дефект шарика}%
	\label{fig3}%
\end{figure}%


\begin{figure}[h]%
	\centering%
	\includegraphics*[width=6.13in, height=4.0in,  keepaspectratio=false]{Figure 4.png}%
	\caption{Разбиение спектра на интервалы и нахождение максимальных значений каждого интервала}%
	\label{fig4}%
\end{figure}%

Наш метод выделения признаков основан на этой особенности сигналов. Мы разделяем ось частот на равные интервалы и внутри каждого ищем частоту с макимальным значением амплитуды. Спектр сигнала при этом предварительно сглажен.

Извлекаемые признаки представляют собой два вектора: первый содержит значения частот с максимальной амплитудой в каждом частичном интервале, второй вектор содержит значения этих амплитуд, другими словами, признаки представляют собой пары $\left(x,y\right)\ \in \ R^{N\times 2}$, где $N$ — количество интервалов.

Мы поделили весь интервал частот в спектре от 0 до 10240 Гц (половина частоты дискретизации) на $N=106$ равных интервалов, каждый из которых имеет длину примерно 100Гц (рис~\ref{fig4}). Таким образом, каждому образцу сигнала из таблицы~\ref{tabOne} мы поставили в соответствие два 106 мерных вектора (значений частот и амплитуд), вычисленные по его спектру.

Чтобы проверить, что извлеченные признаки адекватно представляют объекты различных классов, мы использовали метод t-SNE [4], который позволяет нам отобразить точки из многомерного пространства в трехмерное и посмотреть на их расположение в пространстве.

На на рис.~\ref{fig5} построено такое представление для значений $x$ пар $\left(x,y\right)$, полученное с помощью официального сайта TensorFlow (https://projector.tensorflow.org/). Можно видеть, что 4 типа сигналов сгруппированы в 4 отдельных кластера, что указывает на хорошее качество выделенных признаков и на то, что их можно успешно использовать для решения задачи классификации.
\begin{figure}[h]%
	\centering%
	\includegraphics*[width=4.1in, height=3.0in, keepaspectratio=false]{Figure 5.jpg}%
	\caption{Трехмерная визуализация извлеченных признаков с использованием t-SNE всех типов сигналов IMS}%
	\label{fig5}%
\end{figure}%


\subsection{Модель нейронной сети}

В качестве классификатора, который принимает на вход признаки, извлеченные методом описанным в прошлом разделе, мы использовали нейроннную сеть простой архитектуры. Сеть состоит из двух отдельных модулей, обрабатывающих значения максимальных частот и их амплитуд, соответственно. За входным слоем каждого модуля следуют два скрытых слоя, после чего модули объединяются. Далее следуют еще два скрытых слоя и, наконец, выходной слой из 4 нейронов (по количеству классов), как показано на рис~\ref{fig6}.

При обучении мы традиционно использовали категориальную кросс-энтропию в качестве функции потерь и оптимизатор Адам. Первоначально количество эпох обучения было установлено на 2000, но мы использовали коллбэки Keras: EarlyStopping и ModelCheckpoint, которые на практике уменьшили это число эпох. Эти коллбэки сохраняют лучшую нейронную сеть (с наименьшим значением функции потерь на валидационном наборе). При отстутствии улучшения (уменьшения функции потерь) в течение patience эпох, они останавливают обучение и берут ранее сохраненную лучшую модель. Параметр patience («терпение») был установлен равным 100 эпох.

\begin{figure}[h]%
	\centering%
	\includegraphics*[width=4.0in, height=4.3in, keepaspectratio=false]{Figure 6.jpg}%
	\caption{Схема используемой нейронной сети}%
	\label{fig6}%
\end{figure}%


\section{Результаты экспериментов}
\label{sec-res}

Используем новый предложенный метод для классификации сигналов из набора данных IMS
и, чтобы иметь возможность сравнить с уже существующим результатом в ~\bblcite{Zhang2017}, проведем эксперимент в похожих условиях. Мы также, как и авторы статьи ~\bblcite{Zhang2017}, взяли небольшие отрезки данных (20 секунд) и разделили выборку на обучающие выборки (70\%) и тестовые выборки (30\%) набора.

Для оценки нашей модели используется метрика Accuracy. Accuracy  вычисляется по следующей формуле:
\[Accuracy=\frac{TP+TN}{TP+TN+FP+FN} ,\]

где TP ---истинно положительные образцы (true positive),TN ---истинно отрицательные образцы (true negative), FP ---ложно положительные образцы (false positive), и FN ---ложно отрицательные образцы (false negative).

Авторы статьи ~\bblcite{Zhang2017} решали данную задачу классификации на основе глубокой нейронной сети и, как показал их эксперимент (табл.~\ref{tabtwo}), наилучшие результаты показала сеть с 5 слоями со средней точностью 94.4\% на тестовых данных. В табл.~\ref{tabtwo} приведен также результат, полученный нашим методом, он дает точность 100\% для всех типов сигналов.

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


\begin{table}[h]%
	\centering
	\ttabbox[\FBwidth]{	\caption{Сравнение результатов предложенного решения с [19].}\label{tabtwo}}
	{
\begin{tabular}
{|p{0.5in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|p{0.4in}|} \hline
\multicolumn{1}{|p{0.5in}}{\centering \textbf{Класс} } & \multicolumn{8}{|p{4.5in}|}{ \centering \textbf{  Глубокие нейронные сети,\\ учитывающие временную когерентность}} & \multicolumn{2}{p{0.8in}|}{ \centering \textbf{Наш метод}} \\\cline{2-11}
 & \multicolumn{2}{p{0.8in}|}{3-х слойная модель} & \multicolumn{2}{p{0.8in}|}{4-х слойная модель} & \multicolumn{2}{p{0.8in}|}{5-х слойная модель} & \multicolumn{2}{p{0.8in}|}{6-х слойная модель} &  \multicolumn{1}{p{0.4in}|}{} &\multicolumn{1}{p{0.4in}|}{} \\\cline{2-9}
 & \multicolumn{1}{p{0.4in}|}{\centering прове- рка, \%} & \multicolumn{1}{p{0.4in}|}{\centering тест, \%} & \multicolumn{1}{p{0.4in}|}{прове- рка, \%} & \multicolumn{1}{p{0.4in}|}{тест, \%} & \multicolumn{1}{p{0.4in}|}{\centering прове- рка, \%} & \multicolumn{1}{p{0.4in}|}{\centering тест, \%}& \multicolumn{1}{p{0.4in}|}{\centering прове- рка, \%} & \multicolumn{1}{p{0.4in}|}{\centering тест, \%} & \multicolumn{1}{p{0.4in}|}{\centering прове- рка, \%} & \multicolumn{1}{p{0.4in}|}{\centering тест, \%} \\ \hline


\centering 1 & \centering 98.9 & \centering 98.2 \centering & \centering 99.3& \centering 98.2 & \centering 98.2 & \centering 98.9 & \centering 97.8 & \centering 98.9 & \centering 100  & \multicolumn{1}{p{0,4in}|}{ \centering 100}    \\ \hline
\centering 2 & \centering 92.3 & \centering 92.6 & \centering 91.9 & \centering 91.5 & \centering 93.8 & \centering 91.5 & \centering 90.1 & \centering 90.8 & \centering 100  & \multicolumn{1}{p{0,4in}|}{ \centering 100}    \\ \hline
\centering 3 & \centering 100 & \centering 100 & \centering 100 & \centering 100 & \centering 100 & \centering 100 & \centering 100 & \centering 100 & \centering 100  & \multicolumn{1}{p{0,4in}|}{ \centering 100}    \\ \hline
\centering 4 & \centering 93.0 & \centering 84.9 & \centering 90.8 & \centering 85.7 & \centering 91.9 & \centering 87.1 & \centering 87.9 & \centering 83.5 & \centering 100  & \multicolumn{1}{p{0,4in}|}{ \centering 100}   \\ \hline
\centering среднее &  \centering \textbf{96.6} & \centering 93.9 & \centering95.5 &\centering93.8 &  \centering \textbf{96.0} & \centering \textbf{94.4} & \centering 93.9 & \centering 93.3 &   \centering \textbf{100}  & \multicolumn{1}{p{0,4in}|}{ \centering \textbf{100}}   \\ \hline
\end{tabular}
}
\end{table}





\section*{Заключение}
\label{sec-conclusion}

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

Основная идея предложенного решения заключается в определении новых признаков, которые могут эффективно характеризовать особенности спектров сигналов, отвечающих различным состояниям подшипников качения. Для извлечения признаков мы разделили ось частот на равные интервалы и внутри каждого нашли частоту с максимальным значением амплитуды. Извлеченные признаки представляют собой два вектора: первый содержит значения частот с максимальной амплитудой в каждом частичном интервале, второй вектор содержит значения этих амплитуд. Как показано в работе, данные признаки являются информативными и позволяют хорошо различить сигналы с подшипников разного типа.

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

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

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





%\vspace{1em}
%{\it Благодарность за финансовую поддержку в подготовке статьи необходимо поместить после заключения.}


\insertbiblio

% Сведения об авторах (Ф.И.О.(полностью), степень, звание, кафедра (отдел), вуз, e-mail)
Ибряева Ольга Леонидовна, к.ф.-м.н., доцент, кафедра системного программирования, Южно-Уральский государственный университет (национальный исследовательский университет) (Челябинск, Российская Федерация)

Мохаммад Мохаммад, студент, кафедра системного программирования, Южно-Уральский государственный университет(национальный исследовательский университет) (Челябинск, Российская Федерация)



%\newpage

{\phantom{eee} \hrule \vskip 4 pt}
\classify{} % Код MSC не нужен
%\classify{MSC 00000} % Коды MSC см. http://www.ams.org/msc/

\title{DIAGNOSIS of ROLLING BEARINGS FAULTS USING SPECTRAL-FINGERPRINT \\ and NEURAL NETWORKS}

%\authorsen{}

%\organizations{}

%\emails{\emailsinsert}
\redactiondateen{\rdate}

\begin{abstract}%
The most important components of machine parts are rolling bearings, the condition of which is necessary to control, since possible defects in their design can lead to incorrect operation or general failure of machines. Modern solutions for diagnosing bearing failures typically use complex feature extraction processes, such as their Hilbert spectrum imaging and a further powerful neural network to classify them. In this article, we propose a simple, but nevertheless effective algorithm for solving this problem. To extract features from a signal, we divide the signal spectrum into equal subintervals and find the amplitude maximum and the corresponding frequency value in each of them. In the article, based on the t-SNE method, it is shown that the features selected in this way, despite their small size, well represent different types of signals. At the second stage, the selected features are fed to the input of a simple classifier neural network. The proposed method is computationally simple, both at the stage of feature extraction and at the stage of neural network training. Despite this, the method gives 100\% accuracy for all types of signals on short data from the IMS dataset.

\keywordsen{Diagnostics of malfunctions; rolling bearing; signal spectrum; neural networks; t-SNE }
\end{abstract}

\forcitationen{O.L. Ibryaeva, M.N. Mohammad. Diagnostics of Rolling Bearing Failures using Spectrum Peaks and Neural Networks}

\openaccess

\insertbiblioen

\end{document}
