НЕКОРРЕ́КТНЫХ ЗАДА́Ч ЧИ́СЛЕННЫЕ МЕ́ТОДЫ РЕШЕ́НИЯ
-
Рубрика: Математика
-
-
Скопировать библиографическую ссылку:
Книжная версия:
Электронная версия:
НЕКОРРЕ́КТНЫХ ЗАДА́Ч ЧИ́СЛЕННЫЕ МЕ́ТОДЫ РЕШЕ́НИЯ, приёмы и методы решения задач, для которых не удовлетворяется хотя бы одно из приводимых ниже условий. Задача определения решения уравнения z=R(u), где z – элемент метрического пространства Z (с расстоянием ρZ), по «исходным данным» u из метрич. пространства U (с расстоянием ρU), называется корректно поставленной на паре пространств (Z,U), если: а) для любого u∈U существует решение z∈Z; б) решение определено однозначно; в) задача устойчива на паре (Z,U), т. е. для любого ε>0 существует δ(ε)>0 такое, что для любых u1,u2∈U из неравенства ρU(u1,u2)⩽δ(ε) следует неравенство ρZ(z1,z2)⩽ε, где z1=R(u1),z2=R(u2). Это понятие корректности принадлежит Ж. Адамару (1923), который считал, что всякая математич. задача, соответствующая к.-л. физич. или технич. задаче, должна быть корректной, поскольку в противном случае при сколь угодно малых изменениях исходных данных решения могут сильно отличаться, чему нельзя дать физич. интерпретацию. Однако такая точка зрения, естественная для мн. задач естествознания и техники, не может быть перенесена на все задачи. Например, неустойчивыми являются задачи дифференцирования функций, известных приближенно, численного суммирования рядов Фурье, когда их коэффициенты известны приближенно, задача Коши для уравнения Лапласа.
К некорректным задачам относится широкий класс т. н. обратных задач, возникающих в физике и технике, в частности задачи обработки результатов физич. экспериментов. Пусть величина z – количественная характеристика (функция, вектор) изучаемого явления (объекта), часто о величине z известно, что она принадлежит некоторому подмножеству M метрич. пространства Z. В физич. эксперименте величина z часто недоступна непосредств. измерению, а измеряется лишь её проявление u=Az. Для интерпретации результата измерения необходимо определять z по u, т. е. решать уравнениеAz=u.(1)
Если u∈AM, где AM – образ M при его отображении с помощью оператора A, то решение этого уравнения есть z=A–1u, где A–1 – оператор, обратный оператору A. Так как элемент u часто получают путём измерений, то он обычно бывает известен лишь приближенно; пусть ˜u – его приближенное значение. В этих условиях речь может идти лишь о нахождении приближенного (к z) «решения» уравненияAz=˜u.(2)
Оператор A во многих случаях таков, что обратный ему оператор A–1 не является непрерывным. В этих случаях в качестве приближенного решения уравнения (1) нельзя брать точное решение уравнения (2), т. е. элемент , так как: а) такого решения может не существовать на M, поскольку может не принадлежать множеству AM; б) такое решение, даже если оно существует, может не обладать свойством устойчивости к малым изменениям (поскольку обратный оператор A–1 не является непрерывным) и поэтому не может быть физически интерпретируемым. Задача (2) является некорректной.
Для некорректных задач вида (1), (2) возникают вопросы: а) что понимается под приближенным решением таких задач? и б) каковы алгоритмы построения решений? Эти вопросы были впервые рассмотрены А. Н. Тихоновым (1963).
Метод подбора. В некоторых случаях приближенные решения уравнения (1) находятся методом подбора. Он состоит в том, что из множества M возможных решений, M⊂Z, выбирают элемент ˜z, для которого A˜z приближает правую часть уравнения (1) с заданной точностью, и в качестве искомого приближения берут элемент ˜z. Этот метод применим, когда из неравенства ρU(A˜z,Az)⩽δ следует, что ρZ(˜z,z)⩽ε(δ), где ε(δ)→0 при δ→0. Это имеет место при условии однозначной разрешимости уравнения (1) и при условии, что множество M – компакт. На основе этих условий вводится понятие корректности по Тихонову, называемое также условной корректностью. В применении к уравнению (1) задача называется корректной по Тихонову, если известно, что для точного значения u правой части существует единственное решение z уравнения (1), принадлежащее заданному компакту M. В этом случае оператор A–1 непрерывен на множестве M и если вместо элемента u известен элемент uδ такой, что rU(uδ,u)⩽δ и uδ∈AM, то в качестве приближенного решения уравнения (1) с правой частью u=uδ можно брать элемент zδ=A–1uδ. При этом zδ→z, если δ→0.
Во многих случаях приближенно известная часть ˜u уравнения (2) не принадлежит множеству AM. В этих условиях уравнение (2) не имеет классич. решения и в качестве его приближенного решения берётся обобщённое решение, называемое квазирешением. Элемент ˜z∈М, минимизирующий при данном ˜u функционал ρU(Az,˜u) на множестве M, называется квазирешением уравнения (2) на M. Если M – компакт, то квазирешение существует для любого ˜u∈U, а если ˜u∈AM, то квазирешение ˜z совпадает с классическим (точным) решением уравнения (2). Существование квазирешения гарантируется лишь при условии компактности множества M возможных решений.
Метод регуляризации. Для ряда прикладных задач, приводящих к уравнению (1), характерна ситуация, когда множество M возможных решений не является компактом, оператор A–1 не является непрерывным на AM и изменения правой части этого уравнения, связанные с её приближенным характером, могут выводить её за пределы множества AM. Такие задачи называются существенно некорректными. Разработан подход к решению таких задач, называемый методом регуляризации. В дальнейшем для простоты предполагается, что приближенной является лишь правая часть уравнения (1). В основе подхода лежит понятие регуляризирующего оператора. Оператор R(u,α) из U в Z, зависящий от параметра α, называется регуляризирующим оператором для уравнения Az=u (в окрестности точки u), если он обладает свойствами: 1) существует такое число δ1>0, что оператор R(u,α) определён для любого α>0 и любого uδ∈U, для которого ρU(uδ,u)⩽δ⩽δ1; 2) существует функция α=α(δ) такая, что для любого ε>0 найдётся число δ(ε)⩽δ1 такое, что если uδ∈U и ρU(uδ,u)⩽δ(ε), то ρZ(zδ,z)⩽ε, где zδ=R(uδ,α(δ)). В этом определении не предполагается однозначности оператора R(u,α).
Если ρU(uδ,u)⩽δ , то в качестве приближенного решения уравнения (1) с правой частью uδ можно брать элемент zα(δ)=R(uδ,α(δ)), где α(δ) согласовано с погрешностью исходных данных uδ. Это решение называется регуляризированным решением уравнения Az=uδ, а параметр α называется параметром регуляризации. При δ→0 регуляризированное решение сходится к решению уравнения (1). Т. о., задача нахождения приближенных решений уравнения (1) сводится к нахождению регуляризирующего оператора и к определению параметра регуляризации по дополнит. информации о задаче, напр. по оценке погрешности δ, с которой задаётся правая часть (1).
В основе построения регуляризирующих операторов лежат разл. принципы, среди них – вариационный, при котором используется сглаживающий функционалMα(z,A,u)=ρ2U(Az,u)+αΩ(z),где α>0 – параметр регуляризации, а Ω – регуляризирующий функционал, предназначенный для стабилизации решения. Он выбирается так, чтобы множества {z:Ω(z)⩽c} при всех с>0 были компактными в пространстве Z. Ищется элемент zα, минимизирующий Mα на множестве возможных решений. Параметр α находится по дополнит. информации о задаче.