next up previous
Next: Литература Up: О (НЕ)ЕДИНСТВЕННОСТИ РЕШЕНИЯ Previous: (Не)единственность для совокупности уравнений

3. Исследование прикладных интегральных уравнений

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

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

Метод рентгеноспектрального структурного анализа (РССА). Этот метод предназначен для исследования структуры неупорядоченных материалов. Для однокомпонентных материалов уравнение РССА (см. гл.6 [9]) имеет вид
eqnarray352
где tex2html_wrap_inline1198 и U- известные константы, f(s) и tex2html_wrap_inline1204 - известные функции. Функция tex2html_wrap_inline1206 есть плотность математического числа атомов расположенных на расстоянии r (в ангстремах) от "центрального" атома. Типичные значения параметров a=2, b=8, c=3.5, d=16 определяются как физикой эффекта, так и возможностями эксперимента. Оператор предполагается действующим из tex2html_wrap_inline1212 в tex2html_wrap_inline1214 а искомое решение tex2html_wrap_inline1216 предполагается принадлежащим tex2html_wrap_inline1218 (оператор L есть оператор обратный оператору вложения tex2html_wrap_inline1224) в tex2html_wrap_inline1212).

В [9] п. 1.3 гл.6 подробно описаны модельные расчеты решения уравнения РССА вариационным методом Тихонова с tex2html_wrap_inline1226. На основе расчетов была сформулирована гипотеза, что искомая функция tex2html_wrap_inline1216 отличается от нормального решения уравнения РССА, аппроксимируемого методом Тихонова. Наша цель состоит в том, чтобы подтвердить или опровергнуть эту гипотезу.

Опишем возможную схему оценки tex2html_wrap_inline1230 и r для уравнения РССА (20). Поскольку нам известно точное решение tex2html_wrap_inline1234 то это сразу дает нам шкалу измерения tex2html_wrap_inline1236 Например, будем считать относительную погрешность в изменении решения свыше 300% большой (конечно, такой процент завышен по крайней мере на порядок, он выбран для более наглядного представления результатов), т.е. tex2html_wrap_inline1238

Выберем число r таким образом, чтобы для всех вариаций решения tex2html_wrap_inline1242 величины tex2html_wrap_inline1244 tex2html_wrap_inline1246 были порядка tex2html_wrap_inline1248 Для этого, например, потребуем, чтобы tex2html_wrap_inline1250

Перейдем к оценке точности h задания оператора в уравнении (20). Применительно к уравнению РССА можно утверждать, что h>0, т.к. фаза и амплитуда в уравнении (20) (функции tex2html_wrap_inline1256) заданны приближенно. С другой стороны, по-видимому, эффективной процедуры оценки h нет. Известно, что относительная погрешность задания правой части уравнения (20) при обработке экспериментальных данных есть примерно tex2html_wrap_inline1260 Поэтому потребуем выполнения условия tex2html_wrap_inline1262 т.е. положим tex2html_wrap_inline1264

Далее удобно считать, что собственные вектора из tex2html_wrap_inline1266 нормированы на норму модельного решения: tex2html_wrap_inline1268 Тогда согласно следствию 17 необходимо найти параметр регуляризации tex2html_wrap_inline1270 (такое значение существует) и вычислить величину
eqnarray383
Для того, чтобы tex2html_wrap_inline1272 (случай неединственности), необходимо и достаточно, чтобы выполнялось неравенство tex2html_wrap_inline1274 При этом, в частности, элементы tex2html_wrap_inline1276

Минимум в (21) достигается на собственной функции tex2html_wrap_inline1278 (с нормировкой tex2html_wrap_inline1280) оператора
displaymath1282
отвечающей наименьшему собственному числу этого оператора (для краткости сразу записана матрица, аппроксимирующая оператор tex2html_wrap_inline794). Аппроксимация матрицы tex2html_wrap_inline692 выполнена методом коллокации "на ступеньках" при n=120,m=350. Тогда мы имеем tex2html_wrap_inline1290 Поскольку tex2html_wrap_inline1292 то неравенство tex2html_wrap_inline1294 выполняется и мы имеем для уравнения РССА случай неединственности. Поэтому, не привлекая дополнительной информации, нельзя гарантировать ни для какого метода регуляризации восстановление нужного нам решения.

Отметим, что сделанный вывод о неединственности решения уравнения РССА соответствует интуитивным представлениям об уравнении, имеющем неединственное решение (подробнее см. п.3.3 гл.6 и рис.5,7 в [9]).

Метод дифракции. Другим методом исследования неупорядоченных материалов является метод дифракции. Основное уравнение в этом методе имеет вид
eqnarray406
где функция f(s)- известна. Оператор предполагается действующим из tex2html_wrap_inline1212 в tex2html_wrap_inline1214 а искомое решение предполагается принадлежащим tex2html_wrap_inline1302

В общем случае интервал [c,d] сильно отличается от tex2html_wrap_inline1306 Чтобы проиллюстрировать исследование уравнения, имеющего единственное решение, рассмотрим здесь случай, когда интервал измерения [c,d] близок к интервалу tex2html_wrap_inline1310 Тогда, ввиду связи (22) с прямым tex2html_wrap_inline1312-преобразованием, можно ожидать, что уравнение (22) имеет единственное решение.

Перейдем к исследованию (не)единственности решения уравнения (22), аналогично тому, как это было сделано в предыдущем пункте ( относительную точность задания правой части можно оценить в 1%). Положим tex2html_wrap_inline1314 Поскольку tex2html_wrap_inline1316 то мы имеем h=0.062, r=2.8

Обозначим через tex2html_wrap_inline1320 собственную функцию (с нормировкой tex2html_wrap_inline1322) оператора tex2html_wrap_inline1324 отвечающую наименьшему собственному числу этого оператора
displaymath1326
Матрица tex2html_wrap_inline1328 аппроксимирует оператор tex2html_wrap_inline692 методом коллокации "на ступеньках" при tex2html_wrap_inline1332 есть матрица аппроксимирующая оператор tex2html_wrap_inline1334 обычной второй разностью. Тогда при tex2html_wrap_inline1336 мы имеем, что tex2html_wrap_inline1338 Это значит, что tex2html_wrap_inline1340 (мы взяли tex2html_wrap_inline1342 в "вилку").

Поскольку tex2html_wrap_inline1344 то tex2html_wrap_inline1346 и в рамках заданной точности мы имеем единственное решение.

Как и для уравнения РССА, сделанный вывод о единственности решения дифракционного уравнения соответствует интуитивным представлениям об уравнении, имеющем единственное решение (подробнее см. п.3.4 гл.6 и рис.8 в [9]).

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант N 97-01-00520).

Поступила 30.09.97


next up previous
Next: Литература Up: О (НЕ)ЕДИНСТВЕННОСТИ РЕШЕНИЯ Previous: (Не)единственность для совокупности уравнений