matlab编写魔术公式

魔术公式指的是在一个正方形的格子内填入数值，使得每一行、每一列以及每一条对角线的和都相等。在Matlab中，可以使用循环语句和矩阵运算来编写魔术公式的程序。以下为实现步骤：

1. 定义正方形的长度，例如n=5，表示一个5*5的正方形。
2. 构造一个n*n的矩阵，用0填充每个位置。
3. 定义起始位置，例如从矩阵的中心位置(ceil(n/2),ceil(n/2))开始填数。
4. 根据魔术公式，将首个数字填在起始位置的上一行、左一列处。同时记录下一个填数的位置为(1, n)。
5. 进入循环，从2填到n*n数为止。每次填数都要先检查下一个位置是否超出矩阵范围，若超出则跳出循环。若未超出，则判断该位置是否已有数字，若有则向下移一行。若没有，则在该位置填数，将下一个填数的位置左移一列，上移一行。
6. 输出填好数字的矩阵。

Matlab代码如下：

n = 5;
A = zeros(n,n);
i = ceil(n/2);
j = ceil(n/2);
A(i-1, j) = 1;
j = n;
for k = 2:n*n
if i == 1 && j == n
i = 2;
j = n-1;
elseif i == 1
i = n;
elseif j == n
i = i-1;
j = 1;
elseif A(i-1, j+1) ~= 0
i = i+1;
else
i = i-1;
j = j+1;
end
A(i,j) = k;
end
disp(A);

运行程序后，输出一个5*5的魔术正方形矩阵：

17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9

此矩阵的每行、每列和两条对角线的和均为65。