非线性Volterra积分方程

1、一类第二种非线性Volterra积分方程积分数值解方法1前言微分方程和积分方程都是描述物理问题的重要数学工具，各有优点.相对于 某种情况来说，对于某种物理数学问题，积分方程对于问题的解决比微分方程更 加有优势，使对问题的研究更加趋于简单化，在数学上，利用积分形式讨论存在 性、唯一性往往比较方便，结果也比较完美，所以研究积分方程便得越来越有用， 日益受到重视.积分方程的发展，始终是与数学物理问题的研究息息相关一般认为，从积分 发展的源头可以追溯到国外的数学家克莱茵的著作古今数学思想，该书是被 认为第一个清醒的认为应用积分方程求解的是Abel. Abel分别于1833年和1826 年发表了两篇有关

2、积分方程的文章，但其正式的名称却是由数学家du Bois-Raymond忤次提出的，把该问题的研究正式命名为积分方程。所以最早研 究积分方程的是Abel,他在1823年从力学问题时首先引出了积分方程，并用两 种方法求出了它的解，第一的积分方程便是以Abel命名的方程.该方程的形式 为：上丄力=f(x),该方程称为广义Abel方程，式中a的值在(0, 1)之间.当(x_/)沪时，该式子便成为-g/r = /(x).在此之前，Laplace于1782年所提出 的求Laplace反变换问题，当时这个问题就要求解一个积分方程.但是Fourier 其实已经求出了一类积分方程的反变换，这就说明在早些时候积

3、分方程就已经在 专业性很针对的情况下得到了研究，实际上也说明了 Fourier在研究反变换问题 是就相当于解出了一类积分方程.积分方程的形成基础是有两位数学家 Fredholm和Volterra奠定的，积分方程主要是研究两类相关的方程，山于这两 位数学家的突出贡献，所以这两个方程被命名为Fredholm方程和Volterra方程。 后来乂有德国数学家D. Hilbert进行了重要的研究，并作出了突出的贡献，山于 D. Hilbert领头科学家的研究，所以掀起了一阵研究积分方程的热潮，并出现了 很多重要的成果，后来该理论乂推广到非线性部分。我国在60年代前，积分方 程这部分的理论介绍和相关书本主

4、要靠翻译苏联的相关书籍，那时研究的积分方 程基本是一种模式，即用古典的方法来研究相关的积分方程问题，这样使得问题 的研究变得繁琐、复杂，在内容方面比较单一、狭隘，其至有些理论故意把积分 方程的研究趋向于复杂化。随着数学研究的高速发展，特别是积分方程近年来的 丰富发展，如此单一、刻板的解法已经不能跟上数学研究时代的步伐。在九十年 代我国的数学专家路见可、钟寿国出版了积分方程论，该书选择Z?空间来讨 论古典积分方程，并结合泛函分析的算子理论来分析积分方程的相关问题。最近 出版的比较适合一般读者阅览的积分方程的书有李星岀版的积分方程，该书 从最简单的方法分析研究积分方程的理论问题，并给日后打算研究泛

5、函的读者提 供了基本的实例。山于现代的计算机技术高速发展，对于一些比较复杂，难以求 解的非线性积分方程逐渐采用了比较有效的数值解，常用的方法有逐次逼近法、 Adomian分解法、配置法、haar小波方法、小波-Galerkin方法、泰勒展开等等 一些方法.请预览后下我!现在积分方程的应用广泛，很多问题都可以引出积分方程，并可用积分方程 来解决像弹性弦问题、线性系统响应问题、人口增长问题、等时曲线问题等等. 还有在空气动力学中研究分子运动，对于非均匀流体中悬浮晶粒的布朗位移，导 致了以柯尔莫哥洛夫命名的一类非线性积分方程在确定飞机机翼方面的研究 中，对于气流、升力等问题的计算中也引岀了积分方程的

6、研究现在很多积分方 程方面的研究都取得了不错的进展.2、预备知识积分方程是一个在积分号下出现待求的函数的方程，称作积分方程。含一个未 知函数的线性积分方程的一般形式为：(x)-2p(xO)(yXv = /(x) xea,b 1 我们把积分号中的上下限为常数的积分方程称为Fredholm方程。其中/(小、k(xyy)是已知函数,(刃是未知所要求的函数。般称/(x)为自由项,k(xyy)称为积分方程的核，而兄是积分方程的一个参数，方程的解与几相关 而与Fredholm对应的是Volterra积分方程，与Fredholm的一个不同点是 Volterra积分号下的上下限中的上限是一个变量而不是一个常数

7、。本论文我们 主要研究的是非线性的Volterra积分方程。非线性Volterra积分方程的形式 为卩(x)-/l(k(x,y)F(0(y) dy = /(x)。同样，仅y)为所求。积分方程的解法：一、Fredholm积分方程一般的解法有：有限差分逼近法、逐次逼近法及解核、 泛函修正平均法、Fredholm积分方程退化核解法、退化核近似代替法、待定 系数法。二、八Volterra积分方程的常用解法：有限差分逼近法、逐次逼近法、转化为 常微分方程的初值问题、第二类Volterra积分方程的数值积分解法。第二类线性volterra积分方程与第二类线性Fredholm积分方程的一个很大的差 别是vo

8、lterra方程的解不依赖于参数兄的值，即对于任何的参数几该方程都有 解，而且解具有唯一性。具体的证明过程参见,所以得出定理：如果核k(x,y)与自由函数/(兀)在Volterra积分方程的定域是连续函数，那么无论参数几取何 值该方程都有解，而且解具有唯一性。3、Volterra积分方程3. 1第一类Volterra积分方程3.1.1第一类线性Volterra积分方程形如：卩(刃小= /(x)其中函数k(x,y)、/(x)为已知的，俠x)为所要求的未知函数，这样的方程叫第 一类线性Volterra积分方程.一般来说第一类积分方程由于其不适定性，研究其 解跟第二类积分方程有很大的不同，而且比较复

9、杂，在此主要简单介绍一下.3. 1. 2第一类Volteraa积分方程的一种解法在某种情况下第一类Volterra积分方程通常可以化为第二类volterra积分方 程的解，一般对方程两边求导，当方程的R(x,y)、/(X)可微，且k(x, y) 0 ,就把第一类该方程化为第二类。化为第二类的形式为：*)+1需(恥=需，这样就可以用第二类积分方程的解法来求解.对于一种特殊的第一类Volterra积分方程：Abel方程，Abel方程是Volterra 积分方程的一种特殊情况，其形式为：0(刃d)-dy = f(x)其中当x = y时，该方程出现弱奇性。其解可根据定理：假设Abel积分方程”的自由项

10、/(兀)是连续可微的，而且/(6/) = 0,则它有唯一的(兀_刃at, 口“、 sin 加 df(y) f解即恥)=丁-川穴皿3. 2第二类vol terra积分方程3. 2. 1第二类线性Volterra积分方程第二类方程的形式:卩(x) -k(xyy)(p(y)dy = f(x)其中仅x)是所要求的未知函数，2是已知或是需要讨论的参数，跟Fredholm方程_样k(x.y)是已知的函数，叫Vo It err a方程的核，当k(x.y) =0时，Volterra 方程可以看成特殊形式的Fredholm方程，而且Fredholm方程理论适合于 Volterra方程。Vo It err a方程

11、有自己的特点，例如，Volterr方程没人特征值, 对于任意的自山项它都有解。对于笫一类的Voltrra方程在某种条件下可以转化 为第:类Volteira方程。与第二类线性Fredholm积分方程一样，第二类线性Volterra积分方程也有自己的迭核、解核，其迭核、解核的引出方法跟第二类 线性 Fredholm 一样.请预览后下我!迭核：假设 心兀,y)=(儿，则02(尤)=*2(儿刃/山 (pn W = k”(x,y)/(y)dy ,我们称k(x,y)为 Volterra 方程的 迭核。解核：我们称R(x,y;2)=乞右心匕刃为解核H-1只要知道方程的迭核，就能求得方程的解核，从而求得方程的

12、解。3. 2. 2第二类线性Volerra积分方程的解法：1、逐次逼近法假设方程具有这样一个形式的解fl0(X)= (x) +(p (x)A + + %才=工 0(X)兄r-O如果对于逐次方法来说该方程有解，解次方程一般令% = fw(P = fM + kydy(Pi =/(x) + k(x,y)0yPn =/(x) + k(x,y)%_dy那么，对于上述的级数一定收敛，即对级数jy收敛，可以证明对于任意i-0的参数几方程都有解，依据定理：如果核k(x,y)及自山项/(x)是连续的实函数.那么第二类线性Volterra方程卩 一k(xy y)(p(y)dy = f(x)对于任意的参数兄存在一个

13、唯一的连续解，而且解可以用逐次逼进法求出.迭核：假设込(儿刃=/)“, yW，则(Pi W = k2 (x,df (x) = k(x,y)f(y)dy ,我们称k(x,y)为 Volterra 方程的 迭核。解核：我们称/?(兀;刃=工矿比(圮刃为解核，只要知道方程的迭核，就能 n-1求得方程的解核，从而求得方程的解。3. 2. 3非线性第二类Volterra积分方程非线性第二类volterra积分方程的形式如：卩(x)-几y)F(y) dy = /(x)未知函数为0(x),而/(x)、R(x,刃、F(x)都是已知的。当方程满足一定条件时,可用逐 次逼近法求解。对于第一类非线性Volterra

14、积分方程可以通过转化成第二类非线性 Volterra积分方程求解。具体转化的过程参见。非线性Volterra积分方程的数值解3. 3卷积型Volterra方程的解法3. 3.1第二类卷积型Volterra积分方程的解F(p)1 K(p)1、形如(p(x) = f(x) + Zc(x - y)(p(y)dy称为第_.类卷积型Volterra积分方程， 此类方程一般用Laplace变换来解决。如果方程中的/(x)、y)是足够光滑、 指数阶的函数，那么方程的解也是指数阶的，这样就可以用Laplace变换来解此 方程。设0伙(x)=K(p),莎/(x)=F(p), q(x) = 0(p),通过对方程两

15、边作Laplace 变换，可得(p(p) = F(p) + K(p)(pp),解出 0(P)=当 K(p) H 1 时，0(x) = 0-”-(/?)l-K(p)2、对于第一类Volterra积分方程，即方程 k(x-= f (a)同样对方程两边作Laplace变换，可解得p(p)= 也，所以方程的解为p(x)= K(p)3、非线性Volterra积分方程Laplace变换的解.即方程 p(x) = f(x) +(py(px - y)dy ,设 0(x)=0(p), nrf(x)+F(p),对方程两 边作Laplace变换，可得出卩(切=屛_芒；4)当! 土卍丁“屮)2A2Z存在时，该解就是非

16、线性Vo It err a积分方程的Laplace变换得岀的解.4、积分数值解相关知识4. 1 Newton-Cotes型积分求积公式在这里我们主要讨论的数值计算问题，可以假设/(X)在，闰上可积.在一些时候函数/(X)并不是可积的能用初等函数来表示，所以有的时候并不 能求出该函数的原函数，因此，这里我们来研究用数值方法解函数积分.欲计算积分/(/) = J f(x)W(x)dx ,其中W(x)为权函数，可以假设/(X)在n+1个互异的点：axxx2 - x+l b 的值分别为：/(xj , f(x2),/(G),就可以用/(xj, /(心)，/(畑)的线性组合得出积分的近似解,n+1即/”(

17、/)(/),其中 /”=工人/3)r-1n+1插值求积公式：1(f )=工人/(兀)+ En(f )&(/)是离散误差(兀一兀)心(石)M = (x 一 册)(x X?)(x X”)当我们假设肚列为有限区间，W(x) = l,并将该区间分成n等份，取等距基点为：“=召，4.=.如果令n=2时就可以得到simpson公式.2 - 24.2复合梯形公式我们假设所讨论的积分中函数的定域义为“，在该区间取n+l个互异基点，即 = X = b ,且取步长为h = xi+i -xt =-.在子区间斗+|,兀上H使用体型公式，所以nh “h 91f f(Qdx = - f(x)dx=-f(xi) +f )从

18、而可以得出：f fg = f(a) + fb) + 2/( + 询-炷厂)a/ Lf=ii / /=i舍去JLjj项，于是就得到复合梯形公式：12 /-I/(x)t/x = f(d) + /(b) + 2/(d + 7)5积分方程的数值解方法5.1未知函数展开法在这里，我们将讨论在L2(a,h)中完备的函数系在近似方程解方面的作用， 这些函数系可以是正交的，也可以的非正交的，可以去某个函数系的有限项 当作方程解的近似值.设该函数系为0二，其中该函数系的各个函数是线性 无关的，可令积分方程的解(pg 土臥,把该近似函数代人积分方程： r-1(p(x - k(x, y)(p(y)dy = f(x)

19、, 这 样 就 可 得 到 c 几 q “(x, ydy + f(x),把它整理成：r-1r-1左妙-/ifcj 心y)dy- f(x) = R(x),其中R(x)是残差，如果能是R(x)等 r-1i-1于零，那么方程的近似解就等于该方程的精确解，但是一般来说，要使R(x) 等于零是很难的，一般在R(x)很小的惜况下可以忽略，即得出方程的一种数 值解：C0 =几土cj+，但对残差的不同要求，可以得出r-1r-1不同的解法，一般来说有如下儿种解法：配点法、Galerkin法、最小二乘法 等方法.1、配点法如果要求残差在所选取的基点上满足心石)等于零，其中忑爲是一些互异的 基点，如此便可以得到一下

20、方程组：f 50(无)一几f cf k(无,y)0(yMy = /(无)，(k = ,2,- -ji ),求解该方程组 F-lF-1变可以得到展开系数匕打2、矩量法对于矩量法以下用于Fredholm积分方程，即p(x) -kg y(p(y)dy = /(a).矩量法就是要求残差关于原点到阶的矩为零，即/?(x)xA- = O,可得到如 下的方程组：立0(x)xWy 兄fk(x,y)0(y)心 xk dx = f(x)xk dx,r-1j-I1(A = 1,2,)解此方程变可以得到展开系数匕；二3、Galerkin 法Galerkin法要求残差函数R(x)在平方可积空间即空间la,b与函数0内积

21、 为零，即要求8(劝心=0, 21.2，八所以展开系数可以这样来确定 匕打，取函数系中的前n个函数如(山1.2在d,b上与积分方程 左妙(x)-述qfgy)% = /(忑)两端正交,令0 =c妙于是展开系数满足下列线性方程组:(p(x),(x)=ydy,(x)j + (/(x),.(x)6/, i = 1.2,请预览后下我!其中（/a）,g（x）=f/a）g（x）厶，解该方程组便可以得到展开系数匕打5. 2积分核级数展开法积分核级数展开法乂称退化核近似方法，就是利用某种展开方式把非退化的 核展开成近似退化的核，一般的展开方式有泰勒级数展开、Fourier级数展开、 Ia,b空间内的线性无关的针

22、对为知函数近似展开的函数系等等.如果利用泰勒 展式，那么应该注意保留合适的级数项数，一般来说应该根据积分限的大小来决 定级数项数.也有把未知函数展开求积分方程的未知函数的解.对于用退化核来近似积分方程核的误差有如下估计方法：定理1：设k（x.y）是积分方程核的近似退化核，对于退化核满足条件神*y）-k（x9 y）dth而且以退化核y）为积分核的积分方程的解核Rgy:A）,成立 R（x,y;A）t R则积分方程0（x） = /lk（x,y）0（y）心的解0（X）与用近似退化核代替的积分方程的解0（圮刃，满足/ BA(l + AR)2h 1-|2|/(1 + |2|7?)在式子中,B是f（x）的一

23、个上界.6非线性Volterra积分方程的数值解0（x） +可（卩（刃）心= /（x） xe ,/?,我们假定 /（x）、k（x,y）、F（x）在其定义域上都是连续函数，利用数值求值公式札-怎 A.nK.m F（0,”）= fi心 1,2,.，仏m-1ew是数值积分公式中的权系数，Kim =k（Xi,xm）, m 该方程组是一个n阶的下三角方程组，（p（/） = /（）,有此我们可以顺着方程组 的顶端解出n个数值解，所以我们便得出方程的近视解（P （x） = 2Amk（xi,）F（pm） + f（xi） i = 12儿 wr-l当n趋向于无穷时该解也是趋于精确解6. 1梯形公式我们取h步长，力=迅匸上（n为大于1的正整数），心是x的终点，由梯H形公式儿=，i =a + h,，儿=叫f心y）F（0（y）d）y h壬儿比）弘仇）+心儿）恥（打）+心汕）弘（加）+（儿片了甌片） 俠兀）-h （兀儿）Fe（儿）*（儿 y JFb（x）+心儿-JF0（y”-i）+” （兀儿）He（儿） = /（兀） 心1