求解线性方程组.docxVIP

求解线性方程组 求解线性方程组 求解线性方程组 摘要：自从1946年世界上第一台计算机诞生，计算机在我们的生活作用也越来越重要，它 多样化的功能与日趋简便的操作使得利用计算机求解算数问题成为一种趋势。线性方程组求 解计算量大过程复杂，所以我们小组通过C++软件编写求解线性方程组程序。该程序功能强 大，具有求解线性方程组、三角函数、对数函数、指数函数的功能。而且程序加了新的元素，全数字化控制，使用简单方便。 关键词：高斯消元法；C语言；线性方程组；常用函数求解。 1.1 研究背景 随着人们对计算类软件的需求，求解线性方程组成为常见的应用。线性方程组求解未知 数多，方程复杂，人工计算较难实现。为解决这类问题，本小组通过微软公司的VC6.0软 件设计一款用C语言求解线性方程组的软件，并能进行常用函数的辅助运算。 1.2主要工作 一个软件开发需要大量人力、物力支持，还需要合理有序的工作分配，明确一个软件必 须要有哪些功能，可以添加哪些功能。线性方程组求解程序有求解功能和辅助功能。因此主 要工作是深刻理解高斯消元法原理和初等矩阵，掌握矩阵变换步骤，常用函数的原理；再利 用C语言编写一系列代码，解决方程组求解问题。 1.3本文结构 本文第1部分绪论主要说明了本文的研究背景、主要工作和实现方法。第2部分相关知 识介绍。第3部分系统设计，阐述了需求分析、系统模型设计和系统工作流程设计。第4 部分系统实现介绍了系统各个功能模块，并描述了系统功能图。第5部分软件的调试和验证，进行截图。第六部分总结和心的体会。 以下就线形方程求解过程中和其它的相关知识进行间要介绍。 2.1高斯消元法原理 2.1.1初等矩阵 (1)单位矩阵第i,j两行（列）互换得到的方阵为Pij。将矩阵B的第i，j两行（列)互换 所得矩阵B1，即有PijB=B1 （2）单位矩阵第i行（列）乘以常数k得到初等方阵Di(k),将矩阵B的第i行（列）乘 以k得到矩阵B2，即有B2=Di(k)B。 （3）将单位矩阵的第j行（列）的k倍加到第i行（列）得到初等方阵Tij(k)，矩阵B 的第j行（列）的k倍加到第i行（列）得到矩阵B3，即有B3=Tij(k)B。矩阵B的第i列的k 倍加到第j列得到矩阵B3，即有B3=BTij(k)。 2.1.2高斯消元法 消元法可用来找出下列方程组的解或其解的限制： 2a?b?x?8...X() -3a - b ? 2c ? -11 .(.Y. ) ?2a?b?2c??3...(Z) 这个算法的原理是： 首先，要将X 以下的等式中的a 消除，然后再将Y 以下的等式中的b 消除。这样可 使整个方程组变成一个三角形似的格式。之后再将已得出的答案一个个地代入已被简化的等 式中的未知数中，就可求出其余的答案了。 在刚才的例子中，将3/2 X和Y相加，就可以将Y 中的a 消除了。然后再将X 和Z 相加，就可以将Z 中的a消除。 可以这样写：Y + 3/2 X→ Y ；Z+ X → Z 结果就是：2a + b - c = 8；1/2 b + 1/2c = 1 ；2b + c = 5 现在将 ? 4Y 和Z 相加，就可将Z 中的b 消除：Z + (-4Y) → Z 2a?b?c?8 1/2b ? 1/2c ? 1 这样就完成了整个算法的初步，一个三角形的格式(指:变量的格式而言,上例中的变量各 为(3,2,1个)出现了。 第二步，就是由尾至头地将已知的答案代入其他等式中的未知数。第一个答案就是： 然后就可以将c 代入Y 中，立即就可得出第二个答案： 之后，将c 和b 代入X 之中，最后一个答案就出来了： 3.1需求分析 经过对程序设计题目的分析可知，整个程序的设计实现大致分为三个模块，分别是：输 入方程组，计算方程组，继续运算/退出。要求以菜单的形式调用各功能函数：程序可读性 强，界面新意，输出形式美观。 3.2软件流程设计 打开线性方程求解系统软件，进入主界面： （1) 输入“1”进入线性方程求解系统模块，输入线性方程组相关系求解，按任意键返回主界 （2）输入“2”进入辅助系统模块，可进行三角函数，ex函数，lgx函数，lnx函数求解，以满 足不同 算函数计算的需求； （3）输入“0”退出线性