好文档 - 专业文书写作范文服务资料分享网站

二维积分微分方程初边值问题的Taylor配置解和误差分析

天下 分享 时间: 加入收藏 我要投稿 点赞

二维积分微分方程初边值问题的Taylor配置解和误差分析

赖嘉导,王奇生

【摘 要】摘要:利用Taylor配置方法求解二维Volterra-Fredholm型积分微分方程初边值问题,给出了Taylor配置解的求解格式和误差分析结果,并给出了阐述理论分析结果的数值例子.

【期刊名称】五邑大学学报(自然科学版) 【年(卷),期】2013(000)002 【总页数】9

【关键词】二维积分微分方程;初边值问题;Taylor配置解;误差分析 考虑二维Volterra-Fredholm型积分微分方程问题

其中,为上已知的适当光滑函数,为上的已知核函数,为未知函数,为常数,则微分算子 .

初始条件为: 边界条件为: . (3)

积分微分方程(1)满足初始条件(2),则为积分微分方程初始问题;积分微分方程(1)满足边界条件(3),则为积分微分方程边值问题. 选取配置点,,称

为积分微分方程问题的阶Taylor配置解,其中为待定系数. 本文给出了二维Volterra-Fredholm型积分微分方程初边值问题Taylor配置解的求解格式和误差分析结果.

1 方法介绍

令,,则方程(1)可写成: 对式(5)的两边取阶偏导,得: ,. (6)

1.1 微分部分的矩阵表示 矩阵表示,即令,其中,,, 1.2 积分部分的矩阵表示 在式(6)中, . (8)

用分别替代式(7)、(8)中的,取配置点,则 . (10)

把式(9)、(10)代入式(6),得: 其中, ,

,配置方程组(11)是一个阶线性方程组,其未知量为. 式(11)的矩阵表示如下: . (12) 其中,,,,,. 1.3 初始条件的处理

将式(4)及其导数表示成矩阵形式: ,

.

则初始条件(2)可表示成: , ,

选取个初始点,记为,其中. 将这些点代入上面方程,对应可得到个矩阵方程: ,,,,

即可写成如下形式: (13)

其中系数是常数,与选取的初始点有关. 1.4 边界条件的处理

边界条件(3)分为4种情况,分别表示如下: , ,

选取个初始点,记为,,其中. 将这些点代入上面方程,对应可得到个矩阵方程: ,,,,

即可写成如下形式: (14) 1.5 求解方法 由式(12)得:

令,则,其中,,,. 矩阵方程可表示成: , (16)

由方程(12)确定.

现在考虑增广矩阵(16),用式(13)或式(14)替代式(16)的最后行,得到新的增广矩阵: 或

若秩,则二维Volterra-Fredholm型积分微分方程(1)在初边值条件下有唯一解,此解可表示成: 为阶Taylor配置解.

2 误差分析

对二维积分微分方程问题(1),定义误差函数: 把代入式(1)得: , (19)

其中为扰动项,由式(19)可知:

用式(1)减去式(19),且由式(18)可得: . (21)

构建误差函数的近似函数必须根据阶数的不同,重复计算等式(20),求出再代入式(21),利用与求解式(1)相同的方法可得到误差函数的近似解. 令,则由误差函数的定义得到:

现在的目标是使成立,其中是任意正整数. 令,则随着阶数的不断增大,下列不等式对任意点都成立: .

即随着阶数的增大,误差函数接近于0.

3 数值例子

例1 求解二维积分微分方程初值问题 ,.

初始条件:.

解 取,. 按方程(15)的解法得:

这样,4阶Taylor配置解,恰为原方程的精确解. 例2 求解二维积分微分方程初值问题 初始条件:,这个方程的精确解是. 解 取,. 按方程(15)的解法得: 其中,代入式(4),得到2阶配置解: .

取,按方程(15)的解法得:

其中,代入式(4),得3阶配置解为: .

取,按方程(15)的解法得:

其中,代入式(4),则4阶配置解为:

表1分别给出了时,精确解、数值解及误差函数在相应配置点处的计算结果,同时在图1中给出了误差函数、的图像. 从表1结果可以看出,随着阶数的增加,得到的数值解误差可以任意小. 例3 求解二维积分微分方程边值问题 边界条件:.

二维积分微分方程初边值问题的Taylor配置解和误差分析

二维积分微分方程初边值问题的Taylor配置解和误差分析赖嘉导,王奇生【摘要】摘要:利用Taylor配置方法求解二维Volterra-Fredholm型积分微分方程初边值问题,给出了Taylor配置解的求解格式和误差分析结果,并给出了阐述理论分析结果的数值例子.【期刊名称】五邑大学学报(自然科学版)【年(卷),期】2013(0
推荐度:
点击下载文档文档为doc格式
1jup14icau670et7c26i4qfr01784a016lg
领取福利

微信扫码领取福利

微信扫码分享