201数理方程《数学物理方程》陈才生主编

研究生数理方程期末试题-10-11-1-A-答案

北京交通大学硕士研究生2010-2011学年第一学期 《数学物理方程》期末试题(A 卷) (参考答案) 学院__________ 专业___________ 学号 __________ 姓名____________ 1、( 10分)试证明：圆锥形枢轴的纵振动方程为： 玫［I h .丿&」V h .丿& 其中E是圆锥体的杨氏模量，「是质量密度，h是圆锥的高(如下图所示) 【提示：已知振动过程中，在x处受力大小为ES ，S为x处截面面积。】 ex 【证明】在圆锥体中任取一小段，截面园的半径分别是r1和r2，如图所示。于是，我们有 2、：：u(x dx,t) 2 u(x,t) — 2 u2(x,t) E( D) E( * ) ( A )dx 于 x x t r1 = (h「x)tan : r2= (h _(x dx)) tan : 上式化简后可写成

2 2 ：：U(X,t) 2 :：u(x,t) 2, ；u (x,t) E[(h -x) 卜亠 & -(h -'X) 〔x J - - (h -'X)dx 2 从而有 E ：：[(^x)2；：U(x ,t)H-(^x)2：：u2(x,t) .x :X :t 或成 2 ::[(1「)2汽("]“2(1「)小叩) .x h ：：x h ；：t 其中a^E ，证明完毕。 2、 (20分)考虑横截面为矩形的散热片， 它的一边y=b 处于较高温度U ，其它三边y=0. x = 0和x = a 则处于冷却介质中，因而保持较低的温度 u o 。试求该截面上的稳定温度 分布u(x,y)，即求解以下定解问题: u|y 卫二 %, u|y 生二 U, 0 x a. 【提示：可以令u(x, y)二u 0 v(x, y)，然后再用分离变量方法求解。】 【解】令u(x, y) v(x, y)，则原定解问题变为 Wl x￡=0, V=0, 0cy

数理方程与特殊函数教学大纲

数理方程与特殊函数 课程简介：本课程为电子与通信工程类专业的基础课。学分2，周学时2。本课程由“数学物理方程”与“特殊函数”两大部分组成。“数学物理方程”讲授物理学的一个分支——数学与物理所涉及的偏微分方程。主要介绍物理学中常见的三类偏微分方程及其有关的定解问题和这些问题的几种常用解法。“特殊函数”讲授贝塞尔函数与勒让德多项式，以及如何利用这两种特殊函数来解决数学物理方程的一些定解问题的过程。 教学目的与基本要求：通过数理方程与特殊函数课程的学习，使学生系统的掌握工程数学中数学物理方法的知识和技能，培养学生分析问题解决问题的能力，为后续课程的学习及研究奠定重要的数学基础。本课程的先修课程为：高等数学，复变函数，积分变换 主要教学方法：课堂讲授与课外习题。 第零章预备知识（4学时） 复习先修课程中相关的一些内容，主要包括：二阶线性常微分方程解的结构以及常系数情形解的求法；积分学中的一些重要公式和技巧；傅里叶（Fourier）分析；解析函数的极点及其留数；拉普拉斯（Laplace)变换。 第一章典型方程和定解条件的推导（４学时） 在讨论数学物理方程的求解之前，应建立描述某种物理过程的微分方程，再把一个特定物理现象所具有的具体条件用数学形式表达出来。本章学习的重点和难点是了解数学物

理方程的推导及定解问题的确定过程，学会推导一些简单物理过程的微分方程并能确定某些具体物理现象的定解条件。 第一节基本方程的建立 通过几个不同的物理模型，推导出数学物理方程中的三种典型偏微分方程：波动方程、电磁场方程和热传导方程。 第二节初始条件与边界条件 方程决定了物理规律的数学形式，但具体的物理问题所具有的特定条件也应用数学形式表达出来。用以说明某一具体物理现象的初始状态的条件称为初始条件，用以说明其边界上约束情况的条件称为边界条件。 第三节定解问题的提法 由于每一个物理过程都处在特定的条件之下，所以我们要求出偏微分方程适合某些特定条件的解。初始条件和边界条件都称为定解条件。把某个偏微分方程和相应的定解条件结合在一起，就构成了一个定解问题。 本章习题：3-5题 第二章分离变量法（8学时） 本章主要介绍在求解偏微分方程的定解问题时，如何设法把它们转化为常微分方程来求解。本章学习的重点和难点是掌握分离变量法这一“化繁为简”的典型方法的实质，学会求解常见的定解问题。

数理方程期末考试试题

2013-2014学年度第二学期数理方程（B ）期末考试试题 考后回忆版本 一、求下列偏微分方程的通解),(y x u u =（16分） （1）y x y x u 22=???（2）xy x u y x u y =??+???2二、求下列固有之问题的解。要求明确指出固有值及其所对应的固有函数（10分） ?????=′+∞<<<=+′+′′.0)2(,)0()20(,022y y x y x y x y x λ三、求第一象限}0,0|),{(2 >>∈=y x R y x D 的第一边值问题的Green 函数。（12分） 四、用积分变换法求解下列方程。（12分）???=>+∞<<<=).21(),0(,)(),0(. 1)1,(,0)0,()0,10(,4x x u x x x u t u t u t x u u t xx tt δ?七、用分离变量法求解下列方程。（15分） ?????=<++=++=++0|)1(,1 222222z y x zz yy xx u z y x z u u u 八、求解下列定解问题。（5分） ?????==>+∞<

研究生数理方程期末试题10111A答案

《数学物理方程》期末试题（A 卷） （参考答案） 学院 专业 学号 姓名 1、 （10分）试证明：圆锥形枢轴的纵振动方程为： 其中E 是圆锥体的杨氏模量，ρ是质量密度，h 是圆锥的高（如下图所示）： 【提示：已知振动过程中，在x 处受力大小为u ES x ??，S 为x 处截面面积。】 【证明】在圆锥体中任取一小段，截面园的半径分别是1r 和2r ，如图所示。于是，我们有 上式化简后可写成 从而有 或成 其中2 E a ρ = ，证明完毕。 2、 （20分）考虑横截面为矩形的散热片，它的一边y b =处于较高温度U ，其它三边0y =， 0x =和x a =则处于冷却介质中，因而保持较低的温度0u 。试求该截面上的稳定温度 分布(,)u x y ，即求解以下定解问题： 【提示：可以令0(,)(,)u x y u v x y =+，然后再用分离变量方法求解。】 【解】令0(,)(,)u x y u v x y =+，则原定解问题变为 分离变量：

代入方程得到关于X 和Y 的常微分方程以及关于X 的定解条件： 可以判定，特征值 特征函数 利用特征值n λ可以求得 于是求得特征解 形式解为 由边界条件，有 得到 解得 最后得到原定解问题的解是 3、 （20分）试用行波法求解下列二维半无界问题 【解】方程两端对x 求积分，得 也即 对y 求积分，得 也即 由初始条件得 也即 再取0x =，于是又有 从而得 于是 将这里的()g x 和()h y 代入(,)u x y 的表达式中，即得 4、 （20分）用积分变换法及性质，求解无界弦的自由振动问题： 【提示：可利用逆Fourier 积分变换公式：11 ,||sin []20, ||x at a t F a a x at ωω-?

数理方程第二版 课后习题答案教学教材

数理方程第二版课后 习题答案

第一章曲线论 §1 向量函数 1. 证明本节命题3、命题5中未加证明的结论。 略 2. 求证常向量的微商等于零向量。 证：设，为常向量，因为 所以。证毕 3. 证明 证： 证毕4. 利用向量函数的泰勒公式证明：如果向量在某一区间内所有的点其微商为零，则此向量在该区间上是常向量。 证：设，为定义在区间上的向量函数，因为

在区间上可导当且仅当数量函数，和在区间上可导。所以，，根据数量函数的Lagrange中值定理，有 其中，，介于与之间。从而 上式为向量函数的0阶Taylor公式，其中。如果在区间上处处有，则在区间上处处有 ，从而，于是。证毕 5. 证明具有固定方向的充要条件是。 证：必要性：设具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，于是。充分性：如果，可设，令，其中为某个数量函数，为单位向量，因为，于是 因为，故，从而 为常向量，于是，，即具有固定方向。证毕

6. 证明平行于固定平面的充要条件是。 证：必要性：设平行于固定平面，则存在一个常向量，使得，对此式连续求导，依次可得和，从而，，和共面，因此。 充分性：设，即，其中，如果，根据第5题的结论知，具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，任取一个与垂直的单位常向量，于是作以为法向量过原点的平面，则平行于。如果，则与 不共线，又由可知，，，和共面，于是，其中，为数量函数，令，那么，这说明与共线，从而，根据第5题的结论知，具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，作以为法向量，过原点的平面，则平行于。证毕 §2曲线的概念 1. 求圆柱螺线在点的切线与法平面的方程。 解：，点对应于参数，于是当时，， ，于是切线的方程为：

数理方程期末试题B答案

北 京 交 通 大 学 2007-2008学年第二学期《数理方程与特殊函数》期末考试试卷（B ） （参考答案） 学院_ ____________ 专业___________________ 班级________ ____ 学号_______________ 姓名___________ __ 一、 计算题（共80分，每题16分） 1． 求下列定解问题（15分） 2． 用积分变换法及性质，求解半无界弦的自由振动问题：（15分） 3． 设弦的两端固定于0x =及x l =，弦的出示位移如下图所示。初速度为零，又没有外力 作用。求弦做横向振动时的位移(,)u x t 。 [ 解 ] 问题的定解条件是 由初始条件可得 4． 证明在变换, x at x at ξη=-=+下，波动方程xx tt u a u 2=具有形式解0=n u ξ，并由此求 出波动方程的通解。 5． 用分离变量法解下列定解问题 [ 提示：1) 可以直接给出问题的固有函数，不必推导；2) 利用参数变易法。] [ 解 ] 对应齐次方程的定解问题的固有函数是x l n π sin ，其解可以表示成 把原问题中非齐次项t x t x f l a l π π22sin sin ),(=按照固有函数展开成级数 因此有 利用参数变易法，有 于是 6． 用Bessel 函数法求解下面定解问题 [ 解 ] 用分离变量法求解。令)()(),(t T R t u ρρ=，则可得

以及 设0ρβλn n = 为Bessel 函数)(0x J 的正零点，则问题（II ）的特征值和特征函数分别为 问题（I ）的解为 于是原问题的解是 由初始条件 得到 故 于是最后得到原问题的解是 二、 证明题（共2分，每题10分） 7． 证明平面上的Green 公式 其中C 是区域D 的边界曲线，ds 是弧长微分。 [证明] 设),(),,(y x Q y x p 在D+C 上有一阶连续偏导数，n 为C 的外法线方向，其方向余弦为βαcos ,cos ，则有 再设u,v 在D 内有二阶连续偏导数，在D+C 上有一阶连续偏导数，令 得到 交换u,v ，得到 上面第二式减去第一式，得到 证毕。 8． 证明关于Bessel 函数的等式：

MATLAB编程求解二维泊松方程

%%%% 真解 u=sin(pi*x)*sin(pi*y) %%% %%%% 方程 -Laplace(u)=f %%%%%% %%%% f=2*pi^2*sin(pi*x)*sin(pi*y) %%%%%% %%%%difference code for elliptic equations with constant coefficient %%%%% %clear all %clc N=20; h=1/N; S=h^2; x=0:h:1; y=0:h:1; %%% Stiff matrix A=zeros((N-1)^2,(N-1)^2); for i=1 A(i,i)=4/h^2; A(i,i+1)=-1/h^2; A(i,i+(N-1))=-1/h^2; end for i=N-1 A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,2*i)=-1/h^2; %A(i,i+(N-1))=-1/h^2 end for i=(N-2)*(N-1)+1 A(i,i-(N-1))=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; end for i=(N-1)^2 A(i,i-(N-1))=-1/h^2; A(i,i)=4/h^2; A(i,i-1)=-1/h^2; end for n=2:N-2 i=(N-2)*(N-1)+n; A(i,i-(N-1))=-1/h^2; A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; end for i=2:N-2

A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; A(i,i+(N-1))=-1/h^2; end for m=1:N-3 i=m*(N-1)+1; A(i,i-(N-1))=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; A(i,i+(N-1))=-1/h^2; end for m=2:N-2 i=m*(N-1); A(i,i-(N-1))=-1/h^2; A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,i+(N-1))=-1/h^2; end % for m=1:N-3 % i=m*(N-1)+(N-1); % A(i,i-(N-1))=-1/h^2; % A(i,i-1)=-1/h^2; % A(i,i)=4/h^2; % A(i,i+(N-1))=-1/h^2; % end for m=1:N-3 for n=2:N-2 i=m*(N-1)+n; A(i,i-(N-1))=-1/h^2; A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; A(i,i+(N-1))=-1/h^2; end end %%% Right term F=zeros((N-1)^2,1); for m=0:N-2 for n=1:N-1 i=m*(N-1)+n;

数理方程与特殊函数试卷 3套

2010年6月 一、填空题（20分） 1、微分方程的固有值为 ____________，固有函数为____________。 2、勒让德多项式的母函数为________________________。 3、一长为的均匀直金属杆，x=0端固定，x=l端自由，则纵向震动过程中的边界条件为 ________________________。 4、二阶线性偏微分方程属于____________型方程。 5、微分方程，在条件下的拉氏变换表 达式为____________________________________。 6、埃尔米特多项式的微分表达式为____________________________________。 7、函数是区域内的调和函数，它在上有一阶连续偏导数，则 ____________. 8、定解问题的解为________________________。 9、在第一类奇次边界条件下=____________。 10、=____________，=____________。 二、证明题（10分） 三、建立数学物理方程（10分） 一长为l、截面积为s、密度为、比热容为的均匀细杆，一端保持零度，另一端有恒定的热量q流入，初始温度为试建立热传导方程，写出定界条件（要有必要的步骤）。四、写出下列定解问题的解（35分） 1、

2、 3、 五、将函数展开为广义傅里叶级数（25分） 1、设是的正零点，试将函数展开成的傅里叶贝塞尔级数。 2将函数按埃尔米特多项式展开成级数。 2009年6月 一、填空题（20分） 11、微分方程的固有值为 ____________，固有函数为____________。 12、勒让德多项式的母函数为________________________。 13、一长为的均匀直金属杆，x=0端温度为零，x=l端有恒定的热流流出，则热传导过 程中的边界条件为________________________。 14、二阶线性偏微分方程属于____________型方程。 15、微分方程，在条件下，其拉氏 变换表达式为____________________________________。 16、埃尔米特多项式的微分表达式为____________________________________。 17、函数是区域内的调和函数，它在上有一阶连续偏导数，则 ____________. 18、定解问题的解为 ________________________。 19、在第一类奇次边界条件下=____________。 20、=____________，=____________。 二、证明题（10分）

数理方程试卷及答案2

长沙理工大学考试试卷 ………………………………………………………………………………………………………………… 试卷编号 拟题教研室（或教师）签名 教研室主任签名 ………………………………………………………………………………………………………………… 课程名称（含档次） 数学物理方程与特殊函数 课程代号 专 业 层次（本、专） 本 科 考试方式（开、闭卷） 闭卷 一．判断题：（本题总分25分，每小题5分） 1．二阶线性偏微分方程062242=+++-y x yy xy xx u u u u u 属于椭圆型； （ ） 2．定解问题的适定性包括解的稳定性、解的唯一性和解的存在性； （ ） 3．如果格林函数),(0M M G 已知，且它在Γ+Ω上具有一阶连续偏导数，又若狄利克雷 问题???=Ω∈=?Γ ).,,(|,),,(0z y x f u z y x u 在Γ+Ω上具有一阶连续偏导数的解存在，那么其解可 表示为=)(0M u dS n G z y x f ??Γ??-) ,,(； （ ） 4．设)(x P n 为n 次Legendre 多项式，则0)()(1 1 1050358?-=dx x P x P ； （ ） 5．设)(x J n 为n 阶Bessel 函数，则 [])()(021ax xJ a ax xJ dx d =． （ ） 二．解答题：（本题总分65分） 1．（本小题15分）设有一根长为l 的均匀细杆，它的表面是绝热的，如果它的端点温度为1),0(u t u =，2),(u t l u =，而初始温度为0T ，写出此定解问题． 2．（本小题20分）利用固有函数法求解下面的定解问题 ???????====><<+=. 0),(,0),0(,0)0,(,0)0,(),0,0(cos sin 2t l u t u x u x u t l x l x t A u a u x x t xx tt πω 其中ω,A 是常数． 3．（本小题15分）求出方程xy u u yy xx =+的一个特解． 第 1 页（共 2 页）

矢量分析与数理方程总复习题

矢量分析与场论，数理方程与特殊函数总复习题 矢量和矢性函数 1、 求下列两个矢量的加法、减法、标量积（点乘）和矢量积（叉乘） k j i A 32++= k j i B 654++= 2、 求下列两个矢性函数的加法、减法、标量积（点乘）和矢量积（叉乘） ()k t j t i t t A ++=sin cos ， ()k t j e i t t B t 2++= 3、设k t j i t A 23+-=，k j i B 22+-=，k j t i C -+=3，求() C B A ?? 4、如果 ()k t j t i t t A ++=sin cos ，()k t j e i t t B t 2++= 求 ()dt t A d 和 ()dt t B d 5、如果 ()j i e ???sin cos += ① 求 ()()? ??d e d e =1 ， ② 证明 ()?e ⊥()?1e . 6、如果 ()j i e ???cos sin 1+-= 证明 ()()?? ?e d e d -=1 7、求不定积分 ()? ??d e ， ()? ??d e 1 。 8、计算不定积分 () ? +???d e 122 . 9、求矢量 k j i r -+=22的单位矢量 0r 。 方向导数和梯度 1、求 k j i l 22++= 的方向余弦 2、写出矢径 k z j y i x r ++=的单位矢径0r ，用方向余弦表示0r 3、求矢性函数 () k z j xy i x z y x l 4232,,+-= 的方向余弦 4、求函数2 2 2 z y x u ++=在() 1,0,1M 处沿k j i l 22++=的方向导数 5、求数量场 z y z x u 2 322+= 在点 () 1,0,2-M 处沿 k z j xy i x l 4232+-= 方向的方向导数 6、求下列数量场的梯度 ① 2 2 2 z y x r ++=， ② ??? ? ? ?++=2 221 1z y x r ， ③ 223z xy z x u +-= ③ 3 2 z y x u =， ④ xz yz xy u ++=， ⑥ z y x xy z y x u 623322 2 2 --++++=.

MATLAB编程求解二维泊松方程doc资料

M A T L A B编程求解二维 泊松方程

%%%% 真解 u=sin(pi*x)*sin(pi*y) %%% %%%% 方程 -Laplace(u)=f %%%%%% %%%% f=2*pi^2*sin(pi*x)*sin(pi*y) %%%%%% %%%%difference code for elliptic equations with constant coefficient %%%%% %clear all %clc N=20; h=1/N; S=h^2; x=0:h:1; y=0:h:1; %%% Stiff matrix A=zeros((N-1)^2,(N-1)^2); for i=1 A(i,i)=4/h^2; A(i,i+1)=-1/h^2; A(i,i+(N-1))=-1/h^2; end for i=N-1 A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,2*i)=-1/h^2; %A(i,i+(N-1))=-1/h^2 end for i=(N-2)*(N-1)+1 A(i,i-(N-1))=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; end for i=(N-1)^2 A(i,i-(N-1))=-1/h^2; A(i,i)=4/h^2; A(i,i-1)=-1/h^2; end for n=2:N-2 i=(N-2)*(N-1)+n; A(i,i-(N-1))=-1/h^2; A(i,i-1)=-1/h^2; A(i,i)=4/h^2; A(i,i+1)=-1/h^2; end

2012、11、10、09年电子科技大学研究生数理方程期末试卷

2012、11、10、09年电子科技大学研究生数理方程期末试卷

电子科技大学研究生试卷 （考试时间： 14点 至 16 点 ，共 2小时） 课程名称 数理方程与特殊函数 教师 学时60 学分 3 教学方式 闭卷 考核日期 2012年 12 月 28 日 成绩 考核方式： （学生填写） 1．把方程 22222320u u u x x y y ???++=????化为标准型，指出其 类型，求出其通解. (10分) 2. 设定解问题：(10分) 2000(),0,0,,0(),(),0. tt xx x x l t t t u a u f x x l t u A u B t u x u x x l ?ψ====?-=<<>?? ==>??==≤≤?? 将该定解问题化成可直接分离变量求解的问题(不需要求出解的具体形式)。 学 号 姓 学 院 教 座位 ……………………密……………封……………线……………以……………

第 1页 3． 长为l 的均匀细杆，其侧面与左端保持零度，右端绝热，杆内初始温度分布为()x ?，求杆内温度分布 (,)u x t . (20分) 4．求下面的定解问题：（10分） 22 009,(,0)18,sin 18 t tt xx t t t u u x e x R t u x x u x ==?-=∈>??=++=+??.

第2页 5．求22 cos()a e x d ?τ??+∞-?.（10分） 6． 222 23()(22)(25) s s F s s s s s ++=++++，求Laplace 逆变换1 (())L F s -.（10分）

天津大学研究生课程-数理方程试题

一． 判断题（每题2分）. 1. 2u u x y x y x ??+=???是非线性偏微分方程.( ) 2. 绝对可积函数一定可做Fourier 积分变化.( ) 3. ()(1) 1.n n F x n Legendre F =是次正交多项式， 则 ( ) 4. (,)0xy f x y =的解是调和函数.( ) 5. **12u u 已知,是线性偏微分方程(,)xx yy u u f x y +=的解，则**12u u -是0u ?= 的解.( ) 二． 填空题（每题2分）. 1. ()sin t xx yy u u u xt -+= 是____________型偏微分方程. 2. 内部无热源的半径为R 的圆形薄板,内部稳态温度分布，当边界上温度为()t φ时，试建立方程的定解问题________________________. 3. 2x 的Legendre 正交多项式的分解形式为__________________. 4.某无界弦做自由振动,此弦的初始位移为()x φ,初始速度为()a x φ-,则弦振动规律为______________________________. 5. []()____________.at m L e t s = 三．求解定解问题（12分） 200sin ; 0,0;0. t xx x x x x l t u a u A t u u u ω===-====

四．用积分变换方法求解以下微分方程（每题12分，共24分） （1） 001,0,0; 1,1. xy x y u x y u y u ===>>=+= （2） 00230, 1.t t t y y y e y y =='''+-='== 五．某半无界弦的端点是自由的，初始位移为零，初始速度为cos x ，求弦的自由振动规律。（12分）

数学物理方程与特殊函数期末考试试题卷子2011

XXXXX 大学研究生试卷 （考试时间： 至 ，共 2小时） 课程名称 数理方程与特殊函数 教师 学时60 学分 3 教学方式 闭卷 考核日期 2011年 12 月 28 日 成绩 1．化方程2220xx xy yy x y x u xyu y u xu yu ++++=为标准形. (10分) 2. 把定解问题：(10分) 212(0)(0,)(),(,)() (,0)(),(,0)(),(0) tt xx x x t u a u x l u t h t u l t h t u x x u x x x l ?ψ?=<

3．有一带状的均匀薄板(0x a ≤≤,0y ≤<+∞), 边界0y =上的温度为0u ，其余边界上的温度保持零度，并且当y →+∞时，温度极限为零. 求解板的稳定温度分布. (用分离 变量法求解).(20分) 4．求下面的定解问题：（10分） 090,(,0) 0,sin tt xx t t t u u x R t u u x ==-=∈>??? ==??. 第2页

5．求()2 1,1 (),()0,1 x x F f x f x x ?-≤?=?>??，其中()F ?表示Fourior 变换.（10分） 6．求()2(),()sin(),03 L f t f t t t π =-≥，其中()L ?为Laplace 变换.（10分） 第3页 学 号 姓 名 学 院 教师 座位号 ……………………密……………封……………线……………以……………内……………答……………题……………无……………效……………………

五点差分格式求解泊松方程的第一边值问题(可编辑)

五点差分格式求解泊松方程的第一边值问题（可编辑）五点差分格式求解泊松方程的第一边值问题 摘要:给出了二维泊松方程在单位正方形上的五点差分格式。并运用线性方程组的古典迭代解法??Jacobi迭代求解出在区域上的数值解。最终绘制数值解的图形。 关键字:泊松方程五点差分格式 Jacobi迭代 有限差分法的介绍 有限差分法是求解偏微分方程的主要数值解法之一;其基本思想是把连续问题离散化,即对求解区域做网格剖分,用有限个网格点代替连续区域;其次将微分算子离散化,从而把微分方程组的问题化为线性方程组的求解问题,解方程组就可以得到原问题在离散点上的近似解。 差分法的步骤:1 对求解域做网格剖分 2 插值函数的选择 3 方程组的建立 4 方程组的求解 五点差分格式的构造 二维泊松方程: 在单位正方形上,在正方形边界上的边界条件.在正方形网格上,就是在上离散化,.对于N3如图1所示: 图1 沿方向分别用二阶中心差商代替 2.1 2.2

1、2式相加可得差分方程: 2.3 利用Taylor展式 可得差分算子的截断误差 其中是方程2.3的光滑解。 由于差分方程2.3中只出现在及其四个邻点上的值见图1的中间的粗的点,所以称为五点差分格式。 由边界条件知道,因而2.3式确定了一组具有个未知量的个线性方程。对应的系数矩阵为对称、不可约对角占优,且对角元为正,因而系数矩阵非奇异,且为对称正定阵。 三、方程组的求解 我们已经知道,利用差分方法解椭圆型方程边值问题归结为解大型线性代数方程组的问题。因为差分格式产生的大型线性代数方程组的系数矩阵中非零元素占的比例小,分布很有规律性。而且通过数值线性代数的学习,我们知道对于大型的稀疏矩阵来说,迭代法是比较好的选择,其程序实现比较简单,迭代过程能自动校正计算过程中的偶然误差,要求计算机的存储相对较少。 本文采用了线性方程组古典迭代解法??Jacobi迭代求解由五点差分格式得到的线性方程组。以下对Jacobi迭代作简要的介绍: 给定3.1 令3.2 其中 3.3 那么3.1可以写成,3.4 其中.若给定初始向量,并代入3.4的右端,就可以计算出一个新的向量,即,

数理方程期末试题B答案

数理方程期末试题B答 案 集团文件版本号：（M928-T898-M248-WU2669-I2896-DQ586-M1988）

北 京 交 通 大 学 2007-2008学年第二学期《数理方程与特殊函数》期末考试试卷 （B ） （参考答案） 学院_ ____________ 专业___________________ 班级________ ____ 学号_______________ 姓名___________ __ 一、 计算题（共80分，每题16分） 1． 求下列定解问题（15分） 2． 用积分变换法及性质，求解半无界弦的自由振动问题：（15分） 3．设弦的两端固定于0x =及x l =，弦的出示位移如下图所示。初速度为 零，又没有外力作用。求弦做横向振动时的位移(,)u x t 。 [ 解 ] 问题的定解条件是 由初始条件可得 4． 证明在变换, x at x at ξη=-=+下，波动方程xx tt u a u 2=具有形式解0=n u ξ， 并由此求出波动方程的通解。 5． 用分离变量法解下列定解问题

[ 提示：1) 可以直接给出问题的固有函数，不必推导；2) 利用参数变易法。] [ 解 ] 对应齐次方程的定解问题的固有函数是x l n πsin ，其解可以表示成 把原问题中非齐次项t x t x f l a l π π22sin sin ),(=按照固有函数展开成级数 因此有 利用参数变易法，有 于是 6． 用Bessel 函数法求解下面定解问题 [ 解 ] 用分离变量法求解。令)()(),(t T R t u ρρ=，则可得 以及 设0ρβλn n =为Bessel 函数)(0x J 的正零点，则问题（II ）的特征值和特征函数分别为 问题（I ）的解为 于是原问题的解是 由初始条件 得到 故 于是最后得到原问题的解是 二、 证明题（共2分，每题10分） 7．证明平面上的Green 公式 其中C 是区域D 的边界曲线，ds 是弧长微分。

二维泊松方程很基础详细的求解过程

Topic 2: Elliptic Partial Di?erential Equations Lecture 2-4: Poisson’s Equation: Multigrid Methods Wednesday, February 3, 2010 Contents 1 Multigrid Methods 2 Multigrid method for Poisson’s equation in 2-D 3 Simple V ?cycle algorithm 4 Restricting the Residual to a Coarser Lattice 2 3 5 7

1 MULTIGRID METHODS 5 Prolongation of the Correction to the Finer Lattice 6 Cell-centered and Vertex-centered Grids and Coarsenings 7 Boundary points 8 Restriction and Prolongation Operators 9 Improvements and More Complicated Multigrid Algorithms 8 8 11 11 15 1 Multigrid Methods The multigrid method provides algorithms which can be used to accelerate the rate of convergence of iterative methods, such as Jacobi or Gauss-Seidel, for solving elliptic partial di?erential equations. Iterative methods start with an approximate guess for the solution to the di?erential equation. In each iteration, the di?erence between the approximate solution and the exact solution is made smaller. One can analyze this di?erence or error into components of di?erent wavelengths, for example by using Fourier analysis. In general the error will have components of many di?erent wavelengths: there will be

数理方程第二版_课后习题答案

第一章曲线论 §1 向量函数 1. 证明本节命题3、命题5中未加证明的结论。 略 2. 求证常向量的微商等于零向量。 证：设，为常向量，因为 所以。证毕3. 证明 证： 证毕4. 利用向量函数的泰勒公式证明：如果向量在某一区间所有的点其微商为零，则此向量在该区间上是常向量。 证：设，为定义在区间上的向量函数，因为在区间上可导当且仅当数量函数，和在区间上可导。所以，

，根据数量函数的Lagrange中值定理，有 其中，，介于与之间。从而 上式为向量函数的0阶Taylor公式，其中。如果在区间上处处有，则在区间上处处有 ，从而，于是。证毕 5. 证明具有固定方向的充要条件是。 证：必要性：设具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，于是。 充分性：如果，可设，令，其中为某个数量函数，为单位向量，因为，于是 因为，故，从而

为常向量，于是，，即具有固定方向。证毕 6. 证明平行于固定平面的充要条件是。 证：必要性：设平行于固定平面，则存在一个常向量，使得，对此式连续求导，依次可得和，从而，，和共面，因此。 充分性：设，即，其中，如果，根据第5题的结论知，具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，任取一个与垂直的单位常向量，于是作以为法向量过原点的平面，则平行于。如果，则与不共线，又由可知，，，和共面，于是， 其中，为数量函数，令，那么，这说明与共线，从而，根据第5题的结论知，具有固定方向，则可表示为，其中为某个数量函数，为单位常向量，作以为法向量，过原点的平面，则平行于。证毕 §2曲线的概念 1. 求圆柱螺线在点的切线与法平面的方程。 解：，点对应于参数，于是当时，，，于是切线的方程为：

数学物理方法期末考试试题典型汇总

一、 Mathematical methods for physics 二、 单项选择题（每小题2分） 1．齐次边界条件0),(),0(==t u t u x x π的本征函数是_______。 A) 3,2,1 sin =n nx B) ,2,1,0 cos =n nx C) 2,1,0 )2 1 sin(=+n x n D) 2,1,0 )2 1 cos(=+n x n 2．描述无源空间静电势满足的方程是________。 A) 波动方程 B)热传导方程 C) Poisson 方程 D)Laplace 方程 3．半径为R 的圆形膜，边缘固定，其定解问题是??? ? ? ????====?-??===)(| ),(|0|0),(),(0t 02 22 2ρψρ?ρρρt t R u u u t u a t t u 其解的形式为∑∞ == 1 0)()(),(m m m k J t T t u ρρ，下列哪一个结论是错误的______。 A) )()()()(2 2 2 2t T k a t T dt d t T m m m m -=满足方程 B ）圆形膜固有振动模式是)sin(0t ak m 和)cos(0 t ak m C ）0 m k 是零阶Bessel 函数的第m 个零点。 D ）)()(00ρρm m k J R =满足方程0)(2 202=+'+''R k R R m ρρρ 4．)(5x P 是下列哪一个方程的解_________。 A ）0202)1(2=+'-''-y y x y x B ）0252)1(2=+'-''-y y x y x C ）0302)1(2=+'-''-y y x y x D ）052)1(2=+'-''-y y x y x 5．根据整数阶Bessel 函数的递推公式，下列结论哪一个是正确的________。 A ）)(2)()(120x J x J x J '=- B ）)()()(111x J x x J x xJ '=+

研究生数理方程与特殊函数考题2014

科技大学研究生试卷 （考试时间： 至 ，共 2小时） 课程名称 数理方程与特殊函数 教师 学时60 学分 3 教学方式 闭卷 考核日期 2014年 12 月 日 成绩 考核方式： （学生填写） 1．化简方程22222 (,)(,)(,) 1280u x y u x y u x y x x y y ???++=????并求其通解. (10分) 2. 设有一长度为L 的均匀细棒，其侧面和两端均绝热，初始温度分布为已知。(1)求以后时刻的温度分布；(2)证明：当初始温度分布为常数时，以后时刻的温度分布也必为常数. (20分) 第 1页 3．求解定解问题：(15分) 学 号 姓 名 学 院 教师 座位号 ……………………密……………封……………线……………以……………内……………答……………题……………无……………效……………………

200000 (0,0),t xx x x l t u a u x l t q u u u k u u ===?=<<>? ? ==?? ?=?,00,,,a u k q 均为常数. 4．求函数()() 2 1 ()13f s s s =+- 的Laplace 逆变换.（10分） 第2页 5．求下面的定解问题：（15分） 号 效……………………

2 00,(,0) ,sin tt xx t t t u a u x at x R t u x u x ==?-=+∈>?? ==??. 6．求3()J x dx ? .（10分） 第3页 7．写出平面第一象限的Dirichlets 问题对应的Green 函数及其定解问题.（10分）