Метод _прогонки_MatLab (Отчет)

Посмотреть архив целиком








Лабораторная работа «Решение системы линейных алгебраических уравнений с трехдиагональной матрицей методом прогонки»

Студента Мыльникова Юрия

группы МТ10-32































Теоретическая часть



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



clc,clear all

A=[8 -2 0 0; -1 6 -2 0 ;0 2 10 -4 ; 0 0 -1 6]

b=[6 3 8 5]'

x = A\b

[n,n]=size(A);


P(1)=1;

Q(1)=1;

for i=2:n-1

zn=-A(i,i)-A(i,i-1)*P(i-1);

P(i)=A(i,i+1)/zn;

Q(i)=(A(i,i-1)*Q(i-1)-b(i))/zn;

end

Q(n)=(A(n,n-1)*Q(n-1)-b(n-1))/A(n,n)-A(n,n-1)*P(n-1);

x(n)=Q(n);

for i=n-1:-1:1

x(i)=P(i)*x(i+1)+Q(i);

end;



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








A =

8 -2 0 0

-1 6 -2 0

0 2 10 -4

0 0 -1 6

b =

6

3

8

5

x =

1

1

1

1

2014 год




Случайные файлы

Файл
115375.rtf
176608.rtf
166127.rtf
14616.rtf
169319.rtf