matlab magic

最新推荐文章于 2021-03-17 23:51:10 发布

weixin_34267123

最新推荐文章于 2021-03-17 23:51:10 发布

阅读量1k

点赞数

文章标签： matlab python

原文链接：https://my.oschina.net/donngchao/blog/522100

版权

2019独角兽企业重金招聘Python工程师标准>>>

for n = 8:11
    subplot(2,2,n-7)
    surf(magic(n))
    title(num2str(n))
    axis off
   view(30,45)
    axis tight
end

其中，生成幻方的matlab程序如下

function M = magicV2(n)

%magicV2  magicV2 square.

%   magicV2(N) is an N-by-N matrix constructed from the integers

%   1 through N^2 with equal row, column, and diagonal sums.

%   Produces valid magicV2 squares for all N > 0 except N = 2.

 

%   Copyright 1984-2002 The MathWorks, Inc.

%   $Revision: 5.15 $  $Date: 2002/04/15 03:44:23 $

 

% Historically, MATLAB's magicV2 was a built-in function.

% This M-file uses a new algorithm to generate the same matrices.

 

n = floor(real(double(n(1))));

 

% Odd order.

if mod(n,2) == 1

   [J,I] = meshgrid(1:n);

   A = mod(I+J-(n+3)/2,n);

   B = mod(I+2*J-2,n);

   M = n*A + B + 1;

 

% Doubly even order.

elseif mod(n,4) == 0

   [J,I] = meshgrid(1:n);

   K = fix(mod(I,4)/2) == fix(mod(J,4)/2);

   M = reshape(1:n*n,n,n)';

   M(K) = n*n+1 - M(K);

 

% Singly even order.

else

   p = n/2;

   M = magicV2(p);

   M = [M M+2*p^2; M+3*p^2 M+p^2];

   if n == 2, return, end

   i = (1:p)';

   k = (n-2)/4;

   j = [1:k (n-k+2):n];

   M([i; i+p],j) = M([i+p; i],j);

   i = k+1;

   j = [1 i];

   M([i; i+p],j) = M([i+p; i],j);

end

不同的矩阵的秩产生不同的幻方矩阵。

转载于:https://my.oschina.net/donngchao/blog/522100