求解大型稀疏线性方程组选主元高斯-约当法.docVIP

沈阳航空航天大学 数值分析课程设计 （设计程序） 题目 求解大型稀疏线性方程组的选主元高斯－约当法 班级 / 学号 学 生 姓 名 指 导 教 师 刘勇进 沈阳航空航天大学 课 程 设 计 任 务 书 课 程 名 称 数值分析课程设计 院（系） 理学院 专业 信息与计算科学 班级 学号 2012041401001 姓名 高莹 课程设计题目 求解大型稀疏线性方程组的选主元高斯－约当法 课程设计时间: 2013 年 12月 16日至 2013 年 12月 30日 课程设计的要求及内容： [要求] 学习态度要认真，要积极参与课程设计，锻炼独立思考能力； 严格遵守上机时间安排； 按照MATLAB编程训练的任务要求来编写程序； 根据任务书来完成课程设计论文； 报告书写格式要求按照沈阳航空航天大学“课程设计撰写规范. 要求：a/ 在同一图形窗口画出图形； b/ 将图形窗口分成个小图形区域。 . 及 要求：a/ 分别考虑； b/ . 用数值方法计算 ， 其中n分别取100, 200, 1000。 求由方程所确定的隐函数的导数。 求由参数方程确定的函数的二阶导数。 作出函数的曲线在区间上的两条水平切线。 求解大型稀疏线性方程组的选主元高斯－约当法 要求：1）掌握求解大型稀疏线性方程组的选主元高斯－约当法的基本原理； 2）编写求解大型稀疏线性方程组的选主元高斯－约当法的通用程序； 3）利用所编写的程序求解方程组，其中 ， 指导教师 年 月 日 负责教师 年 月 日 学生签字 年 月 日 沈阳航空航天大学 课 程 设 计 成 绩 评 定 单 课 程 名 称 数值分析课程设计 院（系） 理学院 专业 信息与计算科学 课程设计题目 求解大型稀疏线性方程组的选主元高斯－约当法 学号 2012041401001 姓名 高莹 指导教师评语： 课程设计成绩 指导教师签字 年 月 日  目 录 一、正文 1 二、总结 13 参考文献 14 附 录 15 正文 一·通识题 1、作出下列函数的图形： （1）. 要求：a/ 在同一图形窗口画出图形； b/ 将图形窗口分成个小图形区域。 题目分析： 这道题是利用matlab绘图画出的图像，其中n的取值为2,3,4,5，采用plot函数绘制，x取值范围从0到1，所以选取适当节点，这样可以尽量减少误差。其次，指数函数要用点乘运算，即 y=x.^n。 程序运行结果如图1-1： 图1-1 源程序： x = 0:0.0001:1; n=2;i=1; for n=2:1:5 y = x.^n; subplot(4,1,i); i=i+1; plot(x,y); title(y=x^n); xlabel(x);ylabel(y); end (2) . 题目分析： 利用matlab绘图画出的图像，采用plot函数绘制，x取值为-1到1，所以选取适当节点，可以使画出来的图像更接近真实图像，同样，数与x之间需要用点乘。 程序运行结果如图1-2： 图1-2 源程序： x = -1:0.000001:1; y = 1./(1+25.*x.^2); plot(x,y); title(y=1/(1+25*x^2),x from -1 to 1); xlabel(x); ylabel(y); （3） 题目分析： 利用matlab绘图画出y=arctan(x)的图像，采用plot函数绘制，x的范围为-50到50，选取适当结点，可以使图像更加精准。根据题目要求，编写出程序。 程序运行结果如图1-3-1： 图1-3.1 源程序： x = -50:0.0001:50; y = atan(x); plot(x,y); title(y=atan(x)); text(1,pi/4,y=atan(x)); xlabel(x); y