IT ФЭЛ, кафедра МНЭ, 2 сем, Matlab / ЛР8 / lab8

m=input('Введите кол-во строк: ');

n=input('Введите кол-во столбцов/неизвестных: ');

A=zeros(m,n+1);

T=zeros(m,n);

for i=1:m

for j=1:(n+1)

A(i,j)=input('Введите элемент расширенной матрицы A: ');

end

end

disp(A)

for i=1:m

for j=1:n

T(i,j)=input('Введите элемент основной матрицы A: ');

end

end

disp(T)

U=rank(A); % ранг основной матрицы A

disp('Ранг расширенной матрицы A: ')

disp(U)

P=rank(T); % ранг расширенной матрицы A

disp('Ранг основной матрицы A: ')

disp(P)

if P~=U % по теореме Кронекера-Капелли

disp('Система несовместна, а значит не имеет решения')

else

disp('Решение СЛАУ: ')

for k=1:m

for j=(k+1):(n+1)

A(k,j)=A(k,j)/A(k,k);

end

for i=(k+1):m

for j=(k+1):(n+1)

A(i,j)=A(i,j)-A(k,j)*A(i,k);

end

end

end

x(m)=A(m,n+1);

disp(x(m))

for k=(n-1):-1:1

r=0;

for j=m:-1:(k+1)

r=r+A(k,j)*x(j);

end

x(k)=A(k,n+1)-r;

disp(x(k))

end

end