抛物型最优控制问题的全离散Crank-Nicolson有限元法.pdf
《抛物型最优控制问题的全离散Crank-Nicolson有限元法.pdf》由会员分享,可在线阅读,更多相关《抛物型最优控制问题的全离散Crank-Nicolson有限元法.pdf(11页珍藏版)》请在咨信网上搜索。
1、第41卷第2期2024年3月新疆大学学报(自然科学版中英文)Journal of Xinjiang University(Natural Science Edition in Chinese and English)Vol.41,No.2Mar.,2024抛物型最优控制问题的全离散Crank-Nicolson有限元法张馨丹1,赵建平1,侯延仁1,2(1.新疆大学 数学与系统科学学院,新疆 乌鲁木齐 830017;2.西安交通大学 数学与统计学院,陕西 西安 710049)摘要:针对具有积分控制约束的抛物型最优控制问题,提出了一种基于 Crank-Nicolson 格式的全离散有限元法.使用分段
2、线性有限元对状态进行空间离散,采用 Crank-Nicolson 格式进行时间离散,对控制变量采用分段线性近似,从而得到离散的最优性系统.证明了状态变量、伴随状态变量和控制变量的误差估计,并通过数值算例验证理论结果.关键词:抛物型最优控制问题;Crank-Nicolson 格式;最优性系统DOI:10.13568/ki.651094.651316.2023.05.18.0001中图分类号:O241.82文献标识码:A文章编号:2096-7675(2024)02-0196-010引文格式:张馨丹,赵建平,侯延仁 抛物型最优控制问题的全离散 Crank-Nicolson 有限元法J 新疆大学学报(
3、自然科学版中英文),2024,41(2):196-205+245英文引文格式:ZHANG Xindan,ZHAO Jianping,HOU Yanren Fully discrete finite element method based on theCrank-Nicolson scheme for parabolic optimal control problemJ Journal of Xinjiang University(Natural ScienceEdition in Chinese and English),2024,41(2):196-205+245Fully Discret
4、e Finite Element Method Based on the Crank-NicolsonScheme for Parabolic Optimal Control ProblemZHANG Xindan1,ZHAO Jianping1,HOU Yanren1,2(1.School of Mathematics and System Sciences,Xinjiang University,Urumqi Xinjiang 830017,China;2.School of Mathematics and Statistics,Xian Jiaotong University,Xian
5、Shaanxi 710049,China)Abstract:Aiming at parabolic optimal control problem with integral control constraints,a fully discrete finiteelement method based on Crank-Nicolson scheme is proposed.The state is discretized by piecewise linear finiteelements for the space discretization,Crank-Nicolson scheme
6、for time discretization,and the control variables areapproximated by piecewise linear function approximation,so as to obtain a discrete optimal system.The errorestimates of state variables,adjoint state variables and control variables are proved,and the theoretical results areverified by numerical e
7、xamples.Key words:parabolic optimal control problem;Crank-Nicolson scheme;optimality system0引 言偏微分方程最优控制问题在实际工程领域有着广泛的应用1,例如,大气污染控制、癌症治疗过程中的热处理、低温超导激光能量的爆破、数值天气预测中的资料同化、飞机机翼阻力优化、最优形状设计、石油开采过程优化、流体控制和混凝土大坝的最优温度控制等,都涉及求解偏微分方程描述的最优控制问题.目前对偏微分方程最优控制问题的数值求解思路主要有两种2:第一种是先优化再离散,即先利用Lagrange收稿日期:2023-05-18基
8、金项目:国家自然科学基金“沙丘长时间移动演变行为的动力学模型及其数值模拟研究”(61962056);新疆维吾尔自治区自然科学基金“带界面的抛物型最优控制问题高效数值算法研究”(2022D01C409)作者简介:张馨丹(1996),女,硕士生,从事偏微分方程的数值解法研究,E-mail:通讯作者:赵建平(1981),女,博士,教授,从事偏微分方程的数值解法研究,E-mail:第2期张馨丹,等:抛物型最优控制问题的全离散Crank-Nicolson有限元法197乘子方法推导出最优控制问题的最优性条件(状态方程、伴随方程、变分不等式),再利用数值方法进行离散,使之变成有限维的数值计算问题;第二种是先
9、离散再优化,即先利用数值方法对最优控制问题进行离散,得到有限维的优化问题,然后再利用优化算法进行求解.不论选择哪种方法,数值离散都是不可缺少的,因此好的离散方法对于求出最优解是至关重要的.对于具有积分约束的椭圆型最优控制问题已有大量研究36,但对具有积分约束的抛物型最优控制问题,用Crank-Nicolson有限元方法进行数值求解的研究工作很少.Tang等7研究了一类具有积分约束的抛物型方程的二次最优控制问题,针对最优控制问题,构造了一个全离散的有限元格式,其中:空间离散为有限元离散,时间离散为向后Euler方法.Sun等8研究了具有积分约束的抛物型最优控制问题的向后Euler格式的自适应有限
10、元逼近.王世杰等9考虑了带积分约束的抛物型最优控制问题,对状态变量和伴随状态变量用线性连续函数离散,而控制变量使用分片常数离散,并得到最优的收敛阶.本文主要讨论具有积分控制约束的抛物型最优控制问题的有限元方法.对于状态变量的离散,在空间上采用标准分段线性有限元,在时间上采用Crank-Nicolson格式,控制变量采用分段线性离散,并得到控制变量、状态变量以及伴随状态变量关于时间和空间均为二阶收敛.最后给出了数值算例,验证理论结果.1预备知识设Rn(n=2)凸区域,且具有Lipschitz连续边界.Sobolev空间定义如下:Wm,p=uLp():DuLp(),|m.其中范数记作kkm,p,定
11、义如下:kv km,p,=X|mZ|Dv|pdx1p,p1,),半范数记作|m,p,定义如下:|v|m,p,=X|=mZ|Dv|pdx1p,p1,),其中=(1,n)Nn0为多重指标,|=nPi=1i,D表示阶数为|的微分算子.特别的,当p=2时,记Hm()=Wm,p()且范数kv km,=kv km,2,.设Lr(0,T;Wm,p()为从0,T到Wm,p()范数上所有Lr可积函数构成的Banach空间,其范数定义如下:kv kLr(0,T;Wm,p()=ZT0kv krm,p,dt1r,r1,),定义如下的内积形式:a(u,v)=Zuv dx,u,v H10(),(u,v)=Zuv dx,u
12、,v L2().2模型描述及全离散格式本节主要考虑以下抛物方程分布最优控制问题:minuUadJ(y,u)=12kyydk2L2(0,T;L2()+2kuk2L2(0,T;L2()(1)满足状态方程约束yty=f+Bu,(x,t)Ty=0,(x,t)Ty(,0)=y0,x(2)198新疆大学学报(自然科学版中英文)2024年其中R2凸多边形区域,T=(0,T,T=(0,T.设B是L2(0,T;L2()到L2(0,T;L2()的有界线性算子,积分型控制约束集Uad如下表示:Uad=uL2(0,T;L2():Zudx0.最优控制问题(1)(2)可以如下表示:minuUad12kyydk2L2(0,
13、T;L2()+2kuk2L2(0,T;L2()(yt,w)+a(y,w)=(f+Bu,w),wH10()y(,0)=y0(3)由文献1可知,最优控制问题(3)存在唯一解(y,u)L2(0,T;H10()H2()H1(0,T;H1()Uad.设(y,u)是最优控制问题(3)的解,当且仅当存在伴随状态pL2(0,T;H10()H2()H1(0,T;L2(),使得(y,p,u)满足下面的最优性条件:(yt,w)+a(y,w)=(f+Bu,w),wH10(),t(0,T,y(,0)=y0(4)(pt,q)+a(p,q)=(yyd,q),q H10(),t(0,T,p(,T)=0(5)(u+Bp,vu)
14、0,v Uad,t(0,T(6)其中B是B的伴随算子.变分不等式(6)可表示为u=1nmax(0,Bp)Bpo(7)其中Bp=(RBpdx)/(R1dx)代表Bp在上的积分平均.下面考虑最优控制问题(3)的有限元方法.设Th表示区域上的正规三角形剖分,使得=KThK.有限元空间Vh定义如下:Vh=v C()flflv|KP1(K),K Th,其中P1(K)是所有次不大于1的多项式函数空间.设V0h=VhH10().对于时间离散采用Crank-Nicolson格式.设t=T/N为均匀时间步长,节点用tn=nt(0nN)表示.设vi=v(x,ti),dtvi=(vivi1)/t和|v|=NXn=1
15、tkvn+vn12k20,12.定义Ritz投影Rh:H10()V0h满足如下条件:对任意uH10(),有(Rhu,vh)=(u,vh),vhV0h,和L2正交投影4Qh:UadUad,h满足如下条件:对任意v Uad(vQhv,vh)=0,vhUad,h,接下来给出近似性质:kwRhwk0,Ch2kwk2,第2期张馨丹,等:抛物型最优控制问题的全离散Crank-Nicolson有限元法199kvQhv kr,p,Chr+tkv kt,p,0r,t2.对最优系统(4)(6)建立如下的Crank-Nicolson有限元全离散格式:YnhYn1ht,wh+aYnh+Yn1h2,wh=fn+BUnh
16、+fn1+BUn1h2,wh,Y0h(x)=Rhy0,whV0h,n=1,2,N(8)PnhPn1ht,qh+aPnh+Pn1h2,qh=Ynhynd+Yn1hyn1d2,qh,PNh(x)=0,qhV0h,n=N,N 1,1(9)Unh+Un1h2+BPnh+Pn1h2,vUnh+Un1h20,v Uad,h,n=1,2,N(10)其中Uad,h=v Uadflflv|KP1(K),K Th.3收敛性分析为了得到主要结论,引入如下的辅助变量,对于任意v Uad,设(y(v),p(v)是下列方程的解:yt(v),w+ay(v),w=f+Bv,w,wH10(),t(0,T,y(v)(,0)=y0
17、(11)pt(v),q+ap(v),q=y(v)yd,q,q H10(),t(0,T,p(v)(,T)=0(12)对于任意v Uad,设(yh(v),ph(v)是下列方程的解:yht(v),wh+ayh(v),wh=f+Bv,wh,whV0h,t(0,T,yh(v)(,0)=Rhy0(13)pht(v),qh+aph(v),qh=yh(v)yd,qh,qhV0h,t(0,T,ph(v)(,T)=0(14)对于任意v Uad,设(ynh(v),pnh(v)满足下列方程组:ynh(v)yn1h(v)t,wh+aynh(v)+yn1h(v)2,wh=fn+Bvn+fn1+Bvn12,wh,y0h(v
18、)(x)=Rhy0,whV0h,n=1,2,N(15)pn1h(v)pnh(v)t,qh+apnh(v)+pn1h(v)2,qh=ynh(v)ynd+yn1h(v)yn1d2,qh,pNh(v)=0,qhV0h,n=N,N 1,1(16)定理1设y(u)是方程(11)的解,yh(u)是方程(13)的解,则有如下不等式成立:ky(u)yh(u)k20,Ch4ky k22,+ZT0kyt(u)k21,dt.证明 由方程(11)和(13),可得如下方程yt(u)yht(u),wh+ay(u)yh(u),wh=0.利用Ritz投影的性质,上述方程可以表示为yt(u)Rhyt(u),wh+Rhyt(u)
19、yht(u),wh+aRhy(u)yh(u),wh=0(17)设e=Rhy(u)yh(u),在方程(17)中取wh=Rhy(u)yh(u),利用Young不等式和H older不等式,可得ddtkek20,+C1kek20,C12kek21,+C2kyt(u)Rhyt(u)k21,.200新疆大学学报(自然科学版中英文)2024年利用Rh的性质,则有ddtkek20,+C12kek20,C2h4kyt(u)k21,+C12kek20,(18)对方程(18)在时间上积分,并利用Gronwall引理得kek20,+Zt0kek20,dtCh4ZT0kyt(u)k21,dt(19)利用三角不等式和R
20、h的性质,得ky(u)yh(u)k20,ky(u)Rhy(u)k20,+kRhy(u)yh(u)k20,Ch4ky(u)k22,+Ch4ZT0kyt(u)k21,dt.定理2设yh(u)和ynh(u)分别是方程(13)和(15)的解.当h和t充分小时,存在常数C使得kyh(u)(tn)ynh(u)k20,C(t)4ZT0kyhttt(u)k20,dt(20)证明 在方程(13)中分别取t=tn和t=tn1,则得到如下的方程yht(u)(tn)+yht(u)(tn1)2,wh+ayh(u)(tn)+yh(u)(tn1)2,wh=fn+Bun+fn1+Bun12,wh(21)将方程(15)和(21
21、)相减,得yht(u)(tn)+yht(u)(tn1)2dtynh(u),wh+ayh(u)(tn)+yh(u)(tn1)yn1h(u)ynh(u)2,wh=0(22)设enh=yh(u)(tn)ynh(u),在方程(22)中取wh=enh=(enh+en1h)/2,得12tkenhk20,ken1hk20,+a enh,enhC?dtyh(u)(tn)yht(u)(tn)+yht(u)(tn1)2?0,kenhk0,+ken1hk0,C1(t)1?Ztntn1|(ttn1)(tnt)yhttt(u)|dt?0,kenhk0,+ken1hk0,C2(t)32Ztntn1kyhttt(u)k20
22、,dt12kenhk0,+ken1hk0,C()(t)3Ztntn1kyhttt(u)k20,dt+kenhk20,+ken1hk20,.将上式方程整理得kenhk20,ken1hk20,+C1tk enhk20,C()(t)4Ztntn1kyhttt(u)k20,dt+tkenhk20,+ken1hk20,(23)对方程(23)从n=1到k求和,则有kekhk20,+C1tkXn=1k enhk20,C()(t)4ZT0kyhttt(u)k20,dt+2tkXn=1kenhk20,(24)利用离散Gronwall引理10,则有kekhk20,+c1tkXn=1k enhk20,C()(t)4
23、ZT0kyhttt(u)k20,dt(25)结合定理1和定理2,可以得到以下yn(u)和ynh(u)的误差估计.第2期张馨丹,等:抛物型最优控制问题的全离散Crank-Nicolson有限元法201定理3设y(u)和ynh分别是方程(11)和(15)的解,当h和t充分小时,存在常数C使得kyn(u)ynh(u)k20,C(t)4ZT0kyhttt(u)k20,dt+C1h4ky(u)k22,+ZT0kyt(u)k21,dt(26)定理4在定理3条件成立下,设p(u),ph(u)和pnh(u)分别是方程(12),(14)和(16)的解,则有kpn(u)pnh(u)k20,C(t)4ZT0kpht
24、tt(u)k20,dt+tZT0kyhttt(u)k20,dt+C1h4ky k22,+ZT0kyt(u)k21,dt+kp(u)k22,+ZT0kpt(u)k21,dt.证明 与定理1和定理2的证明类似.定理5设(ynh(u),pnh(u)和(Ynh,Pnh)分别是方程(15)(16)和(8)(9)的解,则kynh(u)Ynhk1,C|uUh|,kpnh(u)Pnhk1,C|uUh|.证明 将方程(8)减去(15)可得dt(Yihyih(u),wh+aYihyih(u)+Yi1hyi1h(u)2,wh=BUihui+Ui1hui12,wh(27)设i=Yihyih(u),在方程(27)中取w
25、h=(ii1)/t,则有dti,dti+ai+i12,ii1t=BUihui+Ui1hui12,ii1t,即dti,dti+12ta(i,i)a(i1,i1)=BUihui+Ui1hui12,dti.利用算子B的有界性以及Young不等式,则有BUihui+Ui1hui12,dtiC()?Uihui+Ui1hui12?20,+kdtik20,(28)对方程(28)从i=1到n求和,从而得到knk21,CNXi=1t?Uihui+Ui1hui12?20,=C|uUh|2.设i=Pihpih(u),类似的a(i1,i1)a(i,i)+Ct?Yihyih(u)+Yi1hyi1h(u)2?20,.将上
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 抛物型 最优 控制 问题 离散 Crank Nicolson 有限元
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。