本科优秀数学本科毕业论文

合集下载

数学本科毕业论文傅里叶级数与傅里叶变换的关系与应用

数学本科毕业论文傅里叶级数与傅里叶变换的关系与应用

本科生毕业论文(申请学士学位)论文题目傅里叶级数与傅里叶变换的关系与应用作者姓名刘军专业名称数学与应用数学指导教师许志才/ 张玲2014年6月学生:(签字)学号:2012220146论文答辩日期:2014年x月xx日指导教师:(签字)目录摘要: 0关键词 0Abstract 01绪论 (1)2傅里叶级数的概念 (1)2.1周期函数 (2)2.2傅里叶级数的定义 (2)3 傅里叶变换的概念及性质 (10)3.1傅里叶变换的概念 (10)3.2傅立叶变换的性质 (11)4傅里叶变换与傅里叶级数之间的区别与联系 (12)5傅里叶级数和傅里叶变换的应用 (12)5.1傅里叶级数的应用 (12)5.2傅里叶变换的应用 (13)参考文献 (14)傅里叶级数与傅里叶变换的关系与应用摘要:傅里叶级数是对周期性现象做数学上的分析,而傅里叶变换则可以看作傅里叶级数的极限形式,它也可以看作是对周期现象进行数学上的分析。

除此之外,傅里叶变换还是处理信号领域的一种很重要的算法。

傅里叶变换是一种分析信号的方法,它可分析信号的成分,也可用这些成分合成信号。

很多波形可以作为信号的成分,例如余弦波,方波,锯齿波等等,傅里叶变换作为信号的成分。

在电子类学科,物理学科,信号处理学科等众多领域都有着广泛的应用。

傅里叶级数针对的是周期性函数,傅里叶变换针对的是非周期性函数,它们在本质上都是一种把信号表示成复正选信号的叠加,存在相似的特性。

关键词:傅里叶级数;傅里叶变换;周期性Fourier series And Fourier TransformsAbstract: Fourier series is made mathematical analysis to cyclical phenomenon, and Fourier transform can be seen as the limit form of Fourier series, it also can be regarded as a mathematical analysis of cycle phenomenon. In addition, the Fourier transform is a kind of very important in the field of signal processing algorithms.Fourier transform is a method of signal analysis, it can analyze signal component, also can use these ingredients synthetic signal. Many waveform can be used as a signal of ingredients, such as cosine wave, square wave, sawtooth wave, etc., the Fourier transformas a signal of composition. In electronics disciplines, physics, signal processing disciplines etc many fields have a wide range of applications.Fourier series is for periodic function, Fourier transform for is a periodic function, they are in essence a kind of papers said the signal into a complex signal superposition, similar features.Key words: Fourier series; Fourier Transform; Periodic1绪论傅里叶级数是法国数学家J.-B.-J.傅里叶在研究偏微分方程的边值问题时提出来的,从而极大的推动了偏微分方程理论的发展,在数学物理以及工程中都具有重要的应用。

浅谈回归分析在葡萄酒等级评估的应用 数学本科毕业论文

浅谈回归分析在葡萄酒等级评估的应用 数学本科毕业论文

浅谈回归分析在葡萄酒等级评估的应用数学本科毕业论文摘要:回归分析是统计学研究中一种重要的分析工具,在葡萄酒等级评估中的应用越来越受到关注。

本文以葡萄酒等级评估为例,具体分析了回归分析在该领域的应用及其优缺点。

通过分析,本文认为,回归分析在葡萄酒等级评估中应用广泛,能够有效的评估葡萄酒的质量等级。

但是,其缺点也不能忽略,比如需要对数据进行清洗、处理等。

同时,回归分析能够结合其他的分析方法,如主成分分析、聚类分析等,进一步提高数据分析的准确性和有效性。

关键词:回归分析;葡萄酒;等级评估;数据分析1. 引言随着葡萄酒消费市场的不断发展,葡萄酒质量等级的评估变得越来越重要。

葡萄酒质量等级评估最初依靠专业人员的品尝鉴定,但由于人口味的差异和评价标准不统一等原因,直接品尝鉴定的评估方法逐渐受到质疑。

为了解决这一问题,研究人员开始采用一些量化的方法进行评估,其中回归分析是最常用的方法之一。

本文将以葡萄酒等级评估为例,详细探讨回归分析在该领域的应用及其优缺点。

2. 回归分析在葡萄酒等级评估中的应用回归分析是一种通过找到一条直线或曲线来描述两个或多个变量之间关系的统计学方法。

在葡萄酒等级评估中,回归分析可以被用来寻找葡萄酒品质和其他因素之间的关系,比如葡萄品种、产地、酿造方法、陈年时间等。

通常情况下,需要收集大量的葡萄酒品质和其他变量的数据,然后使用统计软件进行数据分析。

回归分析的结果可以帮助酿酒师和酿酒商预测葡萄酒的质量等级,从而更好地满足市场需求。

回归分析也可以用来比较不同酿造方法、葡萄品种和陈年时间对葡萄酒品质的影响,从而选择最佳的酿造方式或陈年时间。

3. 回归分析在葡萄酒等级评估中的优缺点3.1 优点(1)快速准确:回归分析是一种非常快速和准确的方法,可以快速处理大量的数据并生成可靠的结果。

(2)可解释性强:回归分析可以清晰地显示不同变量之间的关系,帮助分析人员更好地理解数据。

(3)可预测性强:回归分析可以帮助预测未来的数据趋势,从而更好地满足市场需求。

本科毕业论文_多项式方程的判别式与求根公式

本科毕业论文_多项式方程的判别式与求根公式

本科毕业论⽂_多项式⽅程的判别式与求根公式东莞理⼯学院本科毕业论⽂(2015届)题⽬: 多项式⽅程的判别式与求根公式学⽣姓名: 姚培基学号: 201141410230院(系):计算机学院专业班级: 信息与计算科学(2)班指导教师:起⽌时间: 2015年1⽉—2015年5⽉多项式⽅程的判别式与求根公式摘要: 近代数学史甚⾄能说是⼀部求解多项式⽅程的历史。

对于⾼次⽅程的数值根求解法,⼈们从很早就开始并⼀直探求这样的问题。

⽽且在古代,很多⼈都想出了⼀个办法来解决各种各样的多项式⽅程。

如卡尔⽶诺的《⼤术》,贾宪的《黄帝九章算法细草》,秦九韶的《数书九章》等等。

在⽬前,有关问题求解多项式⽅程根的在⼯程实践中占有举⾜轻重的地位。

如在⼈类的⽣活过程中,经济建设和科学技术的发展过程中,计算⼀直起着⾮常重要的作⽤。

当⼈们在进⾏科学或者⼯程计算时,求解多项式⽅程组更是⾮常容易遇到的问题之⼀。

许多领域如⾃然⽣活和⼯程科学最终都可以归结为求解多项式⽅程组的问题。

这个时候⼈们就通常需要处理求解代数⽅程组的问题,如果当项较简单或变元较少时,计算过程就好相对来说简单⼀些;但是当项⾮常复杂或变元⾮常多的时候,那么其求解的过程中往往会遇到⽐较多的困难。

对多项式⽅程的判别式和求根公式的研究,在理论研究和实际⼯程计算中,具有⼗分重要的意义。

关键词: 多项式; 判别式; 求根公式; MATLABDiscriminant and seek the root of polynomial equationsAbstract: the modern mathematics that would become a history of polynomial equation solution. People long ago began to explore the problem of high order equation of numerical method. But in ancient times, many people have been developed to solve all kinds of method of polynomial equations. Such as "chapter nine of the yellow emperor algorithm fine grass" of jia xian, chiu-shao the number of book chapter nine, Carl mino "big operation" and so on.In nowadays, polynomial equation for the root problem has a pivotal position in the engineering practice. As in human life, economic construction and development of science and technology in the process of calculation is always plays a very important role. In science and engineering calculation, to solve the polynomial equations is one of the most common problems in the natural life and the computing problem in the field of engineering science and many other eventually all boils down to solving the polynomial equations. At this time often need to deal with algebraic equations to solve the problem, if the argument or a simpler, less calculation process is relatively simple; And when the argument is very more or when the item is very complex, itssolving process is often more difficult.The discriminant and seek the root of polynomial equations, in theoretical research and practical engineering calculation, have very important significance.Key words: polynomial; The discriminant. Root formula; MATLAB⽬录⼀、引⾔ (1)(⼀)⼀元⼆次⽅程的判别式和求根与韦达定理 (1)⼆、⼀元多次多项式 (8)(⼀)代数基本定理 (9)(⼆)域论基础 (10)(三)多项式⽅程的判别式 (11)(四)⽜顿恒等式 (12)(五)关于⼀元五次⽅程 (19)三、总结与展望 (20)参考⽂献 (23)致谢 (25)⼀、引⾔在⼈类研究数学的历史长河中,追溯到公元9世纪的波斯,数学家、天⽂学家及地理学家花拉⼦⽶作为第⼀⼈给出了⼀元⼆次⽅程的⼀般解法。

数学与应用数学专业毕业论文范文

数学与应用数学专业毕业论文范文

如何写数学与应用数学专业的论文我是一位大一的学生,导员老师为了虽然我没写过论文,但还是想提点建议,楼主不妨考虑一下。

作为大一学生,限于学识和能力,要写作的所谓“专业论文”,不会要求达到毕业论文那样高的水平,只要对所学过某一方面的知识和方法作一个较为系统的整理就可以了。

鉴于此,下面就楼主所提到的四门课程各拟一题,仅供参考: 1.数学分析:极限的求法; 2.高等代数:行列式的计算方法; 3.空间解析几何:仿射变换及其应用; 4.高等几何:高等几何在平面几何证题中的应用。

应用数学专业毕业论文先修课程:数学与应用数学专业主要课程、教育类课程等适用专业:数学与应用数学(本科、师范)一、目的培养和提高学生综合运用所学知识分析、解决问题的能力(包括数学理论研究和应用研究的能力、教学研究能力、文献检索、科技论文的写作能力)。

使学生获得科学、教学研究方法的初步训练。

培养学生的独立研究能力和重视开发学生的创新能力。

两名或两名以上学生选做同一课题论文时,各人的内容应有较大区别。

学生选定课题后,应填写《毕业论文任务书》,经指导教师同意,方可进行论文工作。

四、毕业论文成绩评定 1.学生毕业论文成绩的评定采取指导教师和毕业论文答辩小组分别单独评分,按比例综合评定,最后由毕业论文答辩委员会综合平衡审定。

2.成绩分5个等级:优秀、良好、中等、及格、不及格。

毕业生毕业论文统一格式要求一、论文用纸:B5纸打印。

二、论文标题: 1、主标题:用小二号黑体字,置于首页第一行,居中。

2、正文采用四级标题,分别以“一、(一)、1、(1)”标明。

其中一级标题用黑体字,二级标题用楷体,三、四级标题与正文字体相同。

三、论文正文: 1、字体:用四号仿宋体。

2、段落:行距为24磅。

3、页码:居中。

四、年级、专业与姓名:四号宋体,置于主标题与正文之间,居中,上下各空一行。

五、注释:如有注释,皆在正文之后注明。

数学与应用数学大学导论课论文怎么写(一)题名(Title,Topic)题名又称题目或标题。

数学应用数学本科毕业论文参考题目

数学应用数学本科毕业论文参考题目

数学应用数学本科毕业论文参考题目1、对数学教学中分层教学的体验和看法2、数学教学中培养学生创新能力浅论3、案例分析:由《立方根》的情景引入所想到的4、浅谈七年级学生数学学习习惯培养策略5、初中数学作业中出现的错误问题及策略6、课堂上如何培养初中学生解决问题的能力7、如何培养初中生的数学探究能力8、浅谈中学数学的函数学习9、现代教育技术在数学教学中的运用10、初中数学课堂情境探究式教学模式的运用探索11、数学教学中良好个性心理品质的培养12、浅谈初中数学学困生的成因及转化策略13、绝对值不等式的解法14、论初中生在数学教学中的数学体验15、主题式教学在初中数学中的应用16、试论如何提高学生代数运算能力17、在数学概念教学中实施“局部探究”的实践18、例谈中学教学的“教与学”19、运用发现法教学,培养学生创新能力20、实践自主学习,促进自主发展1、新课导入环节存在的问题及成因分析2、数学教学目标制定应考虑的几对辩证关系3、提高分层教学实效促进全体学生发展4、初中生数学问题解决观的现状及其分析5、化归思想在数学教学中的应用6、初中生数学学习方式和学习负担的调查分析7、运用数学建模思想提高中学数学教育质量8、人教版和华师版反比例函数编排的比较与探讨9、有效教学的灵魂是以生为本--切线长定理教学案例与分析10、浅谈初中生数学建模能力的培养11、如何培养农村初中学生的数学学习兴趣12、基于学生几何认知水平的教学目标设计探讨13、中考复习导学案设计的实践与思考14、数形结合话三角--三角函数在中考试题中的应用举例15、数形结合在初中数学解题中的应用16、对新课程数学教学中初三复习课的几点思考17、初中数学学案教学教师适应性调查研究18、微课程在初中数学课堂中的功能性研究19、加强初中数学思想方法教学的策略20、试分析新课改下中学数学教学的有效模式1、数学概念教学中有效提问的量化研究2、大、中学数学教学衔接问题的研究综述3、高中数学课程标准下选修课“数学史选讲”教学研究4、普通高中数学课程标准与教学大纲课程编制的对比研究5、新课标下大学概率统计教学与中学数学教学内容的衔接探讨6、让数学文化走进课堂7、高中学生数学建模能力与数学学业成绩关系的调查与分析8、高等数学与新课标下高中数学教学内容对接的研究9、高一数学教学中如何解决好初高中衔接问题10、浅析高中数学生成性课堂的构建策略11、论数学文化视角下的中学数学课堂教学12、高等数学与高中数学衔接改革的研究13、高考数学应用题的特点与启示14、浅谈高中数学导学案教学中存在的问题及对策15、数学课程发展的趋势与思考16、浅议向量在高考数学中的应用17、《函数的概念》教学设计中存在的问题及其解决--兼评网上教学设计18、实施分组分层教学,提高课堂教学效率19、培养反思思维习惯促进创新能力提高感谢您的阅读,祝您生活愉快。

非线性方程求解的不动点算法及研究数学专业设计大学本科毕业论文

非线性方程求解的不动点算法及研究数学专业设计大学本科毕业论文

长沙学院CHANGSHA UNIVERSITY 本科生毕业论文摘要非线性方程在工程实践、经济学信息安全和动力学等方面的大量实际问题中有着极为广泛的应用,而不动点迭代算法作为数学研究的一个新方向,是求解非线性方程问题的一个最基本而又重要的方法.本文主要介绍了非线性方程求解的不动点算法及其研究,首先,综述了非线性方程求解的不动点算法的研究背景、并阐述了本文的主要工作以及介绍了误差、有限差等基本知识;然后,详细介绍了不动点迭代算法的基本思想、在什么条件下方程存在不动点的收敛定理、不动点的收敛阶定理和Atiken加速公式;最后,考虑到方程可能会不满足不动点迭代收敛定理的两个条件的情况提出了反函数法、牛顿迭代法、Steffensen 迭代法和松弛法这四中处理方法.关键词:非线性方程,不动点原理,迭代法ABSTRACTA large number of practical problems of nonlinear equations in engineering practice,economics of information security and other the dynamics has a very wide range of applications.As a new direction in the study of mathematics,fixed point iterative algorithm is a basic and important methods to solving nonlinear equations problem.This paper describes the solving nonlinear equations fixed point algorithm and research. First, the research background of solving nonlinear equations fixed point algorithm and the main word are introduced, the basic knowledge of errors,finite difference are introduced ; Second, the fixed point iterative basic idea, algorithm convergence and convergence rate and the aitken formula are detailed; Last, inverse function method, the newton iterative method,Steffensen iterative method and the relaxation method are proposed when the equation dose not satisfy the fixed point iteration convergence conditions.Keywords:Nonlinear Equation, Fixed Point Theorem, Iterative Method目录摘要 (I)ABSTRACT (I)第1章绪论 (1)1.1 研究背景 (1)1.2 预备知识 (2)1.2.1 误差 (2)1.2.2 有限差 (3)第2章非线性方程求解的不动点迭代算法 (5)2.1不动点迭代算法的基本思想 (6)2.2 不动点迭代算法的收敛性 (7)2.3 不动点迭代算法的收敛速度 (11)2.4 加速不动点迭代算法及其收敛性 (12)第3章非收敛不动点迭代格式的几类处理方法与比较 (14)3.1 非收敛不动点迭代格式的几类处理方法 (15)3.1.1 反函数法 (15)3.1.2 牛顿迭代法 (15)3.1.3 Steffensen迭代法 (15)3.1.4 松弛法 (16)3.2 数值实例 (17)结论 (21)参考文献 (23)附录 (24)致谢 (35)第1章绪论1.1 研究背景非线性数值解的问题是现代数学的主要研究课题之一,这不仅是由于科学技术发展的需要,而且也是由于计算技术的高速发展提供了解决这类问题的可能,利用计算机解决非线性问题时,最终总是将其化成为有限维非线性问题,或称为非线性代数问题.对于求解非线性方程,无论从理论上还是从计算机上,都比解线性问题要复杂的多,一般的非线性方程是很难求出精确解的,往往只能求出近似解、数值解,而长期以来,人们为了得到满足条件的近似值,许多计算工作者致力于研究求解非线性方程的有效方法,尤其是计算机出现后函数方程求根的数值解法得到了蓬勃发展,十七世纪,微积分出现时,Newton和Halley发明了各自的新的数学工具去解非线性方程,十八世纪,随着微积分的快速蓬勃发展,Euler和Lagrange分别找到了一个无穷级数来表示方程解,并以各自的名字来命名,十九世纪,人们开始注重问题分析的严密性,柯西建立了优级数技巧,该技巧不断的被以后的事实证明对于研究方程近似解序列的收敛性是很有成效的,在分析严密性发展的时代,Ostrowski对Newton迭代法的收敛性问题规定了一个合理的假设和一个令人满意的解法,在软件分析完善的年代,Kantorovich把Newton迭代法和Ostrowski的结果推广到Banach空间,从而使许多用硬分析去做非常棘手的有关问题被轻轻松松地推论中得到了令人满意的解决,等等,总之,这些方法不断地被后人完善,但在目前,实际问题中可能还需要求方程的负根,求非线性方程(组)的迭代法,求微分方程迭代法等等,迭代方法还需要更深入的研究,同时意味着迭代法的发展空间将会更广阔.本文将着重介绍求解非线性方程的不动点算法,其中文献[3]是由王则柯先生于1988年总结的单纯不动点算法,他简述了不动点在非线性方程数值解、微分方程初值问题、边值问题、分支问题等许多应用问题方面的十多年的发展,以及对单值连续映射的不动点或零点问题进行了讨论,在文献[4]中,许炎先生简单的阐述了国内外有关不动点理论的发展状况,并主要讨论了L-Lipschitz映射的不动点迭代逼近定理,[3][4]这两篇文献都总结出了不动点问题的研究和解决在实际问题中起到了至关重要的作用,这一系列的文献还有[5][6][7][8],而秦小龙先生在文献[9]中介绍了迭代法的发展情况以及相关定理,为本篇论文提供了大量的基础信息,王公俊先生在文献[10]中分别介绍了常用的求解非线性方程的方法以及收敛性,在文献[11]中,张卷美主要研究了一类不动点迭代法的求解,在迭代格式不满足迭代条件的情况下,运用的几种处理方法,并且用C语言编程上机进行了计算,对迭代收敛结果进行了分析和比较,为本文提供了大量的信息,另外,本文还借鉴了2本不同出版社的《数值分析》教材的大量内容.本文主要介绍了非线性方程求解的不动点算法及其应用,第一章为绪论部分,主要介绍了为什么要研究本文的一些原因、目的,以及价值,也准备了一些预备知识作为对正文的补充;第二章介绍迭代法与不动点的相关思想原理、定理以及迭代法的收敛条件,是本文的一个主要章节和工作重心,并且举出了几个实例来辅助证明了运用不动点迭代法求解非线性方程的方便以及准确性;第三章作为对第二章节的一个完善,非常具有实用性,主要讨论了非收敛不动点迭代格式的几类处理方法,并通过数值实例给予了证明.1.2 预备知识1.2.1 误差误差的来源有多个方面,主要有模型误差、观测误差、截断误差、舍入误差等.例1.1 可微函数)(x f 用泰勒(Taylor)多项式 ,!)0(...!2)0(!1)0()0()()(2n n n x n f x f x f f x p +''+'+= 近似代替,则数值方法的截断误差是 ,)!1()()()()(1)1(+++=-=n n n n x n f x p x f x R ξ ξ在0与x 之间.也就是说,截断误差就是近似值与精确值之间的误差.例1.2 用3.14159近似代替π,表示舍入误差..0000026.014159.3⋅⋅⋅=-=πR同样,可以定义舍入误差是指由于计算机字长有限在表示时产生的误差.定义1.1[1] 设x 为准确值,*x 为x 的一个近似值,称x x e -=**为近似值的绝对误差,简称误差.然而,在实际中,人们是无法准确计算出误差*e 的精确值的,一般是根据需要估计出误差的绝对值不超过某正数*ε,也就是误差绝对值的一个上界,*ε叫做近似值的误差限,它总是正数.对于一般情形,**||ε≤-x x ,即,****εε+≤≤-x x x (1.1)这个不等式有时也表示为.**ε±=x x (1.2)误差的大小有时还不能完全表示近似值的好坏,例如,有两个量110±=x ,51000±=y ,则.5,1000;1,10****====y x y x εε虽然*y ε是*x ε的5倍,但是%5.0**=y y ε比%10**=xx ε小得多,这就说明了*y 近似y 的程度比*x 近似x 的程度要好得多,因此,除了需要考虑误差的大小之外,还应该考虑准确值本身的大小.我们把近似值的误差*e 与准确值x 的比值 ,**xx x x e -= (1.3) 称为近似值*x 的相对误差,记作*r e .在实际计算中,由于真值x 总是不知道的,通常取 ,*****xx x x e e r -== (1.4) 作为*x 的相对误差,条件是***xe e r =较小,此时 ,)(1)()()()(**2*****2*******x e x e e x x e x x x x e x e x e -=-=-=- (1.5) 是*r e 的平方项级,故可忽略不计.相对误差也可正可负,它的绝对值上界叫做相对误差限,记作*r ε,即 .||***x r εε= (1.6) 根据定义,上例中 %10||**=x x ε与%5.0||**=y y ε分别为x 与y 的相对误差限,很显然*y 近似y 的程度比*x 近似x 的程度好得多.在实际运算中,为了避免误差危害,数值计算中通常不采取数值不稳定算法,在设计算法是应该尽量避免误差危害,防止有效数字损失,通常要避免两个相近数字相减和用绝对值很小的数做除数,还要注意运算次序和减少运算次数.1.2.2 有限差定义1.2[2] 分别称),()()(x f h x f x f -+=∆ (1.7) ),()()(h x f x f x f --=∇ (1.8) ⎪⎭⎫ ⎝⎛--⎪⎭⎫ ⎝⎛+=22)(h x f h x f x f δ (1.9)为函数)(x f 在点x 的一阶向前差分,一阶向后差分和一阶中心差分,或者分别简称为一阶前差,一阶后差,一阶中心差,统称为(一阶)有限差,其中)0(>h 表自变量的有限增量,称为步长,∇∆,和δ分别成为(一阶)前差算子、(一阶)后差算子和(一阶)中差算子,统称为(一阶)有限差算,仿此,可以定义高阶有限差,例如,二阶前差记作)(2x f ∆,定义为[]).()()()(2x f h x f x f x f ∆-+∆=∆∆=∆ (1.10) 于是,有).()(2)2()(2x f h x f h x f x f ++-+=∆ (1.11) n 阶前差记作)(x f n ∆,定义为[]).()()()(111x f h x f x f x f n n n n ---∆-+∆=∆∆=∆ (1.12) 同样,二阶后差)(2x f ∇和n 阶后差)(x f n ∇分别定义为[])()()()(2h x f x f x f x f -∇-∇=∇∇=∇ (1.13)和[]).()()()(111h x f x f x f x f n n n n -∇-∇=∇∇=∇--- (1.14) 二阶中心差 )(2x f δ 和n 阶中心差)(x f n δ分别定义为[],22)()(2⎪⎭⎫⎝⎛--⎪⎭⎫ ⎝⎛+==h x f h x f x f x f δδδδδ (1.15)和[].22)()(111⎪⎭⎫⎝⎛--⎪⎭⎫⎝⎛+==---h x f h x f x f x f n n n n δδδδδ (1.16)我们规定0()()f x f x ∆=, 0()()f x f x ∇=, 0()()f x f x δ=.有限差有下列一下性质:(1)常数的有限差恒为零.(2)有限差算子为线性算子,即对任意的实数α,β恒有()),()()()(x g x f x g x f ∆+∆=+∆βαβα (1.17) ()),()()()(x g x f x g x f ∇+∇=+∇βαβα (1.18)()).()()()(x g x f x g x f βδαδβαδ+=+ (1.19)(3)用函数值表示高阶有限差:()),)((1)(0h i n x f C x f ni in i n -+-=∆∑= (1.20)()),(1)(0ih x f C x f n i in i n --=∇∑= (1.21)()),)2((1)(0h i hx f C x f n i i n i n -+-=∑=δ (1.22)其中 .!)1()1(i i n n n C i n +-⋅⋅⋅-= (4)用有限差表示函数值 .)()(0∑=∆=+n i i i nx f C nh x f (1.23)第2章 非线性方程求解的不动点迭代算法2.1不动点迭代算法的基本思想首先讨论解非线性方程)(x g x = (2.1) 的问题. 方程(2.1)的解又称为函数g 的不动点. 为求g 的不动点,选取一个初始值0x ,令⋅⋅⋅==-,2,1),(1k x g x k k (2.2) 已产生序列}{k x . 这一类迭代法称为不动点迭代. )(x g 又被称为迭代函数, 很显然,若迭代序列}{k x 收敛,即有,lim p x k k =∞→ (2.3) 且)(x g 连续,则p 是g 的一个不动点.例2.1[2] 方程042)(23=-+=x x x f 在区间[]2,1中有唯一跟. 我们可以用不同的方法将它化为方程:(1);42)(231+--==x x x x g x(2);22)(212⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-==x x x g x (3);22)(2133⎪⎪⎭⎫ ⎝⎛-==x x g x (4);212)(214⎪⎭⎫ ⎝⎛+==x x g x (5),4342)(2235xx x x x x g x +-+-== 等等.取初始值5.10=x ,分别用式(2.2)的迭代格式计算,结果如下表.表2.1 例2.1迭代公式计算结果从表2.1中可以看到,选取迭代函数为)(4x g ,)(5x g ,分别12次和4次,得到方程的近似根1.130395435.若选取)(3x g 作为迭代函数,则k 为奇数时迭代子序列单调增加,k 为偶数时迭代子序列单调减小,迭代120次得到近似根1.130395436. 若选取)(1x g 作为迭代函数,则迭代序列不收敛, 若选取)(2x g 为迭代函数,则出现了负数开方,因而无法继续进行迭代.2.2 不动点迭代算法的收敛性通过例2.1,可以总结出,对于同一个非线性方程的求解问题,在转化为迭代方程时应该要使其解的迭代次数达到最小,且得到的解应该最精确. 在选择迭代函数)(x g 的基本原则是,首先必须保证不动点迭代序列⋅⋅⋅⋅⋅⋅,,,21k x x x 在)(x g 的定义中,以使迭代过程不至于中断;其次要求迭代序列}{k x 收敛且尽可能收敛得快.定理2.1[2] 假设)(x g 为定义在有限区间[]b a ,上的一个函数,它满足以下条件 (1)对任意[]b a x ,∈有[];,)(b a x g ∈ (2.4) (2))(x g 的导数)(x g '在[]b a ,上有界,且存在正数1<L 使得对一切[]b a x ,∈有 ,1|)(|<≤'L x g (2.5) 那么对于任意初始值[]b a x ,0∈由不动点迭代(2.2)产生的序列都收敛于g 在[]b a ,的唯一不动点p ,并且有误差估计式|,|1||01x x LL e kk --≤,1≥k (2.6) 其中p x e k k -=.证明 首先证明g 的不动点存在且唯一. 令 ).()(x g x x h -= (2.7)据条件(1),0)()(≤-=a g a a h .0)()(≥-=b g b b h又据条件(2),在)(x g '上存在,因此)(x g 在[]b a ,上连续,从而)(x h 在[]b a ,上也连续,因此方程0)(=x h 在[]b a ,上至少有一个跟.现假设方程0)(=x h 在[]b a ,上有两个根q p q p ≠,,,则由Lagrange 中值定理知,在p 与q 之间存在ξ使得|,||)(||))((||)()(|||q p g q p g q g p g q p -'=-'=-=-ξξ 再由(2.5).|||||||)(|q p q p L q p g -<-≤-'ξ这就得到矛盾式:.||||q p q p -<- 因此q p =,即0)(=x h 在[]b a ,中的根是唯一的.其次证明由不动点迭代格式(2.2)产生的序列}{k x 是收敛于p 的.根据定理条件(1)[]b a x k ,∈,⋅⋅⋅=,2,1,0k ,因此不动点迭代过程不会中断.由(2.5)式有).()(1p g x g p x k k -=-- (2.8) 应用Lagrange 中值定理,并根据(2.5)式有|||||)(||)()(|||111p x L p x g p g x g p x k k k k -≤-'=-=----ξ.||0p x L k -≤⋅⋅⋅≤ (2.9) 因为10<<L ,所以,0||lim ||lim 0=-≤-∞→∞→p x L p x k k k k即.lim p x k k =∞→ (2.10)最后,推导估计式(2.6).应用收敛性的证明过程,有|||)()(|||111-++-+++++-≤-=-j k j k j k j k j k j k x x L x g x g x x|,|01x x L j k -≤⋅⋅⋅≤+ (2.11) 于是()∑∑-=+++-=++++-≤-=-11101||m j j k j k m j j k j k k m k x x x xx x||1)1(||010110x x LL L x x Lm k m j jk ---=-≤∑-=+.||101x x LL k--≤ (2.12)在上式中令∞→m ,得.1||||01x x L L p x e kk k --≤-= (2.13) (2.6)式得证.例2.2[2] 讨论例2.1中不动点迭代⋅⋅⋅=⎪⎪⎭⎫ ⎝⎛-==--,2,1,22)(213113k xx g x k k k (2.14) 的收敛性. 为使解的近似值的误差不超过810-,试确定迭代次数.解 迭代法(2.14)的迭代函数为.22)(2133⎪⎪⎭⎫⎝⎛-=x x g )(3x g 的定义域为]4,(3-∞.取初始值5.10=x ,由不动点迭代(2.21)得559016994.01=x ,因此取区间[][]5.1,5.0,=b a .由于,02243)(22133<⎪⎪⎭⎫ ⎝⎛--='-x xx g [],5.1,5.0∈x 因此)(3x g 在[]5.1,5.0上单调减小. 而[],559.0)5.1()(min 335.1,5.0≈=∈g x g x[],399.1)5.0()(max 335.1,5.0≈=∈g x g x于是,当[]5.1,5.0∈x 时,[]5.1,5.0)(3∈x g ,但,04432243)(232333<⎪⎭⎫⎝⎛+-⎪⎪⎭⎫ ⎝⎛--=''-x x x xx g [],5.1,5.0∈x )(3x g '在[]5.1,5.0上单调减小,因此 [][]{}3330.5,1.50.5,1.5max |()|max |(0.5)|,|(1.5)|x x g x g g ∈∈'''=.019.3)5.1(3≈'=g 因此,定理2.1的条件(2)不成立.从表2.1看到,取133074649.130=x 作为初始值0x ,128116321.131=x 作为1x .当[]3031,x x x ∈时,[]303132,31,x x x x ∈从而[]30313,)(x x x g ∈.又由于[]313033,|()|max |()|x x x g x g x ∈''≤{}331330max |()|,|()|g x g x ''= ,1853541077.0)(303<=≈'=L x g因此定理2.1的条件成立.故迭代过程收敛[]3031,x x 中任意取初值,为使解p 的近似值k x 的误差不超过810-,根据误差估计式(2.6)|,|1||01x x L L p x kk --≤- 只要.10||1801-<--x x L L k因此,k 应取为810||lg10lg1lg x x L k L ---->853541077.0lg 146458923.0128116321.1133074649.1lg 8⎪⎭⎫⎝⎛---≈.69977.137≈取138=k .于是迭代138+30=168次必可使近似解的误差不超过810-. 实际上,从表2.1中可以看到,只要迭代110次便可达到所要求解的精度.(2.6)式右端是最大可能的误差界. 就本例来说,估计的迭代次数偏大了.2.3 不动点迭代算法的收敛速度定理2.2[2] 在定理2.1的假设条件下,再设函数)(x g 在区间[]b a , 上)2(≥m 次连续可微,且在方程(2.1)的跟p 处,0)()(=p g j ,1,,1-⋅⋅⋅=m j ,0)()(≠p g m (2.15) 则不动点迭代为m 阶收敛.证明 据定理2.1知,方程(2.1)在[]b a ,上有唯一根p .且对任意初始值[]b a x ,0∈,不动点迭代序列{}k x 收敛于p 由于),()()()(11p g e p g p g x g p x e k k k k -+=-=-=++ (2.16) 据Taylor 公式和定理条件有()mkk k m m k m k k k e e p g m e p g m e p g e p g e )(!1)()!1(1)(!21)()(1121θ++-+⋅⋅⋅+''+'=--+ ,)(!1)(mk k k m e e p g m θ+=其中10<<k θ. 易知,对于充分大的k ,若 01≠-k e ,则 ),1,(0⋅⋅⋅+=≠k k i e k ,从而()).(!1lim 1p g m e e m m k k k =+∞→ (2.17)故证得不动点迭代为m 阶收敛.关于不动点的迭代,还有下面的局部收敛定理.定理2.3[2] 设p 是方程(2.1)的一个根,)(x g 在p 的某领域内m 次连续可微,且 ,0)()(=p g j ,1,,1-⋅⋅⋅=m j ,0)()(≠p g m ),2(≥m则当初始值0x 充分接近p 时(存在正数r ,对一切[]r p r p x +-∈,0),不动点迭代序列{}k x 都收敛于p ,且收敛阶数为m .证明 由于假设()x g '在p 的某领域内连续且()0='p g ,因此必存在0>r 使得对一切[]r p r p x +-∈, 有.1|)(|<≤'L x g 又据Lagrange 中值定理,有),)(()()(p x g p g x g -'=-ξ ξ在x 与p 之间,从而,|||||)(||)()(|r p x p x g p g x g ≤-<-'=-ξ 即.|)()(|r p g x g <- (2.18) 因此当[]r p r p x +-∈,时,[]r p r p x g +-∈,)(.据定理2.2和定理2.3知,对于任意初始值[]r p r p x +-∈,0,不动点迭代收敛,且收敛阶数为m .2.4 加速不动点迭代算法及其收敛性一个收敛的迭代过程将产生一个收敛序列⋅⋅⋅⋅⋅⋅,,,,21n x x x ,如p x n n =∞→lim .这样,只要迭代足够多次,即n 充分大时,如m n =,则可取m x p ≈.但若迭代过程收敛缓慢,则会使计算量变得很大,因此需要加速收敛过程.假设一个序列{}n x :⋅⋅⋅⋅⋅⋅,,,,21n x x x ,线性收敛于p (收敛缓慢),即有λ=--+∞→p x px n n n 1lim ().0≠λ (2.19) 于是当n 足够大时,有,121px px p x p x n n n n --≈---++ 即),)(()()221p x p x p x n n n --≈-++亦即.)(22222121p p x x x x p p x x n n n n n n ++-≈+-++++ (2.20) 解得nn n n n n x x x x x x p +--=++++12212222221121222n n n n n n n n n n n nx x x x x x x x x x x x ++++++-+--=-+.2)(1221nn n n n n x x x x x x +---=+++ 定义⋅⋅⋅=+---=++++,2,1,0,2)(~12211n x x x x x x x nn n n n n n , (2.21)(2.21)称为Aitken 加速公式(方法).Aitken 加速方法得到的序列{}n x ~:⋅⋅⋅⋅⋅⋅,~,,~,~21n x x x 较原来的序列{}n x 更快地收敛于p . 有下面的定理.定理 2.4[2] 设序列}{n x 是线性收敛于p 的,并且对于所有足够大的整数n 有0))((1≠--+n n x p x ,则由Aitken 加速方法(2.21)产生的序列{}n x ~有.0~lim 1=--+∞→p x px n n n (2.22) 证明 由假设序列}{n x 线性收敛于p ,即有,lim 1λ=--+∞→p x px nn n ,0≠λ.记,1λ---=+px px q n n n (2.23) 则有 0lim =∞→n n q ,0lim 1=+∞→n n q . 据(2.21)式,()⎥⎦⎤⎢⎣⎡-+----=--++++p x x x x x x p x p x p x nn n n n n n n n 1221121~2121()1()(2)n n n n n n x x x p x x x +++-=---+21221212111[()]()1()[2()]111..21n n n n n n n n n n n n n n n x p x p x p x p x p x p x p x p x p x p x p x p x p x p x p ++++++++----=-----+-⎡⎤-=--⎢⎥-⎡⎤---⎣⎦-+⎢⎥---⎣⎦ .1)(2))((1)1(112++-++⋅-+-=+λλλλn n n n q q q q (2.24)因此有.012)1(1~lim 221=+---=--+∞→λλλp x p x nn n在绪论中有讲到一阶前差:,1n n n x x x -=∆+ ⋅⋅⋅=,2,1,0n 二阶前差:,2)(122n n n n x x x x x +-=∆∆=∆++ .,2,1,0⋅⋅⋅=n 于是,Aitken 加速公式(2.21)可改写成,)(~221n n n n x x x x ∆∆-=+ .,2,1,0⋅⋅⋅=n (2.25)由于这个缘故,Aitken 加速方法又称为Aitken 2∆加速方法.例2.3[2] 设nx n 1cos =,则1lim =∞→n n x . 由于,111cos 111cos lim 11lim 1=--+=--∞→+∞→nn x x n n n n 因此序列{}nx 收敛于1. 由序列{}nx 应用Aitken 加速方法计算得{}nx ~的开头几项列表如下(表2.2).{}n x ~确实比{}n x 更快的收敛于1.第3章 非收敛不动点迭代格式的几类处理方法与比较在第2章中主要介绍了求解非线性方程的不动点迭代法,其要求是迭代函数要满足收敛定理假定条件,而在现实生活中,明确满足这些条件的迭代函数是很少见的,本章对于迭代函数不满足收敛条件的情况,提出了几类处理方法.3.1 非收敛不动点迭代格式的几类处理方法一个方程的迭代格式不是唯一的,且迭代也不都是收敛的,其收敛性取决于迭代函数)(x g 和初值0x ,关于不动点迭代函数的收敛性,上一章已经进行了讨论, 但假若[]b a x ,∈时,1)(>≥'L x g ,就不满足定理2.1的条件(2)了,于是下面分别介绍了反函数法、牛顿迭代法、Steffensen 迭代法和松弛法这四中处理方法.3.1.1 反函数法因为)(x g x =,有[])(1)(1x g x g '='-,则当[]b a x ,∈时,[]11)(1<≤'-Lx g ,所以方程)(x g x =可写成等价形式)(1x g x -=,从而构造迭代格式)(11k k x g x -+=, ),1,0(⋅⋅⋅=k (3.1) 很明显,)(11k k x g x -+=满足收敛条件.对于)(x g 简单情况, 其反函数)(1x g -容易得到.3.1.2 牛顿迭代法对于迭代格式)(x g x =的情形,采用Newton 迭代格式有 ,)(1)()()(1)(1k k k k k k k k k x g x g x x g x g x g x x x '-'-='---=+ ),1,0(⋅⋅⋅=k (3.2)3.1.3 Steffensen 迭代法根据Aitken 加速算法,对迭代格式)(1k k x g x =+,),1,0(⋅⋅⋅=k ,进行如下修改:),(k k x g y = ),(k k y g z =[]kk k k k k k k k k k k k x x g x g g x g x g g x g g x y z y z z x +---=+---=+)(2))(()())(())((2)(221 (3.3)其中⋅⋅⋅=,1,0k .3.1.4 松弛法将)(x g x =化成等价形式)()1(x wg x w x +-= , 称w 为松弛因子, 从而构造迭代格式),()1(1k k k x wg x w x +-=+ (3.4)其迭代函数为)()1()(x wg x w x g +-= . 记)(min minx g g bx a '='≤≤,)(max max x g g bx a '='≤≤,得到如下结论:(1)当1)(>≥'L x g 时,w 取01)(2max <<-'-w x g 时,)()1(1k k k x wg x w x +-=+迭代收敛;(2)当1)(-<-≤'L x g 时,w 取)(120minx g w '-<< 时,)()1(1k k k x wg x w x +-=+迭代收敛;(3)当1)(<≤'L x g 时,w 取)(1)(11minminx g x g w '-'+<< 时,迭代格式)()1(1k k k x wg x w x +-=+比迭代格式)(1k k x g x =+收敛快. 推导如下:(1)当1)(>≥'L x g 时,由01)(2max<<-'-w x g 得到2)(max->-'w x g w ,其迭代函数为)()1()(x wg x w x f +-=. 因为()()()()()()()max1111111f x w wg x w wg x f x w wg x w g x '''=-+≥-+>-'''=-+=+-<所以有1)(<'x g ,从而)()1(1k k k x wg x w x +-=+迭代收敛.(2)当1)(-<-≤'L x g 时, 由)(120min x g w '-<<得到2)(min->'+-x g w w . 因为()()()()()()()min111,1111f x w wg x w wg x f x w wg x w g x '''=-+≥-+>-'''=-+=+-<所以有1)(<'x f , 从而)()1(1k k k x wg x w x +-=+迭代收敛.(3)当1)(<≤'L x g 时, w 取)(1)(11min minx g x g w '-'+<<,由1>w 得到[]0)(1)1(<'--x g w ,()()()()()1(1)1f x w wg x w g x g x g x '''''=-+=--+<⎡⎤⎣⎦ 由)(1)(1minminx g x g w '-'+<得到0)()(1min min>'+-'+x g w w x g .()()min()11f x w wg x w wg x '''=-+≥-+()()()()minmin 1w wg x g x g x g x ''''≥-++->-所以有)()(x g x f '<', 从而迭代格式)()1(1k k k x wg x w x +-=+ 比迭代格式)(1k k x g x =+收敛快.3.2 数值实例通过以上四种方法都可以解决非收敛不动点迭代格式的问题,现对上述四种给出几个不满足不动点迭代收敛定理的实例,并对结果进行分析和比较. 例3.1 求方程033=--x x 在区间[]2,1内的根,要求精度为510-.解 对于方程033=--x x ,将它化为33-=x x ,所以3)(3-=x x g ,则当[]2,1∈x 时,13)(2>='x x g ,不满足定理2.1的条件(2),因此不能由(2.2)的迭代格式计算.下面分别用反函数方法、牛顿(Newton )迭代法、Steffensen 迭代法、松弛法对迭代函数进行修改,得到相应新的迭代函数,并用C 语言编程上机计算. (1)反函数法:迭代格式为),(11k k x g x -+= 即.)3(311+=+k k x x 取初值5.10=x ,运用程序见附录1.(2)牛顿(Newton )迭代法:迭代格式为,)(1)()()(1)(1k k k k k k k k x g x g x x g x g x g x x x '-'-='---=+即.133231)3(23231-+=----=+k k k k k k k x x x x x x x 取初值5.10=x ,运用计算程序见附录二; (3) Steffensen 迭代法:迭代格式为),(k k x g y = ),(k k y g z = [][].)(2))(()())(())((2221kk k k k k kk k k k k k x x g x g g x g x g g x g g x y z y z z x +---=+---=+即,33-=k k x y,3)3(33--=k k x z[].)3(23)3()3(3)3(3)3(3332333331kk k k kk k x x x x xx x +-----------=+ 取初值5.10=x ,运用如下程序可以得到结果: (4)松弛法:迭代格式为),()1(1k k k x wg x w x +-=+ 即),3()1(31-+-=+k k k x w x w x当[]2,1∈x 时,13)(2>='x x g ,且3)(min='x g ,12)(max ='x g ,所以w 的取值范围为01122>>--w ,现取5.10=x ,15.0-=w ,运用C 语言编程可得到起结果. 以上这四种方法的计算结果见表(3.1),本例中以上四种方法都是收敛的,因此这四种方法均可以解决不满足收敛条件的不动点迭代收敛问题,同时本例中变换后的Newton 迭代法收敛的最快.例 3.2 求方程0124=--x x 在区间[]2,1内的根,要求精度为510-.解 对于方程0124=--x x ,将它化为21214-=x x ,所以2121)(4-=x x g ,则当[]2,1∈x 时,14)(3>='x x g ,因此不满足不动点迭代收敛条件,为求此次方程的解,下面同样分别用本章介绍的四种方法求解此方程. (1)反函数法:迭代格式为),(11k k x g x -+= 将方程变为迭代格式为().12411+=+k k x x取初值5.10=x ,运行附录5的相应程序即可得计算结果. (2)牛顿(Newton )迭代法:迭代格式为,)(1)()()(1)(1k k k k k k k k x g x g x x g x g x g x x x '-'-='---=+代人例题中的数据.12212321212134341-+=-⎪⎭⎫ ⎝⎛---=+k k k k k k k x x x x x x x 取初值5.10=x ,运行附录6的程序即可的计算结果. (3)Steffensen 迭代法:迭代格式为),(k k x g y =),(k k y g z = [][].)(2))(()())(())((2221kk k k k k kk k k k k k x x g x g g x g x g g x g g x y z y z z x +---=+---=+代入例题中的数据有,21214-=k k x y ,2121212144-⎪⎭⎫ ⎝⎛-=k k x z.2121221212121212121212121212121214442444441k k k k k k k x x x x x x x +⎪⎭⎫⎝⎛---⎪⎭⎫ ⎝⎛-⎥⎥⎦⎤⎢⎢⎣⎡⎪⎭⎫ ⎝⎛---⎪⎭⎫ ⎝⎛---⎪⎭⎫ ⎝⎛-=+ 取初值5.10=x ,运行附录7即可算得计算结果. (4)松弛法:迭代格式为),()1(1k k k x wg x w x +-=+ 代入例题中的数据有.2121)1(41⎪⎭⎫ ⎝⎛-+-=+x w x w x k k当[]2,1∈x 时,14)(3>='x x g ,13224)(3max>=⨯='x g ,所以w 取值在01322<<--w ,现取05.0-=w ,初值5.10=x ,运行附录8的程序即可得到计算结果.以上这四种方法的计算结果见表(3.2),本例中以上四种方法都是收敛的,因此这四种方法均可以解决不满足收敛条件的不动点迭代收敛问题,同时本例中变换后的Newton 迭代法收敛的最快.例3.3 求方程032=-+x e x 在区间[]1,0内的根,要求精度为510-.解 将方程化为等价形式x e x 23-=,那么此时x e x g 23)(-=.当[]1,0∈x 时,12)(>-='x e x g ,因此不满足不动点迭代收敛条件.按下面这四种方法处理可以得到近似解.(1)反函数法:首先由反函数处理方法可得到迭代格式,223ln 1⎪⎭⎫⎝⎛-=+k k x x取初值5.00=x ,运用程序见附录9.(2)牛顿(Newton )迭代法:由牛顿迭代法得到迭代格式,21321kk x x k k k e e x x x +-+-=+ 取初值5.00=x ,运用程序见附录10.(3)Steffensen 迭代法:由Steffensen 迭代法得到迭代格式,23k x k e y -= ,2322⎪⎭⎫ ⎝⎛--=k x e k e z()()[](),)23(223232323223231kx e e e k x e e e e x kkx kx kx +------=---+ 取初值5.00=x ,运用程序见附录11. (4)松弛法:由松弛法得到迭代格式为(),23)1(1x k k e w x w x -+-=+当[]1,0∈x 时,122)(-<-≤-='x e x g ,e x g 2)(min-=',所以w 取ew 2120+<<之间的值,现取2.0=w ,初值5.00=x ,运用程序见附录12.以上这四种方法的计算结果见表(3.2),本例中以上四种方法都是收敛的,因此这四种方法均可以解决不满足收敛条件定理的不动点迭代收敛问题,同时本例中变换后的Newton 迭代法收敛的最快.结 论非线性代数问题的解法是现代计算数学的一个重要研究课题,而不动点迭代算法是求解非线性方程近似根的一个重要方法.本文通过搜集大量资料了解了非线性方程求解的不动点迭代算法的的研究背景,以及研究价值,然后通过介绍不动点迭代法的基本思想,对不动点迭代法的收敛性、收敛速度以及加速不动点迭代算法的收敛性进行了研究,并结合实例经行了比较分析,对于不满足收敛条件的情况,本文通过翻阅大量资料和文献,归纳总结出了四种处理方法,分别为反函数法、牛顿(Newton)迭代法、Steffensen 迭代法、松弛法,使得不满足收敛不动点迭代格式的非线性方程的求解得到了解决,同样也给出了相关实例进行了比较和验证.具体来说,对于一般的非线性方程,只要满足第2章定理2.1中的条件(1)和条件(2),那么对于任意的初始值[]b a x ,0∈,则由不动点迭代⋅⋅⋅==-,2,1),(1k x g x k k 产生的迭代序列都收敛于g 在[]b a ,的唯一不动点p ,那么要考虑的是光是收敛还不能很好的解决一个迭代效率问题,于是本文还致力研究了不动点迭代法的收敛速度,以及加速不动点迭代法.而对于一般不满足第2章定理2.1中的条件(1)或者条件(2)的情况,那么有不动点迭代算法产生的迭代序列不是收敛的,就不能求出函数的近似解,那么本文通过阅读大量的资料以及文献,总结出了四种比较好用的处理方法,并且通过3个实例,可以发现这几种方法不仅能得到收敛的迭代序列,而且收敛的速度也比较快,通过分析比较这四种方法,牛顿迭代法的迭代效果最好.这也是本文的亮点所在.由于诸多条件限制,本文也有很多不足之处,比如说没有足够多的实例来充分的证明第2章的定理2.2与定理2.3,并且对于第3章给出的3个实例中的精度要求较低,有待于继续研究,若有条件,可以更深层次的研究非线性方程组的不动点算法及其应用.参考文献[1] 李庆扬,王能超,易大义.数值分析[M].北京:清华大学出版社,2008,32(2):3-8[2] 林成森.数值分析[M].北京:科学出版社,2007,18(1):133-135,255-265[3] 王则柯.单纯不动点算法[J].广州:中山大学数学系,1988,12(12):113-127[4] 许炎.Banach空间中的不动点迭代逼近[J].苏州:苏州科技学院数理学院,2010,13(2):1-44[5] 李素红.非线性算子的不动点迭代逼近[J].天津工业大学学报,2006,16(1):51-58[6] 孙俊萍,徐裕生.非线性算子的不动点定理[J].长春大学学报,2000,10(5):30-31[7] 宋娜娜.几类非线性算子的不动点定理[J].长春工程学院学报(自然科学版),2003,4(3):1-4[8] 段华贵.几类非线性算子的不动点定理及应用[J].哈尔滨师范大学自然科学学报,2004,20(4):31-33[9] 秦小龙.非线性算子方程的迭代算法[J].科学技术与工程研究简报,2010,10(13):3169-3170[10] 王公俊.非线性方程的迭代法研究[J].河北大学学报(自然科学版),2000,20(3):228-231[11] 张卷美.一类不动点迭代法的求解[J].燕山大学学报,1998,22(2):140-142[12] 高尚.不动点迭代法的一点注记[J].大学数学.2003,19(4):30-37[13] 代璐璐.非线性方程组的迭代解法[J].合肥工业大学学报,2012,5(4):1-57[14] Halikias, Galanis, Karcanias, Milonidis.Nearest common root of polynomials,approximate greatest common divisor and the structured singular value[J].IMA Journal of Mathematical Control & Information,2013,30(4):423-442附录附录1(反函数处理法):%main()为主函数%用途:用反函数法求解非线性方程在非收敛不动点迭代格式情况下的解%格式:solut (double x,int k),solut为被调用函数,x为返回输出的值,即为迭代函数产生的迭代序列,k为在一定精度下的迭代次数#include <stdio.h>#include <math.h>double solut (double x,int k){double f;int i;for(i=0;i<k;i++){f=pow(x+3,1.0/3.0);x=f;}return(x);}main( ){double x0=1.5;int i,k;printf("\n input k(k>0):");scanf("%d",&k);for(i=1;i<=k;i++){printf("\n x=%6.5lf\n",solut(x0,i));}}附录2(牛顿(Newton)迭代法):%main()为主函数%用途:用牛顿(Newton)迭代法求解非线性方程在非收敛不动点迭代格式情况下的解%格式:solut (double x,int k),solut为被调用函数,x为返回输出的值,即为迭代函数产生的迭代序列k为在一定精度下的迭代次数#include <stdio.h>#include <math.h>double solut(double x,int k){double f;int i;for(i=0;i<k;i++){f=x-(pow(x,3)-x-3)/(3*pow(x,2)-1);x=f;}return(x);}main( ){double x0=1.5;int i,k;printf("\n input k(k>0):");scanf("%d",&k);for(i=1;i<=k;i++){printf("\n x=%6.5lf\n",solut(x0,i));}}附录3(Steffensen迭代法):%main()为主函数%用途:用Steffensen迭代法求解非线性方程在非收敛不动点迭代格式情况下的解%格式:solut (double x,int k),solut为被调用函数,x为返回输出的值,即为迭代函数产生的迭代序列,k为在一定精度下的迭代次数#include <stdio.h>#include <math.h>double solut(double x ,int k){double f,f1,f2;Int i;for(i=0;i<k;i++){f1=pow(x,3)-3;f2=pow(f1,3)-3;f=f2-(pow(f2-f1,2)/(f2-2*f1+x);x=f;}return(x);}main( ){double x0=1.5;int i,k;printf(“\n input k(k>0):”);scanf(“%d”,&k);for(i=1;i<=k;i++){printf(“\n x=%6.5lf\n”,solut(x0,i);}}附录4(松弛法):%main()为主函数%用途:用松弛法求解非线性方程在非收敛不动点迭代格式情况下的解%格式:solut (double x,int k),solut为被调用函数,x为返回输出的值,即为迭代函数产生的迭代序列,k为在一定精度下的迭代次数#include <stdio.h>。

切比雪夫不等式及其应用(论文)

切比雪夫不等式及其应用(论文)

第一章 绪论概率论是一门研究随机现象数量规律的科学,是近代数学的重要组成部分。

随机现象在自然界和人类生活中无处不在,因而大多数的应用研究,无论是在工业、农业、经济、军事和科学技术中,其本质都是现实过程中的大量随机作用的影响。

这个观点强有力地推动了概率论的飞速发展,使其理论与方法被广泛地应用于各个行业。

而概率论极限理论的创立更使其锦上添花,以至在近代数学中异军突起。

历史上第一个极限定理是属于雅各布·伯努利,后人称之为“大数定律”。

因其遗著《猜度术》于1713 年出版,故概率史家称1713 年为伯努利大数定律创立年。

伯努利大数定律给出了频率估计概率的理论依据,同时开创了概率论中极限理论的先河,标志着概率论成为独立的数学分支。

1837年,泊松对大数定律提出一个较宽松的条件,进而得到泊松大数定律。

之后,由于有些数学家过分强调概率论在伦理学中的应用,又加上概率论自身基础不牢固,大多数数学家往往把概率论看作是有争议的课题,排除在精密科学之外。

切比雪夫正是在概率论门庭冷落的年代从事其研究的。

切比雪夫在1866年发表的论文《论均值》中,提出了著名的切比雪夫大数定律。

该论文给出如下三个定理[1]:定理1.1:若以 ,,,c b a 表示 ,,,z y x 的数学期望,用 ,,,111c b a 表示相应的平方,,,222z y x 的数学期望,则对任何α, +++z y x 落在----+++++++222111c b a c b a c b a α和----+++-+++222111c b a c b a c b a α之间的的概率总小于211α-定理1.2:若以 ,,,c b a 表示 ,,,z y x 的数学期望,用 ,,,111c b a 表示相应的平方,,,222z y x 的数学期望,则不论t 取何值,N 个量 ,,,z y x 的算术平均值和他们相应的数学期望的算术平均值的差不超过Nc b a N c b a t+++-+++2221111 的概率对任何t 都将大于Nt 21-。

浅谈高中数学中最值问题的常用解题方法(数学本科毕业论文)

浅谈高中数学中最值问题的常用解题方法(数学本科毕业论文)

福建师范大学现代远程教育毕业论文题 目: 浅谈高中数学中最值问题的常用解题方法学习中心: 灌 云 奥 鹏 专 业: 数学及应用数学 年 级(入学批次): 201103 学 号: ************ 学生姓名: * * 导师姓名: 严 晓 明2013 年 3月 15 日装 订 线浅谈高中数学中最值问题的常用解题方法201103896627 刘明 指导老师:严晓明摘要: 最值问题是中学数学的重要题型之一。

以最值问题为载体,可以考查中学数学的几乎所有知识点,可以考查分类讨论、数形结合、转化与化归等诸多数学思想和方法,还可以考查学生的思维能力、实践和创新能力。

解决最值问题,从方法上来说,它常用到函数的单调性、二次函数的性质、数形结合法、均值不等式法、导数法、换元法等等。

本文就高中数学的要求,结合一些典型试题进行分析和探讨,说明其解题的思考方法和一般的技能与技巧。

关键词:高中数学 最值 解题方法1、引言在日常生活及科学实验中,常常遇到“最好”、“最省”、“最大”、“最小”、“最低”等问题。

例如质量最好,用料最省,效益最高,成本最低,利润最大,投入最小等等,这类问题在数学上常常归结为求函数的最大值或最小值问题,也就是最值问题.最值问题是一类综合性较强的问题,其题型多样,解法灵活,在高中数学中,最值问题涉及面广,像函数(三角函数,二次函数,指对函数,幂函数),不等式,向量,解析几何,立体几何,圆锥曲线中都能找到最值问题,在高考中,常以一些基础题,小综合的中档题或一些难题的形式出现,是历年高考重点考查的知识点之一,几乎每年的高考试题中都有出现。

2、最大(小)值及其几何意义一般地,设)(x f y =的定义域为A ,如果存在A x ∈0,使得对于任意的A x ∈,都有)()(0x f x f ≤,那么称)(0x f 为)(x f y =的最大值,记为)(0max x f y =;如果存在A x ∈0,使得对于任意的A x ∈,都有)()(0x f x f ≥,那么称)(0x f 为)(x f y =的最小值,记为)(0min x f y =.其几何意义是:函数图象上最高(低)点的纵坐标。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

本科优秀数学本科毕业论文Last revision on 21 December 2020***大学2014 届本科毕业论文论文题目:行列式的计算及应用学生姓名: ***所在院系:数学科学学院所学专业:数学与应用数学(金融方向)导师姓名: ***完成时间: ***年***月***日行列式的计算及应用摘要在高等代数这门课程里,行列式是最基本而又重要的内容之一,同时也是数学研究中的重要的工具之一,在线性代数、数学分析、解析几何等众多课程理论中以及实际问题中许也发挥着重要作用,了解如何计算和应用行列式显得尤为重要。

本文首先阐述行列式的基本理论,在此研究的基础上介绍了降阶法,归纳法,化三角形法等几种常见的且有一定技巧的解行列式的方法,并列举了相关的例子,更直观地了解解行列式方法的精髓。

另外,本文又介绍了行列式在解析几何、代数及其他课程当中的应用,进一步加深了对行列式的理解。

最后本文又列举实例阐述行列式在实际当中的应用,实现了行列式的理论与实际相结合。

研究行列式的计算方法及其应用可以提高对行列式的认识,有利于把行列式的研究推向深入。

通过这一系列的方法可以进一步提升对行列式的认识,为以后学习奠定了基础。

关键词:行列式,因式分解,化三角形法, 归纳法,加边法,Matlab软件Determinant calculation and applicationAbstractThis course in advanced algebra, the determinant is one of the most basic and important content, while many math curriculum theory is one of the important research tools, linear algebra, mathematical analysis, analytic geometry, etc. as well as practical problems also plays an important role in understanding how to calculate and apply the determinant is particularly important.This paper first describes the basic theory of determinants, based on this study describes the reduction method, induction techniques and a certain common determinant of several methods of solution method, the method of the triangle, and cited relevant examples, more intuitive understanding of the essence of the solution determinant method. In addition, thispaper describes the determinant in analytic geometry, algebra and other courses which further deepened the understanding of the determinants. Finally, they provide examples described determinant application in practice to achieve a theoretical and practical determinant combined. Research determinant calculation method and its application can improve the understanding of the determinant, is conducive to deepen the study of determinants. You can further enhance the understanding of the determinants through this series of methods, laid the foundation for future learning.Keywords: determinants, factorization of a triangle, induction, plus side method, Matlab software目录1. 行列式的定义及性质 行列式的定义排列[1]在任意一个排列中,若前面的数大于后面的数,则它们就叫做一个逆序,在任意一个排列中,逆序的总数就叫做这个排列的逆序数.定义[1]n 阶行列式就相当于全部不同行、列的n 个元素的乘积nnj j j a a a 2121 (1-1-1)的代数和,这里n j j j 21是n ,,2,1 的一个排列,每一项(1-1-1)都按下列规则带有符号:当n j j j 21是偶排列时,(1-1-1)是正值,当n j j j 21是奇排列时,(1-1-1)是负值.这一定义可以表述为n nn nj j j j j j j j j nnn n nna a a a a a a a a a a a D21212121)(212222111211)1(∑-==τ, (1-1-2) 这里∑nj j j 21表示对所有n 级排列求和.由于行列指标的地位是对称的,所以为了决定每一项的符号,我们也可以把每一项按照列指标排起来,所以定义又可以表述为n i i i i i i i i i nnn n nnn n a a a a a a a a a a a a D21)(212222111211212121)1(∑-==τ. (1-1-3)行列式的相关性质记 nnn n n n a a a a a a a a a D 212222111211=,nnn nn n a a a a a aa a a D 212221212111'=,则行列式'D 叫做行列式D 的转置行列式.性质1 行列式和它的转置行列式是相等的[2]. 即D D ='. 证明:记D 中的一般项n 个元素的乘积是它处于D 的不同行和不同列,所以它也处于'D 的不同行和不同列,在'D 中应是所以它也是'D 中的一项.反之, 'D 的每一项也是D 的一项,即D 和'D 有相同的项.再由上面(1-2)和(1-3)可知这两项的符号也相同,所以D D ='.性质2 nnn n in i i nnn n n in i i n a a a a a a a a a k a a a ka ka ka a a a212111211212111211=. 证明:inin i i i i nnn n in i i n A ka A ka A ka a a a ka ka ka a a a +++=2211212111211性质3 如果行列式的某行(列)的元素都为两个数之和[2],如nnn n nn n a a a c b c b c b a a a D 21221111211+++=,那么行列式D 就等于下列两个行列式的和: 可以参照性质2的证明得出结论.性质4 对换行列式中任意两行的位置,行列式值相反.即若设 则.1D D -=证明:记D 中的一般项中的n 个元素的乘积是它在D 中处于不同行、不同列,因而在1D 中也处于不同行、不同的列,所以它也是1D 的一项.反之,1D 中的每一项也是D 中的一项,所以D 和1D 有相同的项,且对应的项绝对值相同.现在看该项的符号:它在D 中的符号为由于1D 是由交换D 的i 、k 两行而得到的,所以行标的n 级排列n k i 12变为n 级排列n k i 12,而列标的n 级排列并没有发生变化.因此D 和1D 中每一对相应的项绝对值相等,符号相反,即.1D D -=性质5 如果行列式中任有两行元素完全相同,那么行列式为零.证明:设该行列式为D ,交换D 相同的那两行,由性质4可得D D -=,故.0=D性质6 如若行列式中任有两行或者两列元素相互对应成比例,则行列式为零.证明:设n 阶行列式中第i 行的各个元素为第j 行的对应元素的k 倍,由性质2,可以把k 提到行列式外,然后相乘.则剩下的行列式的第i 行与第j 行两行相同,再由性质5,最后得到行列式为零.性质7 把任意一行的倍数加到另一行,行列式的值不改变.nnn n kn k k in i i n a a a a a a a a a a a a 21212111211=.2. 行列式的计算方法 几种特殊行列式的结果三角行列式nn nn nna a a a a a a a a 221122*********=(上三角行列式).nn nnn n a a a a a a a a a2211212221110=(下三角行列式). 对角行列式nn nna a a a a a22112211000=. 定义法例1 用定义法证明.00000000002121215432154321=e e d d c c b b b b b a a a a a 证明:行列式的一般项可表成.5432154321j j j j j a a a a a 列标543,,j j j 只能在5,4,3,2,1中取不同的值,故543,,j j j 三个下标中至少有一个要取5,4,3中的一个数,则任意一项里至少有一个0为因子,故任一项必为零,即原行列式的值为零.利用行列式的性质计算例2 一个n 阶行列式ij n a D =的元素都满足n j i a a ji ij ,,2,1,, =-=, 那么n D 叫做反对称行列式,证明:奇数阶的反对称行列式的值等于0.证明:由ji ij a a -=知ii ii a a -=,即n i a ii ,,2,1,0 ==所以行列式n D 可写为0000321323132231211312 n n nn n nn a a a a a a a a a a a a D ------=,再由行列式的性质2,'A A =得到0000000321323132231211312321323132231211312nnnnn n n nnn nn n a a a a a a a a a a a a a a a a a a a a a a a a D ------=------=n n nn n n nn n D a a a a a a a a a a a a )1(0000)1(321323132231211312-=-------= ,当n 为奇数时,得n n D D -=,因而得到0=n D .降阶法例3 计算)2(≥n n 级行列式xy y x y x y xd 000000000000=. 解:按第一列展开得到原式阶阶)1(1)1(000000000)1(0000000000000-+--⨯+=n n n y xy y x y y x yx y x y x x)2()1()1(≥-+=+n y x n n n .归纳法形如行列式叫做n 阶范德蒙(Vandermonde )行列式.下面证明,对每一个)2(≥n n ,n 阶范德蒙行列式就等于n a a a ,,,21 这n 个数的所有可能的差)1(n i j a a j i ≤<≤-的乘积.用数学归纳法证明范德蒙德行列式 我们对n 作归纳法. (1)当2=n 时,122111a a a a -=,结果是对的.(2)设对于1-n 级的范德蒙行列式,结论是成立的,先来看n 级的情况.在中,第n 行减第1-n 行的1a 倍,第1-n 行减第2-n 行的1a 倍,即由下而上逐次地从每一行减它上一行的1a 倍,得到22322223223211312111)())((------=n nn n n n n a a a a a a a a a a a a a a a. 最后面这个行列式是1-n 级范德蒙德行列式,再由归纳法假设,它的值就是)1(n i j a a j i ≤<≤-;而所有带有1a 的差即为上式最后等式行列式的前面.所以,结论对n 级范德蒙德行列式也是成立的.由数学归纳法,证明了结论.用连乘号,这个结果可以简写为∏≤<≤-----==ni j j i n nn n n nn n a a a a a a a a a a a a a a D 1113121122322213211111)(. (2-5-1)递推法给定一个递推关系式,再给定某一个较低阶初始行列式的值,就可递推求得所给n 阶行列式的值,运用这种方法计算的方法就叫做递推法。

相关文档
最新文档