Определитель через рекурсию

@luivilla · Регистрация: 01.03.2015

Студворк — интернет-сервис помощи студентам

Требуется написать функцию для вычисления определителя. Я написала рекурсивный алгоритм с разложением по первому столбцу. Не работает. Точнее, выдает все время 0, кроме размера матрицы 1х1. Где я накосячила?

Код

function d1 = opred(A)
 [s,t] = size(A);
 if s==1 
  d1=A(1,1);
 else
  d1=0;
  for i=1:s
    B=A;
    B(:,1)=[];
    B(i,:)=[];
    d1=d1+opred(B)*A(1,i)*(-1)^(i+1);
  end
 end
end

Зосима · 12.10.2015, 09:36

хм... у меня считает, но не верно

Matlab M

>> A = randn(3)
 
A =
 
    0.5377    0.8622   -0.4336
    1.8339    0.3188    0.3426
   -2.2588   -1.3077    3.5784
 
>> opred(A)
 
ans =
 
   -1.9943
 
>> det(A)
 
ans =
 
   -4.7434

@luivilla · 12.10.2015, 15:32 **[ТС]**

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

@kRosis · 12.10.2015, 17:14

Час промучался, пока на листочке не посчитал :\ А ошибка плевая - перепутаны строки и столбцы.

Matlab M

function d = opred(A)
[s,t] = size(A);
if s == 1
    d = A(1,1);
% elseif s == 2
%     d = A(1,1)*A(2,2)-A(1,2)*A(2,1);
else
    d = 0;
    for i = 1:s
        B = A;
        B(1,:) = [];
        B(:,i) = [];
        d = d+opred(B)*A(1,i)*(-1)^(i+1);
    end;
end;
end

@luivilla · 12.10.2015, 19:20 **[ТС]**

мда, я просто дятел: сама же написала про столбец, а умножала на строку.

@luivilla 3 / 3 / 0 Регистрация: 01.03.2015 Сообщений: 122
		1
	Определитель через рекурсию 10.10.2015, 12:41. Показов 2846. Ответов 4 Метки нет (Все метки) Требуется написать функцию для вычисления определителя. Я написала рекурсивный алгоритм с разложением по первому столбцу. Не работает. Точнее, выдает все время 0, кроме размера матрицы 1х1. Где я накосячила? Код function d1 = opred(A) [s,t] = size(A); if s==1 d1=A(1,1); else d1=0; for i=1:s B=A; B(:,1)=[]; B(i,:)=[]; d1=d1+opred(B)A(1,i)(-1)^(i+1); end end end 0

Зосима 5221 / 3552 / 372 Регистрация: 02.04.2012 Сообщений: 6,458 Записей в блоге: 17
	12.10.2015, 09:36	2
	хм... у меня считает, но не верно Matlab M >> A = randn(3) A = 0.5377 0.8622 -0.4336 1.8339 0.3188 0.3426 -2.2588 -1.3077 3.5784 >> opred(A) ans = -1.9943 >> det(A) ans = -4.7434 0

@luivilla 3 / 3 / 0 Регистрация: 01.03.2015 Сообщений: 122
	12.10.2015, 15:32 [ТС]	3
	Да, сорри, тоже считает, но тем не менее, не верно. Где может быть ошибка. Просто не хочется переделывать алгоритм на принципиально другой. По идее, этот же должен работать. 0

@luivilla 3 / 3 / 0 Регистрация: 01.03.2015 Сообщений: 122
	12.10.2015, 19:20 [ТС]	5
	мда, я просто дятел: сама же написала про столбец, а умножала на строку. 0

Опции темы

Определитель через рекурсию

Решение