为了正常的体验网站,请在浏览器设置里面开启Javascript功能!
首页 > 数值分析论文

数值分析论文

2019-11-18 3页 doc 191KB 6阅读

用户头像

is_704284

暂无简介

举报
数值分析论文《数值分析与科学计算概述》研究第一章对象描述一、数值分析与科学计算的概念科学计算即数值计算,科学计算是指应用计算机处理科学研究和工程技术中所遇到的数学计算。在现代科学和工程技术中,经常会遇到大量复杂的数学计算问题,这些问题用一般的计算工具来解决非常困难,而用计算机来处理却非常容易。科学计算是一门工具性、方法性、边缘性的学科,发展迅速,它与理论研究和科学实验成为现代科学发展的三种主要手段,它们相辅相成又互相独立,在实际应用中导出的数学模型其完备形式往往不能方便地求出精确解,于是只能转化为简化模型求其数值解,如将复杂的非线性模型忽...
数值分析论文
《数值分析与科学计算概述》研究第一章对象描述一、数值分析与科学计算的概念科学计算即数值计算,科学计算是指应用计算机处理科学研究和工程技术中所遇到的数学计算。在现代科学和工程技术中,经常会遇到大量复杂的数学计算问题,这些问题用一般的计算工具来解决非常困难,而用计算机来处理却非常容易。科学计算是一门工具性、性、边缘性的学科,发展迅速,它与理论研究和科学实验成为现代科学发展的三种主要手段,它们相辅相成又互相独立,在实际应用中导出的数学模型其完备形式往往不能方便地求出精确解,于是只能转化为简化模型求其数值解,如将复杂的非线性模型忽略一些因素而简化为可以求出精确解的线性模型,但这样做往往不能满足近似程度的要求,因此使用数值方法直接求解做较少简化的模型,可以得到满足近似程度要求的结果,使科学计算发挥更大的作用。自然科学规律通常用各种类型的数学方程式表达,科学计算的目的就是寻找这些方程式的数值解。这种计算涉及庞大的运算量,简单的计算工具难以胜任。在计算机出现之前,科学研究和工程主要依靠实验或试验提供数据,计算仅处于辅助地位。计算机的迅速发展,使越来越多的复杂计算成为可能。利用计算机进行科学计算带来了巨大的经济效益,同时也使科学技术本身发生了根本变化:传统的科学技术只包括理论和试验两个组成部分,使用计算机后,计算已成为同等重要的第三个组成部分。数值分析也称计算方法,它与计算工具发展密切相关。是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科,是数学的一个分支,它以数字计算机求解数学问题的理论和方法为研究对象。为计算数学的主体部分。在电子计算机出现以前,计算工具只有算盘,算图,算表和手摇及电动计算机。计算方法只能计算规模较小的问题。数值分析的任务是研究求解各类数学问题的数值方法和有关理论的学科。数值分析的过程为构造算法、使用算法、分析算法。数值分析是研究数值问题的算法,概括起来有四点:第一,面向计算机,要根据计算机的特点提供切实可行的计算方法。即算法只能包括加、减、乘、除运算和逻辑运算,这些计算是计算机直接能处理的运算。第二,有可靠的理论分析,能任意逼近并达到精度要求,对近似算法要保证收敛性和数值稳定性,还要对误差进行分析,这些都建立在相应数学理论的基础上。第三,要有好的计算复杂性,时间复杂性好是指节省计算时间,空间复杂性好是指节省存储空间,这也是建立算法要研究的问题,它关系到算法能否在计算机上实现。第四,要有数值实验,即任何一个算法除了从理论上要满足上述三点外,还要通过数值试验证明是行之有效的。二、数值计算的相关概念与理论误差来源与分类所谓误差就是测得值与被测量的真值之间的差,可以表示为:误差=测得值-真值误差来源:(一)测量装置误差(二)环境误差(三)方法误差(四)人员误差误差分类:将无限位字长的精确数处理成有限位字长近似数的处理方法称为摄入方法,称为舍入误差。用计算机解决科学计算问题首先要建立数学模型,它是对被描述的实际问题进行描述、简化而得到的,因而是近似的。我们把数学模型与实际问题之间出现的这种误差称为模型误差。在数学模型中往往还有一些根据观测得到的物理量,如温度、长度、电压等,这些参量显然包含误差,这种由观测产生的误差称为观测误差,在数值分析中不讨论这种误差,数值分析只研究用数值方法求解数学模型产生的误差。当数学模型不能得到精确解时,通常要用数值方法求它的近似解,其近似解与精确解之间的误差称为截断误差或方法误差。绝对误差与相对误差某量值的测得值和真值之差为绝对误差。通常简称为误差,即:绝对误差=测得值-真值设x*为准确值x的一个近似值,称e二x-x*为近似值x*的绝对误差(AbsoluteError),简称误差。当同时,则称e为x*的一个绝对误差限。设x*为准确值x的近似值,称绝对误差与准确值之比为近似值x*的相对误差(RelativeError),记为er,即ex-x*e==rxx若leJ-£r,则称er为X*的一个相对误差限。由于绝对误差可能为正值或负值,因此相对误差也可能为正值或负值。对于相同的被测量,绝对误差可以评定其测量精度的高低,但对于不同的被测量以及不同的物理量,绝对误差就难以评定其测量精度的高低,而采用相对误差来评定较为确切。例1已知近似数x*有两位有效数字,试求其相对误差限。解利用有效数字与相对误差的关系。可知n=2且这里ai是1到9之间的数字。(x卜2aX10-n+1<12X1X10-2+1=5%有效数字含有误差的任何近似数,如果其绝对误差界是最末位的半个单位,那么从这个近似数左方起第一个非零的数字,称为第一位有效数字。如果近似值x*的误差限是0.5X10-n,则称x*准确到小数点后第n位,并从第一个非零数字到这一位的所有数字均称为有效数字。22例2问3.142,3.141,7分别作为兀的近似值各具有几位有效数字?因为兀=3.14159265…22x1=3.142,x2=3.141,X3=7兀一x=3.14159…一3.142=—0.00040…矢口2x10-4<一x1<2X10-3因而xi具有4位有效数字。由兀—兀2=3.14159…—3.141=—0.00059…知因而12X10-3<1:兀一x21I配方
、化学反应速度、污染物扩散等都可以作为数值分析的模型实例。通过选取一些贴近事实的问题,使得数值分析与工程实际有机结合,实现数值分析与科学研究零距离对接。五、数学软件在众多的数学软件中,就软件的内核而言,可分为两大类。一类是数值计算型软件,如Matlab、Xmath、Gauss等,这类软件有较强的计算能力和可视化功能,运行效率较高。另外一类是数学分析型软件,如Mathematica、Maple等,它们处理符号运算的功能特别强,可得到解析符号解和任意精度解。MATLAB是MATrixLABoratory的缩写,即矩阵实验室,它整合了非线性方程组、数值积分、三次样条函数、曲线拟合、最优化、常微分方程和绘图工具等功能,但它主要是以EISPACK和LINPACK子程序为基础。MATLAB目前是由C和汇编语言编写的,它的基本结构是执行矩阵运算,是一个对求解线性方程特别有用的功能强大的自包容系统。Mathematica的原始系统是由美国物理学家StephenWolfram领导的一个小组开发来进行量子力学研究的,软件开发的成功促使StephenWolfram推出系统的不同版本。目前国内最常见的是3.0版本和4.0版本。迄今为止,Mathematica已经广泛的应用到工程、应用数学、计算机科学、财经、生物、药学、生命科学以及太空科学等领域。数以万计的、科学、期刊杂志、图书资料、计算机绘图等都是Mathematica的杰作。Mathematica是一个交互式、集成化的计算机软件系统。它的主要功能包括三个方面:符号演算、数值计算、图形功能。4.0版本对Mathematica系统做出了重大的改进,尤其是在处理大容量数据方面,比以前版本增强了许多。另外,在代数、语言和用户界面等方面,也做出了很大的改进,使其运行效率更好,执行速度更快。第二章算法研究一、插值法的定义目前,主要的插值法有Lagrange插值、Newton插值、Hermite插值、三次样条插值。设函数y二f(x)在区间[a,b]上有定义,且已知在点a100Abs[%-g[100]]//N差商(DividedDifference)的定义设有函数f(X),X,X,X,为一系列互不相等的点,称012(i丰j)X-Xijf[x,x]-f(叮-f(X)ij为f(x)关于点x,x的一阶差商(也称均差)。ij一般的,称f[x,x,…,x]-f[x,x,…,x]TOC\o"1-5"\h\zf[x,x,…,x]=0_1k-1+2k-01kx-x0k为f(x)关于点x,x,…,x的k阶差商。01k6.差商的性质(1)各阶差商均具有线性性质,即若p[x,x,…,x],V[x,x,…,x]分别为01k01k申(x),屮(x)关于点x,x,…,x的k阶差商,且f(x)二a^(x)+b屮(x),则有01kf[x,x,…,x]二ap[x,x,…,x]+b屮[x,x,…,x]01k01k01k(2)k阶差商f[x,x,…,x]可表示成f(x),f(x),•••,/(x)的线性组合,即01k01k沪f(x)f[x,x,…,x]=Z^―01kw(x)i=0k+1i其中的w'(x)=(x一x)。n+1iijj=0j知(3)若f(x)是n次多项式,则一阶差商f[x,x]是n-1次多项式。j7.Newton插值多项式N(x)=f(x)+(x-x)f[x,x]+(x-x)(x-x)f[x,x,x]+—TOC\o"1-5"\h\zn00001012+(x—x)(x—x)…(x—x)f[x,x,…,x]01n-101n称为Newton插值多项式。插值误差为R(x)=(x—x)(x—x)…(x—x)f[x,x,x,…,x]n01n01n=w(x)f[x,x,x,…,x]n+101n例2、已知f(x)的函数表如下,求向前差分表,并写出Newton向前插值公式。x0.20.40.60.81.01.2f(x)0.560.791.682.542.873.44解构造向前差分表的Mathematica程序为:data1={{0.2,0.56},{0.4,0.79},{0.6,1.68},{0.8,2.54},{1.0,2.87},{1.2,3.44}}deltaF[p_]=Module[{myD},myD[{a1_,a2_}]:=a1[[2]]-a2[[2]];myD[{a1_,r__,a3_}]:=myD[{r,a3}]-myD[{a1,r}];myD[p]]Table[Map[deltaF,Partition[data1,i,1]],{i,2,6}]//TableForm//Chop运行结果:-0.23-0.89-0.86-0.33-0.57TOC\o"1-5"\h\z-0.660.030.53-0.240.690.5-0.77-0.19-1.27-1.08求出Newton向前插值多项式expr=0.56-0.23t-0.66*t(t-1)/2+0.69*t(t-1)(t-2)/2-0.19*t(t-1)(t-2)(t-3)/2-1.08*t(t-1)(t-2)(t-3)(t-4)/2//Expand%/.t->10x8.Hermite插值多项式设已知函数y=f(x)在n+1互异的节点aSetTitle(Win_Title);//设置画图环境cclientDC.dlgDC(this)//定义变量intItn=100000,I,NUM=9,RandomN=0;D[]=[0.2,0.5,0.7];doubleNewX,NewY,XO,YO,/定义插值点X[]=[0,240,330,600],Y[]=[0,120,240,90],A[9],B,C[9],D[9],E[9],F[9];//定义数据计算出剪切系数For(i=1;i<=NUM;i++){B=X[NUM]-X[0];A[i]=(X(i)-X[i-1])/BE[i]=(X[NUM]*X[i-1]-X[0]*X[i])/BC[i]=(Y[i]-Y[i-1]-D[i]*(Y[NUM]-Y[0]))/BF[i]=(X[NUM]*Y[i-1]-X[0]*Y[i]-D[i]*(X[NUM]*Y[0]-x[0]*Y[NUM]))/B;}dlgDC.MoveTo(50,50);dlgDC.LineTo(50,500);dlgDC.LineTo(750,500);i=0;XO=0;YO=0;//迭代画图while(i0,令x二x+Ax个校正量,称为增量,于是x3-a=0转化为(x+Ax)3=a艮卩X3+3x2Ax+3x(Ax)2+(Ax)3=a000由于Ax是小量,若省略高阶项(心)2,则得xJ12ax—3x4TOC\o"1-5"\h\zAx沁一f±o—x3+3x2Ax+3x(Ax)2ua26x000,艮01J12ax一3x4x=x+Axux±oo-02o6x二x于是o1这里xi不是轟,但它是真值x皿的进一步计算,重复以上过程可得到迭112ax一3x4x=x±kJ代公式k+12k6xk,k=°丄2,…,limx=x*厂(i它可以逐项的求得*S…'若—k,则x*=3a,容易证明序列{»}对任何xo>0均收敛,且收敛很快。11•判断下列命题的正确性:(1)解对数据的微小变化高度敏感是病态的。(正确)(2)高精度运算可以改善问题的病态性。(正确)(3)无论问题是否病态,只要算法稳定都能得到很好的近似值。(错)(4)用一个稳定的算法计算良态问题一定会得到好的近似值。(错)(5)用一个收敛的迭代法计算良态问题一定会得到好的近似值。(错)(6)两个近似数相减必然会使有效数字损失。(错)(7)计算机将1000个数量级不同的数相加,不管次序如何结果都是一样的(错)二、我的课题作业13.f(x)=ln(x-、;x2_1),求f(30)的值,若开平方用6位函数表,问求对数时误差有多大?若改用另一等价公式一ln(x*x2一1)=ln(x一、'x2一1)计算,求对数时误差有多大?解:设y=x—.x2-1,则f(x)二上£(f*)由公式:s(f(x*))uf'(x*).£(x*)得:当x=30时,「=J-•X:-1.二门二;由题意及公式:x-x*|<1X10m-n+1当取四位有效值时,得到:y*-y<^x10-4,y*=0.0167,2>所以:y*-y£(f*)uy*-x10-4<「-0.0167u0.3x10-2。若用等价公式:ln(x-Jx2-1)=-ln(x+Jx2+1)计算,则y*=59.9833此时:y*-y£(f*)u
/
本文档为【数值分析论文】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。 本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。 网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。

历史搜索

    清空历史搜索