Развитие и применение методов регуляризации для обработки экспериментальных данных мёссбауэровской спектроскопии тема автореферата и диссертации по физике, 01.04.01 ВАК РФ

Введение.

Глава 1. Методы восстановления функции плотности вероятности распределения параметров сверхтонкого взаимодействия по экспериментальным данным мёссбауэровской спектроскопии.

1.1. Математическое описание мёссбауэровских спектров и определение функции распределения параметров сверхтонкого взаимодействия.

1.2. Методы определения функции распределения.

- параметризации (Винц, Варрет, Филлипс, Шарон и др.),

- Фурье разложение с фильтрацией (Виндоу, Винц, Дибар, Надь, Стоун, Келлер и др.),

- функций штрафа (Бранд, JTe Кер, Вивел, Тумэй и др.),

- реставрации изображения (Русаков, Николаев и др.),

- линеаризации (Соколов, Силин, Афанасьев, Jle Кер и др.),

- максимума энтропии или статистической регуляризации (Бранд, Ле Кер, Жуковский, Ранкор и др.),

- регуляризации (Хессе, Рубарч, Вивел, Белозерский и др.)

1.3. Обобщенный алгоритм восстановления функции распределения на основе метода регуляризации Тихонова.

1.4. Проблемы существующих методов восстановления функции распределения и перспективы решения.

1.5. Выводы и постановка задачи исследований.

Глава 2. Регулярный алгоритм решения обратной задачи мёссбауэровской спектроскопии.

2.1. Анализ задачи на разрешимость.

2.1.1. Существование и единственность решения.

2.1.2. Практическое значение.

2.2. Модифицированный метод коррекции параметров ядра интегрального уравнения.

2.2.1. Теоретическое обоснование модификации метода.

2.2.2. Алгоритм модифицированного двухэтапного процесса.

2.2.3. Условия и теорема сходимости метода.

2.2.4. Точность определения параметров.

2.2.5. Примеры использования модифицированного метода.

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

2.4. Выводы.

Глава 3. Оценка погрешности функции распределения в равномерной метрике.

3.1. Методика оценки погрешности решения.

3.1.1. Моделирование ошибки.

3.1.2. Оценка статистической ошибки с помощью сингулярного разложения

3.2. Оценка погрешности метода.

3.3. Практическое значение оценки погрешности функции распределения:.

- интерпретация результатов,

- выбор параметра регуляризации,

- рекомендации по проведению эксперимента

3.4. Выводы.

Глава 4. Обработка мёссбауэровских спектров локально неоднородных систем.

4.1. Восстановление заданных на разных интервалах нескольких функций распределения параметров СТВ.

4.2. Модификация обобщенного регулярного алгоритма с целью восстановления нескольких распределений и оценка их погрешности.

4.3. Практическое значение.

4.4. Выводы.

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

В формировании МС принимает участие большое число резонансных ядер атомов, в общем случае не являющихся химически и кристаллографически эквивалентными. Неэквивалентные позиции резонансных ядер различаются параметрами сверхтонкого взаимодействия (СТВ). В случае небольшого числа неэквивалентных позиций и сильно различающихся параметров для описания экспериментального спектра используется дискретное представление. Задача нахождения неизвестных мёссбауэровских параметров для каждой позиции и вероятностей этих позиций решается моделированием формы огибающей спектра и подгонки ее к экспериментальной кривой [2].

В случае неразрешенных и слаборазрешенных спектров трудно предложить обоснованную модель гамма-резонансного спектра, либо отдать предпочтение той или иной модели. Поэтому для них используют непрерывное описание, основанное на линейной суперпозиции непрерывно распределенных парциальных составляющих спектра, различающихся значением одного доминирующего параметра СТВ. В этом случае задача расшифровки МС состоит в определении функции плотности вероятности распределения доминирующего параметра СТВ по экспериментальным данным y(v) и описывается интегральным уравнением Фредгольма 1 рода: ь

Ар(х) = Jk(x,v)p(x)dx = y(v), (0.1) а где y(v) - интенсивность резонансного поглощения, как функция относительной скорости v; К(х, v) - функция, задающая форму элементарных составляющих; р(х) - функция плотности вероятности распределения искомого параметра х; [а, Ь] - интервал существования непрерывного распределения. Такие задачи относятся к классу некорректных задач и характеризуются тем, что сколь угодно малые изменения исходных данных могут приводить к большим изменениям решений [3]. Поскольку вместо точной правой части y°(v) уравнения (0.1) известно приближенное значение y(v), найденное из эксперимента, то для получения устойчивого решения необходимо использовать методы регуляризации [4-7].

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

1) р(х) - непрерывно дифференцируемая функция, что соответствует непрерывному изменению значений искомого параметра СТВ;

2) р(а)=0 или — (а) = 0, р(Ь)=0 или — (Ь) = 0 - граничные условия; dx dx xmax

3) j"p(x)dx = 1 - нормированность; (0.2) xmin

4) p(x)>0 на [a, b] - неотрицательность.

Первые попытки восстановления функции распределения параметров СТВ из МС были предприняты в начале 60-х годов Филлипсом и Тумэйем [8, 9]. Искомую функцию представляли в виде конечного набора дискретных значений или в виде непрерывной функции распределения одного из параметров СТВ, в качестве которого обычно использовались сверхтонкое магнитное поле Н, изомерный сдвиг 5СД или квадрупольное расщепление Q. Последующее развитие методов определения функции распределения параметров СТВ из МС шло по двум направлениям:

1. Методы, основанные на упрощенных моделях функции распределения параметров СТВ, когда эта функция параметризировалась с помощью конечного числа параметров, а затем использовался метод наименьших квадратов [10-18];

2. Методы определения функции распределения параметров СТВ, как регуляризован-ного решения обратной задачи мёссбауэровской спектроскопии [19-28].

Современные методы обработки экспериментальных данных мёссбауэровской спектроскопии, разработанные с использованием теории решения некорректных задач [4-7], позволяют математически корректно восстанавливать функцию распределения параметров СТВ [28-39]. В них реализована возможность учета дополнительной информации о свойствах искомой функции и тем самым получение решения, имеющего физический смысл. Современными методами можно обрабатывать широкий класс спектров разнообразных материалов и систем. Однако опыт использования методов обработки экспериментальных данных мёссбауэровской спектроскопии в последние годы и анализ их возможностей показали несовершенство существующих алгоритмов и необходимость дальнейших исследований.

Самый главный недостаток существующих алгоритмов решения обратной задачи мёссбауэровской спектроскопии состоит в том, что в них либо не определяется погрешность восстанавливаемых параметров, либо указываются точечные оценки для функции распределения параметров СТВ, параметризованной конечным числом подгоночных параметров. Очевидно, что без знания «коридора ошибок» решения, т.е. без оценки погрешности решения в равномерной метрике, трудно говорить о достоверности как получаемых результатов, так и их физической интерпретации. В литературе отсутствуют работы, в которых была бы изложена методика определения полной погрешности (статистическая ошибка и погрешность метода) функции плотности вероятности распределения параметров СТВ в равномерной метрике. Проблема определения погрешности функции распределения параметров СТВ связана с неустойчивостью решения обратной задачи мёссбауэровской спектроскопии. Классические методы получения статистических оценок погрешности искомых параметров, такие как метод наименьших квадратов и метод максимального правдоподобия, неприемлемы для данных задач ввиду их неустойчивости: при стремлении шага дискретизации к нулю статистическая оценка погрешности стремится к бесконечности. Методы регуляризации в свою очередь дают смещенную оценку статистической ошибки решения [40], т.е. математическое ожидание ошибки решения не равно нулю и, следовательно, необходимо учитывать погрешность метода регуляризации.

Другой вопрос, имеющий большое значение в математической обработке и не освещенный в методах решения обратной задачи мёссбауэровской спектроскопии - это возможность определения искомой функции распределения параметров СТВ из экспериментальных данных в рамках выбранной физической модели и ее математического описания, т.е. анализ задачи на разрешимость. Исходные данные в эксперименте всегда известны с погрешностью, поэтому неустойчивость некорректных задач приводит к практической (не)единственности решения и необходимости сформулировать принципы отбора искомого решения среди множества возможных решений.

В ситуации, когда погрешность уравнения от неточного задания параметров ядра K(x,v) интегрального уравнения (0.1) превышает статистическую погрешность экспериментальных данных, восстанавливаемые функции распределения параметров СТВ могут значительно отличаться от истинных [41-44]. Это связано с тем, что значения параметров, входящих в ядро интегрального уравнения, такие как изомерный сдвиг 5СД, квадрупольное расщепление Q и т.д. либо неизвестны, либо часть параметров известна точно, а остальные — приближенно. В результате, необходимо решать обратную задачу мёссбауэровской спектроскопии в условиях приближенно заданного оператора. В частности, в работе [41], было предложено два алгоритма коррекции параметров ядра интегрального уравнения (0.1), основанных на методе Гаусса-Ньютона - двухэтажный итерационный процесс и процесс с проектированием. Эффективность этих методов неоднократно демонстрировалась при обработке МС [4144]. Однако оба метода имеют свои недостатки: ограничение скорости сходимости двухэтап-ного итерационного процесса сверху некоторым значением и необходимость вычисления проекционного оператора на каждом шаге итерации, соответственно. Кроме того, в обоих методах не рассматривается вопрос точности определения значений параметров. Между тем, существуют теоретические предпосылки [45, 46] создания метода коррекции параметров, имеющего большую скорость сходимости и более экономичного.

Для преобразования решения обратной задачи мёссбауэровской спектроскопии с интервалами отрицательных значений в неотрицательное (ограничение 4), (0.2)) традиционно используются различные итерационные процедуры [19, 23, 24, 30-32], которые не всегда приводят к удовлетворительному результату. Функции распределения параметров СТВ, имеющие незначительную по площади долю отрицательных значений, после применения этих алгоритмов иногда могут восстанавливать экспериментальные данные в пределах требуемой точности. В тех случаях, когда функция имеет участки относительно больших отрицательных значений или протяженные интервалы осциллирующих около нуля значений, применение этих итерационных алгоритмов приводит к существенным искажениям восстанавливаемых МС. Необходим эффективный алгоритм получения неотрицательного решения, преобразующий функцию распределения параметров СТВ адаптивно физическим свойствам исследуемого образца и позволяющий восстанавливать спектры в пределах требуемой точности.

В последнее время большой научный интерес направлен на исследования веществ с локальной неоднородностью, которая характеризуется изменением от позиции к позиции окружения атомов одного сорта и, как следствие, параметров СТВ его ядра [34, 35]. Мес-сбауэровские спектры локально неоднородных систем представляют собой совокупность большого числа парциальных спектров и их математическая обработка вызывает особые трудности. Восстановление из таких спектров функции распределения параметров СТВ методами, не созданными для обработки спектров этих систем, не учитывающими специфику задачи, не гарантирует получение физически значимого результата. Например, когда в спектре присутствуют две составляющие с узкими и широкими линиями, которым соответствуют узкое и широкое распределения параметра СТВ, аппроксимация МС одним распределением с одинаковой гладкостью на всем интервале определения параметра может привести к неадекватному описанию МС. Иногда из физической модели следует, что функция распределения параметров СТВ состоит из двух частей, расположенных на отдаленных друг от друга интервалах значений искомого параметра. Восстановление функции на всем интервале определения параметра приводит к физически необоснованным составляющим мёссбауэровского спектра. Таким образом, для обработки спектров локально неоднородных систем необходим физически и математически адаптированный подход.

Одним из современных методов обработки МС, наиболее подробно представленных в литературе и имеющий полное математическое обоснование, является обобщенный регулярный алгоритм восстановления функции распределения параметров СТВ [30, 39, 41-44]. Основанный на вариационном методе Тихонова, обобщенный регулярный алгоритм обеспечивает получение устойчивого решения, удовлетворяющего всем априорным ограничениям (0.2). Он успешно применяется для обработки спектров, характерных как для упорядоченных систем, так и для локально неоднородных. Несмотря на имеющиеся недостатки данный алгоритм представляется наиболее подходящим в качестве базового для развития методов обработки экспериментальных данных мёссбауэровской спектроскопии.

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

В соответствии с поставленной целью решались следующие задачи:

1. Исследование обратной задачи мёссбауэровской спектроскопии на существование и единственность решения.

2. Разработка методики оценки погрешности функции распределения параметров СТВ в равномерной метрике, т.е. получение «коридора ошибок».

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

4. Разработка итерационного алгоритма определения неотрицательной функции распределения параметров СТВ, восстанавливающей МС в пределах требуемой точности.

5. Развитие обобщенного регулярного алгоритма с целью обработки МС локально неоднородных систем, восстановление которых предполагает использование нескольких функций распределения параметров СТВ.

Работа выполнена в Физико-техническом институте Уральского Отделения РАН в соответствии с планами научно-исследовательских работ по темам: «Исследование микроскопических механизмов и кинетики образования метастабильных фаз и нанокристаллических разупорядоченных структур на основе железа при механическом сплавлении» (№ гос. регистрации 01.9.90 002472), «Исследование связи микро- и макроскопических магнитных свойств с локальной атомной структурой, формой и размером частиц в деформированных, компактированных и отожженных сплавах на основе железа» (№ гос. регистрации 01.9.90 002471) и гранта РФФИ 97-01-00520 «Некорректные задачи на классах негладких функций: алгоритмы и приложения».

Научная новизна работы. Для обеспечения достоверных результатов обработки МС впервые разработана методика оценки погрешности функции распределения параметров СТВ, полученной вариационным методом Тихонова, в равномерной метрике.

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

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

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

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

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

Положения, выносимые на защиту:

Схема проверки обратной задачи мёссбауэровской спектроскопии на существование и единственность решения.

Модифицированный двухэтапный метод коррекции параметров ядра интегрального уравнения 1 рода с конечномерной нелинейностью.

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

Методика оценки погрешности функции распределения параметров СТВ, полученной вариационным методом Тихонова, в равномерной метрике и определение параметра регуляризации.

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

Личный вклад автора. Диссертация является самостоятельной работой, обобщившей результаты, полученные лично автором и в соавторстве. Проведено обоснование модификации двухэтапного метода коррекции параметров. Эффективность метода продемонстрирована на примерах обработки модельных и экспериментальных МС, которые выполнены лично автором. Автором совместно с соавторами разработана методика равномерной оценки погрешности решения обратной задачи мёссбауэровской спектроскопии, которая апробирована автором на модельных и реальных задачах. Целесообразность и эффективность предложенного автором итерационного алгоритма получения неотрицательного решения многократно проверены самим автором на примерах функций распределений параметров СТВ, имеющих разнообразные (по протяженности, интенсивности, повторяемости и т.п.) области отрицательных значений. Автором были проанализированы проблемы восстановления функции распределения параметров СТВ из МС, парциальные составляющие которых имеют существенно отличные значения параметров, и предложен способ их решения. Лично автором в рамках научно-исследовательских работ лаборатории выполнена полная математическая обработка МС, полученных для:

• неупорядоченных кристаллических сплавов Fe+15,25, 33at.% Si;

• серии механически сплавленных в течение 2ч., 4ч., 8ч. и 16ч. при ТШМ=300К смеси порошков Fe(75)C(25);

• серии смеси порошков Fe(68)C(32) после механического сплавления в течение 1ч., 2ч., 4ч., 8ч., 12ч. при ТИЗМ=77К;

• выплавленного сплава Fe+8aT.% Sn после механического измельчения при ТИЗМ=300К;

• серии неупорядоченного кристаллического сплава Fe75Si25 после термической обработки.

Идеи и способы развития регулярных методов обработки экспериментальных данных по диссертационной работе сформулированы научными руководителями Агеевым A.JI. и Ворониной Е.В. Основные положения и выводы диссертационной работы сформулированы автором.

Апробация работы и публикации. Результаты работы докладывались и обсуждались на семинарах и конференциях:

• International Colloquium "Mossbauer Spectroscopy in Materials Science - MSMS96",

Lednice, Czech Republic, 1996;

• VII Всероссийская школа-семинар «Современные проблемы математического моделирования», Ростов-на-Дону, 1997;

• Всероссийская научная конференция «Алгоритмический анализ некорректных задач» - ААНЗ, Екатеринбург, 1998,2001;

• Международная научная конференция "Эффект Мёссбауэра: магнетизм, материаловедение, гамма-оптика", Казань, 2000;

• Всероссийская конференция молодых ученых, Ижевск, 2001;

• International Symposium Metastable, Mechanically Alloyed and Nanocrystalline Materials, University of Michigan, Ann Arbor, USA, 2001.

Основное содержание работы изложено в 7 статьях и 7 тезисах докладов (ссылки [44, 58-70] в списке литературы).

Структура диссертации. Работа состоит из введения, четырех глав, заключения, списка литературы.

4.4. Выводы

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

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

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

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

Заключение

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

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

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

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

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

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

Разработанные алгоритмы использовались при обработке мёссбауэровских спектров различных типов: неупорядоченных кристаллических сплавов Fe+15, 25, 33at.% Si; серии механически сплавленных в течение 2ч., 4ч., 8ч. и 16ч. при ТИЗМ=300К смеси порошков Fe(75)C(25); серии смеси порошков Fe(68)C(32) после механического сплавления в течение 1ч., 2ч., 4ч., 8ч., 12ч. при ТИЗМ=77К; выплавленного сплава Fe+8aT.% Sn после механического измельчения при ТИЗМ=300К; серии неупорядоченного кристаллического сплава Fe75Si25 после термической обработки.

1. Вертхейм Г.К. Эффект Мёссбауэра. -М.: Мир. -1966. -172С.

2. Литвинов B.C., Каракишев С.Д., Овчинников В.В. Ядерная гамма-резонансная спектроскопия сплавов. -М.: Металлургия. -1982. -144С.

3. Hadamard J. Sur les problemes aux derivees partielles et leur signification physique// Bull.Univ.Princeton. -1902. -Y.13. -P.49-52.

4. Тихонов A.H., Арсенин В.Я. Методы решения некорректных задач. -М.: Наука. -1986. -288С.

5. Морозов В.А. Регулярные методы решения некорректно поставленных задач. -М.: Наука. -1987. -240С.

6. Иванов В.К., Васин В.В., Танана В.П. Теория линейных некорректных задач и ее приложения. -М.: Наука. -1978. -206С.

7. Ортега Дж., Рейнболдт В. Итерационные методы решения нелинейных систем уравнений со многими неизвестными. -М.: Мир. -1975. -558С.

8. Phillips D.L. A technique for the numerical solution of certain integral equations of the first kind//J. Ass. Comput. Mach. -1962. -V.9. -P.84-97.

9. Twomey S. On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature// J. Ass. Comput. Mach. -1963. -V.10. -P.97-101.

10. Window B. Hyperfine field distributions from Mossbauer spectra// J. Phys. E: Sci. Instrum. -1971. -V.4. -P.401-402.

11. Varret F., Gerard A. and Imbert P. Magnetic field distribution analysis of the broadened Mossbauer spectra of zine feirite// Phys. Stat. Sol.(b). -1971. -V.43. -P.723-730.

12. Sharon Т.Е., Tsuei С.С. Magnetism in amorphous Fe-Pd-P alloys// Phys. Rev. -1972. -V.B5. -P. 1047-1064.

13. Billard L., Chamberod A. On dissymmetry of Mossbauer spectra in iron-nickel alloys// Solid State Commun. -1975. -V.17. -P. 113-118.

14. Tsuei C.C., Lilienthal H. Magnetization distribution in an amorphous ferromagnetic// Phys. Rev. B. -1976. -V.13. -P.4899-4906.

15. Logan J., Sun E. A Mossbauer study of amorphous iron-phosphorus alloys// J. Non-Crystalline Solid. -1976. -V.20. -P.285-298.

16. Vincze I. Evalution of complex Mossbauer spectra in amorphous and cristallne ferromagnets// Solid State Commun. -1978. -Y.25. -P.689-693.

17. Stone A.J. Improvement of resolution of Mossbauer and other spectra by Fourier transform techniques// Chem. Phys. Lett. -1970. -Y.6. -№4. -P.331-335.

18. Dibar Ure M.C., Flinn P.A. Technique for the removal of the "blackness" distortion of Mossbauer spectra// In "Mossbauer effect methodology". New York. Plenum Press. -1971. -V.7. -P.245-262.

19. Nagy D.L., Kulsar K. Resolution of Mossbauer spectroscopy// Proc. Int. Conf. Mossbauer spectrometry. Dresden. -1971. -Y.2. -P.618-624.

20. Hesse J. and Riibartsch A. Model independent evaluation of overlapped Mossbauer spectra// J. Phys. E.: Sci. Instrum. -1974. -V.7. -P.526-532.

21. Mongin P., Marchal G., Piecuch M. and Janor C. Mossbauer spectra analysis in amorphous system studies// J. Phys. E: Sci. Instrum. -1976. -V.9. -№12. -P.l 101-1105.

22. Belozerskii G.N., Khimich Yu.P. and Gitsovich Y.N. Parameters of hyperfine interaction from Mossbauer spectra of polycrystals in an externally applied magnetic field// Phys. Stat. Solidi(b). -1977. -V.79. -P. 125-129.

23. Le Caer G. and Dubois J.M. Evaluation of hyperfine parameter distributions from overlapped Mossbauer spectra of amorphous alloys// J. Phys. E: Sci. Instrum. -1979. -V.12. -P.1083-1090.

24. Wivel C. and Morup S. Improved computational procedure for evaluation of overlapping hyperfine parameter distributions in Mossbauer spectra// J. Phys. E: Sci. Instrum. -1980. -Y.14. -P.605-610.

25. Keller H. Evaluation of hyperfine field distributions from Mossbauer spectra using Window's Fourier method// J. Appl. Phys. -1981. -V.52. -№8. -P.5268-5273.

26. Овчинников В.В., Первухин Н.А. Сверхтонкая структура ЯГР спектров и анализ упорядочения сплавов методом восстановления функций плотности// ФТТ. -1981. -Т.23. -№10. -С.3095-3098.

27. Vincze I. Fourier evaluation of broad Mossbauer spectra// J. Nucl. Instr.and Meth. -1982. -V.199. -P.247-262.

28. Nikolov S. Method for decomposition of Mossbauer spectra of amorphous magnetic materials// Solid State Commun. -1983. -V.48. -№9. -P.761-763.

29. Белозерский Г.Н., Семенов Н.Ф., Чекашова Г.П. Изучение сверхтонких взаимодействий в аморфных сплавах, содержащих хром, с помощью мёссбауэровской спектроскопии// ФММ. -1983. -Т.55. -№5. -С.923-929.

30. Воронина Е.В., Ершов Н.В., Агеев A.JI., Бабанов Ю.А. Регулярный алгоритм решения задачи в ЯГР-спектроскопии. Свердловск. -1986. -27С. -Деп. рук. ВИНИТИ, №5704-В.

31. Rancourt D.G. and Ping J.Y. Voigt-based methods for arbitrary-shape static hyperfine parameter distributions in Mossbauer spectroscopy// J. Nucl. Instr. and Meth. B. -1991. -V.58. -P.85-97.

32. Nagy D.L. and Rohlich U. An overview on model-independent data reduction methods in Mossbauer spectroscopy// Hyperfine Interactions. -1991. -V.66. -P. 105 -126.

33. Kaptas D. and Vincze I. Evalution of the iron hyperfine field distribution in amorphous Fe90Zr10// Hyperfine Interactions. -1990. -V.55. -P.987-992.

34. Русаков B.C. Восстановление функции распределения сверхтонких параметров мёссбауэровских спектров локально неоднородных систем. -Алматы. -2000. —431 С.

35. Русаков B.C. Мёссбауэровская спектроскопия локально неоднородных систем // Изв. АН СССР, Сер. физическая. -1988. -Т.52. -В.9. -С.1783-1786.

36. Математические задачи обработки эксперимента// Сборник под. ред. В.Я.Галкина, П.Н.Заикина. -М.: Изд-во МГУ. -1984. -232С.

37. Горьков В.П., Меченов А.С. Численные методы решения обратных задач мёссбауэровской спектроскопии// В сб. Численные методы решения обратных задач матем. физики. Изд-во МГУ. -1988. -С.

38. Cieslak J. and Dubiel S.M. Harmonic analysis of Mossbauer spectra// J. Nucl. Instr. and Meth. B. -1995. -V.95. -P. 131-140.

39. Voronina E.Y., Ershov N.Y., Ageev A.L., Babanov Yu.A. Regular algorithm for the solution of the inverse problem in Mossbauer spectroscopy// Phys. Stat. Sol.(b). -1990. -V.160. -P.624-634.

40. Худсон Д. Статистика для физиков// Пер. с англ. под ред. Лейкина Е.М. -М.: Мир. -1970. -296С.

41. Ageev A.L. and Voronina E.V. Algorithm for parameter correction for solving the inverse problem of Mossbauer spectroscopy// J. Nucl. Instr. and Meth. B. -1996. -V.108. -P.417-424.

42. Агеев A.JI. Решение уравнений 1 рода с конечномерной нелинейностью// Изв. Вузов. Математика. -1997. -№>3. -С.68-72.

43. Агеев А.Л., Антонова Т.В., Воронина Е.В. Методы уточнения параметров при решении интегральных уравнений 1 рода// Матем. моделирование. -1996. -№12. -С.110-124.

44. Ageev A.L., Nemtsova O.M., Voronina E.V., E.P.Yelsukov Algorithm of parameters correction for solving the problem of Mossbauer spectroscopy// Chechoslovak Journal of Physics. -1997. -V.47. -№5. -P.547-552.

45. Агеев A.JI. Регуляризованный спектральный анализ и решение уравнений 1рода// Изв. Вузов. Математика. -1995. -№11. -С.3-16.

46. Годунов С.К., Антонов А.Г., Кирилюк О.П., Костин В.И. Гарантированная точность решения систем линейных уравнений в евклидовых пространствах// Новосибирск.: Наука. -Сиб. отделение. -1988. -456С.

47. Шпинель B.C. Резонанс гамма-лучей в кристаллах// М.: Наука. -1969. -407С.

48. Иркаев С.М., Кузьмин Р.Н., Опаленкао А.А. Ядерный гамма-резонанс// М.: Изд-во Моск. Ун-та. -1970. -235С.

49. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. -М.: Наука. -1989.-624С.

50. Жуковский Е.Л., Плискин С.Ю. Соотношение неопределенности и регуляризация систем линейных алгебраических уравнений//ДАН СССР. -1982. -Т.262. -№6. -С.1301-1304.

51. Жуковский E.JI. Метод наименьших квадратов для вырожденных и плохо обусловленных систем линейных алгебраических уравнений// ЖВМ и МФ. -1977. -Т. 17. -№4. -С. 814-826.

52. Котельников В.А. Теория потенциальной помехоустойчивости. М.: Госэнергоиздат. -1956. -243С.

53. Соколов С.Н., Силин И.Н. Нахождение минимума функционалов методом линеаризации// Препринт ОИЯИ-Д-810. -1961. -31 ОС.

54. Васин В.В., Агеев А.Л. Некорректные задачи с априорной информацией. Екатеринбург: Наука. -1996. -260С.

55. Агеев A.Jl. Регуляризованный спектральный анализ и решение уравнений 1 рода с конечномерной нелинейностью// Дис. доктора физ.-мат. наук. -Свердловск. -1997. -198С.

56. Дутчак Я.И. Рентгенография жидких кристаллов// Львов. -Высш. шк. -1977. -215С.

57. Немцова О.М. Об одном методе решения обратной задачи ядерного гамма-резонанса// Современные проблемы мат. моделирования, VII Всерос. школа-семинар. -1997. -С.95-99.

58. Немцова О.М. О решении уравнений 1 рода с конечномерной нелинейностью// Деп. рук. ВИНИТИ N2720-B00 Москва. -2000. -13С.

59. Дорофеев Г.А., Елсуков Е.П., Фомин В.М., Коныгин Г.Н., Загайнов А.В., Немцова О.М. Фазообразование в системе Fe(68)/C(32) при механическом легировании// Физика, химия обработки материалов. -2001. -№5. -С.71-78.

60. O.M.Nemtsova, A.L.Ageev, E.V.Voronina The estimation of the error of the hyperfine interaction parameter distribution from Mossbauer spectra// J. Nucl. Instr.and Meth.B. -2002. -V. 187. -№1. -P. 132-136.

61. Yelsukov E.P., Fomin F.M., Nemtsova O.M., Shklyaev D.A. Evolution of structure in the disordered nanocrystalline Fe75Si25 alloy tinder heat treatment// Mater. Sci. Forum. -2002. -V.386-388. -P.181-186.

62. Немцова O.M., Агеев А.Л. О решении операторных уравнений 1 рода с конечномерной нелинейностью// Тезисы докладов Всероссийской научной конференции ААНЗ, Екатеринбург. -1998. -С. 185-186.

63. Немцова О.М. Математическая обработка экспериментальных данных мёссбауэровской спектроскопии// Тезисы докладов на Всероссийской конференции молодых ученых. Ижевск.-2001.-С.22.

64. Немцова О.М., Воронина Е.В. Оценка ошибки решения обратной задачи мёссбауэровской спектроскопии// Тезисы докладов Всероссийской научной конференции ААНЗ. Екатеринбург. -2001. -С.298-299.

65. Елсуков Е.П., Дорофеев Г.А., Ульянов A.JL, Немцова О.М., Порсев В Е. Твердофазные реакции в системе Fe(68)Ge(32) при механическом сплавлении// ФММ. -2002. -Т.94. -№4. (принята в печать).

66. Гапоненко Ю.Л. О единственности решения интегрального уравнения// Вестн. МГУ. Сер. Вычисл. Матем. И кибернетика. -1977. -№4. -С.14-17.

67. Лаврентьев М.М., Романов В.Г., Шишатский С.П. Некорректные задачи математической физики и анализа//М.: Наука. -1980. -257С.

68. Регуляризующие алгоритмы и априорная информация/ А.Н.Тихонов, А.В.Гончарский, В.В.Степанов, А.Г.Ягола. -М.: Наука. -1983. -200С.

69. Като Т. Теория возмущений линейных операторов. М.: Мир. -1972. -С.98.

70. Свешников А.Л. Основы теории ошибок// Л.: ЛГУ. -1972. -123С.

71. Гантмахер Ф.Р. Теория матриц// М.: Наука. -1988. -552С.