数值分析课程设计

1、课 程 设 计 报 告课程名称 课题名称 专 业 班 级 学 号 姓 名 指导教师 年 月 日湖 南 工 程 学 院课 程 设 计 任 务 书课程名称 数值分析 课 题 专业班级 学生姓名 学 号 指导老师 审 批 任务书下达日期 2009 年 5 月 4 日任务完成日期 2009 年 5 月 18日一、设计内容与设计要求1设计内容：对课程计算方法中的常见算法进行综合设计或应用（具体课题题目见后面的供选题目）。2设计要求： 课程设计报告正文内容a. 问题的描述及算法设计；b. 算法的流程图（要求画出模块图）；c. 算法的理论依据及其推导；d. 相关的数值结果（通过程序调试，；e. 数值计算结果

2、的分析；f. 附件（所有程序的原代码，要求对程序写出必要的注释）。 书写格式a要求用A4纸打印成册b正文格式：一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。c正文的内容:正文总字数要求在3000字左右（不含程序原代码）。d封面格式如下页。 考核方式指导老师负责验收程序的运行结果，并结合学生的工作态度、实际动手能力、创新精神和设计报告等进行综合考评，并按优秀、良好、中等、及格和不及格五个等级给出每位同学的课程设计成绩。具体考核标准包含以下几个部分：a平时出勤 （占10%）b系统需求分析、功能设计、数据结构设计及程序总体结构合理与否（占10%）c程序能否完整、准确地运行

3、，个人能否独立、熟练地调试程序（占40%）d设计报告（占30%）注意：不得抄袭他人的报告（或给他人抄袭），一旦发现，成绩为零分。e独立完成情况（占10%）。 课程验收要求a判定算法设计的合理性，运行相关程序，获得正确的数值结果。b回答有关问题。c提交课程设计报告。d提交软盘（源程序、设计报告文档）。e依内容的创新程度，完善程序情况及对程序讲解情况打分。三、进度安排1、 班级： 信息与计算科学：0601、0602、06032、 主讲教师：聂 存 云3、 辅导教师：聂 存 云上机时间安排：第 12 周 星期一 8时：30分11时：30分星期三 8时：30分11时：30分星期五 8时：30分11时：

4、30分第 13 周 星期三 8时：30分11时：30分星期五 8时：30分11时：30分数理系课程设计评分表课程名称： 项 目评 价设计方案的合理性与创造性设计与调试结果设计说明书的质量答辩陈述与回答问题情况课程设计周表现情况综合成绩教师签名： 日 期： 计算方法课程设计供选课题1. 线性代数系统的求解设计（供5人选：学号1-5）一、设计问题：数值求解下面的微分方程。-u=f(x u(0=u(1=0 （1）二、设计内容与方案（1） 对上述问题进行离散，得到相应的线性系统（2） 采用高斯消去（追赶法）和一种迭代方法（Jacobi或Gausiseidd）求解线性方程组（3） 利用Matlab描绘u

5、(x的图像（对数值解结果）（4） 针对（1）并进行一定理论分析问题 (1 对应的每一个精确解，可供一人进行课程设计。1人 u=x2(x-12 1人 u=sin(x 1人 u=xsin(x 1人 u=x3(x-12 1人 u= x2sin(x 2. 二维椭圆问题的离散求解（供4人选：学号6-9）一、设计问题： y数值求解下面的微分方程。 -uxx-uyy=f(x,y1u|=0 0 1 x设精确解为u=sin(x*sin(y (2人u=x2(x-12 （2人）二、设计内容与方案区域划分节点排序（自然排序）得到离散的线性代数方程组采用Gauss-Seidel迭代法求解3. 非线性问题的求解设计（供4

6、人选：学号10-13）一、设计问题 x2-3x+2-ex=0 3x12-x22=0 x2+xy+y=3 （3） 3x1x22-x13-1=0 x2-4x+3=0x(0=(1,1T二、 设计内容与方案设计各种线性收敛的迭代方法求解，然后采用Stiffensen加速方法计算用牛顿迭代方法求解（），并与1的结果进行比较。采用牛顿法迭代求解（）采用f(xk+f(xk(x-xk+(f“（xk）/2!*(x-xk2=0方法求解三、设计要求1、 每一种方法均必须输出设计的结果2、 给出算法流程3、 算法描述注：设计题目为（1）（2）为一组（2人）；（1）（3）为一组（2人）每组中的内容(1(2(3 和(1(

7、2(4各供1人4.初值问题的数值求解设计（供4人选：学号14-17）一、设计问题 (1 (2 y= -50y+50x2+2x y=-1/x2-y/x,0x1 1x2 y(0=1/3 y(1=1准确解为 y(x=（1/3）*e-50x+x2二、设计内容与方案1 采用Euler方法求解2 采用改进的Euler方法求解3 采用梯形公式法求解4 四阶Runge-Kutta方法求解三、设计要求1、 得到各剖分节点处的精确解、近似解、误差2、 进行一定的理论分析每一问题可供2人，分别为 （1）（2）(3 供1人 ；（1）（2）(4 供1人5. Runge现象及其修正方法设计（供2人选：学号18-19）一、

8、内容：插值被插值函数 (供不同同学选择f(x=1/(1+x2 x-5,5f(x=1/(1+5x2 x-1,1二、要求采用Lagrange插值：将-5,5分成10等分 x=Xi=-5+(i-1h, i=1、2.、n得到L10(x,验证Range现象，利用Matlab绘出插值函数L10(xd的图像采用分段线性插值：将-5,5分成10等分，x=(5-(-5/10=1在第i个子区间上有：xi-1,xiLi(x= - 计算函数值L1i(0.5;利用Matlab绘图分段线性插值函数的图像 从理论上给出（2）的插值误差估计结果6. 对称矩阵的条件数的求解设计（供4人选：学号20-23）一、 求矩阵A的二条件

9、数问题1 A= 问题2 A=二、设计内容：1 采用幂法求出A的 .2 采用反幂法求出A的.3计算A的条件数 A2* A-12=cond2（A）=/.（精度要求为10-6）三、设计要求1、 求出A2。2、 并进行一定的理论分析。每个问题可供2人设计内容1，3 供1人； 设计内容2 供1人7. 数值积分方法的加速和自适应算法设计（供10人选：学号24-33）一、设计内容1、 数值积分加速收敛方法。2、 自适应选取求积步长二、设计问题x21、 I=xdx2、 I=3、 -1dx4、 I= 其中令 =15、 I= 三、设计方案1、 简单梯形/Simpson公式求积2、 复化梯形/Simpson公式求积

10、3、 对复化梯形公式的结果，进行Romberg算法得到数值积分的加速4、对复化梯形公式（3点）根据给定的精度=10-6，设计自适应选取积分步长5、采用高斯公式求积6、从理论上给出2、3中的误差估计，并列出相应的精确值、近似值、误差。设计方案中的梯形公式, Simpson公式各供1人，每个设计题1人1 实验中，考虑如下情形 (i a=1 和 b=0, (ii a=1 和 b=1 ,(iii a=1 和 b=-1 ,(iv a=1 和 b=0 .取f=0.2 实验中可取a=0,b=1, 分别利用如下方法，验算相应的精度：(i 复化梯形(ii 复化辛普生公式(iii 对复化梯形构造自适应算法(iv

11、对复化辛普生公式构造自适应算法(v 对于二重、三重积分情形考虑复化中矩形公式，。3（i） 其中 a(x=1.0（ii） 其中（iii） 其中.（iv） 其中4利用(i 向前Euler方法 求解 (ii 向后Euler方法 求解 (iii 改进的Euler方法 求解 (iv 考虑如下迭代方法（二阶R-K方法）求解5求解的根（1） 简单迭代方法 ，加速技术1： ，加速技术2：Aitken 加速（2） 单点和双点弦截法，尝试加速因子，（3） 二分法 在1.5,2的根。（4） 牛顿法 修正的牛顿法 m为重根数。6、情形（i）情形（ii）采用五点差分格式（两种）含网格剖分（两种以上）、形成离散格式、数值求解、数值结果分析。7考虑区域的网格剖分技术与局部加密技术。（给出矩形区域、圆形区域、正规和扭曲、局部加密）8对问题若（i）u为调和函数，验算如下格式为六阶格式。（ii）u为一般函数，验算如下格式为二阶格式。（见专业英语书P50）(iii u为一般函数，验算如下格式为二阶格式。（见专业英语书P50）(iv u为一般函数，验算如下格式为二阶格式。9对AU=f 构造迭代格式（Gauss消去法、Gauss-Seidel迭代、Jacobi迭代、追赶法）其中 10. （i）利用SOR迭代方法求解