一次同余式组的Matlab程序解法毕业论文.docVIP

本科生毕业论文 一次同余式组的Matlab程序解法 二级学院 ： 数学与计算科学学院 专 业 ： 数学与应用数学 年 级 ： 学 号 ： 作者姓名 ： 指导教师 ： 完成日期 ：  目录 １．引言 1 ２．一次同余式组的中国剩余定理解法 1 ３．一次同余式组的矩阵解法 3 ４．一次同余式组的逐步逼近解法 4 参考文献 6 附录 7 附录1　孙子解法Matlab程序 7 附录2　矩阵解法Matlab程序 10 附录3　逐步逼近解法Matlab程序 13 一次同余式组的Matlab程序解法 摘 要: 本文主要介绍一次同余式组的解法及其Matlab程序解，即中国剩余定理解法、矩阵解法和逐步逼近解法. 关键词: 中国剩余定理；矩阵解法；逐步逼近解法 Matlab program Methods for Congruences of Order One Abstract：In this paper, the main purpose is to introduce Matlab program methods for congruences of order one. These methods include Chinese Remainder Theorem, matrix and successive approximation solution． Key words：Chinese Remainder Theorem; matrix method; successive approximation １．引言 在世界数学的初期发展史上，中国古代的数学成就，尤其是同余方程方面的成就占有相当领先的地位．其中，一次同余式组是数论的重要组成部分，其求解更为引人注目.本论文介绍了一次同余式组的三种求解方法及Matlab求解程序．任意给出已知同余式，只要模间互素，程序可直接给出其特解．在我们知道其特解时，便可以直接得到一次同余式的一般解．用Matlab程序对一次同余式机械求解，直接快速，在通信密钥中，有着重要应用． ２．一次同余式组的中国剩余定理解法 定理1[1](大衍求一术) 一次同余式 　　　　　　　　　　　　　　　　　　(1) 的解，即不定方程 　　　　　　　　　　　　　　　　　　　(2) 中的值．(2)的解为 ，， 其中． ， ， ，是对应用Euclid算法所得的商． 定理2[2] (孙子定理) 设且则，同余方程组 （3） 都有解，且其解在模意义下是唯一的．同余方程组(3)的全部解为 ，　　　　　　　　　　　　　　(4) 其中，，，(叫做乘率)是同余方程 　　　　　　　　　　　　　　　　　(5) 的解. 算法步骤　已知模和余数，，求一次同余式组(3)的特解主要步骤是 步1．模间互素检查． 步2．求与． 步3．求对应的乘率． 步4．给出一个特解． 上述过程的Matlab程序见附录1． 例1　今有物不知数，三三数之剩二，五五数之剩三，七七数之剩二，问物几何． 解　设是所求物数，依题意得一次同余式组 执行程序 B=[2 3 2];N=[3 5 7];x=linear_congruence_equation_sunzi(B,N) 得到解 x=23． ３．一次同余式组的矩阵解法 定理3[3] 设同余式组(3)，如(5)，构造阶矩阵 ． 对A施行初等行变换，则矩阵A的某一行总可以变成(1,)的形式，则同余式组(3)的解为． 证明　因为，所以，存在整数使得 ， (6) 对应地， ． 　 (7) 给定，当时，，从而，由(6)式得有，(7)两边取模得 ， 从而满足同余式组(3)的每一个同余式，再由孙子定理的唯一性知，同余式组(3)有唯一解． 算法步骤　已知模和余数，，求一次同余式组(3) 的特解主要步骤是: 步1．模间互素检查． 步2．求与． 步3．把A的第一列中含有0的项剔除． 步4．对A以第１列元素由小到大进行排序． 步5．对A施行初等行变化(变换中必须也只能出现整数)，使得A的第一行为(1,)形式.． 步6．是特解． 上述过程的Matlab程序见附录2． 例2　求相邻四个整数，它们依次可被及整除． 解 设这四个相邻整数是依题意得 执行程序 B=[1 0 -1 -2];N=[4 9 25 49]; f=linear_congruence_equation_matrix(B,N) 得到一次同余式组的一个特解 f = -588051． ４．一次同余式组