%-----------MikTex------------------
\documentclass[12pt]{article} % это понятно
\usepackage{a4wide} % аналог fullpage.sty
%\usepackage{a4} % стандартные значения полей страницы
\usepackage[cp1251]{inputenc} % подключение кодировки Windows
%\usepackage[cp866]{inputenc} % подключение кодировки MS DOS
\usepackage[russian]{babel} % подключение русификации (переносы)
\usepackage {graphicx} % подключение графики
\usepackage{epsfig}                                                 %%
\usepackage{amssymb}                                                %%
\usepackage{amsmath}                                                %%
\usepackage{texdraw}                                                %%
\usepackage{rotating}                                               %%
\usepackage{amsfonts}                                               %%
%-----------MikTex------------------
\usepackage{demopack}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand{\rank}{\mathop{\rm rank}\nolimits}                      %%
\newcommand{\ie}{т.$\!$~е. }                                          %%
\newcommand{\td}{т.$\!$~д. }                                          %%
\newcommand{\tp}{т.$\!$~п.}                                          %%
\newcommand{\dr}{др. }                                          %%
\newcommand{\smr}{см.~рис.~}                                        %%
\newcommand{\ris}{рис.~}                                            %%
\newcommand{\tab}{табл.~}                                           %%
\newcommand{\smt}{см.~табл.~}  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand{\abs}[1]{\mathop{\left|#1\right|}\nolimits}  %%
\newcommand{\E}[1]{\mathop{{\rm \bf E}\left\{#1\right\}}\nolimits}  %%

\newcommand{\sgn}[1]{\mathop{\rm sign \ }[#1]}                      %%
\newcommand{\hm}[1]{#1\nobreak\discretionary{}{\hbox{\ensuremath{#1}%%
}}{}}%%продолжение команды дублирования знаков при переносе формулы%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%
\emergencystretch=4pt
\hfuzz=0.5pt
\tolerance=400
%\pagestyle{empty}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%*** SRART ***%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}
%\begin{center}
{\noindent
\begin{tabular}{@{}p{2.8cm}@{}@{}p{7.6cm}p{6.14cm}@{}}
\raisebox{2.7mm}{\begin{tabular}{@{}p{2.8cm}@{}}
\centerline%
{\includegraphics[scale=0.85]{a4}}\\
\raisebox{-0.4mm}{\centerline{\footnotesize 08.12.2009}}
\end{tabular}}
&
\begin{tabular}{@{}p{7.6cm}@{}}\textbf{\large Ульяновский}\\ 
\textbf{\large государственный} 
\textbf{\large университет}\\
ул. Л. Толстого 42 \\ 432970 Ульяновск, Россия \end{tabular}
&
\begin{tabular}{@{}p{6.14cm}@{}} {\small каф. Информационные\hfill  технологии} \\ \hline
{\small тел: \hfill +7 (8422) 32-1029}\\
%Fax: \hfill +7 (8422) 41-2340\\
{\texttt {\small http://staff.ulsu.ru/semoushin/}} \\
{\small e-mail:}\hfill  {\small \texttt{innokentiyvsem@gmail.com}} \end{tabular}\\\hline
\end{tabular}%\hrule
}\\[2ex]
%\end{center}
%{Проректору по научной работе\\ и информационным технологиям\\
%д.ф.-м.н, проф. В. Н. Голованову\\
%д.т.н, проф. И. В. Семушина}

\begin{center}
%\vspace{0.5cm}
\noindent
\begin{tabular}{p{14.5cm}}
%\vspace{1pt}
\centerline{{\bf ОТЗЫВ}}
\centerline{официального оппонента диссертации В.~Е.~Зотеева}\\
\centerline{ЛИНЕЙНО-ПАРАМЕТРИЧЕСКИЕ ДИСКРЕТНЫЕ МОДЕЛИ В ФОРМЕ}
\centerline{РАЗНОСТНЫХ УРАВНЕНИЙ В ЗАДАЧАХ ИДЕНТИФИКАЦИИ}
\centerline{ДИССИПАТИВНЫХ МЕХАНИЧЕСКИХ СИСТЕМ}\\
\centerline{на соискание ученой степени доктора технических наук}
\centerline{по специальности 05.13.18---<<Математическое моделирование,} 
\centerline{численные методы и комплексы программ>>}\\
\hline
%\\[-2ex]
\end{tabular}
\end{center}
\bigskip
%\vspace{0.5cm}


\medskip

На оппонирование представлена диссертационная работа на 400 листах, включающая 126 рисунков и 52 таблицы. 
Работа разделена на Введение, шесть глав основной части, Заключение и Список используемых источников и литературы (206 наименований, из которых 68 --- ссылки на труды автора). Приложение  (58 страниц) содержит документы о внедрении / использовании результатов работы и листинги разработанного пакета прикладных программ.\\[-1.5ex]

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{A}     %{Актуальность}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[1.] Актуальность
\end{description}
Работа посвящена актуальной теме. Она находится в русле приоритетных, фундаментальных направлений совершенствования численных методов и реализуемых на их основе математических информационных технологий для успешного развития отечественного машиностроения. 
На любом этапе жизненного цикла сложной механической системы необходим надежный контроль ее технического состояния. 
Математически, задачи контроля оказываются достаточно сложными. 
Возникнув в технике, они решаются средствами математики ---  теории обнаружения внезапных нарушений, теории обнаружения моментов нарушений, теории диагностики нарушений или теории идентификации моделей систем. 
Практически все эти теории являются <<model-based>>, \ie основанными на моделях, чаще всего --- параметрических моделях. 
Построение модели --- определение ее структуры, а затем и ее параметрическое <<наполнение>> --- служит отправной точкой упомянутых теорий. 
По сути дела, то, как построена модель, предопределяет или во многом определяет то, каким будет метод решения возникающих задач.\\[-1.5ex]

%\textbf{1.} 
\emph{Какие же препятствия существуют на данный момент в области построения моделей процессов в механических системах?}

Первое и основное препятствие --- нелинейность моделей, далее --- их многомерность, многосвязность и иногда --- распределенность. 
Из этого набора автор данной диссертационной работы выбирает первое  препятствие --- нелинейность, \ie его объект исследования --- нелинейные диссипативные механические системы (ДМС), сводимые при их модельном (теоретическом) описании к системам скалярных дифференциальных уравнений четырех видов относительно модельной координаты $\tilde y(t)$:
\begin{equation}\label{e:kulon}
m\tilde y^{\prime\prime}(t)+b\frac{\tilde y^{\prime}(t)}{\abs{\tilde y^{\prime}(t)}}+c\tilde y(t)=P(t),
\quad \mbox{\rm{---кулоново трение}}
\end{equation}
\begin{equation}\label{e:lin-vyaz}
m\tilde y^{\prime\prime}(t)+b\tilde y^{\prime}(t)+c\tilde y(t)=P(t),
\quad \mbox{\rm{---линейно-вязкое трение}}
\end{equation}
\begin{equation}\label{e:turbo}
m\tilde y^{\prime\prime}(t)+b\tilde y^{\prime}(t)\abs{\tilde y^{\prime}(t)}+c\tilde y(t)=P(t),
\quad \mbox{\rm{---турбулентное трение}}
\end{equation}
\begin{equation}\label{e:hister}
m\tilde y^{\prime\prime}(t)+b\sgn{\tilde y^{\prime}(t)}\cdot\frac{a^n}{\pi}\sqrt{1-\frac{\tilde y^2(t)}{a^2}}+c\tilde y(t)=P(t),
\quad \mbox{\rm{---гистерезисное трение}}
\end{equation}
((\ref{e:lin-vyaz}) --- единственный случай, когда модель линейная). 
Отмечу, что в диссертации эти модели записаны на сс.~22--23 относительно не  $\tilde y(t)$, а  $y(t)$, хотя в дальнейшем (начиная со с.~145) под $\tilde y_k$ понимаются именно значения отклика модели (ЛПДМ --- линейно-параметрической дискретной модели), формируемые по решениям $\tilde y(t)$  уравнений (\ref{e:kulon})--(\ref{e:hister}). При этом  $y_k$ обозначает измеренное значение отклика $y(t)$ реальной ДМС, а $\varepsilon_k$ --- отклонение $y_k$ от $\tilde y_k$, так что $y_k=\tilde y_k+\varepsilon_k$.

Основная, базовая часть исследований автора заключается в построении ЛПДМ для различных типов ДМС по ее отклику  $\tilde y(t)$  на тестовые воздействия $P(t)$ нескольких типов и для различных областей съема отклика ДМС: 
\begin{itemize}
\item если отклик снимается во временн\'ой области --- по виброграмме (осциллограмме, сейсмограмме и \tp), по кривой разгона (переходная характеристика ДМС), по кривой ползучести (неупругого деформируемого материала) или по огибающей колебаний ДМС, то ЛПДМ, естественно, строится в той же области,   --- основной материал глав 2, 3, 4 и 5, 
\item если отклик снимается в частотной области --- по экспериментально построенной амплитудно-частотной характеристике системы, то ЛПДМ строится в той же области, ---  п.~5.3.
\end{itemize}
Независимо, в какой области фиксируются отклики и, соответственно,  строится ЛПДМ, принцип  $y_k=\tilde y_k+\varepsilon_k$ сохраняется, а именно: $\varepsilon_k$ обозначает расхождение (невязку) между откликом $y_k$ реальной системы (ДМС) и откликом $\tilde y_k$ ее модели (ЛПДМ), которую автор предлагает строить и показывает, как это делать, --- применяя различные аппроксимации, для каждого конкретного приложения с точностью до некоторого вектора параметров $\lambda$.
При этом задача приводится к подгонке параметров $\lambda$: по экспериментальным замерам  $y_k$ и  $\tilde y_k$ ($k=1,\ldots, N$, $N$ называется объемом выборки), полученным при одинаковом (для ДМС и ЛПДМ) плане эксперимента $F$ и сведенным, соответственно, в векторы $b$ и  $F\lambda$, найти численным методом такие значения $\lambda^\star$ параметров, которые минимизируют квадрат евклидовой нормы эквивалентного расхождения $\eta=b-F\lambda$. 
Таким образом, автор применяет методологию метода наименьших квадратов (МНК) к широкому кругу задач, возникающих при технической диагностике сложных нелинейных ДМС. \\[-1.5ex]

%\textbf{2.} 
\emph{Почему эта методология не применялась в этой области ранее? Казалось бы, это должно было произойти раньше? Ответ на это любой может увидеть, если возьмется за такую работу, --- здесь существуют большие препятствия!}  \\[-1.5ex]


Главные препятствия, которые преодолевает автор:
\begin{itemize}
\item
Во-первых, неочевидно, как --- (это очень непросто!) --- построить необходимые ЛПДМ. Автор это делает весьма изобретательно. 
Даже там, где, казалось бы, это невозможно сделать в силу сложных нелинейных зависимостей, он делает это математически скрупулезно, обосновывая вынужденные аппроксимации. 
Все выкладки, необходимые для этого, он приводит в полном объеме, методично и строго, преследуя одну и ту же цель --- разработать \emph{научный подход} к построению линейно-параметрических дискретных моделей по тем модельным зависимостям --- аналитическим решениям уравнений (\ref{e:kulon})--(\ref{e:hister}), которые изначально являются континуальными (не являются дискретными, \ie не связывают между собой последовательные отсчеты функций в этих зависимостях).
\item
Во-вторых, в базовом выражении $b=F\lambda+\eta$ для задачи метода обыкновенных наименьших квадратов (МОНК --- Ordinary Least Squares), невязка $\eta$ не должна быть функцией от модельного параметра $\lambda$, но здесь она каждый раз сложным образом зависит от  $\lambda$: $\eta= P_{\lambda}\varepsilon$, 
где $P_{\lambda}$ --- $N\times N$-матрица, зависящая от  $\lambda$, а 
$\varepsilon$ --- вектор значений $\varepsilon_k$.
Для  преодоления этого препятствия автор разработал оригинальный и очень эффективный численный метод. Метод заключается в \emph{итеративном}  применении МОНК: 
на каждой $n$-й итерации для вычисления очередной оценки $\hat\lambda^{(n)}$ используются предыдущее МНК-решение $\hat\lambda^{(n-1)}$.
\item
В-третьих,  матрица $F$ плана эксперимента оказывается зависима от экспериментальных значений $y_k$, \ie фактически от погрешностей $\varepsilon_k$. (Выражения для этих зависимостей систематизированы в 
табл.~3.1--3.4 на с.~132--134 диссертации и далее раскрываются по тексту в многочисленных конкретных примерах).
Это тоже серьезное препятствие для применения МОНК. Строго, традиционный МОНК применить нельзя. 
Однако автор доказывает, что и в этом случае предложенный им итерационный метод  МОНК \emph{сходится}   за небольшое число (единицы) итераций достаточно точно к неизвестному значению вектор-параметра $\lambda$, --- без применения более сложного метода полных наименьших квадратов (Total Least Squares). 
Эти доказательства даны в форме теорем 3.1--3.2 и следствий 3.1--3.4 на с.~179--184, дополнены явным вычислением матриц, необходимых для реализации этого численного метода --- на с.~185--193, и результатами численно-аналитического исследования  на с.~193--205, 
показывающими уменьшение ошибок смещения в оценках параметров $\lambda$ модели на 2--3 порядка по сравнению с традиционной (неитерационной) процедурой МОНК, --- табл.~3.11--3.21 и рис.~3.9--3.21.
\end{itemize}

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{О/Д} %{Обоснованность}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[2.] Обоснованность и достоверность
\end{description}
%\textbf{3.}
\emph{Насколько детально автор обосновывает  свои главные научные результаты и насколько эффективно он их  излагает?}

Автор диссертации применяет единую многошаговую схему получения и демонстрации своих результатов:
\begin{enumerate}
\item 
Дает мотивировку --- отмечает важность решения рассматриваемой задачи.
\item 
Проводит предварительный анализ  существующих методов решения и отмечает найденные недостатки.
\item 
Строит ЛПДМ для рассматриваемой задачи --- демонстрирует всю технику получения рекуррентных соотношений, связывающих последовательные отсчеты  $\tilde y_k$ теоретической (модельной) зависимости $\tilde y(t)$.
\item 
Заменяет в этих ЛПДМ  модельные отсчеты $\tilde y_k$ экспериментальными отсчетами $y_k$, получаемыми от реальной ДМС, по формуле $y_k=\tilde y_k+\varepsilon_k$, тем самым вводя неизбежные расхождение $\varepsilon_k$ между $y_k$ и $\tilde y_k$ как вектор $\varepsilon = (\varepsilon_1 \bigm| \ldots  \bigm| \varepsilon_N)^{T}$.
\item 
Представляет построенную таким образом ЛПДМ в стандартной форме
\begin{eqnarray}
b &=& F\lambda+\eta \label{e:ls1}\\
\eta &=& P_{\lambda}\varepsilon, \quad \det{P_{\lambda}}\ne 0 \label{e:ls2}
\end{eqnarray}
алгебраической задачи решения переопределенной несовместной СЛАУ --- системы линейных алгебраических уравнений $b = F\lambda$ в смысле критерия МОНК
\begin{equation}\label{e:ls3}
\| \varepsilon  \|^2= \varepsilon^T\varepsilon =  \| P_\lambda^{-1}b-P_\lambda^{-1}F\lambda  \|^2 
{\Longrightarrow} \min\limits_\lambda
\end{equation}
\ie численного отыскания $\hat\lambda= \arg\min\limits_\lambda \| \varepsilon  \|^2$.
\item
Вводит статистическую интерпретацию модели (\ref{e:ls1})--(\ref{e:ls2}), трактуя невязки $\varepsilon$ как случайные векторы с нормальным распределением и двумя моментами
\[
\E{\varepsilon}=\boldsymbol{0}, \qquad 
\E{\varepsilon\varepsilon^T}=I\sigma_{\varepsilon}^2.
\]
\item
Конкретизирует вид и наполнение всех матриц, входящих в задачу 
(\ref{e:ls1}), (\ref{e:ls2}) и (\ref{e:ls3}), для каждой из многочисленных подзадач, рассмотренных в рамках данного единого подхода.
\item
Конкретизирует все детали разработанного итерационного численного метода, решающего поэтапно последовательность задач МОНК.
\item
Реализует этот метод в виде компьютерного эксперимента, чтобы получить конкретные численные результаты --- оценки $\hat\lambda$. 
Для этого автор планирует эксперимент в пакете прикладных программ, созданном в рамках данного диссертационного исследования в среде Microsoft Visual Basic 6.0.
\item
Подтверждает, что результаты применения его метода дают значительно более точные оценки по сравнению с другими методами. Для сравнения автор проводит (для тех же условий) вычислительный эксперимент по известным методам, заложенным в приложение Microsoft Excel системы MS Office XP.
\item
Вычисляет --- на основании оценок $\hat\lambda$, полученных от своего оригинального метода, --- оценки технических характеристик (собственной частоты, декремента затуханий и \tp) исследуемой ДМС, пользуясь выведенными заранее аналитическими зависимостями этих характеристик от точных $\lambda$.
\item
Оценивает погрешности проведенных таким образом вычислений, сравнивая предлагаемый метод и известные методы по достигаемой точности.
\item
Завершает обоснование своего вклада в решение задачи выводами о преимуществах разработанного численного метода и тем подтверждает ценность данного научного подхода к построению и практическому использованию 
линейно-параметрических дискретных моделей в форме
разностных уравнений в задачах идентификации
диссипативных механических систем.
\end{enumerate}
%\\[-1.5ex]

Такое детальное обоснование дается автором для каждого конкретного применения предлагаемого подхода, рассмотренного в работе. 
Но и в целом вся диссертация выстроена в соответствии с требованиями, предъявляемыми к научным текстам:
\begin{itemize}
\item 
Глава 1 формулирует проблему повышения точности, быстродействия и расширения функциональных возможностей методов определения характеристик диссипативных механических систем. Этим продиктована внутренняя логика представления всего материала исследования, его целей и задач в тексте диссертации.
\item 
Глава 2 разрабатывает математическое описание диссипативных механических систем в форме линейно-параметрических дискретных моделей для ДМС различного типа и при различных типовых внешних воздействиях. Разработаны  основополагающие принципы построения  линейно-параметрических дискретных моделей и даны оценки точности строящихся ЛПДМ.
\item 
Глава 3 разрабатывает и исследует свойства нового численного метода определения характеристик ДМС на основе введенных ЛПДМ. 
Математически строго обоснована сходимость предлагаемого метода.
\item 
Глава 4 занимается исследованием --- повышением точности и численной устойчивости предложенного численного метода. 
Учитываются погрешности разного происхождения: методические и инструментальные (от квантования измерений по времени, от округления при вычислениях, подробно --- в табл.~4.1, c.~213--215). 
Для улучшения численной обусловленности СЛАУ (нормальной системы уравнений 
$F^TF\hat\lambda=F^Tb$) разработан (п.~4.3, cc.~236--247) структурный метод 
<<прореживания>> отсчетов экспериментальных наблюдений, доказана теорема~4.2
об оптимальном значении $l_0=\left[ \frac{T}{4\tau} \right]$ параметра прореживания $l$ в формуле $\tau_l=l\tau$, гарантирующем наибольшее  уменьшение числа обусловленности матрицы  $F$, а значит и $F^TF$ (для случая модели (\ref{e:lin-vyaz})), где $\tau_l$ --- период взятия отсчетов измерений после прореживания и $\tau$ ---  этот период до прореживания, $T$ --- период колебаний в ДМС.
.
\item 
Глава 5 изучает вопросы построения линейно-параметрических дискретных моделей и применения в них разработанного численного метода идентификации для задач из различных областей науки и техники. Обсуждаются: 
задача неупругого реологического деформирования материалов и элементов конструкций, 
задача определения параметров передаточной функции систем автоматического управления, 
задача определения некоторых параметров САУ по амплитудно-частотной характеристике,
задача определения демпфирующих свойств ДМС с различными диссипативными силами и
результаты практического применения разработанных методов в научно-технических проектах --- при оценке технического состояния силовых элементов шасси
самолета, в цифровом осциллографировании (упоминается использование ЛПДМ в серийных приборах) и при обнаружении некачественной сборки деталей прессованием (описана работа в пакете прикладных программ, созданном в среде Delphi). 
\item 
Глава 6 завершает основную часть работы представлением практических разработок автора --- 
программного обеспечения, реализованного в среде MS Visual Basic 6.0, и пяти специализированных устройств определения динамических характеристик ДМС на основе ЛПДМ.
\end{itemize}

Таким образом, проведенный анализ диссертационной работы В.~Е.~Зотеева, автореферата этой работы и имеющихся научных публикаций автора по теме диссертации позволяет мне
с уверенностью констатировать следующее:
\begin{itemize}
\item 
\textbf{Актуальность избранной темы} несомненна. Она диктуется необходимостью надежного технического диагностирования сложных объектов машиностроения, в том числе, и в целях безопасности, с применением современных компьютерных средств.
\item 
\textbf{Обоснованность и достоверность} положений, выносимых на защиту, обеспечены, а именно: 

\emph{Теоретические основы и принципы построения ЛПДМ} даны в полном объеме.

\emph{Линейно-параметрические дискретные модели} построены для ДМС с различными показателями нелинейности диссипативных сил.

\emph{Новые структурные соотношения} во временн\'ой области между отсчетами  колебаний в ДМС и динамическими характеристиками выведены математически строго.

\emph{Численный метод} определения характеристик ДМС на основе ЛПДМ построен, проанализирован теоретически и испытан в вычислительных экспериментах многократно.

\emph{Структурные методы повышения численной устойчивости} данного метода доказаны теоретически и проверены практически.

\emph{Численный метод определения параметров кривой ползучести} на основе ЛПДМ, построен и испытан на реальных деформируемых материалах.

\emph{Линейно-параметрические дискретные модели для огибающей амплитуд колебаний} в ДМС и численный метод определения диссипативных характеристик ДМС на их основе математически построены и экспериментально подтверждены

\end{itemize}
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{Н}          %{Новизна}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>

\begin{description}
\item[3.] Новизна
\end{description}
\textbf{Новизна} проведенных исследований и полученных результатов (на момент их опубликования) имеется. Оригинальными (впервые полученными  в данной работе) являются следующие результаты:
\begin{itemize}

\item \emph{новый научный подход} к решению задачи определения характеристик ДМС, который заключается в совокупности десяти этапов, характеризованных выше в пп.~3--12 на с.~4--5 этого отзыва,

\item \emph{основы теории и техники построения ЛПДМ} --- новых рекуррентных соотношений, связывающих последовательные отсчеты  $\tilde y_k$ теоретической (модельной) зависимости $\tilde y(t)$.

\item \emph{ЛПДМ для различных типов нелинейности диссипативной силы}, новые как по структуре, так и по найденным связям их параметров с динамическими характеристиками исследуемой системы,

\item \emph{итерационный численный метод}, решающий поэтапно последовательность задач МОНК, с доказательством сходимости метода (новое применение известного, классического метода обыкновенных наименьших квадратов),

\item \emph{новые аналитические соотношения}, выражающие оценку погрешности результатов определения динамических характеристик исследуемой системы через финальные МНК-оценки параметров ЛПДМ,

\item \emph{новые структурные методы} повышения численной устойчивости МНК-оценок параметров ЛПДМ при плохой обусловленности задачи МНК,

\item \emph{новые ЛПДМ и численные алгоритмы} определения параметров кривой ползучести при деформировании материалов и элементов конструкций,

\item \emph{новые ЛПДМ для огибающей колебаний ДМС} с диссипативными силами общего вида --- пропорциональными $n$-й степени скорости движения ($n=0, 1, 2, \ldots $),

\item \emph{новый численный алгоритм} определения диссипативных характеристик
систем типа (\ref{e:lin-vyaz}) по ЛПДМ, построенной для амплитудно-частотной характеристики ДМС,

\item \emph{новые специализированные устройства} для измерения различных диссипативных характеристик --- декремента затухания и \dr

\end{itemize}

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{З/С} %{Значимость/Соот}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[4.] Значимость
\end{description}
\begin{itemize}
\item 
\textbf{Значимость результатов работы} для науки и практики. 
Отмеченный выше \emph{новый научный подход} применим для широкого класса ДМС, но его значение для науки не ограничивается этим классом. 
Подход принципиально применим для параметрической идентификации характеристик систем различной природы: электротехнических, биологических, химико-технологических, экономических и \dr 
\item 
\emph{Методологическая база}, созданная автором, предоставляет исследователям и специ\-алистам-практикам, занимающимся определением динамических характеристик процессов, практический инструментарий и руководство по внедрению этих методов в новые аппаратные или программные средства.
\end{itemize}

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{С} %{Соответствие}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[5.] Соответствие критериям
\end{description}
\begin{itemize}
\item 
\textbf{Соответствие диссертации критериям}, установленным Положением о порядке присуждения ученых степеней (Постановление Правительства РФ от 30.01.2002 №74). Диссертация соответствует этим критериям в полной мере.

\item 
Соблюдены необходимые \textbf{принципы соответствия}: соответствие целей и задач, соответствие --- автореферата и диссертации, содержания диссертации и содержания опубликованных работ, а также темы диссертации и научной специальности.

\item 
\textbf{Строгость авторской аргументации} в тексте автореферата, в тексте публикаций и в самой диссертации, академический стиль изложения, отличное полиграфическое качество оформления, практическое отсутствие опечаток или грамматических ошибок свидетельствуют о высокой требовательности к себе и о высокой культуре автора как научного работника.

\item 
\textbf{Автореферат диссертации соответствует ее тексту}. Несоответствие обнаружено в данных по публикациям автора по теме данного исследования, а именно:
\begin{center}
\begin{tabular}{c@{\quad}|c@{\quad}|c@{\quad}|c}
 {\small Труды} & {\small по списку литературы}  & {\small по списку трудов} & {\small стр.~16 диссертации и}  \\
  &{\small в диссертации} & {\small в автореферате} &{\small стр.~11 автореферата}  \\\hline
 {\small статьи (ВАК)}&36 &36 &36 \\
 {\small в трудах конф.}&19 &26 &25 \\
 {\small в сб. трудов}&7 &7 &46 \\
 {\small монографии}&1 &1 &1 \\
 {\small АС СССР}&5 &5 &5 \\\hline
 {\small Всего}&68 &75 &113 \\\hline
\end{tabular}
\end{center}
\end{itemize}
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{О/И}%{Опубликование/Исп}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[6.] Опубликование / Использование результатов
\end{description}
\begin{itemize}
\item 
\textbf{Основные результаты диссертационного исследования опубликованы в печати.} Отмеченное выше несоответствие количественных данных не является существенным: научные положения диссертации полностью отражены в  ведущих  рецензируемых научных журналах и изданиях, утвержденных ВАК.

\item 
\textbf{Сведения о практическом использовании результатов диссертации} приведены в достаточном объеме --- пять актов о внедрении.  Дальнейшее внедрение / использование результатов этой работы целесообразно проводить в организациях и учреждениях, занимающихся разработкой и исследованием сложных изделий отечественного машиностроения. В их числе: Самарский государственный технический университет, Самарский государственный аэрокосмический университет 
им. акад. С.~П.~Королева, Санкт-Петербургский государственный университет аэрокосмического приборостроения, Ульяновский государственный технический университет, Московский государственный технический университет имени Н.~Э.~Баумана, Санкт-Петербургский государственный политехнический университет, Государственный технологический университет Московский институт стали и сплавов и \dr

\item 
В работе явно прослеживается \textbf{внутреннее единство} новых научных результатов, которое свидетельствует о том, что все результаты работы получены автором лично.

\end{itemize}

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{Критика}    %{Критика}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[7.] Критика
\end{description}
Вместе с отмеченными достоинствами этой солидной работы нельзя не отметить и \textbf{некоторые недостатки}, не являющиеся, впрочем, существенными для общей положительной оценки данного исследования. 
По существу, заслуживает внимания лишь один недостаток, и он заключается в не вполне верной трактовке задачи (\ref{e:ls1}), (\ref{e:ls2}) и (\ref{e:ls3}) метода наименьших квадратов. Как известно, этот метод как алгебраическая процедура отыскания  $\hat\lambda= \arg\min\limits_\lambda \| \varepsilon  \|^2$ был предложен Лежандром в 1805 году и позднее подтвержден как статистическая процедура Гауссом в 1809 году (хотя рукопись Гаусса уже была доступна в 1806 году). Гаусс в работе 1809 года объявил, что он уже использовал МНК как алгебраическую процедуре с 1795 года, чем вызвал большое раздражение Лежандра\footnote{см. Stigler, S.~M. (1977), An Attack on Gauss, Published by Legendre in 1820, \emph{Historia Mathematica, 4}, pp.~31--35.}
Различие между нестатистической (чисто алгебраической) и статистической задачами иллюстрируется \ris\ref{f:modest}. 
Применительно к формулам  (\ref{e:ls1}), (\ref{e:ls2}) и (\ref{e:ls3}) на 
\ris\ref{f:modest} следует считать, что $z=P_{\lambda}^{-1}b$ и $A=P_{\lambda}^{-1}F$, тогда  $z=A\lambda+\varepsilon$. В диссертационной работе решается нестатистическая задача подгонки модели под реальную систему
в контексте \ris\ref{f:modest}$(a)$. Решение $\hat\lambda= \arg\min\limits_\lambda \| \varepsilon  \|^2= (A^TA)^{-1}A^Tz$. 
После его отыскания оно интерпретируется статистически. Для этого невязка
$\varepsilon$ трактуется как случайный гауссов вектор с характеристиками   в контексте \ris\ref{f:modest}$(b)$. 
Однако это не означает, что автор тем самым превращает разностные уравнения ЛПДМ в стохастические разностные уравнения. 
Он ведь не получает (да и не обсуждает) вероятностное выборочное пространство, на котором были бы заданы реализации (ансамбль выборочных функций) решения этих уравнений и соответствующая вероятностная мера. 
Поэтому фразы о ЛПДМ в форме стохастических уравнений включены неоправданно.
Дается простая статистическая интерпретация алгебраического  решения переопределенной СЛАУ полного столбцового ранга. 
Слова о статистической обработке включены также неоправданно, поскольку обработка производится по \emph{одной} реализации экспериментальных отсчетов процесса, а не по многим реализациям (выборкам) из генеральной совокупности (ансамбля) реализаций.
Термины <<выборка, объем выборки $N$>> в работе имеют другой смысл: $N$ означает число отсчетов, взятых от одной реализации процесса.

Статистическая интерпретация решения нормальных уравнений по  \ris\ref{f:modest} 
в целом должна включать следующие пункты:
\begin{enumerate}
\item 
Критерий 
$\min\limits_{\lambda} (z-A\lambda)^T(z-A\lambda)$ не является статистическим критерием. Решение, отвечающее ему, удовлетворяет нормальным уравнениям
$(A^TA)\hat \lambda= A^Tz$ и, когда $\rank A=n$, единственно: 
$\hat\lambda=(A^TA)^{-1}A^Tz$.
\item Статистическая интерпретация этого решения (\ris\ref{f:modest}$(b)$) имеет следующий смысл:

1) $\lambda$ есть вектор (случайный или нет --- это зависит от полноты интерпретации), который надо оценивать по наблюдениям  $z=A\lambda+\varepsilon$,

2) $A$ есть матрица наблюдений (в регрессионном анализе --- матрица регрессоров, в теории планирования эксперимента --- матрица плана эксперимента), которая показывает, какие элементы неизвестного вектора 
$\lambda$ и в каких комбинациях включены в вектор наблюдаемых значений --- в вектор $z$,

3)  $z$ есть значение измеряемого вектора $A\lambda$, искаженное присутствием аддитивной случайно помехи $\varepsilon$, 

4) $\varepsilon$ есть случайный вектор погрешностей измерений,

5) $\Lambda =A^TA$ есть информационная матрица решения $\hat\lambda=(A^TA)^{-1}A^Tz$ при данных определениях, а его ковариационная матрица $P$ равна $\Lambda^{-1}$.

\item При условии гауссовых распределений оценка максимума правдоподобия вектора $\lambda$, как и оценка по среднеквадратическому критерию (\ris\ref{f:modest}(b)), определяется той же формулой $\hat\lambda=(A^TA)^{-1}A^Tz$.

\item При этом условии оценка максимума апостериорной вероятности параметра $\lambda$ дается формулами фильтра Калмана: 
$\hat\lambda=\bar\lambda +K(z-A\bar\lambda)$, $K=\bar P A^T(A\bar P A^T +R)^{-1}$, где $\bar\lambda$ -- безусловное (априорное) среднее значение случайного  $\lambda$, а  $\bar P$ -- его безусловная (априорная) ковариация,
$R$ --- ковариация вектора $\varepsilon$.
\end{enumerate}

Последний из отмеченных пунктов указывает на то, что переопределенную СЛАУ, получаемую в этой работе, можно численно решать более эффективно, 
чем отмечается в работе (с разложением Холецкого или ортогональным разложением матрицы $\Lambda =A^TA$) на с.~211. 
А именно: <<расщепляя>> большую СЛАУ на априорную и текущую порции уравнений, можно решать систему последовательно по алгоритму калмановской фильтрации, беря, правда, его не в том исходном виде, как показано выше, а используя современные факторизованные и совмещенные алгоритмы, обладающие высокой численной устойчивостью. Однако, это замечание можно отнести к пожеланиям для продолжения работ в этом направлении.

\begin{figure}%[h]
\input{modest.txt}
\caption{Алгебраически эквивалентные задачи, решаемые по методу наименьших квадратов значений невязки $\varepsilon$ или по критерию минимума среднего квадрата погрешности  $e$: ($a$) --- оптимальное моделирование неизвестной системы по экспериментальным условиям $A$ и данным $z$; ($b$) --- оптимальное оценивание неизвестного вектора по наблюдениям вектора $A\lambda$ в присутствии аддитивных случайных помех $\varepsilon$ с характеристиками $\E{\varepsilon}=0$  и $\E{\varepsilon \varepsilon^T}=I$}\label{f:modest}
\end{figure}

%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
%\section{Заключение}   %{Вывод}%%%%%%%=======>>%%%%%%%=====>>%%%%%======>>
%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>%%%%%%%=======>>
\begin{description}
\item[8.] Заключение
\end{description}
Таким образом, диссертационная работа В.~Е.~Зотеева дает эффективное решение крупной научной проблемы, имеющей важное народнохозяйственное значение, --- 
\emph{проблемы идентификации нелинейных диссипативных механических систем в процессе их эксплуатации или прочностных промышленных испытаний}.
Работа содержит новые научные результаты, имеет общетеоретическую значимость для науки и представляет практическую ценность для специалистов в области конструирования, производства, испытаний и эксплуатации сложных изделий отечественного машиностроения. 
Работа полностью отвечает требованиям, предъявляемым ВАК к диссертациям на соискание ученой степени доктора технических наук по специальности 05.13.18 --- <<Математическое моделирование, численные методы и комплексы программ>>.\\

\vspace{0.5cm}

\noindent\begin{tabular}{p{13.3cm}}
%\noindent\includegraphics{my2}
\end{tabular}

\noindent И.~В. Семушин, д-р техн. наук\\
профессор\\



\noindent Подпись И.~В. Семушина заверяю:\\

\hbox to \textwidth{Секретарь Ученого Совета\hrulefill О. А. Литвинко}


\end{document}