数值线性代数课程设计_代数

《数值代数》课程设计

评分标准 (1)

交作业的内容 (1)

附录1 论文结构撰写规范 (2)

附录2 (2)

参考论文1 (2)

参考论文2 (13)

参考论文3 (16)

参考论文4 (21)

1. 2-3两天查资料;

2. 1-2天论文构思,列出提纲;

3. 2-3确定计算方法和计算程序的详细内容和设计步骤；

写课程设计报告；

论文内容: 问题背景介绍, 建立正确的数学模型, 求解模型的数学原理, 计算过程,编写计算程序, 解释计算结果.

论文字数: 1800-2200字. 注: 字数是指正文的字数.

标题、摘要、关键字、目录、图表说明、参考文献、程序等不算入字数统计.

论文格式: 5号字体,A4版面, 上下左右各留2cm,单倍行距,页面设置选为无网格; 编页码;

图形尺寸：长度不超过4cm.一行可以放若干张图. 文稿中数据保留3位小数点.

（1）要有勤于思考、刻苦钻研的学习精神和严肃认真、一丝不苟、

有错必改、精益求精的工作态度，对有抄袭他人设计论文或找他人代设计、

代做论文等行为的弄虚作假者一律按不及格记成绩。

（2）要敢于创新，勇于实践，注意培养创新意识。

（3）掌握课程的基本理论和基本知识，概念清楚，设计计算正确，

结构设计合理，实验数据可靠，软件程序运行良好，绘图符合标准，

论文撰写规范，答辩中回答问题正确。

附录1 论文结构撰写规范

（1）题目、院系、班级、学生姓名。

（2）摘要

摘要是论文内容的简短陈述，一般不超过150字。

关键词应为反映论文主题内容的通用技术词汇，一般为4个左右，一定要在摘要中出现。

（3）正文

正文应按目录中编排的章节依次撰写，要求计算正确，论述清楚，文字简练通顺，

插图简明，书写整洁。文中图、表及公式不能徒手绘制和书写。

（4）参考文献（资料）

参考文献必须是学生在课程设计中真正阅读过和运用过的，文献按照在正文中的出现顺序排列。

各类文献的书写格式如下：

a．图书类的参考文献

序号作者名·书名·（版次）·出版单位，出版年：引用部分起止页码。

b．翻译图书类的参考文献

序号作者名·书名·译者·（版次）出版单位，出版年：引用部分起止页码。

c．期刊类的参考文献

序号作者名·文集名·期刊名·年，卷（期）：引用部分起止页码。

附录2

数值代数课程设计题

参考论文1

《数值线性代数》课程设计

信息与计算科学0440501106 蔡云

摘要：

Gauss消去法是目前求解中小规模线性方程组（即阶数不要太高，例如不超过1000）常用的方法，它一般用于系数矩阵稠密（即矩阵的绝大多数元素都是非零的）而没有任何特殊结构的线性方程组。选主元Gauss消去法弥补了不选主元的Gaauss消去法的不足，但是也付出了昂贵的代价。平方根是求解对称正定线性方程组最常用的方法之一。

关键词：Gauss 消去法，选主元Gauss 消去法，平方根法

正文：

一，问题背景：

在古代数学和近代数学数值计算和工程应用中，求解线性方程组都是重要的课题。西方的Gauss 消元法在求解未知数多的大型线形方程组则是在20世纪中叶电子计算机问世后才成为现实。这次实验将充分运用各种数值计算方法，及计算机的强大数值计算功能来解决所遇到的问题。

二，问题提出：

1.先用你所熟悉的计算机语言讲不选主元和列主元Gauss消去法编写成通用的子程序，然后用你所编写的程序求解下面的84阶方程组

最后，将你的计算结果与方程组的精确解进行比较，并就次谈谈你对Gauss 消去法的看法。

2．先用你所熟悉的计算机语言将平方根法和改进的平方根法编写成通用的子程序，然后用你编写的程序求解对称正定方程组Ax=b,其中 （1）b 随机的选取，系数矩阵为100阶矩阵

1011101110111011101110⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥

⎢⎥

⎢⎥⎣⎦

；

（2）系数矩阵为40阶Hilbert 矩阵，即系数矩阵A 的第i 行第j 列元素为1

1

i j a i j =+-，向量的

第i 个分量为11

1

n

i j b i j ==

+-∑。 3．用第1题的程序求解第2题的两个方程组并比较所有的计算结果，然后评论各个方法的优

劣。

三，模型建立与数学原理： （1）Gauss 消去法

Gauss 消去法即是把一个给定的矩阵分解为一个下三角L 和一个上三角U 的乘积。为此我们找到了一个最简单的下三角矩阵T k k k L I l e =-.其中1,(0,,0,,,)k k k nk l l l =+

即

1,,1111k k k

n k L l l +⎡⎤

⎢⎥⎢⎥⎢⎥=⎢

⎥-⎢⎥⎢⎥

⎢

⎥-⎢⎥⎣⎦

这种类型的下三角称作Gauss 变换。而称向量k l 为Gauss 向量。对于一个给定的向量1(,,)T n n x x x R =∈ ，我们有111,(,,,,,)T k k k k k k n k n

k L x x x x x l x x l ++=-- 只要取,1,,i ik k x

l i k n x ==+ ，便有1(,,,0,,0)T k k L x x x = 。于是将Gauss 变换作用于一个矩阵，

就将对应的列向量相应项变为0。一般的阶矩阵A,在一定条件下，我们也可以计算n-1个Gauss 变

123828384617861158611586115861158614x x x x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥

=⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦

⎣⎦