1009010122付琳数值研究课设.docVIP

成 绩 评 定 表 学生姓名 付琳 班级学号 1009010122 专 业 信息与计算科学 课程设计题目 数值分析算法案例 评 语 组长签字： 成绩 日期 20 年 月 日 课程设计任务书 学 院 理学院 专 业 信息与计算科学 学生姓名 付琳 班级学号 1009010122 课程设计题目 数值分析算法案例 实践教学要求与任务: 要求：格式以学校毕业论文格式要求为准，不准粘贴图片，尤其公式。 对每个试验，要求有：实验基本原理，实验目的，实验内容及数据来源和实验结论。 以班级为单位统一装订封皮。 6月25日，十八周的周二交论文 每人至少四个实验，最少15页 任务（实验项目）： 线性方程组数值解法 参考题目：（1) 列主元Gauss消去法（2）LU分解法 插值法和数据拟合 参考题目：（1）Lagrange插值（2）Newton插值（3）最小二乘拟合 数值积分 参考题目：（1）复化Simpon积分（2）变步长的梯形积分公式（3）龙贝格求积公式 常微分方程数值解 Runge-Kutta方法 数值方法实际应用 用数值方法解决实际问题（自选） 工作计划与进度安排: 线性方程组数值解法 （4学时） 插值法和数据拟合 （4学时） 数值积分 （4学时） 常微分方程数值解 （4学时） 数值方法实际应用 （4学时） 答辩 （4学时） 指导教师： 201 年 月 日 专业负责人： 201 年 月 日 学院教学副院长： 201 年 月 日 摘 要 计算机与计算数学的发展以及他们在工程及科学技术问题中的广泛应用，使得数值分析（计算方法）课程对高等院校理工科学生来说越来越重要。在学习的过程中我们需要掌握数值分析的理论基础、相关理论知识、MATLAB操作，至于编程相关只需要了解与应用，编程实现了线性方程组直接三角形求解（LU分解），利用Newton插值法实现对数进行拟合并得到相应的差商表（Newton插值法）以及公式，并且使用复化Simpson以及复化梯形公式编程实现数值积分问题的解决过程，同时采用Runge-Kutta法解决了常微分方程的数值解问题。 此外，数值分析是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科，是数学的一个分支，它以数字计算机求解数学问题的理论和方法为研究对象为计算数学的主体部分。主要内容插值法，函数逼近，曲线拟和，数值积分，数值微分，解线性方程组的直接方法，解线性方程组的迭代法，非线性方程求根，常微分方程的数值解法。 1、1 1 1、2实验基本原理 1 1、3 Matlab程序 2 1、4 数据来源 3 1、5 实验结论 5 实验二、插值法与数值拟合 6 2、1 实验目的 6 2、2 实验基本原理 6 2、3 Matlab程序 7 2、4 数据来源 7 2、5 实验结论 8 实验三、数值积分 9 3、1 实验目的 9 3、2 实验基本原理 9 3、3 Matlab程序 9 3、4 数据来源 10 3、5 实验结论 11 实验四、常微分方程数值解 12 4、1 实验目的 12 4、2 实验基本原理 12 4、3 Matlab程序 12 4、4 数据来源 13 4、5 实验结论 14 实验五、天气预报问题 15 5、1 实验目的 15 5、2 实验原理 15 5、3 Matlaba程序 15 5、3 数据来源 16 5、5 实验结论 19 总 结 20 参考文献 21 实验一、线性方程组数值解法 1、1for k=1:n-1; for i=k:n; if(abs(A(i,k))==max(abs(A(k:n,k)))) P(k)=i; temp=A(k,:); A(k,:)=A(i,:); A(i,:)=temp; end end for j=k+1:n; A(j,k)=A(j,k)/A(k,k); A(j,k+1:n+1)=A(j,k+1:n+1)-A(j,k)*A(k,k+1:n+1); end end P(n)=n; L(1,1)=1; L(2:n,1)=A(2:n,1); L(1,2:n)=0; U(1,1)=A(1,1); U(2:n,1)=0; U(1,2:n