数值分析法求正弦余弦积分函数.pdf

天津职业技术师范大学 课 程 设 计 任 务 书 理 学院 数学 1403 班 学生 张群 课程设计课题： 用数值积分法计算正弦积分函数和余弦积分函数 一、课程设计工作日自 2016 年 7 月 4 日至 2016 年 7 月 5 日 二、同组学生： 无 三、课程设计任务要求（包括课题来源、类型、目的和意义、基本要求、完成时 间、主要参考资料等）： 课题来源： 教师自拟 类型： 理论研究 目的和意义：培养学生对数值分析中主要算法的应用能力， 探索相关算法之 间的内在联系。 基本要求： 根据数值分析课程所学的知识，应用 Matlab 软件编写程序，完 成对算法及其内在原理的实验研究。 完成时间： 参考资料： 《数值分析》 李庆扬等 清华大学出版社 指导教师签字： 教研室主任签字： 一、问题叙述 用数值积分法计算正弦积分函数和余弦积分函数 提示： 正弦积分，余弦 si(x) x sint dt 函数 ci (x) x cost dt 0 t t 要求：（ 1）编写函数，对任意给定的 x，可求出值。 （2）使用尽可能少的正余弦函数的调用， 计算更精确的值。 （用 多种方法，创新方法） 二、问题分析 1 、数值积分基本原理：用数值分析求解积分的数值方法有很多， 如简单的梯形法、矩形法、辛普森（ Simpson）法、牛顿 - 科斯特 （Newton-Cotes ）法等都是常用的方法。 它们的基本思想都是将整个 积分区间 [a ，b] 分成 n 个子区间 [x i ,x i+1 ] ，i=1 ，2，, ， n，其中 x 1=a， xn+1=b。这样求定积分问题就分解为求和问题。 2、本题要求用数值积分法计算正弦积分函数和余弦函数积分，那么 应该从编写函数的入手，建立 function 的 m文件，通过对函数的调 用，对任意跟定的 x 值，求出积分函数的值。 三、数值积分法求解问题 1、 梯形公式、矩形公式 首先，积分中值定理告诉我们，在积分区间 [a，b] 内存在一点 ，成 b 立 f （x ） dx= （b-a）f （ ），就是说，底为 b-a 而高为 f （ ）的矩 a 形面积恰等于所求区边梯形的面积。 如果我们用两端点 “高度”f （a） 与 f （b）的算术平均值作为平均高度 f （ ）的近似值，这样导出的 b b - a 求积公式 f （x ） dx ≈ [f （a）+f （b）]便是我们熟悉的梯形公式。 a 2 （b - a） 将积分区间 [a,b]n 等分，每个小区间宽度均为