《数值分析课程设计83613》由会员分享,可在线阅读,更多相关《数值分析课程设计83613(17页珍藏版)》请在金锄头文库上搜索。
1、课 程 设 计 报 告课程名称 课题名称 专 业 班 级 学 号 姓 名 指导教师 年 月 日湖 南 工 程 学 院课 程 设 计 任 务 书课程名称 数值分析 课 题 专业班级 学生姓名 学 号 指导老师 审 批 任务书下达日期 2009 年 5 月 4 日任务完成日期 2009 年 5 月 18日一、设计内容与设计要求1设计内容: 对课程计算方法中的常见算法进行综合设计或应用(具体课题题目见后面的供选题目)。2设计要求:l 课程设计报告正文内容a. 问题的描述及算法设计;b. 算法的流程图(要求画出模块图);c. 算法的理论依据及其推导;d. 相关的数值结果(通过程序调试),;e. 数值计
2、算结果的分析;f. 附件(所有程序的原代码,要求对程序写出必要的注释)。l 书写格式a要求用A4纸打印成册b正文格式:一级标题用3号黑体,二级标题用四号宋体加粗,正文用小四号宋体;行距为22。c正文的内容:正文总字数要求在3000字左右(不含程序原代码)。d封面格式如下页。l 考核方式指导老师负责验收程序的运行结果,并结合学生的工作态度、实际动手能力、创新精神和设计报告等进行综合考评,并按优秀、良好、中等、及格和不及格五个等级给出每位同学的课程设计成绩。具体考核标准包含以下几个部分:a平时出勤 (占10%)b系统需求分析、功能设计、数据结构设计及程序总体结构合理与否(占10%)c程序能否完整、
3、准确地运行,个人能否独立、熟练地调试程序(占40%)d设计报告(占30%)注意:不得抄袭他人的报告(或给他人抄袭),一旦发现,成绩为零分。e独立完成情况(占10%)。l 课程验收要求a判定算法设计的合理性,运行相关程序,获得正确的数值结果。b回答有关问题。c提交课程设计报告。d提交软盘(源程序、设计报告文档)。e依内容的创新程度,完善程序情况及对程序讲解情况打分。三、进度安排1、 班级: 信息与计算科学:0601、0602、06032、 主讲教师:聂 存 云3、 辅导教师:聂 存 云上机时间安排:第 12 周 星期一 8时:30分11时:30分 星期三 8时:30分11时:30分 星期五 8时
4、:30分11时: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
5、) 利用Matlab描绘u(x)的图像(对数值解结果)(4) 针对(1)并进行一定理论分析 问题 (1) 对应的每一个精确解,可供一人进行课程设计。 1人 u=x2(x-1)2 1人 u=sin(x) 1人 u=xsin(x) 1人 u=x3(x-1)2 1人 u= x2sin(x) 2. 二维椭圆问题的离散求解(供4人选:学号6-9)一、设计问题: y数值求解下面的微分方程。 -uxx-uyy=f(x,y) 1u|=00 1 x设精确解为u=sin(x)*sin(y) (2人)u=x2(x-1)2 (2人)二、设计内容与方案区域划分节点排序(自然排序)得到离散的线性代数方程组采用Gauss-
6、Seidel迭代法求解3. 非线性问题的求解设计(供4人选:学号10-13)一、设计问题 x2-3x+2-ex=0 3x12-x22=0x2+xy+y=3 (3) 3x1x22-x13-1=0x2-4x+3=0 x(0)=(1,1)T二、 设计内容与方案设计各种线性收敛的迭代方法求解,然后采用Stiffensen加速方法计算用牛顿迭代方法求解(),并与1的结果进行比较。采用牛顿法迭代求解()采用f(xk)+f(xk)(x-xk)+(f“(xk)/2!)*(x-xk)2=0方法求解三、设计要求1、 每一种方法均必须输出设计的结果2、 给出算法流程3、 算法描述注:设计题目为(1)(2)为一组(2
7、人);(1)(3)为一组(2人) 每组中的内容(1)(2)(3) 和(1)(2)(4)各供1人4.初值问题的数值求解设计(供4人选:学号14-17)一、设计问题 (1) (2 )y= -50y+50x2+2xy=-1/x2-y/x,0x11x2y(0)=1/3y(1)=1准确解为 y(x)=(1/3)*e-50x+x2二、设计内容与方案 采用Euler方法求解 采用改进的Euler方法求解 采用梯形公式法求解 四阶Runge-Kutta方法求解三、设计要求1、 得到各剖分节点处的精确解、近似解、误差2、 进行一定的理论分析每一问题可供2人,分别为 (1)(2)(3) 供1人 ;(1)(2)(4
8、) 供1人5. Runge现象及其修正方法设计(供2人选:学号18-19)一、内容:插值被插值函数 (供不同同学选择)f(x)=1/(1+x2) x-5,5f(x)=1/(1+5x2) x-1,1二、要求采用Lagrange插值:将-5,5分成10等分 x=5-(-5)10=1Xi=-5+(i-1)h, i=1、2.、n得到L10(x),验证Range现象,利用Matlab绘出插值函数L10(x)d的图像采用分段线性插值:将-5,5分成10等分,x=(5-(-5)/10=1在第i个子区间上有:xi-1,xiLi(x)=x-xixi-1-xi - x-xi-1xi-xi-1yi计算函数值L1i(
9、0.5);利用Matlab绘图分段线性插值函数的图像从理论上给出(2)的插值误差估计结果6. 对称矩阵的条件数的求解设计(供4人选:学号20-23)一、 求矩阵A的二条件数问题1 A=110.5110.250.50.252问题2 A=2-10-12-1-12-10-12二、设计内容: 1 采用幂法求出A的 max. 2 采用反幂法求出A的min. 3计算A的条件数 A2* A-12=cond2(A)=max/min.(精度要求为10-6)三、设计要求1、 求出A2。2、 并进行一定的理论分析。每个问题可供2人设计内容1,3 供1人; 设计内容2 供1人7. 数值积分方法的加速和自适应算法设计(
10、供10人选:学号24-33)一、设计内容1、 数值积分加速收敛方法。2、 自适应选取求积步长二、设计问题x21、 I=02x2e-x2dxxdx2、 I=234cotxdx3、 23x12-1dx4、 I= 01sinxxdx其中令 sinxx=1=15、 I= 01ln(1+X)dx三、设计方案1、 简单梯形/Simpson公式求积2、 复化梯形/Simpson公式求积3、 对复化梯形公式的结果,进行Romberg算法得到数值积分的加速4、对复化梯形公式(3点)根据给定的精度=10-6,设计自适应选取积分步长5、采用高斯公式求积6、从理论上给出2、3中的误差估计,并列出相应的精确值、近似值、误差。 设计方案中的梯形公式, Simpson公式各供1人,每个设计题1人1 实验中,考虑如下情形 (i) a=1 和 b=0,