嗨,英格丽德,
我猜你的意思是你想计算te和对角线的交换矩阵包含
所有
非零元素。在这种情况下,我想这样的程序:
%定义一个
一个= [0 0.4 0.3 0.3;0 0 0.5 - 0.5;0.5 - 0.3 0 0.2;0.5 0.3 0.2 0);
你写的:
[n, n] =大小(一个)
。这是奇怪的(对我来说),我宁愿这样的程序:
[n m] =大小(一个);% * * *
断言(n = = m,“矩阵必须广场!”)% * * *
%找到所有的排列
B =重塑((:,flipud(烫发(1:4)。”),n, n的阶乘(n));
我需要使用内置的函数
诊断接头
,所以我不得不重命名
诊断接头
来
diag_
%计算对角金额与非零元素
diag_ = 0(1阶乘(n));%诊断接头函数具有相同名称的“阴影”
为i = 1:阶乘(n)
如果任何(诊断接头(B(:,:我))= = 0),继续;结束
diag_ (i) =跟踪(B(:,:我));
结束
diag_”
程序的错误是:
%定义一个
一个= [0 0.4 0.3 0.3;0 0 0.5 - 0.5;0.5 - 0.3 0 0.2;0.5 0.3 0.2 0);
[n, n] =大小(一个);
%找到所有的排列
B =重塑((:,flipud(烫发(1:4)。”),n, n的阶乘(n));
%计算对角金额与非零元素
诊断接头= 0(1阶乘(n));
为i = 1:阶乘(n)
为j = 1: n
如果B (j, j,我)= = 0
考虑这样一种情况:
B (:,:, j) = (1 0 0 0;0 2 0 0;0 0 0 0;0 0 0 0]
。在前两个迭代
B (j, j,我)~ = 0
,所以你分配
我诊断接头(1)
的痕迹
B(:,:,我)
。但在第三次迭代中,
B (3 3 i) = 0
:如果你打破这个循环,它是强制性的重置价值
诊断接头(j)
为0,否则它将存储错误的结果。
顺便说一句,在这个算法的计算跟踪无益地重复很多次。
我诊断接头(1)= 0;%这是失踪
打破
结束
诊断接头(1,)=跟踪(B(:,:我));
结束
结束
诊断接头
当然,两个版本完全相同的结果。
