不具备Nagumo条件的二阶微分方程两点边值问题解的存在性

常微分方程边值问题与不动点定论文

目录 引言 (1) 1预备知识 (2) 定义1.1（奇异Sturm-Liouville边值问题的正解） (2) 引理1.1.1 (2) 定义1.2（凸集的概念) (3) 定义1.3锥的定义 (3) 定义1.4（全连续算子的概念） (3) 1.5 （常微分边值问题的定义） (4) 定义1.6混合单调算子得定义） (4) 2 常微分方程边值问题正解得存在性 (5) 2.1 奇异Sturm-Liouville常微分边值问题的正解存在在 (5) 子 (8) 2.2 一类二阶边值问题的存在性 (9) 3一类混合单调算子应用 (11) 3.1一类混合单调算子的存在唯一性?........................ 错误!未定义书签。 3.2 求常微分边值问题的例题 (13) 结束语 (15) 参考文献 (15) 致 (16)

常微分方程边值问题与不动点定 （数学与统计学院 11级数学与应用数学2班）指导教师：攀峰 引言 从历史上看在有了微积分这个概念以后，紧接着出现了常微分方程。发展初期是属于“求通解”得时代，当人们从初期的热潮中结束要从维尔证明了卡帝方程中是一定不会存在一般性的初等解的时候开始的，并且柯西紧接着又提出了初值问题，常微分方程开始从重视“求通解”转向重视“求定解”的历史时代。 大学我们都学习了常微分方程这门学科，如果要研究它的定解问题，我们首先就会知道是常微分方程的初值问题。然而，在科学技术、生产实际问题中，我们还是提出了另一类定解问题-边值问题。对于常微分方程边值问题，伟大的科学家最早在解决二阶线性微分方程时，提出了分离变量法。[]1.在牛顿时期，科学家们已经提出过常微分的边值问题，牛顿也对常微分边值问题进行过研究，并且在1666年10月牛顿已经在这个领域取得了很大的成就，但是由于种种原因当时并没有整理成论文，所以没有及时出版。但在1687年他终于把在常微分方程上研究的成果发表了，虽然不是在数学著作中，却是他的一本力学著作中（《自然哲学的数学原理》）。 在微积分刚创立时期，雅克.伯努利来自瑞士的科学家提出了远著文明的问题-悬链线问题，紧着的地二年著名数学家莱布尼兹就给出了正确的解答，通过对绳子上个点受力分析，建立了以下方程 这个方程满足的定解条件是y（a）=α；y(b)=β.这是一个典型的常微分方程的边值问题。从这开始，常微分边值问题已经是科学家研究微分方程是不可或缺的工具，我就简单列举几个例子：（比如种族的生态系统；梁的非线性震动）等。对于怎么研究它，

常微分方程边值问题的数值解法

第8章 常微分方程边值问题的数值解法 引 言 第7章介绍了求解常微分方程初值问题的常用的数值方法；本章将介绍常微分方程的边值问题的数值方法。 只含边界条件(boundary-value condition)作为定解条件的常微分方程求解问题称为常微分方程的边值问题(boundary-value problem). 为简明起见，我们以二阶边值问题为 则边值问题(8.1.1)有唯一解。 推论 若线性边值问题 ()()()()()(),, (),()y x p x y x q x y x f x a x b y a y b αβ'''=++≤≤?? ==? (8.1.2) 满足 (1) (),()p x q x 和()f x 在[,]a b 上连续； (2) 在[,]a b 上, ()0q x >, 则边值问题(8.1.1)有唯一解。 求边值问题的近似解，有三类基本方法： (1) 差分法(difference method)，也就是用差商代替微分方程及边界条件中的导数，最终化为代数方程求解； (2) 有限元法(finite element method)；

(3) 把边值问题转化为初值问题，然后用求初值问题的方法求解。 差分法 8.2.1 一类特殊类型二阶线性常微分方程的边值问题的差分法 设二阶线性常微分方程的边值问题为 (8.2.1)(8.2.2) ()()()(),,(),(), y x q x y x f x a x b y a y b αβ''-=<

偏微分方程边值问题的数值解法论文

求解偏微分方程的边值问题 本实验学习使用MATLAB 的图形用户命令pdetool 来求解偏微分方程的边值问题。这个工具是用有限元方法来求解的，而且采用三角元。我们用个例题来说明它的用法。 一、MATLAB 支持的偏微分方程类型 考虑平面有界区域D 上的二阶椭圆型PDE 边值问题： ()c u u f α-??+=g (1.1) 其中 (1) , (2) a,f D c x y ?????=? ????? 是上的已知函数(3)是标量或22的函数方阵 未知函数为(,) (,)u x y x y D ∈。它的边界条件分为三类： （1）Direchlet 条件： hu f = (1.2) （2）Neumann 条件： ()n c u qu g ?+=g (1.3) （3）混合边界条件：在边界D ?上部分为Direchlet 条件，另外部分为Neumann 条件。 其中,,,,h r q g c 是定义在边界D ?的已知函数，另外c 也可以是一个2*2的函数矩阵，n 是沿边界的外法线的单位向量。 在使用pdetool 时要向它提供这些已知参数。 二、例题 例题1 用pdetool 求解 22D 1 D: 10u x y u ??-?=+≤??=?? (1.4)

解：首先在MATLAB 的工作命令行中键入pdetool ，按回牟键确定，于是出现PDE Toolbox 窗口，选Genenic Scalar模式． ( l ）画区域圆 单击椭圆工具按钮，大致在（0，0）位置单击鼠标右键，拖拉鼠标到适当位置松开。为了保证所绘制的圆是标准的单位园，在所绘园上双击，打开 Object Dialog 对话框，精确地输入

偏微分方程数值习题解答

李微分方程数值解习题解答 1-1 如果0)0(' =?，则称0x 是)(x J 的 驻点（或稳定点）.矩阵A 对称（不必正定），求证0x 是)(x J 的驻点的充要条件是：0x 是方程组 b Ax =的解 证明:由)(λ?的定义与内积的性线性性质,得 ),()),((2 1 )()(0000x x b x x x x A x x J λλλλλ?+-++=+= ),(2 ),()(2 00x Ax x b Ax x J λλ+ -+= ),(),()(0'x Ax x b Ax λλ?+-= 必要性:由0)0(' =?,得,对于任何n R x ∈,有 0),(0=-x b Ax , 由线性代数结论知, b Ax b Ax ==-00,0 充分性: 由b Ax =0,对于任何n R x ∈, 0|),(),()0(00'=+-==λλ?x Ax x b Ax

即0x 是)(x J 的驻点. §1-2 补充: 证明)(x f 的不同的广义导数几乎处处相等. 证明:设)(2I L f ∈,)(,221I L g g ∈为)(x f 的广义导数,由广义导数的定义可知,对于任意 )()(0I C x ∞∈?,有 ??-=b a b a dx x x f dx x x g )()()()(' 1?? ??-=b a b a dx x x f dx x x g )()()()('2?? 两式相减,得到 )(0)()(021I C x g g b a ∞ ∈?=-??? 由变分基本引理,21g g -几乎处处为零,即21,g g 几乎处处相等. 补充:证明),(v u a 的连续性条件(1.2.21) 证明: 设'|)(|,|)(|M x q M x p ≤≤,由Schwarz 不等式

常微分方程组(边值)

常微分方程组边值问题解法 打靶法Shooting Method （shooting.m ） %打靶法求常微分方程的边值问题 function [x,a,b,n]=shooting(fun,x0,xn,eps) if nargin<3 eps=1e-3; end x1=x0+rand; [a,b]=ode45(fun,[0,10],[0,x0]'); c0=b(length(b),1); [a,b]=ode45(fun,[0,10],[0,x1]'); c1=b(length(b),1); x2=x1-(c1-xn)*(x1-x0)/(c1-c0); n=1; while (norm(c1-xn)>=eps & norm(x2-x1)>=eps) x0=x1;x1=x2; [a,b]=ode45(fun,[0,10],[0,x0]'); c0=b(length(b),1); [a,b]=ode45(fun,[0,10],[0,x1]'); c1=b(length(b),1) x2=x1-(c1-xn)*(x1-x0)/(c1-c0); n=n+1; end x=x2; 应用打靶法求解下列边值问题： ()()??? ????==- =010004822y y y dx y d 解：将其转化为常微分方程组的初值问题

()????? ? ?????==-==t dx dy y y y dx dy y dx dy x 0011221 048 命令： x0=[0:0.1:10]; y0=32*((cos(5)-1)/sin(5)*sin(x0/2)-cos(x0/2)+1); 真实解 plot(x0,y0,'r') hold on [x,y]=ode45('odebvp',[0,10],[0,2]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,5]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,8]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,10]'); plot(x,y(:,1))

第十一章 常微分方程边值问题的数值解法汇总

第十一章 常微分方程边值问题的数值解法 工程技术与科学实验中提出的大量问题是常微分方程边值问题．本章将研究常微分方程边值问题的数值求解方法.主要介绍三种边界条件下的定解问题和两大类求解边值问题的数值方法,打靶法算法和有限差分方法. 11.1 引言 在很多实际问题中都会遇到求解常微分方程边值问题. 考虑如下形式的二阶常微分方程 ),,(y y x f y '='', b x a <<, (11.1.1) 在如下三种边界条件下的定解问题: 第一种边界条件: α=)(a y , β=)(b y (11.1.2) 第二种边界条件: α=')(a y , β=')(b y (11.1.2) 第三种边界条件: ? ? ?=-'=-'101 0)()()()(b b y b y a a y a y βα, (11.1.13) 其中0 0, ,00000>+≥≥b a b a . 常微分方程边值问题有很多不同解法, 本书仅介绍打靶方法和有限差分方法. 11.2 打靶法 对于二阶非线性边值问题 ()()().,,βα==≤≤'=''b y a y b x a y y x f y ，，， (11.2.1) 打靶法近似于使用初值求解的情况． 我们需要利用一个如下形式问题初值解的序列： ()()v a w a w b x a w w x f w ='=≤≤'='')(,,，，，α， （11.2.2） 引进参数v 以近似原边界值问题的解．选择参数k v v =，以使： ()()β==∞ →b y v b w k k ,lim ， （11.2.3）

其中),(k v x w 定义为初值问题（11.2.2）在k v v =时的解，同时()x y 定义为边值问题（11.2.1）的解． 首先定义参数0v ，沿着如下初值问题解的曲线，可以求出点),(αa 对应的初始正视图 ()()v a w a w b x a w w x f w ='=≤≤'='')(,,，，，α． （11.2.4） 如果),(0v b w 不严格收敛于β，那么我们选择1v 等值以修正近似值，直到),(0v b w 严格逼近β． 为了取得合适的参数k v ，现在假定边值问题（11.2.1）有唯一解，如果),(v x w 定义为初始问题（11.2.2）的解，那么v 可由下式确定： 0),(=-βv b w ． （11.2.5） 由于这是一个非线性方程，我们可以利用Newton 法求解．首先选择初始值0v ，然后由下式生成序列 ),)(()),((111----- =k k k k v b dv dw v b w v v β，此处),(),)(( 11--=k k v b dv dw v b dv dw ， （11.2.6） 同时要求求得),)(( 1-k v b dv dw ，因为),(v b w 的表达式未知，所以求解这个有一点难度；我们只能得到这么一系列的值。 ，，，),(),(),(),(1210-??k v b w v b w v b w v b w 假如我们如下改写初值问题（11.2.2），使其强调解对x 和v 的依赖性 ()()v v a w v a w b x a v x w v x w x f w ='=≤≤'=''),(,),(),,(,，，，α，（11.2.7） 保留初始记号以显式与x 的微分相关．既然要求当k v v =时),)((v b dv dw 的值，那么我们需要求出表达式（11.2.7）关于v 的偏导数．过程如下： )),(),,(,(),(v x w v x w x v f v x v w '??=?''? ),()),(),,(,()),(),,(,(v x v w v x w v x w x w f v x v x w v x w x x f ??'??+??'??= ) ,()),(),,(,(v x v w v x w v x w x w f ?'?''??+ 又因为x 跟v 相互独立，所以当b x a ≤≤上式如下；

常微分方程组(边值)

常微分方程组边值问题解法 打靶法Shooti ng Method (shoot in g.m ) % 丁靶法求常微分方程的边值问题 function [x,a,b ,n]=shooti ng(fu n, xO,x n, eps) if nargin<3 eps=1e-3; end x1=x0+ra nd; [a,b]=ode45(fu n, [0,10],[0,x0]'); c0=b(le ngth(b),1); [a,b]=ode45(fu n, [0,10],[0,x1]'); c1=b(le ngth(b),1); x2=x1-(c1-x n)*(x1-x0)/(c1-c0); n=1; while (no rm(c1-x n)>=eps & no rm(x2-x1)>=eps) x0=x1;x 仁x2; [a,b]=ode45(fu n,[ 0,10],[0,x0]'); cO=b(le ngth(b),1); [a,b]=ode45(fu n,[ 0,10],[0,x1]'); c1= b(le ngth(b),1) x2=x1-(c1-x n)*(x1-x0)/(c1-c0); n=n+1; end x=x2; 应用打靶法求解下列边值问题: y 10 0 解：将其转化为常微分方程组的初值问题

命令: xO=[O:O.1:1O]; y0=32*((cos(5)-1)/si n( 5)*si n(x0/2)-cos(x0/2)+1); plot(xO,yO,'r') hold on [x,y]=ode45('odebvp',[0,10],[0,2]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,5]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,8]'); plot(x,y(:,1)) [x,y]=ode45('odebvp',[0,10],[0,10]'); plot(x,y(:,1)) dy i dx y 2 dy 2 dx y i 0 y 4 y o dy dx X0 真实解 30 ' 12^4567^9 10

几类分数阶微分方程边值问题和初边值问题的

几类分数阶微分方程边值问题和初边值问题的 数值方法 聂宁明 摘要 分数阶微积分实际上已有三百多年的历史,由于缺乏物理、力学背景的支持,它的发展极其缓慢.直到20世纪后期,人们才意识到在石油渗流、地下水污染防治、黏弹性材料、信号与图象处理、控制、量子力学、金融及生命科学等多个领域,分数阶微积分和分数阶微分方程都有重要的应用,从而开始重视对其数值算法的研究. 本文针对几类分数阶常微分方程边值问题和分数阶偏微分方程初边值问题,构造不同的数值方法进行求解并给出相应的误差分析. 第一章介绍分数阶微积分的历史和发展现状,并给出分数阶微积分的基本定义和性质. 第二章考虑Riemann-Liouville分数阶微分方程两点边值问题,先讨论分数阶微分方程两点边值问题解的存在唯一性条件,然后用打靶法对其进行数值求解,并对线性情形,给出误差分析. 第三章介绍如何用三次样条配置法来数值求解分数阶微分方程两点边值问题,并给出误差估计和数值算例. 第四章运用谱方法求解分数阶微分方程的边值问题和初边值问题.第一节用谱方法求解高阶导数是二阶,低阶导数是分数阶的微分方程两点边值问题,分析谱逼近解的收敛性,并通过数值算例验证谱精度.第二节,将谱方法用于求解稳态分数阶对流扩散方程,分析算法的稳定性和收敛性,并通过数值计算验证算法的可行性.第三节以地下水污染问题中抽象出来的空间分数阶扩散方程为例,介绍谱方法求解分数阶偏微分方程初边值问题的过程.对该方程,在空间方向用Galerkin谱方法进行数值逼近,在时间方向用向后Euler差分格式进行离散求解,分析方法的稳定性和收敛性,并通过数值算例验证理论分析结果. 为了进一步说明分数阶微分方程的意义,在第五章中给出了分数阶微分方程在石油渗流问题中的一个应用实例.通过对破裂可形变地层中裂口附近的渗流情况的研究,给出了分数阶微分方程模型的建模过程,并用数值计算的结果说明分数阶微分方程模型在实际中优于整数阶微分方程模型的事实.

微分方程的边值问题

微分方程边值问题的数值方法 本部分内容只介绍二阶常微分方程两点边值问题的的打靶法和差分法。 二阶常微分方程为 (,,),y f x y y a x b '''=≤≤ (1.1) 当(,,)f x y y '关于,y y '为线性时，即(,,)()()()f x y y p x y q x y r x ''=++，此时(1.1)变成线性微分方程 ()()(),y p x y q x y r x a x b '''--=≤≤ (1.2) 对于方程(1.1)或(1.2)，其边界条件有以下3类： 第一类边界条件为 (),()y a y b αβ== (1.3) 当0α=或者0β=时称为齐次的，否则称为非齐次的。 第二类边界条件为 (),()y a y b αβ''== (1.4) 当0α=或者0β=时称为齐次的，否则称为非齐次的。 第三类边界条件为 0101()(),()()y a y a y b y b ααββ''-=+= (1.5) 其中00000,0,0αβαβ≥≥+>，当10α=或者10β=称为齐次的，否则称为非齐次的。微分方程(1.1)或者(1.2)附加上第一类，第二类，第三类边界条件，分别称为第一，第二，第三边值问题。 1 打靶法介绍 下面以非线性方程的第一类边值问题(1.1)、(1.3)为例讨论打靶法，其基本原理是将边值问题转化为相应的初值问题求解。 【原理】假定()y a t '=，这里t 为解()y x 在x a =处的斜率，于是初值问题为 (,,) ()()y f x y y y a y a t α '''=?? =??'=? (1.6) 令z y '=，上述二阶方程转化为一阶方程组

(完整版)二阶常微分方程边值问题的数值解法毕业论文

二阶常微分方程边值问题的数值解法 摘要 求解微分方程数值解的方法是多种多样的，它本身已形成一个独立的研究方向，其要点是对微分方程定解问题进行离散化．本文以研究二阶常微分方程边值问题的数值解法为目标，综合所学相关知识和二阶常微分方程的相关理论，通过对此类方程的数值解法的研究，系统的复习并进一步加深对二阶常微分方成的数值解法的理解，为下一步更加深入的学习和研究奠定基础． 对于二阶常微分方程的边值问题，我们总结了两种常用的数值方法：打靶法和有限差分法．在本文中我们主要探讨关于有限差分法的数值解法．构造差分格式主要有两种途径：基于数值积分的构造方法和基于Taylor展开的构造方法．后一种更为灵活，它在构造差分格式的同时还可以得到关于截断误差的估计．在本文中对差分方法列出了详细的计算步骤和Matlab

程序代码，通过具体的算例对这种方法的优缺点进行了细致的比较．在第一章中，本文将系统地介绍二阶常微分方程和差分法的一些背景材料．在第二章中，本文将通过Taylor展开分别求得二阶常微分方程边值问题数值解的差分格式．在第三章中，在第二章的基础上利用Matlab求解具体算例，并进行误差分析． 关键词：常微分方程,边值问题，差分法，Taylor展开，数值解

The Numerical Solutions of Second-Order Ordinary Differential Equations with the Boundary Value Problems ABSTRACT The numerical solutions for solving differential equations are various. It formed an independent research branch. The key point is the discretization of the definite solution problems of differential equations. The goal of this paper is the numerical methods for solving second-order ordinary differential equations with the boundary value problems. This paper introduces the mathematics knowledge with the theory of finite difference. Through solving the problems, reviewing what have been learned systematically and understanding the ideas and methods of the finite difference method in a deeper layer, we can establish a foundation for the future learning.

Matlab求解常微分方程边值问题的方法

Matlab 求解常微分方程边值问题的方法：bvp4c 函数 常微分方程的边值问题，即boundary value problems ，简称BVP 问题，是指表达形式为 (,)((),())0'=??=?y f x y g y a y b 或(,,)((),(),)0'=??=? y f x y p g y a y b p 的方程组（p 是未知参数），在MA TLAB 中使用积分器bvp4c 来求解。 [命令函数] bvp4c [调用格式] sol=bvp4c(odefun,bcfun,solinit,options,p1,p2,…) sol 为一结构体，sol.x 、sol.y 、sol.yp 分别是所选择的网格点及其对应的y(x)与y'(x)数值; bvp4c 为带边值条件常微分方程积分器的函数命令；odefun 为描述微分方程组的函数文件；bcfun 为计算边界条件g(f(a),f(b),p)=0的函数文件；solinit 为一结构体，solinit.x 与solinit.y 分别是初始网格的有序节点与初始估计值，边界值条件分别对应a=solinit.x(l)和b=solinit.x(end)； options 为bvpset 命令设定的可选函数，可采用系统默认值；p1, p2…为未知参数。 例 求常微分方程0''+=y y 在(0)2=y 与(4)2=-y 时的数值解。 [解题过程] 仍使用常用方法改变方程的形式： 令1=y y ，21'=y y ，则原方程等价于标准形式的方程组1221 ?'=??'=-??y y y y ； 将其写为函数文件twoode.m ； 同时写出边界条件函数对应文件twobc.m ； 分别使用结构solinit 和命令bvp4c 确定y-x 的关系； 作出y-x 的关系曲线图。 [算例代码] solinit =bvpinit(linspace(0,4,5),[1 0]); % linspace(0,4,5)为初始网格，[1,0]为初始估计值 sol=bvp4c(@twoode,@twobc,solinit); % twoode 与twobc 分别为微分方程与边界条件的函数，solinit 为结构 x=linspace(0,4); %确定x 范围 y=deval(sol,x); %确定y 范围 plot(x,y(1,:)); %画出y-x 的图形 %定义twoode 函数（下述代码另存为工作目录下的twoode.m 文件） function dydx= twoode(x,y) %微分方程函数的定义 dydx =[y(2) -abs(y(1))]; %定义twobc 函数（下述代码另存为工作目录下的twobc.m 文件） function res= twobc(ya,yb); %边界条件函数的定义 res=[ya(1);yb(1)+2];

二阶常微分方程边值问题的数值解法

摘要 本文主要研究二阶常微分方程边值问题的数值解法。对线性边值问题，我们总结了两类常用的数值方法，即打靶法和有限差分方法，对每种方法都列出了详细的计算步骤和Matlab程序代码，通过具体的算例对这两类方法的优缺点进行了细致的比较。 关键字：常微分方程边值问题；打靶法；差分法；

ABSTRACT This article mainly discusses the numerical methods for solving Second-Order boundary value problems for Ordinary Differential Equations. On the one hand, we review two types of commonly used numerical methods for linear boundary value problems, i.e. shooting method and finite difference method. For each method, we give both the exact calculating steps , we compare the advantages and disadvantages in detail of these two methods through a specific numerical example. Key words：Boundary-Value Problems for Ordinary Differential Equations；Shooting Method；Finite Difference Method；

常微分方程边值问题的数值解法

常微分方程论文 题目：常微分方程边值问题的数值解法 组长：数学132文洲 组员：数学131王琦 数学132姚瑶 信息132郭斌 院（系）：理学院 指导教师：岳宗敏 时间：2015年6月9日

常微分方程边值问题的数值解法 摘要：作为一类定解问题，补充条件由以自变量取某些值时，未知函数及其导数的值而定，称其为边值条件。许多物理和数学问题都归结为边值问题。本文介绍边值问题的待定常数法和格林函数。 关键词：边值问题 待定常数法 格林函数 Abstract: as a kind of definite solution problems, the supplementary conditions by took the certain values in the independent variable, the value of the unknown function and its derivative, referred to as boundary value conditions. Many physical and mathematical problems boil down to boundary value problems. In this paper, the boundary value problem of the method of undetermined constants and green's function. Keywords: boundary value problem Method of undetermined constants Green's function 11.1 引言 在很多实际问题中都会遇到求解常微分方程边值问题. 考虑如下形式的二阶常微分方程 ),,(y y x f y '='', b x a <<, (11.1.1) 在如下三种边界条件下的定解问题: 第一种边界条件: α=)a (y , β=)(b y (11.1.2) 第二种边界条件: α=')(a y , β=')(b y (11.1.2) 第三种边界条件: ?? ?=-'=-'1 01 0)()()()(b b y b y a a y a y βα, (11.1.13) 其中0 0, ,00000>+≥≥b a b a .

椭圆型偏微分方程边值问题的一种数值解

椭圆型偏微分方程边值问题的一种数值解 为了解不规则区域上的椭圆型偏微分方程边值问题, 首先要对区域进行剖分，这样做使得在整个解题过程中进行了两次边值问题的求解。在学习中得到启发看到了一个方法，它将区域剖分的问题及求解的问题结合起来进行, 使整个求解过程得到简化这个方法求得的是未知函数的一组等值线,这在某些物理问题中是方便的。 （1） 其中Ω是区域；Γ 1、 Γ 2、 Γ 3、 Γ4Ω的边界。且Γ 1、 Γ3相对，Γ 2、 Γ4相对。 公式的系数分别是Ω上的连续函数。φ1φ2是单调函数但可以不连续。u 0，u n 是常数。又设d>0,c<=0,u n >u 0.特殊的，Γ1、Γ2、Γ3、Γ4中至多有两个可以退化为一点。为了求解上式，引入辅助问题 （2） 00:;m m v v v v <其中、是常数且 34??、是单调函数, 也可以不连续， 034m v v ??、、、可按解题方便来选取作变换 (3) 变换（3）区域Ω变为Ω`由椭圆型方程的性质可见(3)是可逆的。 设（3）的逆变换是 （4） 变换（3）将（1）（2）中的方程变为

（5） （6） 其中： ，易见仍有即式（3）和（6）是一个拟线性椭圆型方程组。设曲线的几何方程分别是 解下面四组联立方程 并分别记它们的解为 于是（3）将（1）（2）、中的边界条件变为 （7）现将方程（5）（6）加上边界条件（7）称为问题（1`）向题(1`), 虽然方程复杂, 但定解区域是矩形,用差分法离散, 迭代法求解是很方便的。(1`) 的解形如(4).将u 视为常数, v是参数, (4)就是u的等值线的参数方程。 参考文献 1、刘家琦。应用求解拉普拉斯方程的边值问题建立有限元网格。计算数学1988，5（1）：1~9 2、李子才。具有奇点的Laplace方程边值问题的原始能量有限元结合法。计算数学，1980，2（4）：319~328