创建6阶幻方矩阵A MATLAB,用Matlab构造奇数阶幻方矩阵

第30卷 第1期 高 师 理 科 学 刊 Vol. 30 No.1 2010 年 1 月 Journal of Science of Teachers′College and University Jan. 2010 文章编号：1007-9831(2010)01-0015-03 用 Matlab 构造奇数阶幻方矩阵 陈豫眉，谭代伦 (西华师范大学 数学与信息学院，四川 南充 637002) 摘要：幻方矩阵因其诸多神奇性质而受到广泛关注．介绍一种利用矩阵加法构造奇数阶幻方的方法，并用 Matlab 程序将其实现．在此程序基础上，给出一种改进算法，可以生成同阶新的奇数阶幻方矩阵． 关键词：奇数阶幻方；矩阵加法；Matlab 中图分类号：O151.21 文献标识码：A n阶幻方是指由 2 n 个数构成的方阵，要求其每行、每列与 2 条对角线上元素之和相等．此元素之和称 为幻和，并记为 n H ，易知 2 )1( 2 + = nn Hn ．通常，方阵的元素由 2 , ,2 ,1 n" 的正整数组成，而且每个元素 仅能使用一次．这样的幻方通常称为古典幻方[1]．构造奇数阶幻方的方法有调数法、先定对角线法、凸十字形补空法与对角圆筒填写法等[2-4]． 1 构造法 矩阵加法构造幻方是指先按某种规则构造几个 nn× 矩阵，然后把这几个矩阵加起来即成n阶幻方．按照参考文献[2]中方法，对于n阶奇数阶幻方 n M ，将幻方各边中点依次连接而得到一个斜放的正方形(每 边上有 2 1+n 个数)，叫做这个n阶幻方的菱形块，记为 D．在菱形块上共有 2 12 +n 个数，这个n 阶奇数阶 幻方 n M 的其余 2 12 −n 个数分布在 4 个角上，每个角上有 8 12 −n 个数，都呈三角形，将左上角、右上角、 左下角、右下角各块依次记为 4321 , , , C CC C ． 为了构造奇数阶幻方 n M ，需要先构造 2 个 nn × 矩阵 A 和 B ．先构造 A 矩阵，令 A 中第i行第 j 列元 素为 2 )1( )1()1(),( 2− +−−+= n jninjia ．矩阵B 由分块矩阵 4321 , , , , C CC CD 构成．将 4321 , , , , C CC CD 中所有元素分别取为 1 ,1 ,1 ,1 ,1 22 −−−−−−− nnnn ．则奇数阶幻方 BAM += n ．例如：当 5=n 时， ⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎜ ⎝ ⎛ −−−−− −−−−− −−−−− −−− − = 6612626 611126 11111 241114 2424144 B ，于是，可得矩阵 收稿日期：2009-09-13 基金项目：四川省教育厅青年基金项目(07ZB087，07ZB154)；西华师范大学校极科研项目(2008-2010) 作者简介：陈豫眉(1971-)，女，四川眉山人，副教授，硕士，从事微分方程数值解及数学建模研究．E-mail：xhshuxue@163.com 16 高 师 理 科 学 刊 第 30 卷 ⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎜ ⎝ ⎛ = ⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎜ ⎝ ⎛ −−−−− −−−−− −−−−− −−− − + ⎟ ⎟ ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎜ ⎜ ⎝ ⎛ −− =+= 12162548 61519232 59131721 24371120 182211014 6612626 611126 11111 241114 2424144 1822263034 1216202428 61