基于matlab的数值计算中的优化技术_毕业论文.doc

《基于matlab的数值计算中的优化技术_毕业论文.doc》由会员分享，可在线阅读，更多相关《基于matlab的数值计算中的优化技术_毕业论文.doc（22页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、-Hefei University 毕业论文（设计）BACHELOR DISSERTATION论文题目： 基于matlab的数值计算中的优化技术 学位类别： 理 学 学 位 学科专业： 信 息 与 计 算 科 学 基于matlab的数值计算中的优化技术中文摘要优化是人们寻求的目标，数值计算中优化技术采用的好，能从时间与空间上得到巨大的好处。一个算法除了正确外，还要空间能存贮程序数据，且运行时间短，因而算法优化技术就很重要，一个程序无法存贮到计算机内存中，或运行时慢得无法等候是没有任何实际意义的。由于数值计算的优化技术有很多方面，为此选用数值积分进行说明，数值积分计算有很多的方法，用不同的方法所

2、计算的积分精确度不同，所需要的时间也不同，通过一些实例的分析，对优化技术进行归纳与总结。本论文是基于matlab在数值计算中的优化技术，优化技术是算法设计的重要而关键的课题，本论文选取数值分析中的一些著名的优化技术进行讨论，并在matlab中加以实现，通过tic、tuc、cputime等函数的使用对其进入深入分析。关键字：数值计算；优化技术；matlab；数值积分Numerical optimization technique based on MATLABABSTRACTOptimization is that people seek target and numerical optimiz

3、ation technique used is good. It can get huge benefits from the time and space. An algorithm not only in order to right, but also space can store data, and short running time. So the algorithm optimization technique is very important and a program cannot be stored in the computer memory, or run slow

4、 and practical of significance. The optimization technology of numerical calculation is used in many respects. The numerical integration is described and numerical integral calculation has many methods. Integral accuracy is calculated by different methods and the time required is different also, and

5、 through the analysis of some examples, and summarizes the optimization technique.This paper is to optimize the technology based on MATLAB in the numerical calculation. Optimization technique is the key task in the algorithm design, this paper selects some well-known optimization techniques in numer

6、ical analysis are discussed, and implemented in MATLAB, through in-depth analysis using tic, tuc and cputime and other functions of the entry.KEYWORD: numerical calculation; optimization techniques;MATLAB;numerical integration第一章前言1第二章 数值积分的计算22.1 数值求积公式的构造22.1.1求积公式的推导22.1.2几个低次牛顿-科特斯求积公式42.2 复化求积公

7、式62.2.1复化梯形求积公式62.2.2复化辛浦生求积公式62.2.3复化科特斯求积公式72.3 高精度数值积分算法72.3.1 龙贝格求积公式8第三章 线性方程组的求解103.1线性方程组的介绍103.2 线性方程组的迭代法113.2.1 Jacobi迭代法113.2.2 Gauss-Seidel迭代法123.2.3 SOR迭代法13第四章 各种求积公式的MATLAB编程实现与应用144.1 对数值积分运行结果及其分析144.1.1 数值积分运行结果144.1.1 数值积分运行结果分析154.1 线性方程组运行结果及其分析154.2.1 线性方程组运行结果154.2.2 线性方程组运行结果

8、分析16附录17参考文献23致 谢24【精品文档】第 19 页第一章 前言数值计算是有效使用数字计算机求数学问题近似解的方法与过程，以及由相关理论构成的学科。数值计算主要研究如何利用计算机更好的解决各种数学问题，包括连续系统离散化和离散形方程的求解，并考虑误差、收敛性和稳定性等问题。从数学类型分，数值运算的研究领域包括数值逼近、数值微分和数值积分、数值代数、最优化方法、常微分方程数值解法、积分方程数值解法、偏微分方程数值解法、计算几何、计算概率统计等。随着计算机的广泛应用和发展，许多计算领域的问题，如计算物理、计算力学、计算化学、计算经济学等都可归结为数值计算问题。优化是人们寻求的目标，数值计

9、算中优化技术采用的好，能从时间与空间上得到巨大的好处。一个算法除了正确外，还要空间能存贮程序数据，且运行时间短，因而算法 优化技术就很重要，一个程序无法存贮到计算机内存中，或运行时慢得无法等候是没有任何实际意义的。 在研究基于matlab在数值计算中的优化技术有很多方面求数值积分就是具有代表性的一点。 求某函数的定积分时，在多数情况下，被积函数的原函数很难用初等函数表达出来，因此能够借助微积分学的牛顿-莱布尼兹公式计算定积分的机会是不多的。另外，许多实际问题中的被积函数往往是列表函数或其他形式的非连续函数，对这类函数的定积分，也不能用不定积分方法求解。由于以上原因，数值积分的理论与方法一直是计

10、算数学研究的基本课题。通过这个课题的研究，我们将会更好地掌握运用数值积分算法求特殊积分函数的定积分的一些基本方法、理论基础；并且通过matlab软件编程的实现，得出计算数值积分的最优化的方法，并应用于实际生活中。第二章 数值积分的计算2.1 数值求积公式的构造人们根据积分的定义得出Newton-Leibniz求定积分的公式,但是这些公式并不是能求出所有式子的积分,而是针对许多特殊的例子,但是有许多都是球不出来的,如: ,等。所以采用积分的几何意义来设计出积分公式从而求出近似值。2.1.1求积公式的推导建立数值积分公式的途径比较多，其中最常用的优两种：（1）对于连续函数，优积分中值定理：其中是被

11、积函数在积分区间上的平均值。因此，如果能给出求平均值的一种近似方法，相应地就可以得到计算定积分的一种数值方法。（2）先用某种简单函数近似逼近，然后在区间的积分值近似表示在区间上的定积分，即取一般情况下，我们可以取为前面介绍的插值多项式或拟合多项式进行近似计算。若取为插值多项式，则相应得到的数字微分公式就是插值型求积公式把区间等分，其分点为、，过这个节点，可以构造一个次插值多项式：其中，用代替被积分函数则有其中。上式叫做牛顿-科茨公式，使用牛顿-科茨公式关键是计算系数，用变量替换，于是而 这样 记 则 这时是不依赖于函数和区间的常数，可以事先计算出来，叫做牛顿-科茨系数。 表2.1牛顿-科茨系数

12、2.1.2几个低次牛顿-科特斯求积公式1、梯形求积公式定义2.1 在牛顿-科特斯求积公式中，如果取，用一次多项式代替被积函数，即用梯形面积代替曲边梯形的面积，则有其中，,查表可得代入上式得出称式为梯形求积公式根据牛顿-科特斯求积公式的误差理论，梯形求积公式的误差估计为是被积函数二阶导数在点的取值，2、辛浦生求积公式定义2.2 在牛顿-科特斯求积公式中，如果取，用二次多项式代替被积函数，即曲边用抛物线代替，则有其中，,查表可得，代入上式得出称式为辛浦生求积公式，也称抛物线求积公式。辛浦生求积公式的误差估计式3、科特斯求积公式定义 2.4 在牛顿-科特斯求积公式中，如果取时，牛顿科特斯公式为称式为

13、科特斯求积公式。同理可求得其误差估计式2.2 复化求积公式2.2.1复化梯形求积公式在上一节求积分的过程只是求粗约的近似值，所以应根据积分的可加性，可以将区间分为许多部分使得积分值更加接近精确值，从而优化了梯形积分公式，辛普生积分公式和科特斯积分公式，这就是复化求积分公式的思想。定义 2.5 将积分区间进行等分，记为,在每个小区间 上用梯形公式求和，得若将所得的近似值记为，整理得称式为复化梯形公式。记为 复化梯形公式的截断误差2.2.2复化辛浦生求积公式在辛普生积分公式上加以复化可以得到复化辛普生积分公式。定义2.6 将积分区间分成等分，分点为， 在每个小区间上。用Simpson公式求积分，得

14、到式就称为复化辛浦生求积公式。记为 如果， 则由Simpson插值余项公式可得复化公式的截断误差为2.2.3复化科特斯求积公式定义2.7 将积分区间等分为个子区间，每个子区间的中点,，子区间长度， 在每个子区间上用科特斯公式求和，得 式就称为复化科特斯求积公式，式中类似地可以推出复化科特斯公式的截断误差为2.3 高精度数值积分算法求积分时，复化积分公式采用逐步分段法，是一种比较有效的方法，但是它也存在许多的弊端，它收敛于积分真值的速度缓慢，从而人们在复化求积分公式上进行改进。在求积分时步长的大小 也会影响到积分的效果，步长太长积分值就不会太精确，步长太短则会增加许多的运算量。运算的时间也会随之

15、增加。在计算器中编程计算结果要花费大量的时间。以下采用变步长的计算方法，从而避免了这一点。2.3.1 龙贝格求积公式梯形法的算法简单，单精度低，收敛的速度缓慢。由此引出了龙贝格公式。由梯形的递推法可以看出，将积分区间等分时，用复化梯形公式计算的结果作为积分的近似值，其误差近似值为。可以设想，如果用这个误差作为的一种补偿，即将 作为积分的近似值，可望提高其精确度。直接根据复化求积公式，不难验证这说明，将区间对分前后两次复化梯形公式的值，作线性组合恰好等于复化辛浦生公式的值，它比更接近于近似值。同样，用于作线性组合会得到比更精确的值，且通过直接验证可得用与作线性组合，又可得到比更精确的值，通常记为

16、,即式就称为龙贝格求积公式。上述用若干个积分近似值推算出更为精确的积分近似值的方法，称为外推方法。我们将序列,和分别称为梯形序列、辛浦生序列、科特斯序列和龙贝格序列。由龙贝格序列当然还可以继续进行外推，得到新的求积序列。在积分区间逐次分半的过程中，利用外推法算式将粗糙近似值逐步加工成越来越精确的近似值。也就是说，将收敛速度缓慢的梯形序列逐步加工成收敛速度越来越快的序列。根据这种原理设计的计算积分近似值的方法称为龙贝格积分方法，又称为数值积分逐次分半加速收敛法。利用龙贝格序列求积的算法称为龙贝格算法。这种算法具有占用内存少、精确度高的优点。因此，成为实际中常用的求积方法。在优化技术方面的考虑龙贝

17、格方法比较合适的选择，在第四章会用程序进行讨论说明。第三章 线性方程组的求解上一章讲述了有关于优化技术在数值积分计算方法的应用，为了更加体现优化技术的应用本章将会讨论优化技术在线性方程组中的应用。实际中，存在大量的解线性方程组的问题。很多数值方法到最后也会涉及到线性方程组的求解问题：如样条插值的M和m关系式，曲线拟合的法方程，方程组的Newton迭代等问题。求解线性方程组有很多的方法，如gauss消去法，按比例主元消去法，用Cholesky分解解线性方程组，平方根法和追赶法等等。3.1线性方程组的介绍一般地设阶线性方程组为表示成矩阵形式其中为系数矩阵高斯消元法是按照消元和回代两个过程。高斯消元

18、法的改进为高斯主元消元法，并且主元消元法主要有列主元，按比例主元和全主元。高斯消元法的基本思想：首先将A化为上三角阵，再回代求解 得第二步 得 类似下去我们有 第步 n1步以后，我们可以得到变换后的矩阵为：从此再回代可以 解出线性方程组的解。但是高斯解线性方程组一般都是针对中小型的，一下介绍几种线性方程组的迭代法，从而求解线性方程组的近似解，利用优化技术判断哪个方法最优。3.2 线性方程组的迭代法3.2.1 Jacobi迭代法设阶线性方程组的系数非奇异（nonsigular），且。将方程组改为任取，将各分量代入上式的右边得将代入上式的右边，得以此类推，可得由此可得向量序列。称由迭代式建立的迭代

19、法为Jacobi迭代法。3.2.2 Gauss-Seidel迭代法在Jacobi迭代法中，每次迭代计算时用的是前一次迭代的全部分量。实际上，在计算分量时，最新的分量，已经算出，但没有被利用，而且，如果Jacobi迭代收敛，最新算出的分量一般比的精度更高。因此，可以对Jacobi迭代法加以改进，即在迭代过程中，每个分量计算出来之后，计算一下分量时就利用最新计算出的近似结果，具体地，即用新分量去替换右端的各项中，可得新的迭代公式此式称为Gauss-Seidel迭代法。3.2.3 SOR迭代法为了提高收敛速度，对Gauss-Seidel迭代法进一步用Gauss-Seidel迭代公式计算得到第个近似解

20、将前一步迭代值与Gauss-Seidel迭代值做加权平均，即其中是参数，整理得此式称为松弛迭代法，其中参数为松弛因子当时，式称为超松弛法；当时，式称为低松弛法；当时，式就是Gauss-Seidel。一般称这些方法为SOR方法。第四章 各种求积公式的MATLAB编程实现与应用4.1 对数值积分运行结果及其分析 不同的方法计算出来的积分求值的结果不同，所以一下针对进行求解，运用不同的方法得到近似值，余项以及运行程序所需要的时间，在精度不同的情况下，再进行分析判断出哪个方法最优。4.1.1 数值积分运行结果表 精度为的积分公式运行结果方法近似值时间梯形求积公式11.4685924646452720.

21、029394761280958辛普生求积公式11.4685924646452720.020265210920251科特斯求积公式11.4685924646452720.012255210320241复化梯形求积公式81.4641040335416940.141131524566874复化辛普生求积公式21.464109199415270.12400000000000复化科特斯求积公式21.464109199415270.12400000000000龙贝格积分公式21.4641062763466080.002670表 精度为的积分公式运行结果方法近似值时间梯形求积公式11.46859246464

22、52720.029394761280958辛普生求积公式11.4685924646452720.020265210920251科特斯求积公式11.4685924646452720.012255210320241复化梯形求积公式81.4641040335416940.148061234217113复化辛普生求积公式81.4641040335416940.143152625817205复化科特斯求积公式41.464101802284210.28100000000000龙贝格积分公式4表 精度为的积分公式运行结果方法近似值时间梯形求积公式11.4685924646452720.02939476128

23、0958辛普生求积公式11.4685924646452720.020265210920251科特斯求积公式11.4685924646452720.012255210320241复化梯形求积公式161.4641017695399430.277173394958087复化辛普生求积公式161.4641017695399430.283381996592803复化科特斯求积公式81.464101618606490.56100000000000龙贝格积分公式表 精度为的积分公式运行结果方法近似值时间梯形求积公式11.4685924646452720.029394761280958辛普生求积公式11.46

24、85924646452720.020265210920251科特斯求积公式11.4685924646452720.012255210320241复化梯形求积公式321.4641016248413570.567891012749313复化辛普生求积公式321.4641016248413570.563100795217624复化科特斯求积公式81.464101618606490.59300000000000龙贝格积分公式表 精度为的积分公式运行结果方法近似值时间梯形求积公式11.4685924646452720.029394761280958辛普生求积公式11.4685924646452720.0

25、20265210920251科特斯求积公式11.4685924646452720.012255210320241复化梯形求积公式641.4641016157450771.129131380037762复化辛普生求积公式321.4641016157450771.153032881910490复化科特斯求积公式161.464101615194791.17000000000000龙贝格积分公式4.1.1 数值积分运行结果分析根据matlab积分的命令的准确值为I=1.464101615137754。从运行的结果可以看出梯形积分公式和辛普生积分公式运行的结果都差不多，但是科特斯公式运行的结果较为准确，

26、但是运行的时间相比科特斯积分公式的时间比较长。相同的精度但是科特斯所分割的节点少。这是因为梯形公式、辛普生公式是低精度公式，但对被积函数的光滑性要求不高，他对对被积分光滑性较差的积分很有效。特别是梯形积分公式对被积分函数式周期函数积分时，效果更加突出。高阶科特斯求积分公式稳定性差，收敛较慢。从而，为了提高收敛速度建立的复化梯形积分公式，复化辛普生积分公式。但是相比之下，龙贝格积分公式是算法简单，是一个很好的加速方法。从优化技术的角度看，一般都是选择龙贝格积分方法。4.1 线性方程组运行结果及其分析4.2.1 线性方程组运行结果线性方程组表方程组进行运行结果迭代方法次数精确解近似解Jacobi迭

27、代法71.080303030303031.064935064935071.269047619047621.080042900000001.064748410000001.26874770000000Gauss-Seide迭代法51.080303030303031.064935064935071.269047619047621.080271541369921.064922680419551.26903884435789SOR方法51.080303030303031.064935064935071.269047619047621.080271541369921.064922680419551.26

28、9038844357894.2.2 线性方程组运行结果分析从运算的结果分析SOR运行的结果与精确值相比较为准确，而且迭代的次数少。在日常生活中，迭代法常用的优Jacobi迭代法,Gauss-Seidel迭代法，SOR方法中Jacobi迭代法简单，并具有很好的串行算法，很适合并行计算，但收敛速度较慢。Gauss-Seidel迭代法是典型的串行算法，在Jacobi迭代法与Gauss-Seidel迭代法同时收敛的条件下，后者比前者收敛的快，但两种迭代收敛收敛域互不相容，不能互相替代。SOR方法是一种应用极为广泛的方法，但选取最佳松弛因子比较困难，常通过试算来确定最佳松弛因子。各种方法都有其利弊，但是

29、在计算方程组时要求精确度很高时往往选用SOR方法。附录1、复化梯形积分公式的源程序function out1,out2,out3=fhtx(a,b,f) format long; clc; if nargin=3 wc=0.5*10(-6); end% 控制输入参数结束 disp(fhtx) tic; n=1; s1=(subs(f,x,a)+4*subs(f,x,(a+b)/2)+subs(f,x,b)*(b-a)/n/6; n=2; s2=s1+1; while abs(s2-s1)wc s1=s2; n=n*2; s=0; h=(b-a)/n; for i=1:n s=s+4*subs(

30、f,x,a+(i-1/2)*(h); end; for i=1:n-1 s=s+2*subs(f,x,a+i*h); end; s=(subs(f,x,a)+subs(f,x,b)+s)*h/6; s2=s; end; time=toc; jsz=s; % 控制输出参数结束 if nargout=1 out1=jsz; end if nargout=2 out1=jsz; out2=n; end if nargout=3 out1=jsz; out2=n; out3=time; end2、复化辛普生源程序function out1,out2,out3=xps(a,b,f) format lon

31、g; clc; if nargin=3 wc=0.5*10(-6); end%此处填复化辛普生公式 tic; n=1;s1=(subs(f,x,a)+4*subs(f,x,(a+b)/2)+subs(f,x,b)*(b-a)/n/6; n=2;s2=s1+1; while abs(s2-s1)wc s1=s2; n=n*2; s=0; h=(b-a)/n; for i=1:n s=s+4*subs(f,x,a+(i-1/2)*(h); end; for i=1:n-1 s=s+2*subs(f,x,a+i*h); end; s=(subs(f,x,a)+subs(f,x,b)+s)*h/6;

32、s2=s; end; jsz=s2; time=toc; % 控制输出参数开始 if nargout=1 out1=jsz; end if nargout=2 out1=jsz; out2=n; end if nargout=3 out1=jsz; out2=n; out3=time; end 3、复化牛顿科特斯源程序function jsz,n,time=kts(a,b,f) format long; clc; fd=zeros(1,5); wc=0.5*10(-3);%此处填复化科特斯公式 tic; wqja=a;wqjb=b; for i=1:5 fd(i)=wqja+(i-1)*(wq

33、jb-wqja)/4; end; s1=(7*subs(f,x,fd(1)+32*subs(f,x,fd(2)+12*subs(f,x,fd(3)+32*subs(f,x,fd(4)+7*subs(f,x,fd(5)*(wqjb-wqja)/90; n=2; h1=(b-a)/n; s2=0; for i=1:n wqja=a+(i-1)*h1;wqjb=a+i*h1; for i=1:5 fd(i)=wqja+(i-1)*(wqjb-wqja)/4; end s2=s2+(7*subs(f,x,fd(1)+32*subs(f,x,fd(2)+12*subs(f,x,fd(3)+32*subs

34、(f,x,fd(4)+7*subs(f,x,fd(5)*(wqjb-wqja)/90; end while abs(s2-s1)wc s1=s2; n=n*2;h1=(b-a)/n; s2=0; for i=1:n wqja=a+(i-1)*h1;wqjb=a+i*h1; for i=1:5 fd(i)=wqja+(i-1)*(wqjb-wqja)/4; end s2=s2+(7*subs(f,x,fd(1)+32*subs(f,x,fd(2)+12*subs(f,x,fd(3)+32*subs(f,x,fd(4)+7*subs(f,x,fd(5)*(wqjb-wqja)/90; end en

35、d jsz=s2; time=toc; end4、龙贝格积分公式的源程序function t=Romberg(fname,a,b,e)%Romberg法求函数的积分%fname是被积函数,a是积分上限,b是积分下限,e为精度ticif nargine i=i+1;h=h/2; T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:b-h/2+0.001*h)*h/2; for j=1:i T(i+1,j+1)=4j*T(i+1,j)/(4j-1)-T(i,j)/(4j-1); endend%Tt=T(i+1,j+1);toc5、jacobi迭代法求方程组%P:范数

36、的名称，P=1,2，inf%error：近似值x的误差%maxl迭代的最大次数function Jacobi(A,b,X0,P,error,maxl)n n=size(A);X=zeros(n,1);for k=1:maxl for j=1:n X(j)=(b(j)-A(j,1:j-1,j+1:n)*X0(1:j-1,j+1:n)/A(j,j); end X errX=norm(X-X0,P); X0=X;X1=Ab; if(errX=error) disp(请注意Jacobi迭代次数已经超过最大迭代次数maxl )end6、GS迭代法求方程组%P:范数的名称，P=1,2，inf%error：

37、近似值x的误差%maxl迭代的最大次数function GS(A,b,X0,P,error,maxl)n n=size(A);X=zeros(n,1);for k=1:maxl for j=1:n; XX=0; for i=1:n if ij XX=XX+A(j,i)*X0(i); end end X(j)=(b(j)-XX)/A(j,j); end X errX=norm(X-X0,p); X0=X;X1=Ab; if(errX=error)disp(请注意Jacobi迭代次数已经超过最大迭代次数maxl) end7、SOR迭代法求方程组 %P:范数的名称，P=1,2，inf%error：近

38、似值x的误差%maxl迭代的最大次数%w：松弛因子function SOR(A,b,X0,P,error,maxl)n=size(A);X=zeros(n,1);for k=1:maxl for j=1:n; XX=0; for i=1:n if ij XX=XX+A(j,i)*X0(i); end end X(j)=(1-w)*X0(j)+w*(b(j)-XX)/A(j,j); end errX=norm(X-X0,P); X0=X;X1=Ab; if(errX=error)disp(请注意Jacobi迭代次数已经超过最大迭代次数maxl ) end 参考文献1朱晓临.数值分析M，中国科技大

39、学出版社，20102杜延松，沈艳君.数值分析及实验M，科学出版社，2006.3王沫然. Matlab6.0功能科学计算M，北京，电子工业出版社，2001，166-191.4薛莲,将金生. 数值计算方法M，电子工业出版社，2007.5楼顺天,于卫,闫华梁. Matlab程序设计语言M，西安，西安电子科技大学出版社，1997.6薛定宇,陈阳泉.高等应用数学问题的Matlab求解M,北京，清华大学出版社,2004.7李庆扬,王能超,易大义.数值分析M,北京，清华大学出版社,1988.8肖筱南,赵来军,党林立. 现代数值计算方法M, 北京，北京大学出版社2003.9封建湖,车刚明,聂玉峰.数值分析原理M. 北京，科学出版社2001.10孙海东. EXCEL在数值计算方法方面的应用J,软件世界,1999,No.11,137.致 谢毕业论文致谢词模版本论文是在导老师的悉心指导下完成的。导师渊博的专业知识，严谨的治学态度，精益求精的工作作风，诲人不倦