41900

ИТЕРАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ СИСТЕМ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ

Лабораторная работа

Математика и математический анализ

Цель работы: научиться решать системы линейных алгебраических уравнений СЛАУ методом простых итераций МПИ и методом Зейделя с помощью ЭВМ. Изучить метод простых итераций и метод Зейделя для решения СЛАУ. Сравнить скорости сходимости метода простых итераций и метода Зейделя. Построить рабочие формулы МПИ и метода Зейделя для численного решения системы.

Русский

2013-10-26

244.14 KB

29 чел.

ЛАБОРАТОРНАЯ РАБОТА №5.

ИТЕРАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ

СИСТЕМ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ

Вариант №5.

Выполнил:

Студент группы 24275

Кожевников Е.И.

Проверил:

Доцент

Горбунов Д.В.

Цель работы: научиться решать системы линейных алгебраических уравнений (СЛАУ) методом простых итераций (МПИ) и методом Зейделя с помощью ЭВМ.

Содержание работы:

1. Изучить метод простых итераций и метод Зейделя для решения СЛАУ.

2. На конкретном примере усвоить порядок решения СЛАУ с помощью ЭВМ указанными методами.

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

Задание.

1. Аналитически решить СЛАУ вида:

2. Построить рабочие формулы МПИ и метода Зейделя для численного решения системы.

3. Составить программу на любом языке программирования, реализующую построенные итерационные процессы.

Решение.

 

Аналитический метод.

Аналитическим решением системы являются значения:

Метод простых итераций. 

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

Разделив каждое уравнение системы на соответствующий диагональный коэффициент, сформируем столбец в левой части и перенесем остальные слагаемые в правую часть и получим рабочие формулы МПИ вида:

 

где k=0,1,2…

 

Начальное приближение обычно выбирают равным столбцу свободных членов преобразованной системы . Процесс заканчивается при одновременном выполнении трех условий:

, , .

В этом случае значения являются приближенными значениями решения СЛАУ.

Метод Зейделя. Более быструю скорость сходимости имеет метод Зейделя, в котором найденное -е приближение сразу же используется для получения -го приближения последующих координат (Рис.1).

Рис.1

Рабочие формулы метода Зейделя запишутся так:

где k=0,1,2…

 

Условия выхода итерационного процесса и выбор начального приближения аналогичны МПИ.

Блок-схема метода простых итераций и метода Зейделя приведена на рисунке 2.

Рис.1 Схема МПИ и метода Зейделя.

Тексты программ:

  1. МПИ:

Program P5;

uses Crt;

var n:integer;

     x0,x,xn,y0,y,yn,z0,z,zn,eps,rx,ry,rz:real;

begin

clrscr;

n:=0; x0:=1; x:=x0; y0:=-0.6; y:=y0; z0:=0.2; z:=z0; eps:=0.001;

writeln ('  n   x(i)   x(i-1)  RAZ(x)   y(i)   y(i-1)  RAZ(y)   z(i)   z(i-1)  RAZ(z)');

repeat

 xn:=1-(1/6)*y+(5/12)*z;

 rx:=x;

 yn:=-0.6+0.6*x+0.1*z;

 ry:=y;

 zn:=0.2-0.2*x+0.2*y;

 rz:=z;

 writeln (n:3, x:8:4, xn:8:4, (xn-x):8:4, y:8:4, yn:8:4, (yn-y):8:4, z:8:4, zn:8:4, (zn-y):8:4);

 x:=xn;

 y:=yn;

 z:=zn;

 n:=n+1;

until (abs(x-rx)<=eps) and (abs(y-ry)<=eps) and (abs(z-rz)<=eps);

readln;

end.

  1. Метод Зейделя:

Program P5; {Zeydel}

uses Crt;

var n:integer;

     x0,x,xn,y0,y,yn,z0,z,zn,eps,rx,ry,rz:real;

begin

clrscr;

n:=0; x0:=1; x:=x0; y0:=-0.6; y:=y0; z0:=0.2; z:=z0; eps:=0.001;

writeln ('  n   x(i)   x(i-1)  RAZ(x)   y(i)   y(i-1)  RAZ(y)   z(i)   z(i-1)  RAZ(z)');

repeat

 xn:=1-(1/6)*y+(5/12)*z;

 rx:=x;

 yn:=-0.6+0.6*xn+0.1*z;

 ry:=y;

 zn:=0.2-0.2*xn+0.2*yn;

 rz:=z;

 writeln (n:3, x:8:4, xn:8:4, (xn-x):8:4, y:8:4, yn:8:4, (yn-y):8:4, z:8:4, zn:8:4, (zn-y):8:4);

 x:=xn;

 y:=yn;

 z:=zn;

 n:=n+1;

until (abs(x-rx)<=eps) and (abs(y-ry)<=eps) and (abs(z-rz)<=eps);

readln;

end.


Результаты отработки программы:

  1.  Рис.2 – программы, работающей по методу простых итераций;
  2.  Рис.3 – программы, работающей по методу Ньютона.

Рис.2. Ответ: х(10)≈0.9999, y(10)≈-0.0002, z(10)≈0.0002.

Рис.2. Ответ: х(5)≈1.0001, y(5)≈0.0000, z(5)≈-0.0000.


 

А также другие работы, которые могут Вас заинтересовать

16249. Конфигурация функции IGMP Snooping 724.66 KB
  Лабораторная работа №1 Конфигурация функции IGMP Snooping 1 Цель работы 1.1Научиться конфигурировать протокол управления групповой multicast рассылкой на коммутаторах Dlink. 2 Литература 2.1 Смирнова Е.В. Пролетарский А.В. Баскаков И.В. Федотов Р.А. Построение комму
16250. Создание профиля многоадресной рассылки 724.41 KB
  Лабораторная работа №2 Создание профиля многоадресной рассылки 1 Цель работы 1.1 Научиться управлять multicast трафиком. 2 Литература 2.1 Смирнова Е.В. Пролетарский А.В. Баскаков И.В. Федотов Р.А. Построение коммутируемых компьютерных сетей: учебное пособие. –...
16252. Настройка QoS 93 KB
  Лабораторная работа №3 Настройка QoS. Цель работы: 2.1 Изучить настройку приоритизации трафика управления полосой пропускания на коммутаторах DLink.Исследовать эффективность работы приоритизации. 2 Литература: 2.1 Смирнова Е.В. Пролетарский А.В. Баскаков И.В.
16253. Организация IP вещания 852.79 KB
  Лабораторная работа №1 Организация IP вещания Цель работы Исследовать структуру сети IP вещания. Научиться настраивать IP стример. Научиться настраивать IP STB. Литература 2.1 Техническое описание приборов. Подготовка к р...
16254. Изучение головной станции кабельного телевидения 2.42 MB
  Лабораторная работа №1 Лабораторная работа №1 Изучение головной станции кабельного телевидения 1 Цель работы: 1.1 Изучить принципы построения головной станции кабельного телевидения. 1.2 Научиться производить настройку головной станции кабельного телевидения. ...
16255. Исследование структуры и принципа действия конвертора телевизионных сигналов 856.5 KB
  Лабораторная работа №2 Исследование структуры и принципа действия конвертора телевизионных сигналов 1 Цель работы: 1.1 Изучить структурную схему конвертора СТ07. 1.2 Научиться производить настройку конвертора СТ07. 2 Литература: 2.1 Джакония В.Е. Телевидение. М.:
16256. Исследование спектра сигнала кабельной сети 991 KB
  Лабораторная работа №4 Исследование спектра сигнала кабельной сети 1 Цель работы: 1.1 Исследовать спектр сигнала Архангельской телевизионной компании. 1.2 Научиться пользоваться анализатором спектра DL4. 1.3 Научиться пользоваться программным обеспечением S.M.A.R.T. ...
16257. Удосконалення бухгалтерської звітності підприємства 473.5 KB
  Предметом дослідження є сукупність теоретичних, методологічних і організаційно економічних питань з удосконалення бухгалтерської звітності в ППА «Коровинці» Недригайлівського району Сумської області