Решение систем, линейных уравнений |
Решение систем, линейных уравнений |
Altair |
7.10.2004 16:01
Сообщение
#1
|
Ищущий истину Группа: Модераторы Сообщений: 4 824 Пол: Мужской Реальное имя: Олег Репутация: 45 |
Решение систем линейных уравнений методом Гаусса
uses crt; -------------------- Помогая друг другу, мы справимся с любыми трудностями!
"Не опускать крылья!" (С) |
volvo |
28.01.2005 14:36
Сообщение
#2
|
Гость |
Численное решение систем линейных уравнений методом Зейделя
Const |
volvo |
21.03.2005 10:46
Сообщение
#3
|
Гость |
Решение систем линейных уравнений методом Крамера
Позволяет решить неоднородную систему n линейных алгебраических уравнений с n неизвестными, если определитель основной матрицы не равен нулю. {$N+} |
hiv |
22.03.2005 18:57
Сообщение
#4
|
Профи Группа: Пользователи Сообщений: 660 Пол: Мужской Реальное имя: Михаил Репутация: 11 |
Решение системы линейных уравнений методом Холецкого
Этот метод менее подвержен накопительным ошибкам в сравнении с методом Гаусса. Поэтому им можно решать системы из многих десятков уравнений! Код program HOLETSKY; {$M 65520, 65520, 65520} const max_eq = 50; type RealType=double; MCvector = array [0..max_eq] of RealType; {вектор} MCmatr = array [0..max_eq,0..max_eq] of RealType; {запись типа матрица} function maxval(x,y:realtype):realtype; begin if x<y then maxval:=y else maxval:=x; end; procedure Holetskiy(N:integer;A0:MCmatr; B0:MCvector; var X:MCvector; var error:RealType); var i,j,m,k,p :byte; s,e,max :RealType; A :MCmatr; begin if (N>max_eq)or(N<2) then begin writeln('Ошибка в передаваемых данных.'); halt(1); end; A:=A0; for i:=0 to N-1 do A[i,N]:=B0[i]; for i:=0 to N-1 do {выбор главного элемента по диагонали} begin max:=0; for j:=i to N-1 do if max<abs(A[i,j]) then begin max:=abs(A[i,j]); p:=j; end; if p<>i then begin for j:=0 to N do A[N,j]:=A[i,j]; for j:=0 to N do A[i,j]:=A[p,j]; for j:=0 to N do A[p,j]:=A[N,j]; end; end; {метoд хoлецкoгo} for j:=1 to N do A[0,j]:=A[0,j]/A[0,0]; for m:=1 to N-1 do begin for i:=m to N-1 do begin s:=0; for k:=0 to m-1 do s:=s+A[i,k]*A[k,m]; A[i,m]:=A[i,m]-s end; for j:=m+1 to N do begin s:=0; for k:=0 to m-1 do s:=s+A[m,k]*A[k,j]; A[m,j]:=(A[m,j]-s)/A[m,m] end; end; X[N-1]:=A[N-1,N]; for i:=N-2 downto 0 do begin s:=0; for k:=i+1 to N-1 do s:=s+A[i,k]*X[k]; X[i]:=A[i,N]-s end; {расчет ошибки} error:=0; for i:=0 to N-1 do begin s:=0; for j:=0 to N-1 do if i=j then s:=s+(A0[i,j]+1)*X[j] else s:=s+A0[i,j]*X[j]; s:=s-B0[i]; error:=maxval(error,abs(X[i]-s)); end; end; var A :MCmatr; B,X :MCvector; err :RealType; n,i,j :integer; begin {of program} writeln('Решение системы линейных уравнений метoдoм Хoлецкoгo'); writeln('с выбoрoм главнoгo элемента.'); write('Введите числo уравнений = ');read(N); for i:=0 to N-1 do begin writeln('Введите кoэффициенты ',i:0,'-го уравнения:'); for j:=0 to N-1 do begin write('A[',i:0,',',j:0,'] = '); readln(A[i,j]); end; write('B[',i:0,'] = '); readln(B[i]); end; Holetskiy(N,A,B,X,err); writeln('Решение уравнения:'); for i:=0 to N-1 do writeln('X[',i:0,'] = ',X[i]); writeln('Точность резутьтата =',err); end. -------------------- Никогда не жадничай. Свои проблемы с любовью дари людям!
|
Текстовая версия | 27.04.2024 19:19 |