АвтоАвтоматизацияАрхитектураАстрономияАудитБиологияБухгалтерияВоенное делоГенетикаГеографияГеологияГосударствоДомДругоеЖурналистика и СМИИзобретательствоИностранные языкиИнформатикаИскусствоИсторияКомпьютерыКулинарияКультураЛексикологияЛитератураЛогикаМаркетингМатематикаМашиностроениеМедицинаМенеджментМеталлы и СваркаМеханикаМузыкаНаселениеОбразованиеОхрана безопасности жизниОхрана ТрудаПедагогикаПолитикаПравоПриборостроениеПрограммированиеПроизводствоПромышленностьПсихологияРадиоРегилияСвязьСоциологияСпортСтандартизацияСтроительствоТехнологииТорговляТуризмФизикаФизиологияФилософияФинансыХимияХозяйствоЦеннообразованиеЧерчениеЭкологияЭконометрикаЭкономикаЭлектроникаЮриспунденкция

Рассмотрим второй и третий случаи

Читайте также:
  1. V2: Культура Российской империи второй половины XIX – начала ХХ вв.
  2. XX съезд КПСС. Процесс политической реабилитации и десталинизации во второй половине 1950 – начале 1960-х гг. и его значение.
  3. Актерское искусство второй половины XIX века
  4. Английские колониальные захваты в Индии во второй половине XVIII в.
  5. Антигоспитальное в области психиатрии движение в мире во второй половине XX века
  6. Археологические исследования второй половины XIX – первой трети XX вв. (с.43)
  7. Белорусские города во второй пол 13 – первой пол 17 вв. Развитие ремесла и торговли.
  8. Билет23. Коренной перелом в ходе Великой Отечественной войны и Второй мировой войны.
  9. Билет25 Классификация систем линейных уравнений по числу решений, ступенчатый вид расширенной матрицы системы в каждом случаи.
  10. Болгарская Православная Церковь во второй половине ХХ – начале XXI вв. Церковный раскол в Болгарии.
  11. БУРЖУАЗНЫЕ РЕФОРМЫ ВТОРОЙ ПОЛОВИНЫ XIX в. И ИХ ЗНАЧЕНИЕ
  12. В младшем школьном возрасте (второй период детства)

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

к курсовой работе

По программированию

на тему:

«Решение системы линейных уравнений»

 

 

Выполнил: Принял:

ст.гр.020603 Навроцкий А.А.

Червоный А.В.

 

 

Минск 2001г.


Содержание

Введение.

1. Анализ существующих методов решения задачи.

2. Описание используемого метода.

3. Анализ результатов.

Вывод.

Список использованной литературы.

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

 


Введение

 

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

В общем виде СЛАУ выглядит так: Ax=b, и при решении данной системы возможно несколько вариантов.

1. Матрица А – квадратная и det(A)≠0.

2. А-квадратная, но det(A)=0.

3. А – прямоугольная.

Рассмотрим второй и третий случаи.

В первом случае часто решают через обратную матрицу x = A-1b.

Было вы не плохо чтобы для наших случаев решение можно было написать в аналогичном виде. Но для таких матриц обратной матрицы не существует. И для решения системы Ax = b под решением понимается вектор х, минимизирующий квадрат длины вектора невязки (r,r) = (Ax-b,Ax-b)→min. Если задача решается неоднозначно(как обычно бывает) то выбирают вектор наименьшей длины(norm(x)→min). При этом решение записывается в виде

x = A+b. Под A+ часто понимают псевдообратную матрицу Мура-Пенроуза X= A+, удовлетворяющую условиям:

· AXA=A;

· XAX=X;

· AX – симметричная;

· XA – симметричная.

Покажем, как получается матрица X:

Ax=b;

ATAx = ATb;

(ATA)-1ATAx = (ATA)-1ATb;

X = (ATA)-1AT;

A=U∑VT;

X = V∑+UT;

σ1+ 0 0 … 0 σ1+ 0 0 … 0

+ = (… …0 … 0) ∑+ = (… …0 … 0)

0 … σN+…0 0 … σN+…0

 

Где σi+ = 1/σ, если σ≠0 и 0, если σ=0. Часто для повышения надежности результатов используют эффективную псевдообратную матрицу A+ε у которой σε+ = 1/σ, если σ>ε иначе 0. Сингулярное разложение матрицы не является единственным и его не просто искать «руками», поэтому будем использовать пакет Matlab и его функцию [U, S, V] = svd(M);

Рассмотрим несколько примеров

Возьмем матрицу A mxn, где m>n

При решении «обычными» методами, мы бы получили, что решений нет.

Получим сингулярное разложение матрицы A

Теперь легко оценить число обусловленностей матрицы A

Cond(A) = 3.053

Получим ∑+ и A+

Получим наконец решение уравнения Ax=b

Проверим наш результат

Оценим погрешность наших вычислений

ε=1-norm(A*x)/norm(b) = 1-3.1543/3.3166 =1- 0.951 = 0.049

Проделаем те же действия но для матрицы, у которой m<n

При решении «обычными» методами, мы бы получили бесконечно много решений.

Получим сингулярное разложение матрицы A

Теперь легко оценить число обусловленностей матрицы A

Cond(A) = 3.053

Получим ∑+ и A+

Получим наконец решение уравнения Ax=b

Проверим наш результат

Оценим погрешность наших вычислений

ε=1-norm(A*x)/norm(b) 0

Проделаем те же действия но для квадратной матрицы 3x3, у которой пределитель = 0

При решении «обычными» методами, мы бы получили, что решений нет.

Получим сингулярное разложение матрицы A

Теперь легко оценить число обусловленностей матрицы A

Cond(A) = 2.99

Получим ∑+ и A+

Получим наконец решение уравнения Ax=b

Проверим наш результат

Оценим погрешность наших вычислений

ε=1-norm(A*x)/norm(b) 1-0.9914 = 0.0086


1 | 2 |

Поиск по сайту:



Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Студалл.Орг (0.012 сек.)