用Matlab构造奇数阶幻方矩阵.pdfVIP

第 3O卷 第 1期 高 师 理 科 学 刊 Vo1．30 No．i 2010钲 1月 JournalofScienceofleachers CollegeandUniversity Jan． 2010 文章编号 ：1007—9831(2010)01—0015—03 用Matlab构造奇数阶幻方矩阵 陈豫眉，谭代伦 (西华师范大学 数学与信息学院，四川 南充 637002) 摘要：幻方矩阵因其诸多神奇性质而受到广泛关注．介绍一种利用矩阵加法构造奇数阶幻方的方 法，并用Matlab程序将其实现．在此程序基础上，给出一种改进算法，可以生成同阶新的奇数阶 幻方矩 阵． 关键词：奇数阶幻方；矩阵加法；Matlab 中图分类号：0151．21 文献标识码：A n阶幻方是指由n个数构成的方阵，要求其每行、每列与2条对角线上元素之和相等．此元素之和称 为幻和，并记为H ，易知日=堕 ． 通常，方阵的元素由1，2，…， 的正整数组成，而且每个元素 仅能使用一次．这样的幻方通常称为古典幻方 “．构造奇数阶幻方的方法有调数法、先定对角线法、凸十 字形补空法-q对角圆筒填写法等拉． 1 构造法 矩阵加法构造幻方是指先按某种规则构造几个n×n矩阵，然后把这几个矩阵加起来即成 l／阶幻方．按 照参考文献[2]dO方法，对于I／：阶奇数阶幻方M ，将幻方各边中点依次连接而得到一个斜放的正方形 (每 边上有 个数)，叫做这个n阶幻方的菱形块，记为D．在菱形块上共有 ． 个数，这个n阶奇数阶 Z 幻方 的其余 ． 个数分布在 4个角上，每个角上有 ，儆 ，都呈三角形，将左上角、右上角、 O 左下角、右下角各块依次记为C，C ，C ，c ． 为了构造奇数阶幻方M ，需要先构造 2个，z×n矩阵A和曰．先构造A矩阵，令A中第 i行第 ‘『列元 ，， 1、2 素为a(i，)=(，z+1)i一(n一1)j+ — ．矩阵B由分块矩阵D，Cl，c2，C3，C构成．将D，c。，c2，c3，c 中所有元素分别取为一1，n一1，n。一1，一n 一1，一Ft一1．则奇数阶幻方M =A+B．例如：当，z=5时， 4 4 —1 24 24 4 — 1 一l 一1 24 B = — 1 — 1 一l 一 1 —1 ， 于是，可得矩阵 — 26 —1 —1 —1 ～6 — 26 —26 —1 —6 ～6 收稿 日期：2009—09—13 基金项目：四川省教育厅青年基金项目 (07ZB087，07ZB154)；西华师范大学校极科研项 日(2008—2010) 作者简介：陈豫眉 (1971一)，女 ，四川眉山人，副教授 ，硕士，从事微分方程数值解及数学建模研究．E-mail l6 高 师 理 科 学 刊 第 30卷 l0 6 2 —2 —61 f4 4 —1 24 5 l6 l2 8 4 0 l l4 一l 一1 —1 = 22 18 l4 l0 6 』+l一1 —1 —1 —1 28 24 20 l6 12Il一26 —1 —1 —1 A + 34 30 26 22 18J I一26 —26 —1 —6 8 显然，矩阵M5为5阶幻方，其幻和等于 65． ，，，， ．．． ．．．．． ．．．． ．．． ．I．_．。．。． ．．．．