.RU

§4. Метод наименьших квадратов - Конспект лекций по методам конечных элементов


§4. Метод наименьших квадратов.

Пусть и для оператора существует ограниченный обратный оператор . Тогда общая задача принимает вид

(4.1)

И алгоритм (1.3) называется методом наименьших квадратов. Формулировка его заключается в следующем:

1) Выбираются базисные функции

2) Приближенное решение ищется в виде



3)Коэффициенты определяются из системы уравнений

(4.2)

которую можно записать в матричной форме

(4.3)

где



Вычислив элементы симметричной матрицы , компоненты вектора и решив систему (4.3) можно определить приближенное решение .

Часто при формулировке данного алгоритма соотношение (4.2) получают из условия минимизации функционала невязки на линейной оболочке системы .

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

(4.4)

Дифференцируя (4.4) по получаем систему уравнений

(4.5)

Т.е. систему (4.2). Решив (4.5), определим коэффициенты , которые будут минимизировать функционал невязки.

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

Доказательство. Пусть . имеет только нулевое решение . Тогда, как известно ранее, функции линейно независимы при любом .Необходимым и достаточным условием линейной независимости системы является неравенство нулю определителя Грамма: Следовательно, определитель системы (4.5) отличен от нуля и она имеет единственное решение при любом .

Рассмотрим теперь вопрос о сходимости.

Теорема 1. Метод наименьших квадратов дает последовательность приближенных решений , сходящихся к точному решению , если последовательность является А-полной и существует константа такая, что для любой функции имеет место соотношение

(4.6)

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

Следовательно, по лемме 1, приближенные решения могут быть построены при любом . Далее, поскольку А-полна, то при заданном можно найти такие что



Это неравенство является верным, если заменить на , определяемые из (4.5), так как тогда левая часть этого неравенства достигает минимума. В этом случае при имеем и из (4.6) следует, что , т.е. при . Теорема доказана.

Следствие. Если условия теоремы 1 выполнены, то



. (4.7)

Из этого следствия вытекает, что если построены по методу наименьших квадратов, то формула (4.7) позволяет судить о погрешности приближенных решений .

Отметим связь между методом наименьших квадратов и методом Ритца. Предположив, что где сопряженный оператор, и подействовав на (4.1) оператором , получаем уравнение

(4.8)

Если оно разрешимо, то его решение минимизирует функцию





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

Преимущества и недостатки метода наименьших квадратов. Пусть, например А – дифференциальный оператор порядка . Решая уравнение (4.1) методом наименьших квадратов, мы фактически решаем методом Ритца уравнение (4.8), оператор которого имеет уже порядок . Еще одна трудность – требование (требование гладкости и удовлетворение краевых условий). Преимущество : кроме обратимости А никаких условий больше не накладываются.

Лекция №№

Пример 1. Пусть задача



(4.9)

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



И решая систему (4.5), находим : имеем систему



и, значит,

.

Следовательно, приближенное решение имеет вид

.

Итак, при выбранных базисных функциях методы Ритца, Бубнова-Галеркина и наименьших квадратов в рассматриваемом примере приводят к одному и тому же результату. Оценим сходимость. Непосредственное вычисление показывает, что

(4.10)

(4.11)

Если для оценки погрешности воспользоваться неравенством (4.7), то Сравнивая этот результат с (4.10) делаем вывод, что оценка, полученная при помощи (4.7), в некоторых случаях может оказаться огрубленной.


§ 5. Обобщенный метод наименьших квадратов.

Рассмотрим задачу (1.1):

(5.1)

и пусть в (1.3) выполнены следующие ограничения: , и - произвольные операторы, причем предполагается существование . Алгоритм приближенного решения задачи в этом случае имеет вид

, (5.2)

(5.3)

и называется обобщенным методом наименьших квадратов.

Сформулируем достаточные условия сходимости к точному решению в этом алгоритме. Введем в новое скалярное произведение



и норму

.

Пополнив по норме получим гильбертово пространство, которое обозначим через . Теперь уравнение (5.1) можно записать в виде

(5.4)

а уравнения для приближенного решения можно представить в эквивалентной форме

(5.5)

Теорема 1. Пусть: 1) уравнение (5.1) однозначно разрешимо; 2) полна в ; 3) оператор вполне непрерывен в . Тогда система приближенных уравнений (5.3) однозначно разрешима при достаточно больших и сходятся к при как в метрике пространства , так и в метрике пространства .

Доказательство теоремы повторяет доказательство теоремы 1 § 3 (метод Бубнова-Галеркина) и остается в качестве задачи.


§ 6. Обобщенный метод моментов.


Пусть в гильбертовом пространстве (в этом параграфе – комплексный) рассматривается уравнение (1.1):

. (6.1)

Предположим, что в общей постановке (1.3) выполнены следующие ограничения: , оператор является -положительно определенным (или положительно определенным в обобщенном смысле), т.е.



Где - постоянные, . Алгоритм (1.3) в этом случае называется обобщенным методом моментов, он состоит в следующем:

  1. выбирается базисная система

  2. приближенное решение ищется в виде

. (6.2)

  1. Коэффициенты определяются из системы уравнений

(6.3)

или, что одно и то же, из системы

(6.4)

Прежде всего отметим, что в силу -положительной определенности оператор обладает ограниченным обратным оператором. Действительно, так как и , то . В силу теорем об обратном операторе существует и . Из условий положительной определенности вытекает также, что число вещественно. Покажем, что

(6.5)

Для этого выберем и запишем тождество



(6.6)

и аналогичное ему, поменяв и местами:



(6.7)

Но вещественно, так что сравнивая слагаемые в (6.6) и (6.7) с учетом свойств скалярного произведения, получаем, что правые части в (6.6) и (6.7) совпадают, т.е. имеем (6.5).

На основании свойств (6.5) и -положительной определенности оператора на можно ввести скалярное произведение

(6.8)

в силу чего D (A) превращается в новое гильбертово пространство HK, которое можно считать полным *). Метрику в HK определим соотношением =(u, u).

Подействовав А на обе части (6.1), получим уравнение

u+Tu = f, f= Af, T = AB.

Определение. Элемент uHK назовем обобщенным решением уравнения (6.1), если он удовлетворяет (6.9) в пространстве HK, т. е. удовлетворяет равенству

(6.10)

Очевидно, что если элемент u удовлетворяет (6.1), то он является также и обобщенным решением (обратное, вообще говоря, неверно).

Теперь систему (6.4) можно записать в виде

j = 1, …, N (6.11)

или

j = 1, …, N

Алгоритм (6.12) можно рассматривать как процесс определения приближенного обобщенного решения u.

Теорема 1. пусть уравнение (6.1) имеет единственное обобщенное решение и оператор T = AB вполне непрерывен в HK. Тогда: 1) существует такое целое N, что при любом N N система (6.12) имеет единственное решение а; 2) приближенные решения u сходятся в HK (u в H) к решению уравнения (6.1).

Доказательство. Обозначим через систему , полученную путем ортогонализации в метрике пространства HK при нормировке = 1. Тогда обобщенное решение u может быть представлено в виде






причем ряд сходится в HK и коэффициенты c удовлетворяют системе



т. е. системе



где = (T,), b= (f, K).

Системы (6.14) и (3.9) имеют одинаковый вид, и дальнейшее доказательство теореме повторяет рассуждение в доказательстве теоремы 1 … 3.

Отметим, что так как здесь T = AB предполагается вполне непрерывным в HK, то



Кроме того,



Следовательно в обобщенном методе моментов можно сделать предположение о принадлежности к HK (а не к D(A)). Поэтому, если в этом случае предположить плотность базиса в HK, то утверждения теоремы 1 остается справедливым. Система для определения коэффициентов абудет записываться здесь в виде (6.12) (а не в (6.4)), а сама задача решения (6.1) рассматривается в обобщенной постановке (6.10).


§ 7. проекционный метод в гильбертовом пространстве


В этом параграфе будет изложен проекционный метод решения уравнения

Lu = Au + Bu = f, f H, (7.1)

Причем он будет рассмотрен лишь для случая гильбертовых пространств. Ниже мы покажем также, что все приведенные раннее алгоритмы можно рассматривать как частные случаи проекционного метода.

Пусть L – вообще говоря, неограниченный оператор, действующий в H и обладающий ограниченным обратным оператором L. Предположим, что множество определения D(L) и множество значения R(L) плотны в H.

Введем в Н линейно независимую систему . Соответствующие N-мерные подпространства, порождаемые , обозначим через M. Предположим, что последовательность предельно плотна в ^ Н. Зададим последовательность проекционных операторов P, каждый из которых отображает Н на соответствующее пространство M. Предположим в дальнейшем, что , N = 1,2, … (Здесь P не обязательно являются ортопроекторами, т. е. от них требуется лишь выполнение свойств, . Отметим, что для приводимых ниже рассуждений полезно рассмотреть определение проектора, ортопроектора и пример 3 из § 3.)

Введем также линейно независимую систему . Пространство, порождаемое , обозначим через НN, а линейную оболочку системы - через LHN. Предположим, что последовательность подпространств предельно плотна в Н и для любого элемента



Приближенное решение задачи (7.1) ищем в виде



Где определяется из уравнения



Теорема 1. Пусть для любого N и любого элемента

> 0,

где постоянная не зависит от N. Тогда при любом N уравнение (7.3) имеет единственное решение .

Доказательство. Пусть . Тогда (7.3) записывается в виде



Обозначим через сужение на , т.е. - оператор с областью определения , для которого при . Согласно (7.4) и теореме об обратном операторе обладает ограниченным обратным оператором , переводящим на LHN, при этом Следовательно, существует единственный элемент удовлетворяющий (7.5). Но тогда - единственный элемент, удовлетворяющий (7.3). Теорема доказана.

Теорема 2. если выполнены условия теоремы 1, то при любом N существует приближенное значение причем невязка стремится к нулю при и справедлива оценка

(7.6)

где

Доказательство. Пусть - решение (7.3) (которое согласно теореме1 существует), .выберем произвольный элемент . Так как то



Отсюда, используя оценку получаем



В силу произвольности элемента , а также условия плотности в Н, получаем сходимость невязки к нулю и оценку (7.6).

Следствие. Если то имеют место сходимость к и при и оценка погрешности



т.е.



Замечание. Пусть А, В, К – операторы, относительно которых справедливы предположения …1. Предположим, что есть базис в Н. Ведем систему линейно независимых функций которая К–полна. В этом случае система может быть принята в качестве нового базиса Н. Пусть в (7.3) принято оно из основных ограничений: является оператором ортогонального проектирования. В этом случае уравнение (7.3) эквивалентно системе уравнений (см. пример 3 из … 3)

,

т. е. система вида

,

которая с точностью до обозначений совпадает с системой (1.3). Таким образом , если в (7.3) оператор ортопроектором (а не просто проектором на ), то алгоритм (1.3) совпадает с (7.3), т. е. (1.3) есть частный случай проекционного метода.


4-poryadok-priema-na-obuchenie-perevod-prinimat-ustav-naprimer-sovet-konferenciya-obshee-sobranie-i-dr.html
4-poryadok-provedeniya-otkritogo-zaprosa-predlozhenij-instrukcii-po-podgotovke-predlozhenij.html
4-poryadok-rassmotreniya-ocenki-i-sopostavleniya-zayavok-razvitie-grazhdanskoj-morskoj-tehniki.html
4-poryadok-vedeniya-bankovskogo-scheta-pravila-polzovaniya-bankovskimi-kartami-dlya-derzhatelej-kart-visa-bank24-ru-oao.html
4-posledovatelnost-vipolneniya-proekta-uchebnoe-posobie-izdatelstvo-tpu-tomsk-2005.html
4-postanovleniya-pravitelstva-ob-utverzhdenii-instrukcii-o-poryadke-podgotovki-i-soglasovaniya-proektov-normativnih.html
  • write.bystrickaya.ru/glava-9vospitanie-i-obuchenie-detejs-tyazhelimi-porazheniyamisluha-i-zreniya.html
  • literature.bystrickaya.ru/doklad-direktora-stranica-3.html
  • bukva.bystrickaya.ru/podrazdel-stalej-na-tipi-i-klassi-prikaz-komiteta-po-nadzoru-za-ohranoj-truda-ukraini-gosnadzorohrantruda.html
  • desk.bystrickaya.ru/perspektivnoe-planirovanie-rebenok-i-okruzhayushij-mir-prirodnoe-okruzhenie-yavleniya-obshestvennoj-zhizni.html
  • grade.bystrickaya.ru/o-molekulah-i-ih-puteshestviyah.html
  • universitet.bystrickaya.ru/tablica-4-zakonomernost-i-sluchajnost-sluchajnaya-izmenchivost-v-tochnih-naukah-v-biologii-i-medicine.html
  • urok.bystrickaya.ru/prilozheniya-plan-doklada-na-sezde-rabochih-poselyan-i-gorcev-chernomorskogo-okruga-26-iyunya-1920-goda-11-1920-goutline-.html
  • college.bystrickaya.ru/23-organizaciya-zashiti-vkr-metodicheskie-ukazaniya-po-vipolneniyu-i-oformleniyu-vipusknoj-kvalifikacionnoj-raboti.html
  • shpargalka.bystrickaya.ru/viigiliya-shestaya-ernst-teodor-amadej-gofman-zolotoj-gorshok-skazka-iz-novih-vremen.html
  • studies.bystrickaya.ru/esli-vi-zabludilis-v-tajge.html
  • knigi.bystrickaya.ru/resheniem-edinstvennogo-akcionera-stranica-2.html
  • znaniya.bystrickaya.ru/razdel-ii-ni-a-takzhe-voprosi-ee-organizacii-podgotovki-kadrov-povishe.html
  • kolledzh.bystrickaya.ru/a-s-gormin-teoreticheskie-osnovi-obucheniya-i-vospitaniya-odarennih-detej-v-paradigme-barernoj-pedagogiki.html
  • learn.bystrickaya.ru/g-l-a-v-a-3-stroevoj-ustav-vooruzhennih-sil-respubliki-tadzhikistan.html
  • portfolio.bystrickaya.ru/otchet-o-rezultatah-samoobsledovaniya-belomorskogo-filiala-petrgu-po-sostoyaniyu-na-01-05-2008-stranica-4.html
  • student.bystrickaya.ru/1-organizacionno-pravovoe-obespechenie-deyatelnosti-obrazovatelnogo-uchrezhdeniya-6-stranica-20.html
  • ucheba.bystrickaya.ru/programma-mezhdunarodnoj-konferencii-mestnoe-samoupravlenie-i-razvitie-territorij.html
  • otsenki.bystrickaya.ru/rol-proizvoditelnosti-truda-v-formirovanii-stran-liderov-v-istorii-proizvodstva.html
  • credit.bystrickaya.ru/ouli-1-molekulali-kinetikali-teoriyani-negizgi-imdari-dene-rilimini.html
  • pisat.bystrickaya.ru/tabl-23-g-kalendar-podvizhnih-cerkovnih-dat-utochnenie-dat-po-prazdnikam.html
  • notebook.bystrickaya.ru/izpitna-programa-za-opredelyane-na-godishna-ocenka-po-blgarski-ezik-i-literatura-zadlzhitelna-podgotovka-za-x.html
  • znaniya.bystrickaya.ru/programma-professionalnoj-perepodgotovki-po-napravleniyu-menedzhment-kachestvom-normativnij-srok-prohozhdeniya.html
  • uchit.bystrickaya.ru/uchebnij-plan-uchebnie-programmi-kursov-povisheniya-kvalifikacii-nalchik.html
  • zadachi.bystrickaya.ru/s-bondarenko-v-kurilskij-neizvestnie-strugackie-pisma-rabochie-dnevniki-1942-1962g-g-stranica-10.html
  • uchebnik.bystrickaya.ru/vopros-4-v-publicisticheskom-stile-russkij-yazik-v-sovremennom-mire.html
  • nauka.bystrickaya.ru/voprosi-deklarirovaniya-i-kontrolya-tamozhennoj-stoimosti-reguliruyutsya-i-nacionalnim-zakonodatelstvom.html
  • uchit.bystrickaya.ru/tairibi-trapeciyani-trt-tamasha-sizii-turali-sekciyasi-matematika-zhetekshs.html
  • notebook.bystrickaya.ru/i-n-denisov-8-fevralya-2001-g.html
  • university.bystrickaya.ru/glava-6-kak-dobivatsya-soglasiya-tri-osnovnih-obyazatelnih-usloviya-i-chetire-sposoba-obespechit-ih-soblyudenie-pri-zaklyuchenii-sdelki.html
  • write.bystrickaya.ru/glava-2-8-obekt-nakazaniya-n-i-faleev-celi-voennogo-nakazanii.html
  • shpargalka.bystrickaya.ru/uchebno-metodicheskij-kompleks-disciplini-pedagogika-vvedenie-v-pedagogicheskuyu-deyatelnost-vvedenie-v-pedagogicheskuyu-deyatelnost.html
  • notebook.bystrickaya.ru/iii-teoriya-pismennih-istochnikov-d-gatri-vvedenie-v-novij-zavet.html
  • studies.bystrickaya.ru/33-subekti-prava-na-sredstva-individualizacii-zashita-prav-intellektualnoj-sobstvennosti-problemi-i-puti-resheniya.html
  • shpargalka.bystrickaya.ru/volshebnaya-palochka-dlya-uchitelya-raboti-pobeditelej-vserossijskogo-konkursa-pedagogicheskih-esse.html
  • holiday.bystrickaya.ru/mmg-4210-matematicheskie-metodi-v-geologii-vostochno-kazahstanskij-gosudarstvennij-tehnicheskij.html
  • © bystrickaya.ru
    Мобильный рефератник - для мобильных людей.