如何用MATLAB计算矩阵的行列式利用矩阵的Laplace展开定理,而非用det命令求得.请教M程序该如何编写?如下,请指教function y=detnn(A)y=0;[n,n]=size(A);if n==3 y=det33(A);endfor i=1:n y=y+(-1)^(i+1)*A(1,i)*detnn(A(2

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/05 03:01:43
如何用MATLAB计算矩阵的行列式利用矩阵的Laplace展开定理,而非用det命令求得.请教M程序该如何编写?如下,请指教function y=detnn(A)y=0;[n,n]=size(A);if n==3    y=det33(A);endfor i=1:n    y=y+(-1)^(i+1)*A(1,i)*detnn(A(2

如何用MATLAB计算矩阵的行列式利用矩阵的Laplace展开定理,而非用det命令求得.请教M程序该如何编写?如下,请指教function y=detnn(A)y=0;[n,n]=size(A);if n==3 y=det33(A);endfor i=1:n y=y+(-1)^(i+1)*A(1,i)*detnn(A(2
如何用MATLAB计算矩阵的行列式
利用矩阵的Laplace展开定理,而非用det命令求得.请教M程序该如何编写?
如下,请指教
function y=detnn(A)
y=0;
[n,n]=size(A);
if n==3
y=det33(A);
end
for i=1:n
y=y+(-1)^(i+1)*A(1,i)*detnn(A(2:n, [1:(i-1) (i+1):n]));
end

如何用MATLAB计算矩阵的行列式利用矩阵的Laplace展开定理,而非用det命令求得.请教M程序该如何编写?如下,请指教function y=detnn(A)y=0;[n,n]=size(A);if n==3 y=det33(A);endfor i=1:n y=y+(-1)^(i+1)*A(1,i)*detnn(A(2
clear
clc
n0=20;
e=1e0;
r0=(1:n0);
p0=poly(r0);
t=1;
for k=14:2:20
p=p0;
p(n0-k+1)=p(n0-k+1)+e;
r(:,k)=roots(p);
subplot(2,2,t)
plot(real(r(:,k)),imag(r(:,k)),'*')
xlabel('根的实部')
ylabel('根的虚部')
grid on
t=t+1;
end