边值问题和本征值.pdfVIP

第第第6节节节 边边边值值值问问问题题题和和和本本本征征征值值值 二阶常微分方程 y00 = f(x, y, y0) (a ≤ x ≤ b)的边值条件有: 第一类: y(a) = α 或 y(b) = β 第二类: y0(a) = α 或 y0(b) = β 0 0 第三类: y (a) ? α0y(a) = α1 或 y (b) + β0y(b) = β1 微分方程加上第一、第二、第三边值条件,分别称为第一、第二、第三边 值问题. 1. 边值问题的直接积分 初值问题 边值问题 自变量是时间变量 自变量是空间变量 在起始点有两个条件 两个边界上各有一个边界条件 解法思路 利用初值问题解法 在任一边界上补充一个猜测的边界条件，按照初值问题来解方程，所得的 解通常不会满足另一端的边界条件，需要改变猜测的边界条件或者对所得 的解加以某些修正，重新解方程，直到找出解为止。 1 例题:电荷分布为ρ(r) = e?r 求所产生的静电势Φ. 8π 静电势Φ满足泊松方程 52Φ = ?4πρ 对于球对称的ρ 和Φ, 这个方程简化为 · ? 1 d dΦ r2 = ?4πρ r2 dr dr 作代换 Φ(r) = r?1φ(r) 上述方程化为 d2φ = ?4πrρ dr2 这种电荷分布的总电荷是 Z Z ∞ Q = ρ(r)d3r = ρ(r)4πr2dr = 1 0 直接积分方程,得到的精确解是 1 φ(r) = 1 ? (r + 2)e?r 2 由此可得静电势Φ = r?1φ. 可以看出这个解的行为特性是：当r = 0时它等于 零，当r大时，φ → 1,它相当于Φ → r?1即单位电荷产生的库仑势。 按照上面讨论,构成如下边值问题 作为以原点为起点的初值问题 2 已知φ(0) = 0 d φ 1 ?r = ? re 用解析解计算 dr2 2 dφ/dr|r=0 φ|r=0 , φ|r→∞ → 1 所用程序如下 function bzwt r=0:0.1:20; s= - 0.5*r.*exp(-r); %电势 exact = 1 - 0.5* (r + 2).* exp(-r); dexact = exact(2)/0.1; [x,y0]=ode45(@bzfun,r,[0,dexact]); figure(1) plot(x,y0(:,1), r,exact,’:’) [x,y1]=ode45(@bzfun,r,[0,dexact+0.01]); m=(y1(end,1)-y1(end-10,1))/(x(end)-x(end-10)); phi=y1(:,1)-m*x; figure(2) plot(x,y1(:,1) ,x,phi,’-.’) function ydot=bzfun(x,y) ydot = [ y(2); - 0.5*x.*exp(-x) ]; 图5.1 (a)解析解的曲线与用精确解