正倒向随机微分方程组的数值解法_赵卫东

第十章 偏微分方程数值解法偏微分方程问题，其求解十分困难。

除少数特殊情况外，绝大多数情况均难以求出精确解。

因此，近似解法就显得更为重要。

本章仅介绍求解各类典型偏微分方程定解问题的差分方法。

§1 差分方法的基本概念1.1 几类偏微分方程的定解问题椭圆型方程：其最典型、最简单的形式是泊松（Poisson ）方程),(2222y x f yu x u u =∂∂+∂∂=∆ 特别地，当0),(≡y x f 时，即为拉普拉斯（Laplace ）方程，又称为调和方程2222=∂∂+∂∂=∆yux u u Poisson 方程的第一边值问题为⎪⎩⎪⎨⎧Ω∂=Γ=Ω∈=∂∂+∂∂Γ∈),(),(),(),(),(2222y x y x u y x y x f y ux u y x ϕ 其中Ω为以Γ为边界的有界区域，Γ为分段光滑曲线，ΓΩ称为定解区域，),(y x f ，),(y x ϕ分别为Ω，Γ上的已知连续函数。

第二类和第三类边界条件可统一表示为),(),(y x u u y x ϕα=⎪⎪⎭⎫ ⎝⎛+∂∂Γ∈n 其中n 为边界Γ的外法线方向。

当0=α时为第二类边界条件， 0≠α时为第三类边界条件。

抛物型方程：其最简单的形式为一维热传导方程220(0)u ua a t x∂∂-=>∂∂ 方程可以有两种不同类型的定解问题：初值问题⎪⎩⎪⎨⎧+∞<<∞-=+∞<<-∞>=∂∂-∂∂x x x u x t x u a tu )()0,(,0022ϕ初边值问题221200,0(,0)()0(0,)(),(,)()0u ua t T x l t x u x x x lu t g t u l t g t t Tϕ⎧∂∂-=<<<<⎪∂∂⎪⎪=≤≤⎨⎪==≤≤⎪⎪⎩其中)(x ϕ，)(1t g ，)(2t g 为已知函数，且满足连接条件)0()(),0()0(21g l g ==ϕϕ边界条件)(),(),(),0(21t g t l u t g t u ==称为第一类边界条件。

Advances in Applied Mathematics 应用数学进展, 2023, 12(4), 1698-1703 Published Online April 2023 in Hans. https:///journal/aam https:///10.12677/aam.2023.124176波浪能发电装置输出功率的微分方程模型马家耀1，吴彬睿1，吕 平2*1杭州师范大学经亨颐教育学院，浙江 杭州 2杭州师范大学数学学院，浙江 杭州收稿日期：2023年3月24日；录用日期：2023年4月18日；发布日期：2023年4月27日摘要波浪能是当今社会一种十分具有前景的清洁可再生能源，提高波浪能发电装置的转化效率尤为重要。

针对波浪能发电装置只做垂荡运动的情况，本文建立基于微分方程的垂荡模型。

以垂荡模型为基础，本文进一步计算波浪能发电装置的平均输出功率，并建立以波浪能发电装置平均输出功率最大为目标的最优阻尼系数模型。

最后，本文通过多重搜索算法，遍历求解得到最大输出功率及对应的最优阻尼系数。

关键词微分方程，多重搜索法，波浪能发电装置，阻尼系数The Differential Equation Model for the Output Power of Wave Energy ConverterJiayao Ma 1, Binrui Wu 1, Ping Lv 2*1Jing Hengyi School of Education, Hangzhou Normal University, Hangzhou Zhejiang 2School of Mathematics, Hangzhou Normal University, Hangzhou ZhejiangReceived: Mar. 24th, 2023; accepted: Apr. 18th, 2023; published: Apr. 27th, 2023AbstractNowadays wave energy as the clean and renewable energy resource has a promising future, so it is particularly significant to enhance the conversion efficiency of wave energy converter. Aiming at the situation that the wave energy converter only does vertical motion, this paper establishes a vertical oscillation model based on differential equation. In terms of this model, the paper further calculates the average output power of wave energy converter, and establishes an optimal damp-ing coefficient model with the goal of the maximum average output power. Finally, the maximum output power and the corresponding optimal damping factor are obtained by traversing the solu-tion through multiple search algorithms.*通讯作者。

数学年刊A辑2021,42(1):75-88DOI:10.16205/ki.cama.2021.0007条件平均场随机微分方程的最优控制问题吴霜1提要作者研究了一个条件平均场随机微分方程的最优控制问题.这种方程和某些部分信息下的随机最优控制问题有关，并且可以看做是平均场随机微分方程的推广.作者以庞特里雅金最大值原理的形式给出最优控制满足的必要和充分条件.此外，文中给出一个线性二次最优控制问题来说明理论结果的应用.关键词条件平均场随机微分方程，随机最大值原理，倒向随机微分方程，线性二次最优控制，黎卡堤方程MR(2000)主题分类93E20,60H10中图法分类0225,0231.3文献标志码A文章编号1000-8314(2021)01-0075-141引言本节先给出在以后章节中出现的符号.令(C,兀｛兀｝0€圧丁，町是一个完备的概率空间，这里的T>0是一个常数,兀=<7(叭s);0W s5尸=斤，而｛VT(t)｝舜o表示一个标准的d维布朗运动.令QtQTt是一个给定的子-代数，它表示控制者在时刻t所能获得的信息.通篇我们以K"表示n维欧式空间，以R nxd表示nxd矩阵空间，以上标r表示向量或矩阵的转置.对给定的欧式空间，我们以〈•,•〉(〔•|)表示内积(范数)，以U 表示一个非空凸子集.此外若一个卯-值,兀-适应的(或弘适应的)随机过程如)满足妙(t)Fdt<+oo,我们记为讽・)e仍(0,T;肿)(或仍(0,T;腐));若一个肿-值,兀-适应的连续随机过程炉(・)满足E[sup⑷⑴鬥<+oo,我们记为炉(・)e S|(0,T;R").研究的初衷始于如下部分信息下的随机控制问题，其中的状态方程是一个受控的随机微分方程(SDE)'dX(t)=雎,X(t),叫X(t)|$],"(t),叫u(t)覘)dt<+a{t,X(t),叫X(t)|$],呃),叫u(t)覘)dW(t),X(0)=x,指标泛函为』(“(•))=町/T lit,X(t),IE[X(t)|$],u(t),IE["(t)|$])dt+/z(X(T))],而w(-)e U ad是容许控制•这里的容许控制集为U ad={"(•)e L^(O,T;R fe)|u(t)G uyt e[0,T]}.(1-1) (1-2) (1.3)本文2019年5月11日收到，2020年11月15日收到修改稿.1■中国石油大学(华东)应用数学系，山东青岛266580.E-mail:namozhuntipusa©76数学年刊A辑42卷然后我们考虑如下的控制问题.CMF问题对于(1.1)和(1.2),寻找iT(J e%山使得丿("*(•))=inf丿("(•))・(1.4)U(-)eWad方程(1.1)是更为一般的平均场型随机微分方程.因其含有条件期望,我们称之为条件平均场随机微分方程(CMF-SDE).显然,CMF-SDE可以看做是对近年来广受学者关注的平均场随机微分方程(MF-SDEs)的推广.MF-SDEs通常用于描述大量个体的群体行为，比如交互作用的粒子系统，其研究工作可以追溯到文[1]并由此引出了所谓的McKean-Vlasov 随机微分方程•从那时起，陆续有很多关于MF-SDEs的研究文献,领域涉及方程理论、最优控制、微分对策以及金融工程等•参见文[2-8],本文组织如下•第二节会给出有关CMF-SDEs的基本结果.第三节规范CMF-SDEs 的最优控制问题并建立最优性条件(最大值原理).第四节讨论了一个条件平均场型的线性二次最优控制问题(LQ问题)来展示所获结果的应用.最后，在第五节对本文做了最后的总结.2基本结果这一节给出CMF-SDEs解的存在唯一性，尽管一个CMF-SDE含有条件期望项，但仍然可以运用不动点理论来证明这一点.现在从下面的CMF-SDE开始.<'dX(t)=b(t,X(t),E[X(t)|$])dt+c(t,X(t),E[X(t)£])dW(t),〔X(0)=z€1R",i•这里&:x[0,T]x R n x R"->R",o-:Q x[0,T]x R"x R n R nXd是给定的可测映射并满足(H2.1)b(-,x,x)和是兀-适应的；存在C>0,使得\b(t,x,x)-b(t,y,y)\+ \cr(t,x,x}-<7(i,y,y)|W C(\x—示|+\y-y\),t e\0,T],x,x,y,y e R n,B卩b和是一致Lipschitz的；(H2.2)&(-,0,0)e^(0,7;®"),ct(-,O,O)e L^(0,T;R nXd).以下定理确保了方程⑵1)解的存在唯一性.定理2.1假设(H2.1)和(H2.2)成立，则CMF-SDE(2.1)有唯一的兀-适应解X(・)e L^(0,T;R n).证为简化论述，对于认)e仍(o,T；R"),使用符号x(t)=叫叹)|羽并引入如下SDE:dX(t)=&(t,a:(t),方(t))dt+(i,a:(i), ))d TV(f),(2-2)、X(0)=x.1期吴霜CMF-SDE的最优控制问题77由于(2.2)有唯一的解X(.)e0纟(0,丁;氏")，我们可以定义一个映射I:0纟(O,T;R")T 仍(0,T;肿)，使得X(・)=弘(.)].现需要证明I是压缩映射.为此，引入Banach空间0纟(0,丁;氏")的等价范数~T丄归(•川0=何/。

微分方程求解-解微分方程微分方程求解求解微分方程：简单地说，就是去微分，将方程化成自变量与因变量关系的方程。

近来做毕业设计遇到微分方程问题，搞懂后，特发此文，来帮广大同学，网友。

1.最简单的例子：——————》求微分方程的通解。

dx解方程是可分离变量的，分离变量后得两端积分：得：从而：又因为。

仍是任意常数，可以记作C 。

非齐次线性方程2y 求方程的通解解：非齐次线性方程。

先求对应的齐次方程的通解。

5，，用常数变易法：把C换成u(x)，即令则有，dx12，代入原方程式中得两端积分，得。

33再代入式即得所求方程通解。

3法二：假设待求的微分方程是:我们可以直接应用下式得到方程的通解，其中，2，代入积分同样可得方程通解5，3232.微分方程的相关概念：(看完后你会懂得各类微分方程)一阶微分方程：或可分离变量的微分方程：一阶微分方程可以化为的形式，解法：得：称为隐式通解。

，即写成的函数，解法：dxxydydududxduy设，则，，分离变量，积分后将代替u，齐次方程：一阶微分方程可以写成即得齐次方程通解。

一阶线性微分方程：当时,为齐次方程，当时，为非齐次方程，，全微分方程：如果中左端是某函数的全微分方程，即：应该是该全微分方程的通解。

二阶微分方程：时为齐次时为非齐次二阶常系数齐次线性微分方程及其解法：，其中p,q为常数；求解步骤：1、写出特征方程：，其中r2，r的系数及常数项恰好是(*)式中的系数；2、求出式的两个根r1,r23、根据r1,r2的不同情况，按下表写出(*)式的通解：，p,q为常数型，为常数；型3.工程中的解法：四阶定步长Runge-Kutta算法其中h 为计算步长，在实际应用中该步长是一个常数，这样由四阶Runge-Kutta算法可以由当前状态变量Xt 的值求解出下状态变量Xt +1 的值亲们，你们满意吗？一阶微分方程的解一阶微分方程的常数变易法的应用探析The exploration of linear ordinary differential equation of first order with method of leadingvariables作者：刘*专业：数学与应用数学指导老师：杜* *完成时间：2016年9月1号摘要常数变易法是作为求解一阶线性方程的解法给出的。