数值分析课程设计积分方法的.docx

  1. 1、本文档共18页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析课程设计积分方法的 实际应用 数值分析课程设计报告 求积公式的实际应用 学 院 数学与统计学院 专 业 信息与计算科学 学 号 姓 名 指导教师 成 绩 教师评语: 指导教师签字: 2018年1月 8日 1 绪论 数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用 计算机求解各种数学问题的数值检索方其理论与软件的实现。 随着计算机和计算方法的飞速发展,几乎所有学科都走向定量化和精确化,从而产 生了一系列计算性的学科分支,如计算物理、计算化学、计算生物学、计算地质学、计 算气象学和计算材料学等,计算数学中的数值计算方法则是解决“计算”问题的桥梁和 工具。我们知道,计算能力是计算工具和计算方法的效率的乘积,提高计算方法的效率 与提高计算机硬件的效率同样重要。科学计算已用到科学技术和社会生活的各个领域 中。 数值计算方法,是一种研究并解决数学问题的数值近似解方法, 是在计算机上使 用的解数学问题的方法,简称计算方法。在科学研究和工程技术中都要用到各种计算方 法。例如,在航天航空、地质勘探、汽车制造、桥梁设计、天气预报和汉字字样设计 中都有计算方法的踪影。计算方法既有数学类课程中理论上的抽象性和严谨性, 又有实 用性和亚丁实验性的技术特征,计算方法是一门理论性和实践性都很强的学科。在 70 年代,大多数学校仅在数学系的计算数学专业和计算机系开设计算方法这门课程。 随着 计算机技术的迅速发展和普及,现在计算方法课程几乎已成为所有理工科学生的必修课 程。计算方法的计算对象是微积分,线性代数,常微分方程中的数学问题。内容包括: 插值和拟合、数值微分和数值积分、求解线性方程组的直接法和迭代法、计算矩阵特征 值和特征向量和常微分方程数值解等问题。 2 Gauss求积公式 2.1基本原理 求积公式 b n f (x)d x Akf(xQ (2.1) a k 0 含有2n+2个待定参数, Xk, Ak(k 0,1,L , n).当Xk为等距节点时得到的插值求积公其 代数精度至少为n次,如果适当选取xk(k 0,1丄,n),有可能使求积公式(2.1)具有2n+1 次代数精度,这类求积公式称为高斯求积公式为具有一般性,研究带权积分Ib f(x) (x)d x,这里(X) 次代数精度,这类求积公式称为高斯求积公式 为具有一般性,研究带权积分I b f(x) (x)d x,这里(X)为权函数,求积公式为 a b f (x) (x)d x a n Ak f(Xk) k 0 (2.2) Ak(k 0,1,L , n)为不依赖于f(x)的求积系数, Xk(k 0,1,L , n)为求积节点,可适当选取 Xk及Ak,使(2.2)具有2n+1次代数精度。 如果求积公式(2.2) 如果求积公式(2.2)具有2n+1次代数精度, 则称其节点Xk(k 0,1,L , n)为高斯点,相 应求积公式(2.2)称为高斯求积公式。 应求积公式(2.2)称为高斯求积公式。 根据定义要使(2.2)式具有2n+1次代数精度,只要对f(x) xm, (m 0,1,L 2n 1 ), 令(2.2)式精确成立,即 nmA n m AkXk k 0 b xm (x)dx m 0,1,L 2n a (2.3) 当给定权函数(x),求出右端积分,则可由(2.3)式解得Xk(k 0,1L , n)及 Ak(k 0,1,L , n) 2.2程序实现 建立gaussl.m文件,写入如下内容: function s=gaussl(a,b ,n) h=(b-a)/n; s=0.0; for m=0:(1* n/2-1) s=s+h*(gaussf(a+h*((1-1/sqrt(3))+2*m))+gaussf(a+h*((1+1/sqrt(3))+2*m))); end 2.3实例分析 例 计算积分o --xlogxdx 解 建立gaussf.m文件以调用gaussl.m文件中的函数,再写入如下内容: function y=gaussf(x) y=sqrt(x)*log(x); 再在命令行中输入: >> s=gaussl(0,1,20) 得出如下结果: -0.4456 3 高斯 - 勒让德求积公式 3.1 基本原理 在高斯求积公式 (2.1)中,若取权函数 (x) 1 ,区间为 -1,1 ,则得公式 1n -1 f(x)dx Ak f(xk). (3.1) -1 k0 由于勒让德多项式是区间 -1,1 上的正交多项式,因此,勒让德多项式 Pn 1(x) 的零点就 是求积公式 (3.1)的高斯点。形如 (3.1)式的高斯公式称为高斯 -勒让德求积公式。 3.2 程序实现 建立 guasslegendre.m 文件,写入如下内容: function [ql,Ak,xk]=gu

文档评论(0)

xiaozu + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档