作业帮如何用MATLAB计算矩阵的行列式
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/05/01 18:28:44
如何用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
利用矩阵的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
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
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