偏微分方程期末试题A卷

偏微分方程期末试题A卷
偏微分方程期末试题A卷

安徽大学20 08 —20 09 学年第 二 学期 《 偏微分方程 》考试试卷(A 卷)

(闭卷 时间120分钟)

院/系 年级 专业 姓名 学号

一、填空题(每小题3分,共15分)

1.对常系数方程x y z u au bu cu du f ?++++=作未知函数的变换 可以将所有一阶微商消失.

2.设:R R Φ→是光滑凸函数,(,)u x t 是热传导放程0t u u -?=的解,则()u Φ是热传导方程 的 (下解;上解;解).

3.上半平面的Green 函数G(x,y)为 ,其中12(,)y y y =为上半平面中某固定点.

4.设函数u 在以曲面Γ为边界的区域Ω内调和,在ΩΓ 上有连续的一阶偏导数,则u

dS n Γ

????= ,其中n 是Γ的外法方向.

5.热传导方程2()0t xx yy u a u u -+=的特征曲面为 .

二、计算题(每小题10分,共40分) 1.求解初值问题

0,(,)(0,)(,0),,t x u bu cu x t R u x g x R ++=∈?∞??=∈? 其中,,b c R ∈都是常数.

2.试用延拓法求解半有界直线上的热传导方程的边值问题:

200

0,0,0,|(),

|0.t xx t x u a u x t u x u ?==?-=>>?

=??=?

3.试求解

2

2

008(),

|,|.tt xx yy zz t t t u u u u t u xy u z ==?-++=??==??

4.写出定解问题:

200

(),0,0,|0,|0,

|().t xx x x l t u a u f x x l t u u u g x ===?-=<<>?

==??=?

解的一般形式.

三、判断分析题(10分)

试判断下面命题是否成立,并说明原因.

在证明Hopf 引理的过程中,我们能够作出一个辅助函数()v x 满足 (a)在球面()R B y ?上0;v =

(b)v 沿球()R B y 的半径方向的方向导数v

ν

??<0;

(c)在整个球()R B y 内下调和.

四、分析计算题(15分)

试判断下列方程

2222

222sin cos cos 0u u u u x x x x x y y y

????---=?????

的类型,并根据标准型求出此方程的通解.

五、证明题(下面两道题请任选一题)(20分)

1.设G 是2R

满足方程0xx yy u u +=的函数u(x,y)在G 上的最大值不会超过它在边界G ?上的最大值.

2.试用能量法(即用格林第一公式法)证明n 维Laplace 方程的第三边值问题

12n u(x)0,x=(x ,x ,,x ),0u u f n σσ?Ω?=∈Ω

??

????

+=> ????

?? 是常数

的解的唯一性,其中Ω为边界光滑的有界区域.

偏微分方程数值解期末试题及标准答案

偏微分方程数值解试题(06B ) 参考答案与评分标准 信息与计算科学专业 一(10分)、设矩阵A 对称,定义)(),(),(2 1)(n R x x b x Ax x J ∈-=,)()(0x x J λλ?+=.若0)0('=?,则称称0x 是)(x J 的驻点(或稳定点).矩阵A 对称(不必正定),求证0x 是)(x J 的驻点的充要条件是:0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点,对于任意的n R x ∈,令 ),(2),()()()(2 000x Ax x b Ax x J x x J λλλλ?+-+=+=, (3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若n R x ∈0满足b Ax =0,则对于任意的x ,)(),(2 1)0()1()(00x J x Ax x x J >+==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的展开式3分, 每问3分,推理逻辑性1分 二(10分)、 对于两点边值问题:?????==∈=+-=0 )(,0)(),()('b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ],[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式:求泛函极小的Ritz 形式和Galerkin 形式的变分方程。 解: 设}0)(),,(|{11=∈=a u b a H u u H E 为求解函数空间,检验函数空间.取),(1b a H v E ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. (3分)

偏微分方程数值解复习题(2011硕士)

偏微分方程数值解期末复习(2011硕士) 一、考题类型 本次试卷共六道题目,题型及其所占比例分别为: 填空题20%;计算题80% 二、按章节复习内容 第一章 知识点:Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等; 要求: 会辨认差分格式, 判断线性多步法的误差和阶; 第二章 知识点:矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和、稳定性等; 要求: 建立椭圆型方程边值问题的差分格式, 极值原理; 第四章 知识点:最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和稳定性等; 要求: 建立抛物型方程边值问题的差分格式, 计算局部截断误差; 第五章 知识点:左偏心格式、右偏心格式、中心格式、LF格式、LW格式、Wendroff 格式、跳蛙格式、特征线、CFL条件等; 要求: 建立双曲型方程边值问题的差分格式, 计算局部截断误差; 第七章 要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式

三 练习题 1、 已知显格式21131()22 n n n n u u h f f +++-=-,试证明格式是相容的,并求它的阶。 P39+P41 2、用Taylor 展开原理构造一元函数一阶导数和二阶导数的数值微分公式。 提示:向前、向后和中心差商与一阶导数间关系,二阶中心差商与二阶导数 之间的关系 课件 3、用数值微分方法或数值积分方法建立椭圆型方程 2222(,),(,),u u f x y x y x y ??--=?∈Ω?? :01,01x y Ω≤≤≤≤ 内点差分格式。 P75+课件 4、构造椭圆型方程边值问题的差分格式. P101 (4)题 5、构建一维热传导方程220,(0)u u Lu a a t x ??=-=>??的数值差分格式(显隐格式等)。 参考P132-135相关知识点 6、设有逼近热传导方程22(0)u u Lu a f a const t x ??≡-==>??的带权双层格式 ()()1111111122(1)2k k j j k k k k k k j j j j j j u u a u u u u u u h θθτ++++-+-+-??=-++--+?? 其中[0,1]θ∈,试求其截断误差。并证明当2 1212h a θτ=-时,截断误差的阶最 高阶为24()O h τ+。 P135+P165+课件 7、传播因子法证明抛物型方程22(0)u u Lu a f a const t x ??≡-==>??的最简显隐和六点CN 格式稳定性。 P156+课件 8、对一阶常系数双曲型方程的初边值问题 0,0,0,0,(,0)(),0,(0,)(),0, u u a t T x a t x u x x x u t t t T φψ???+=<≤<<∞>?????=≤<∞??=≤≤?

(完整版)偏微分方程的MATLAB解法

引言 偏微分方程定解问题有着广泛的应用背景。人们用偏微分方程来描述、解释或者预见各种自然现象,并用于科学和工程技术的各个领域fll。然而,对于广大应用工作者来说,从偏微分方程模型出发,使用有限元法或有限差分法求解都要耗费很大的工作量,才能得到数值解。现在,MATLAB PDEToolbox已实现对于空间二维问题高速、准确的求解过程。 偏微分方程 如果一个微分方程中出现的未知函数只含一个自变量,这个方程叫做常微分方程,也简称微分方程;如果一个微分方程中出现多元函数的偏导数,或者说如果未知函数和几个变量有关,而且方程中出现未知函数对几个变量的导数,那么这种微分方程就是偏微分方程。 常用的方法有变分法和有限差分法。变分法是把定解问题转化成变分问题,再求变分问题的近似解;有限差分法是把定解问题转化成代数方程,然后用计算机进行计算;还有一种更有意义的模拟法,它用另一个物理的问题实验研究来代替所研究某个物理问题的定解。虽然物理现象本质不同,但是抽象地表示在数学上是同一个定解问题,如研究某个不规则形状的物体里的稳定温度分布问题,由于求解比较困难,可作相应的静电场或稳恒电流场实验研究,测定场中各处的电势,从而也解决了所研究的稳定温度场中的温度分布问题。 随着物理科学所研究的现象在广度和深度两方面的扩展,偏微分方程的应用范围更广泛。从数学自身的角度看,偏微分方程的求解促使数学在函数论、变分法、级数展开、常微分方程、代数、微分几何等各方面进行发展。从这个角度说,偏微分方程变成了数学的中心。

一、MATLAB方法简介及应用 1.1 MATLAB简介 MATLAB是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。 1.2 Matlab主要功能 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 数字信号处理 通讯系统设计与仿真 财务与金融工程 1.3 优势特点 1) 高效的数值计算及符号计算功能,能使用户从繁杂的数学运算分析中解脱出来; 2) 具有完备的图形处理功能,实现计算结果和编程的可视化; 3) 友好的用户界面及接近数学表达式的自然化语言,使学者易于学习和掌握; 4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ,

偏微分方程数值解期末试题及答案(内容参考)

偏微分方程数值解试题(06B) 参考答案与评分标准 信息与计算科学专业 一(10分)、设矩阵A 对称,定义)(),(),(2 1 )(n R x x b x Ax x J ∈-= ,)()(0x x J λλ?+=.若0)0('=?,则称称0x 是)(x J 的驻点(或稳定点).矩阵A 对称(不必正定),求证0x 是)(x J 的驻点的充要条件是:0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点,对于任意的n R x ∈,令 ),(2 ),()()()(2 000x Ax x b Ax x J x x J λλλλ?+ -+=+=, (3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有 0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若 n R x ∈0满足 b Ax =0,则对于任意的 x ,)(),(2 1 )0()1()(00x J x Ax x x J >+ ==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的展开式3分, 每问3分,推理逻辑性1分 二(10分)、 对于两点边值问题:????? ==∈=+-=0 )(,0)() ,()(' b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ] ,[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式:求泛函极小的Ritz 形式和 Galerkin 形式的变分方程。 解: 设}0)(),,(|{11 =∈=a u b a H u u H E 为求解函数空间,检验函数空间.取),(1 b a H v E ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. (3分)

偏微分方程数值解法

一、 问题 用有限元方法求下面方程的数值解 2 u u u f t ?-?+=? in (]0,T Ω? 0u = on []0,T ?Ω? ()00,u x u = in Ω 二、 问题分析 第一步 利用Green 公式,求出方程的变分形式 变分形式为:求()()21 00,;u L T H ∈Ω,使得 ()())(2 ,,,,u v u v u v f v t ???+??+= ???? ()10v H ?∈Ω (*) 以及 ()00,u x u =. 第二步 对空间进行离散,得出半离散格式 对区域Ω进行剖分,构造节点基函数,得出有限元子空间:()12,,,h NG V span ???=???,则(*)的Galerkin 逼近为: []0,t T ?∈,求()()1 0,h h u t x V H ∈?Ω,使得 ()()()()() () )(2 ,,,,h h h h h h h d u t v u t v u t v f v dt +??+= h h v V ?∈ (**) 以及()0,0h h u u =,0,h u 为初始条件0u 在h V 中的逼近,设0,h u 为0u 在h V 中的插值. 则0t ?≥,有()()1 N G h i i i u t t ξ? == ∑,0,h u =01 N G i i i ξ?=∑,代人(**)即可得到一常微分方程组. 第三步 进一步对时间进行离散,得到全离散的逼近格式 对 du dt 用差分格式.为此把[]0,T 等分为n 个小区间[]1,i i t t -,其长度1i i T t t t n -?=-= ,n t T =. 这样把求i t 时刻的近似记为i h u ,0 h u 是0u 的近似.这里对(**)采用向后的欧拉格式,即 ()()() () )(2 11 11 1 ,,,,i i i i h h h h h h h i h u u v u v u v f v t ++++-+??+ = ? h h v V ?∈ (***) i=0,1,2…,n-1. 0 h u =0,h u 由于向后欧拉格式为隐式格式且含有非线性项,故相邻两时间步之间采用牛顿迭代,即:

北京理工大学数学专业偏微分方程期末试题2014级A卷(MTH17178)

课程编号:MTH17178 北京理工大学2016-2017学年第一学期 2014级偏微分方程期终考试(A ) 1.(10分)利用特征线方法求解一阶波动方程初值问题:()22,,0,0,t x x u u u x t u x e x -+=∈>???=∈?? 。 2.(10分)利用Fourier 变换方法求解:()() (),,,0,0,t x u bu cu f x t x t u x x x ?--=∈>???=∈?? 。 3.(10分)利用行波法求解:()()()()0,,,0,,0 tt xx u u t x u x x x x u x x x x ?ψ?-=>?-=?。 给出适当的相容性条件。如果?在(],0a -上给定,ψ在[)0,b 上给定,给出其决定区域。 4.(15分)求解初边值问题:()()()20,01,00,0,1,0,0,0,01 t xx x x u a u u x t u t u t t u x A x ?-+=<<>?==>??=<?==∈??=+=≥? 推导边界条件齐次化的公式(不需要解方程)。 6.(13分)对于有界区域()(],0,T Q a b T =?上的热方程()2 ,0t xx u a u c x t u -+=,其中(),c x t 下有界,证明如果(),u x t 在抛物边界上非正,则(),u x t 在T Q 上非正。 7.(15分)考虑波动方程初边值问题[]()()()()[]()()()20,0,,0,0,,0,0,0,0,,,0,0 tt xx t x x u a u x L t u x x u x x x L u t u L t u L t t ?ψσ?-=∈>?==∈??=+=≥?,其中 0σ>,令t 时刻的能量()()()22222011,22 L t x E t u a u dx a u L t σ=++?,证明()E t 守恒,并由此证明相应的一般非齐次方程非齐次初边值问题的解的唯一性。 8.(20分)设() ()1,02,1T T u C Q C Q ∈ 且满足初边值问题()()()()[]()()[] ,,,,0,0,0,,0,0,t xx T x u u f x t x t Q u x x x L u t u L t t T ??-=∈?=∈??==∈?,证明:[]()()()()22220000000,sup ,,,L T L L T L x t T u x t dx dt u x t dx M x dx dt f x t dx ?∈??+≤+??????????,其中M 仅依赖于T 。 提示:Gronwall 不等式:设(][]1 0,0,G C T C T ∈ ,()00G =,且对于任意的[]0,t T ∈,有()()()G t CG t F t '≤+,其中C>0,F 非负单调递增,则有 ()()()()()11,Ct Ct G t C e F t G t e F t -'≤-≤。

偏微分方程期末试题A卷

安徽大学20 08 —20 09 学年第 二 学期 《 偏微分方程 》考试试卷(A 卷) (闭卷 时间120分钟) 院/系 年级 专业 姓名 学号 一、填空题(每小题3分,共15分) 1.对常系数方程x y z u au bu cu du f ?++++=作未知函数的变换 可以将所有一阶微商消失. 2.设:R R Φ→是光滑凸函数,(,)u x t 是热传导放程0t u u -?=的解,则()u Φ是热传导方程 的 (下解;上解;解). 3.上半平面的Green 函数G(x,y)为 ,其中12(,)y y y =为上半平面中某固定点. 4.设函数u 在以曲面Γ为边界的区域Ω内调和,在ΩΓ 上有连续的一阶偏导数,则u dS n Γ ????= ,其中n 是Γ的外法方向. 5.热传导方程2()0t xx yy u a u u -+=的特征曲面为 .

二、计算题(每小题10分,共40分) 1.求解初值问题 0,(,)(0,)(,0),,t x u bu cu x t R u x g x R ++=∈?∞??=∈? 其中,,b c R ∈都是常数. 2.试用延拓法求解半有界直线上的热传导方程的边值问题: 200 0,0,0,|(), |0.t xx t x u a u x t u x u ?==?-=>>? =??=?

3.试求解 2 2 008(), |,|.tt xx yy zz t t t u u u u t u xy u z ==?-++=??==?? 4.写出定解问题: 200 (),0,0,|0,|0, |().t xx x x l t u a u f x x l t u u u g x ===?-=<<>? ==??=? 解的一般形式.

微分方程几种求解方法

第五章 控制系统仿真 §5.2 微分方程求解方法 以一个自由振动系统实例为例进行讨论。 如下图1所示弹簧-阻尼系统,参数如下: M=5 kg, b=1 N.s/m, k=2 N/m, F=1N F 图1 弹簧-阻尼系统 假设初始条件为:00=t 时,将m 拉向右方,忽略小车的摩擦阻力,m x 0)0(= s m x /0)0(=? 求系统的响应。 )用常微分方程的数值求解函数求解包括ode45、 ode23、ode113、ode15s 、ode23s 等。 wffc1.m myfun1.m 一、常微分方程的数值求解函数ode45求解 解:系统方程为 F kx x b x m =++??? 这是一个单变量二阶常微分方程。

将上式写成一个一阶方程组的形式,这是函数ode45调用规定的格式。 令: x x =)1( (位移) )1()2(? ?==x x x (速度) 上式可表示成: ??????--=??????=??? ???????)1(*4.0)2(*2.02.0)2()2()2()1(x x x x x x x && 下面就可以进行程序的编制。 %写出函数文件myfun1.m function xdot=myfun1(t,x) xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)]; % 主程序wffc1.m t=[0 30]; x0=[0;0]; [tt,yy]=ode45(@myfun1,t,x0); plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') hold on plot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-k') legend('位移','速度',’加速度’)

最新偏微分方程期末复习笔记

《偏微分方程》期末考试复习 一、波动方程(双曲型方程)U tt -a 2U xx 二f (x,t) (一)初值问题(柯西问题) < 2 U tt —a U xx = f(x,t) 1、一维情形 Ut t^a (x) (1) 解法(传播波法): 由叠加原理,原初值问题的解可表示为下述初值问题的解之 和, * 2 * 2 U tt —a U xx =o U tt —a U xx = f (x,t) (i) J U t^=

②决定区域:区间[x1,X2】的决定区域为:{(x,t)|捲? at込x込X2-at}

偏微分方程期末考试试题(06)

课程名称:偏微分方程数值解法 课程编号:24014110 适用专业(班级):数学 共1页 命题人:潘晓丽 教研室主任: 第1页 一、(15分)写出三类典型泛定方程并分别说明其名称和特点. 二、(10分)求一维波动方程()()()()()22 222 ,,0,0,,0t u u a x t t x u x x u x x ?ψ???=-∞<<+∞>?????==? 的通解. 三、(15分)写出达朗贝尔公式并利用公式求解 ()()()2,0,,0sin ,0cos tt xx t u a u t x u x x u x x ?=>-∞<<+∞? =?? =? 四、(10分)计算积分()32x J x dx -?. 五、(15分)设1,1≥≥n m ,证明 ()()()dx x p x m dx x p x n m n m n m ??--=++1 111 1 六、(15分)用分离变量法求解 ()()()()()20,0,0,00,,00,0,,0 tt xx t u a u x l t u x u x x u t u l t ?-=<<>? ==?? ==? 七、(10分)解固有值问题()()()''0,''0 y y l x l y l y l λ+=-<

课程名称:偏微分方程数值解法 课程编号:24014110 适用专业(班级):数学 共1页 命题人:潘晓丽 教研室主任: 第1页 一、解:波动方程:()22 2,u a u f t x t ?=?+? 热传导方程: ()2,u a u f t x t ?=?+? 位势方程:()u f x ?= ……………………….5分 其中()12,,,n x x x x = ,a 为常数,(),f t x 及()f x 为已知函数,在波动方程及 热传导方程中,未知函数u 是时间变量t 和空间坐标变量()12,,,n x x x x = 的函数,在位势方程中,未知函数u 是空间坐标变量()12,,,n x x x x = 的函数,而与时间t 无关,三类典型方程均为二阶线性偏微分方程。……………………….15分 二、解:首先判别方程的类型, 20a ?=> ………………………2分 即此方程在整个全平面上都是双曲型的。 特征方程为:()()2 2 20dx a dt -= () ()2 2 200dx a dt dx adt -=?= 特征曲线为1 2 x at c x at c -=??+=? ………………………6分 做变量替换,令x at x at ξη=-??=+?, 由链式法则得 0u ξη= 通解()()()()u f g f x at g x at ξη=+=-++ ……………………….10分

偏微分方程期末考试试题(06)

黑龙江科技学院考试试题 课程名称:偏微分方程数值解法 课程编号:24014110适用专业(班级):数学 命题人:潘晓丽 教研室主任: 、(15分)写出三类典型泛定方程并分别说明其名称和特点 2 2 U 2 U 一、(10分)求一维波动方程 t 2 x 2 ,t 0 的通解 x u x,0 x , u t x,0 三、(15 分) 写出达朗贝尔公式并利用公式求解 u tt a 2 u xx , t 0, x u x,0 sinx U t x,0 cosx 四、(10分)计算积分 x 3 J 2 x dx . 五、(15分)设m 1,n 1,证明 六、(15分)用分离变量法求解 2 u tt a U xx 0, 0 x l,t 0 u x,0 0,u t x,0 x u 0,t 0,u l,t 0 八、(10分)叙述斯图模-刘维尔定理. 黑龙江科技学院考试试题答案 七、(10分)解固有值问题 y'' y 0, y' l y' l 第一套 共1页 第1页 n 1 0x m p n xdx 1 m 1 , m 0 x p n 1 x dx

2 一、解:波动方程:一a2u f t,x t - 热传导方程:汁a2 u f t,x 位势方程:u f x (5) 其中x X j,x2,L ,x n,a为常数,f t,x及f x为已知函数,在波动方程及热传导方程中,未知函数u是时间变量t和空间坐标变量x x1,x2,L ,x n的函数,在位势方程中,未知函数u是空间坐标变量x 为必,L ,人的函数,而与时间t无关,三类典型方程均为二阶线性偏微分方程。 (15) 二、解:首先判别方程的类型, a20 ............. 2 分 即此方程在整个全平面上都是双曲型的。 特征方程为:dx $ a2 dt $ 0 2 2 2 dx a dt 0 dx madt 0 x at 特征曲线为G x at C2 做变量替换,令 x at x at 由链式法则得u 0 通解u f g f x at g x at ....................... .10 ................................ 分

金融工程期末复习题

一、简述题(30分) 1.金融工程包括哪些主要内容? 答:产品与解决方案设计,准确定价与风险管理是金融工程的主要内容P3 2.金融工程的工具都有哪些? 答:基础证券(主要包括股票和债券)和金融衍生产品(远期,期货,互换和期权)P4 3.无套利定价方法有哪些主要特征? 答:a.套利活动在无风险的状态下进行 b.无套利的关键技术是“复制”技术 c.无风险的套利活动从初始现金流看是零投资组合,即开始时套利者不需要任何资金的 投入,在投资期间也不需要任何的维持成本。P16 4.衍生证券定价的基本假设为何? 答:(1)市场不存在摩擦 (2)市场参与者不承担对手风险 (3)市场是完全竞争的 (4)市场参与者厌恶风险,且希望财富越多越好 (5)市场不存在无风险套利机会P20 5.请解释远期与期货的基本区别。 答:a.交易场所不同 b.标准化程度不同 c.违约风险不同 d.合约双方关系不同 e.价格确定方式不同 f.结算方式不同 g.结清方式不同P44 6.金融互换的主要有哪些种类? 答:利率互换与货币互换和其它互换(交叉货币利率互换、基点互换、零息互换、后期确定互换、差额互换、远期互换、股票互换等等)P104 7.二叉树定价方法的基本原理是什么? 答:二叉树图方法用离散的模型模拟资产价格的连续运动,利用均值和方差匹配来确定相关参数,然后从二叉树图的末端开始倒推可以计算出期权价格。P214 8.简要说明股票期权与权证的差别。 答:股本权证与备兑权证的差别主要在于: (1)有无发行环节; (2)有无数量限制; (3)是否影响总股本。 股票期权与股本权证的区别主要在于: (1)有无发行环节 (2)有无数量限制。P162 9.影响期权价格的因素主要有哪些?它们对欧式看涨期权有何影响? 答: 1)标的资产的市场价格(+) 2)期权的协议价格(—) 3)期权的有效期(?) 4)标的资产价格的波动率(+) 5)无风险利率(+) 6)标的资产收益(—) “+”表示对欧式看涨期权正向的影响,“—”表示反向的影响,“?”表示不确定P175 10.蒙特卡罗模拟法的主要优缺点。 答:优点:A.在大多数情况下,人们可以很直接地应用蒙特卡罗模拟法,而无需对期权定价模型有深刻的理解,所用的数学知识也很基本 B.为了获得更精确的答案,只需要进行更多的模拟 C.无需太多工作就可以转换模型。 缺点:A.难以处理提前执行的情形,因此难以为美式期权定价 B.为了达到一定的精确度,一般需要大量的模拟运算P226 11.用蒙特卡罗法确定期权价格的基本过程是什么? 答:由于大部分期权价值等于期权到期回报的期望值的贴现,因此先模拟风险中性世界中标的

偏微分方程求解方法及其比较

偏微分方程求解方法及其比较 发表时间:2008-12-11T09:32:01.530Z 来源:《科海故事博览科教创新》2008年第10期供稿作者:曹海洋吕淑娟王淑芬 [导读] 近些年来,无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高,使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法,它具有无穷阶收敛性.因此,谱方法也就引起人们更多的关注. 摘要:近些年来,无限维动力系统得到了很大的发展.随着对它研究的深入和计算能力的迅速提高,使得与之相关的数值研究越来越被人们关注.谱方法作为一种数值求解偏微分方程的方法,它具有无穷阶收敛性.因此,谱方法也就引起人们更多的关注. 关键词:谱方法;偏微分;收敛;逼近; 1偏微分方程及其谱方法的介绍 偏微分方程主要借助于未知函数及其导数来刻画客观世界的物理量的一般变化规律。理论上,对偏微分方程解法的研究已经有很长的历史了。最初的研究工作主要集中在物理,力学,几何学等方面的具体问题,其经典代表是波动方程,热传导方程和位势方程(调和方程)。通过对这些问题的研究,形成了至今仍然使用的有效方法,例如,分离变量法,fourier变换法等。早期的偏微分方程研究主要集中在理论上,而在实际操作中其研究方法和研究结果都难以得到广泛的应用。求解的主要方法为:有限差分法,有限元法,谱方法。 谱方法起源于Ritz-Galerkin方法,它是以正交多项式(三角多项式,切比雪夫多项式,勒让得多项式等)作为基函数的Galerkin方法、Tau 方法或配置法,它们分别称为谱方法、Tau方法或拟谱方法(配点法),通称为谱方法。谱方法是以正交函数或固有函数为近似函数的计算方法。从函数近似角度看.谱方法可分为Fourier方法.Chebyshev或Legendre方法。前者适用于周期性问题,后两者适用于非周期性问题。而这些方法的基础就是建立空间基函数。 下面介绍几种正交多项式各种节点的取值方法及权重。 1) Chebyshev-Gauss: 2) Chebyshev-Gauss-Radau: x0 =1, 3) Chebyshev-Gauss-Lobatto: x0 =1, xN =1, 4)Legendre-Gauss: xj 是的零点且 5) Legendre-Gauss-Radau: xj 是的N+1个零点且 6) Legendre-Gauss-Lobatto: x0=-1,xN=1其它N-1个点是的零点且 下面介绍谱方法中最重要的Jacobi正交多项式其迭代公式为: 其中: Jacobi正交多项式满足正交性: 而Chebyshev多项式是令时Jacobi多项式的特殊形式,另外Legendre多项式是令时Jacobi多项式的特殊形式。 2 几种典型的谱方法 谱方法是以正交函数或固有函数为近似函数的计算方法。谱近似可以分为函数近似和方程近似两种近似方式。从函数近似角度看.谱方法可分为Fourier方法.Chebyshev或Legendre方法。前者适用于周期性问题,后两者适用于非周期性问题。从方程近似角度看,谱方法可分为在物理空间离散求解的Collocation法、在谱空间进行离散求解的Galerkin法,以及先在物理空间离散求积,再变换到谱空间求解的Pseudo-spectral法。Collocation法适用于非线性问题.Galerkin法适用于线性问题,而Pseudo-spectral法适用于展开方程时的非线性项的处理。谱方法的特点是对光滑函数指数性逼近的谱精度;以较少的网格点得到较高的精度;无相位误差;适合多尺度的波动性问题;计算精度高于其他方法。快速傅立叶变化的提出大大促进了谱方法的发展,迄今已有各种的谱方法计算格式被提出.并被应用于天文学、电磁学、地理学等各种问题的计算。 下面介绍一下应用于各个区域的几种谱方法: 1)以Fourier谱方法为例介绍谱方法解方程的主要过程 以一阶波动方程为例: 其中u(x,t)为方程的解,L是包含u和u关于空间变量的导数的算子,除了方程以有初始条件和适当的边界条件。 故可设其中为试探空间的基函数,ak(t)为展开系数,对于傅立叶谱方法中的共轭有: 其中从而利用其正交性和周期性可以减少工作量,另外再结合边界条件就可以求出来。 2) Galerkin方法是谱方法中十分经典的解偏微分方程的方法,但还有其局限性,而利用Hermite谱方法中依赖时间的权函数对经典的Galerkin方法进行拓展后的新的方法能适用范围扩大了很多。它能很好的应用在微分方程最优控制问题有限元方法的分析中,并且如果能够灵活运用利用Chebyshev方法、Galerkin方法和配置方法,则会形成更强的计算方法。如将Tau方法的思想成功地应用于奇数阶微分方程Petrov-Galerkin谱方法。 3)在无界区域上谱方法和拟谱方法发展了以Hermite函数和Laguerre函数为基函数的正交逼近和插值理论,在这些结果的基础上发展了全空间和半空间上数理方程的谱方法和拟谱方法,从而形成一种新的能更好解决误解区域问题的方法,此种方法被很好的应用于统计物理、量子力学和流体力学中。 4) 我们利用非一致带权Sobolev空间中的Jacobi多项式正交逼近和Jacobi-Gauss型插值理论,提出以Jacobi多项式为基函数的Jacobi谱方法和拟谱方法用来解决一些奇异问题和计算某些特定的无界区域问题。 5)有限谱方法是基于有限点、有限项的局域谱方法。这种方法要求近似函数应具有等同隔网格和非周期性的性质。有限谱方法分为基于非

偏微分方程数值解法答案

1. 课本2p 有证明 2. 课本812,p p 有说明 3. 课本1520,p p 有说明 4. Rit2法,设n u 是u 的n 维子空间,12,...n ???是n u 的一组基底,n u 中的任一元素n u 可 表为1n n i i i u c ?==∑ ,则,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???=== -=-∑∑是12,...n c c c 的二次函数,(,)(,)i j j i a a ????=,令 () 0n j J u c ?=?,从而得到12,...n c c c 满足1 (,)(,),1,2...n i j i j i a c f j n ???===∑,通过解线性方程组,求的i c ,代入1 n n i i i u c ?==∑, 从而得到近似解n u 的过程称为Rit2法 简而言之,Rit2法:为得到偏微分方程的有穷维解,构造了一个近似解,1 n n i i i u c ?== ∑, 利用,11 11()(,)(,)(,)(,)22j n n n n n n i j i j j i j j J u a u u f u a c c c f ???===-=-∑∑确定i c ,求得近似解n u 的过程 Galerkin 法:为求得1 n n i i i u c ? == ∑形式的近似解,在系数i c 使n u 关于n V u ∈,满足(,)(,) n a u V f V =,对任 意 n V u ∈或(取 ,1j V j n ?=≤≤) 1 (,)(,),1,2...n i j i j i a c f j n ???===∑的情况下确定i c ,从而得到近似解1 n n i i i u c ?==∑的过程称 Galerkin 法为 Rit2-Galerkin 法方程: 1 (,)(,)n i j i j i a c f ???==∑ 5. 有限元法:将偏微分方程转化为变分形式,选定单元的形状,对求解域作剖分,进而构 造基函数或单元形状函数,形成有限元空间,将偏微分方程转化成了有限元方程,利用 有效的有限元方程的解法,给出偏微分方程近似解的过程称为有限元法。 6. 解:对求解区间进行网格剖分,节点01......i n a x x x x b =<<<<=得到相邻节点1,i i x x -

偏微分方程考试重点与作业内容

偏微分方程考试重点与作业内容 记号说明: 第一个数字为章,第二个数字为节,第三个数字为小节。例如§3.1.2节就是第三章第一节的第二小节——Poisson公式。 考试重点: 1.§ 2.1节一阶线性方程的特征线解法(要求会用特征线法求解一阶线性偏微分方程, 计算题)灵活运用 2.§2.2.2节解的表达式(要求牢记解的表达式,会运用解的公式求解全平面上的波动方 程或者证明相关问题,计算题或证明题)灵活运用 3.§2.2.5节半无界问题(掌握奇延拓和偶延拓,计算题或证明题) 4.§2.4.1节分离变量法(掌握分离变量的方法,会用分离变量的方法求解波动方程和热 传导方程,计算题)灵活运用 5.§3.1.1节Fourier变换(牢记Fourier变换和Fourier逆变换的公式,会求函数的Fourier 变换和逆变换,判断题、填空题、选择题)灵活运用 6.§3.1.2节Poisson公式(牢记Poisson公式,填空题、选择题、计算题) 7.§3.1.3节广义函数简介(掌握广义函数的相关定义,会求广义函数以及广义函数的导 数,判断题、填空题、选择题)灵活运用 8.§3.3.1节弱极值原理(书上例题和作业题很重要,证明题) 9.§4.1.1节基本解与Green公式(Green公式的应用,判断题、填空题、选择题) 10.§4.2.1节极值原理(书上例题和作业很重要,证明题) 11.§4.2.4节调和函数的性质(熟记调和函数的相关性质,判断题、选择题、填空题) 作业内容: 第一章29页16题(1) 第二章100页第三题;101页10题;102页11题、12题;104页23题的(1)(3)小题; 26题的(2)(3)小题; 第三章161页2题、3题、4题;162页5题;163页9题;166页18题; 第四章212页1题;213页4题;218页25题; 考试说明: 1.考试内容:考试内容出自上述重点章节,主要是书上例题、作业题;要求灵活运用的内 容,会在例题作业题的基础上做微小改动,比例不超过期末考试的30%。 2.考试分数:平时成绩20分(国培计划的自动满分);期末考试80分。 3.考试题型:判断题;选择题;填空题;计算题;证明题 4.考试难度:考虑到大四学生实习的辛苦、找工作的焦虑、考研的烦躁,考试内容会尽量 简单。

求解偏微分方程三种数值方法

数值模拟偏微分方程的三种方法介绍 (有限差分方法、有限元方法、有限体积方法) I.三者简介 有限差分方法(Finite Difference Methods)是数值模拟偏微分方程最早采用的方法,至今仍被广泛使用。该方法包括区域剖分和差商代替导数两个步骤。首先将求解区域划分为差分网格,用有限个网格节点代替连续的求解区域。其次,利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替进行离散,从而建立以网格节点上的值为未知量的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且十分成熟的数值方法。 差商代替导数后的格式称为有限差分格式,从格式的精度来考虑,有一阶格式、二阶格式和高阶格式。从差分的空间离散形式来考虑,有中心格式和迎风格式。对于瞬态方程,考虑时间方向的离散,有显格式、隐格式、交替显隐格式等。目前常见的差分格式,主要是以上几种格式的组合,不同的组合构成不同的差分格式。差分方法主要适用于结构网格,网格的大小一般根据问题模型和Courant 稳定条件来决定。 有限元方法(Finite Element Methods)的基础是虚位移原理和分片多项式插值。该方法的构造过程包括以下三个步骤。首先,利用虚位移原理得到偏微分方程的弱形式,将计算区域划分为有限个互不重叠的单元(三角形、四边形、四面体、六面体等),在每个单元上选择合适的节点作为求解函数的插值点,将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式,得到微分方程的离散形式。利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。 有限元方法有较完善的理论基础,具有求解区域灵活(复杂区域)、单元类型灵活(适于结构网格和非结构网格)、程序代码通用(数值模拟软件多数基于有限元方法)等特点。有限元方法最早应用于结构力学,随着计算机的发展已经渗透到计算物理、流体力学与电磁学等各个数值模拟领域。

偏微分方程数值解复习题(2013硕士)

偏微分方程数值解期末复习(2012硕士) 一、考题类型 本次试卷共六道题目,题型及其所占比例分别为:填空题20%;计算题80% 二、按章节复习内容 第一章 知识点:Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等; 要求: 熟练一元函数的数值微分公式;会辨认差分格式, 计算线性多步法的局部截断误差和阶; 第二章 知识点:矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容、收敛和稳定性等; 要求: 熟练多元函数的数值微分公式;会建立椭圆型方程边值问题的差分格式;计算局部截断误差;了解极值原理讨论格式的收敛性和稳定性; 第四章 知识点:最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式等; 要求: 会建立抛物型方程边值问题的经典差分格式;计算局部截断误差; 会计算格式的传播因子或传播矩阵;会讨论格式的稳定性; 第五章 知识点:依赖区域、左偏心格式、右偏心格式、中心格式、LF格式、LW 格式、Wendroff格式、跳蛙格式、特征线、CFL条件等; 要求: 建立双曲型方程边值问题的差分格式;计算局部截断误差; 会计算格式的传播因子或传播矩阵;讨论格式的稳定性; 第七章 知识点:单元、线性元、线性基、(单元)刚度矩阵、(单元)荷载向量等;

相关文档
最新文档