一维抛物方程的两层网格算法优越性研究_李小珍.pdf
《一维抛物方程的两层网格算法优越性研究_李小珍.pdf》由会员分享,可在线阅读,更多相关《一维抛物方程的两层网格算法优越性研究_李小珍.pdf(5页珍藏版)》请在咨信网上搜索。
1、第 38 卷哈尔滨师范大学自然科学学报Vol 38,No 6 2022第 6 期NATUAL SCIENCE JOUNAL OF HABIN NOMAL UNIVESITY一维抛物方程的两层网格算法优越性研究*李小珍(安徽国防科技职业学院)【摘要】为深入研究一维抛物方程的求解问题,利用 Matlab 软件使用 Jacobi迭代法、Gauss Seidel 迭代法和两层网格算法对一维热传导抛物方程进行求解,分析不同算法的数值解与精确解 研究发现,在取得相同迭代误差的条件下,两层网格算法所需的迭代次数最少,且大大节省了运算时间,提高了运算的效率【关键词】抛物方程;网格算法;优越性;效率【中图分类号
2、】O241【文献标识码】A【文章编号】1000 5617(2022)06 0018 04收稿日期:2022 05 22*基金项目:安徽省教学研究项目“课程思政视域下 高等数学 教学探索与实践”(2021jyxm0232)0引言抛物型方程属于一种偏微分方程,这一类方程往往与人类的生产生活密切相关,可以通过求解偏微分方程来解释一些自然现象和物理意义求解偏微分方程的数值解,目前来说,主要归结为直接法和迭代法这两种常见的方法1 其中,直接法是通过有限差分方式把微分方程中的微分形式用差商来代替,以此将一个抛物型偏微分方程转化为便于研究的代数方程或方程组 但这种方法也有着诸多不便 相较于迭代法来说,它占用
3、内存大,运算速度慢2 通过讨论两种较常见的迭代法,即雅可比迭代法、高斯 赛德尔迭代法剖析网格算法的优势 这两种迭代法均可收敛,但通过计算比较得出,在相同精度下,高斯 赛德尔迭代法所需迭代的次数更少,可节省更多时间34 然而,高斯 赛德尔迭代法的计算量依然存在计算量较大和收敛速度慢的问题5 因此,利用两重网格算法进行抛物方程迭代求解6 两层网格方法的主要思想就是在细网格上迭代,再在粗网格上校正,如果没有达到所需误差精度便一直重复上述两步骤 这种方法运算速度快,占用内存少,一般来说,仅需很少的迭代次数便可达到所需的精度1抛物方程及算法概述1 1抛物型偏微分方程抛物型偏微分方程是偏微分方程的一种 通
4、过判别式的符号来对偏微分方程进行分类,在微分方程1 ut K2ux2=f(x,t)(1)在公式(1)中,x 表示在一维平面中的位置,t 表示时间 当 f(x,t)=0 的情况,即抛物型偏微分方程 一维抛物型微分方程中的典型例子是一维的热传导方程(在此处与扩散方程表达式相同,因为热扩散的方式是一样的)基本形式可表示为5:ut K2ux2=0,0 x 1(2)有边值条件为u(0,t)=u(1,t)=0,并有初值条件u(x,0)=u0(x)一维的热传导方程可模拟杆第 6 期一维抛物方程的两层网格算法优越性研究中的热传导,在这种情况时,u(x,t)表示热流的密度,它与温度成正比;K 表示热传导系数,它
5、取决于材料和它自身导热的程度 为了后文研究和计算方便,(2)式中取热传导系数 K=1,即化为ut2ux2=0,0 x 1(3)1 2迭代法的基本概念(1)迭代法的基本概念对于任意一个给定的线性方程组 Ax=b,其中矩阵 A=(aij)nn,b=(b1,bn)T 假设它有唯一解,可以通过迭代产生序列,并构造以下的向量序列x(k+1)=Mx(k)+g(4)其中,k 表示迭代次数(k=0,1,2,n),M 为n 阶 矩阵,对于给定的线性方程组,用公式(4)逐步代入求近似解的方法称为迭代法5(2)雅可比迭代法3 雅可比迭代法是以普鲁士的著名数学家雅可比(Carl Gustav Jacobi)来命名的
6、假设一个线性方程组为Ax=b,采用Jacobi的方法,则要将A分解为三部分,分别为下三角矩阵 L,对角阵 D,上三角矩阵 U,A=L+U+D,如式(5)所示:也就是说,b=(L+D+U)x,移项得:Dx=(L+U)x+b,解得:x=(D)1(L+U)+(D)1b 对于公式(4),经过如下一系列迭代:(6)得到迭代公式:x(k+1)i=1aij(binj=1aijxkj),i=1,2,n(3)高斯 赛德尔迭代法高斯 赛德尔迭代法是以德国数学家卡尔弗里德里希高斯和路德维希赛德尔来命名的,也称为李伯曼方法或连续位移方法4 先假设一个线性方程组为 Ax=b,同样将系数矩阵 A分解为三部分,对角阵 D,
7、下三角矩阵 L,上三角矩阵 U 不同的是,A=D L U,故高斯 赛德尔迭代法的迭代公式为x(k+1)i=(bii1j=1aijx(k+1)jnj=i+1aijxkj)/aii,(i=1,2,n;k=0,1,2,n)(7)(4)两层网格算法两层网格方法主要思想是通过细网格迭代来减少误差中的高频部分,粗网格校正来减少误差的低频部分6 通过这种方法,可以满足“收敛速度与网格步长 h 无关”这种要求抛物方程用有限差分法将其化作如下的矩阵形式:Au=f,然后用两层网格的方法进一步求解:第一步:使用 Jacobi 迭代(也可用其他迭代如 Gauss Seidel 迭代)迭代 V1次(V1可以自由选取合适
8、值,此处 V1选取 3),记得结果为 un第二步:计算余量 rh,通过等式 rh=fnAn un 并做余量限制 r2h=h*rh第三步:解粗网格方程:A2hE2h=r2h 得到 E2h,其中 A2h 可由之前定义的差值矩阵和约束矩阵得到,为:A2h=hAhIh;第四步:做粗网格校正,并将校正量加到 uh上,即此时的近似解u=uh+Eh 校正量Eh可通过粗网格计算 Eh=IhE2h 得到;第五步:计算此时的误差精度是否达到要求,如达到,则已完成所有步骤,上一步骤中的 u即为所需要的解;如没有达到,则返回第一步中的细网格迭代,迭代方程中的 un用最新得到的 u来代替 如此循环几次,直到达到接近真实
9、值为止91哈尔滨师范大学自然科学学报2022 年 第 38 卷2结果与讨论2 1不同算法的数值解和精确解比较利 用 Matlab 软 件 使 用 Jacobi 迭 代 法、Gauss Seidel迭代法和两层网格算法对一维热传导抛物方程进行求解 当 n=10,k=1 时的数值解和精确解结果如图 1 所示图 1一维热传导抛物方程不同算法的数值解与精确解由图 1 可知,Jacobi 迭代法、Gauss Seidel迭代法和两层网格法的数值解与精确解走势图变化趋势非常接近,未出现多峰等现象,说明三者计算结果的准确度相差不大2 2Jacobi 和 Gauss Seidel 迭代法中 n 与 k 的关系
10、(1)当 k=1 时的迭代次数为剖析Jacobi迭代法和GaussSeidel迭代法中 n 值和 k 值对迭代次数的影响 先固定 k=1,取不同的 n 值,如 10,20,40,80,160,320,取误差值精度为 0 01,比较此时 Jacobi 迭代的次数由 Matlab 软件计算,结果见表 1表 1k=1 时不同算法的迭代次数结果算法n 值10204080160320Jacobi 迭代法6042316898335232139372554206Gauss Seidel 迭代法308116544981762369692277110(2)当 n=160 时的迭代次数然后固定 n=160,取不同
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一维抛物 方程 网格 算法 优越性 研究 李小珍
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。