北京理工大学偏微分方程期末考试

黑龙江科技学院考试试题

第一套

课程名称：偏微分方程数值解法 课程编号：24014110 适用专业（班级）：数学 共1页 命题人：潘晓丽 教研室主任： 第1页

一、（15分）求方程x2uxx 2xyuxy y2uyy xux yuy 0的通解。

y'' y' 3 y 0二、（10分）解固有值问题

y 0 y 1 0

0 x 1

三、（10分）写出二阶线性偏微分方程的分类和标准式。

四、（15分）求半无界弦的自由振动问题

utt a2uxx t 0,x 0

u t,0 0

u 0,x x ,ut 0,x x 五、（15分）计算积分 pn x dx,n为偶数。

01

0

1

六、（15分）将函数f x

2 1

1 x x

按勒让德多项式展开。

x 1

七、（10分）截面为矩形的无穷长棱柱，内部无热源，两相对侧面与外界绝热，另

两相对侧面温度分别保持零度和与高度无关的稳恒分布，求此棱柱的温度分布。 八、（10分）写出Legendre多项式的各种表示。

黑龙江科技学院考试试题答案

第一套

课程名称：偏微分方程数值解法 课程编号：24014110 适用专业

黑龙江科技学院考试试题

第一套

课程名称：偏微分方程数值解法 课程编号：24014110 适用专业（班级）：数学 共1页 命题人：潘晓丽 教研室主任： 第1页

一、（15分）求方程x2uxx 2xyuxy y2uyy xux yuy 0的通解。

y'' y' 3 y 0二、（10分）解固有值问题

y 0 y 1 0

0 x 1

三、（10分）写出二阶线性偏微分方程的分类和标准式。

四、（15分）求半无界弦的自由振动问题

utt a2uxx t 0,x 0

u t,0 0

u 0,x x ,ut 0,x x 五、（15分）计算积分 pn x dx,n为偶数。

01

0

1

六、（15分）将函数f x

2 1

1 x x

按勒让德多项式展开。

x 1

七、（10分）截面为矩形的无穷长棱柱，内部无热源，两相对侧面与外界绝热，另

两相对侧面温度分别保持零度和与高度无关的稳恒分布，求此棱柱的温度分布。 八、（10分）写出Legendre多项式的各种表示。

黑龙江科技学院考试试题答案

第一套

课程名称：偏微分方程数值解法 课程编号：24014110 适用专业

1，2.1 Why are girls not as good at math and science?

Girls don't do as well at maths and science as boys.

There are several reasons for this.The first and most important reason is that they aren't encouraged to play with toys that build up interest in math and science problems and that build skills for problem solving or understanding how things work. Girls are encouraged to play with toys that help foster language and human relations skills. As a result,they can grow up not knowing how an engine works or how to build a

数学与计算科学学院

实 验 报 告

实验项目名称 用Eular方法求解一阶常微分方程数值解 所属课程名称 偏微分方程数值解 实 验 类 型 验证性 实 验 日 期 2015-3-26

班 级 信计12-2班 学 号 201253100215 姓 名 张洪清 成 绩

一、实验概述： 【实验目的】 学会使用显性Eular方法和隐形Eular方法 应用显性Eular方法和隐形Eular方法求解一般一阶常微分方程的近似数值解。 学会用MATLAB解决数学问题。 【实验原理】 1、Eular方法： 一阶线性微分方程初值问题 ?y'?f(x,y),a?x?b??y(

实 验 报 告

课程名称：偏微分方程数值解院 系：专业班级：学 号：学生姓名：指导教师：开课时间： 数学科学系 信计1101 1131120140 张军 沈 林 2013至2014学年第二学期

一、学生撰写要求

按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。

学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。

二、教师评阅与装订要求

1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。

2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零

科目代码 241 243 244 245 333 360 360 611 611 611 613 616 617 617 619 624 625 626 6英语 无参考书目 《标准日本语》初级（上、下）中级（上） 无参考书目 简明法语教程 (上、下) (修订版) 《大学数学简明教程[M]》 数学分析教程（上，下） 数学分析（上，下） 人民教育出版社 科目名称 参考书目 出版社名称 编者 版次 日语 第一版 德语 法语 商务印书馆 清华大学出版社 高等教育出版社 高等教育出版社 机械工业出版社 高等教育出版社 高等教育出版社 孙辉 盛祥耀，陈魁，王飞燕 李忠 方丽萍 陈纪修 於崇华 金路 张润琦 陈一宏 何书元 2005 数学B 数学分析 1 2 2007 2006年6月 2007 数学分析 数学A 微积分(上、下册) 数学A 概率论与数理统计 数学A 线性代数 杨刚、吴惠彬 药理学综合 药理学 李瑞 第六版 基础英语 法学基础（法理学、宪法学） 法学基础（法理学、宪法学） 基础日语 无参考书目 北京大学、高等教育出版社 北大、高教出版社 高等教育出版社 高等教育出版社 高等教育出版 宪法

偏微分方程数值解试题

1、考虑一维的抛物型方程：

?u?2u??2, x?[0,?], 0?t?T?t?x u(x,t)x?0?u0, u(x,t)x???u?u(x,0)??(x)（1）导出时间离散是一阶向前Euler格式，空间离散是二阶精度的差分格式；

（2）讨论（1）中导出的格式的稳定性； （3）若时间离散为二阶精度的蛙跳格式，

?uun?1?un?1 ??tt?tn2?t空间离散是二阶精度的中心差分，问所导出的格式稳定吗？为什么？

2、考虑Poission方程

??2u(x,y)?1, (x,y)???u ?0, in AB and AD?nu(x,y)?0, in BC and CD其中Ω是图1中的梯形。

图1 梯形

使用差分方法来离散该方程。由于梯形的对称性，可以考虑梯形的一半，如图2，

图2 从物理空间到计算区域的几何变换

?，然后在??上使用差分为了求解本问题，采用如下方法：将Ω的一半投影到正方形区域??上用N?N个网格点，空间步长为方法来离散该方程。在计算区域???????1/N(?1) 。

?（带有坐标?,?）（1）引入一个映射T将原区域

为了促进学术交流，加强学术合作，提高国内研究生和年轻教师的业务水平和科研能力，给研究生和青年教师提供学习、交流和了解国内外关于非线性偏微分方程最新研究动态的机会，由香港中文大学数学研究所辛周平教授倡导和组织，国内华南师范大学、

华中师范大学、兰州大学、武汉大学、南京师范大学、首都师范大学、四川师范大学、西北大学、湘潭大学、中山大学等学

校联合主办的“非线性偏微分方程暑假讲习班”已历时八届了。每次参加讲习班学习的研究生和青年教师达100多人，同时也有很多知名专家、学者应邀来讲习班作精彩的学术报告，产生了良好的社会效应和深远的学术影响。第九届“非线性偏微分方程暑假讲习班”将于2011年7月13日至7月28日在中山大学举办。现将有关事宜通知如下：

一、 学术委员会

主 席

辛周平 香港中文大学数学研究所 委 员（排名以汉语拼音为序）

曹道民（中科院数学与系统科学研究院） 陈 化（武汉大学） 邓引斌（华中师范大学）

丁夏畦（中科院数学与系统科学研究院） 范先令（兰州大学）

郭柏灵（北京应用物理与计算数学研究所） 洪家兴（复旦大学） 黄云清（湘潭大学）

江 松（北京应用物理与计算数学研究所） 倪维明（University of Minn

偏微分方程数值习题解答

李微分方程数值解习题解答 1-1 如果 (0) 0，则称x0是J(x)的 驻点（或稳定点）.矩阵A对称（不必正定），求证x0是J(x)的驻点的充要条件是：x0是方程组 Ax b的解

证明:由 ( )的定义与内积的性线性性质,得 ( )

'

J(x0 x)

12

(A(x0 x),x0 x) (b,x0 x)

J(x0) (Ax0 b,x)

'

2

2

(Ax,x)

( ) (Ax0 b,x) (Ax,x)

必要性:由 (0) 0,得,对于任何x R,有 (Ax b,x) 0,

由线性代数结论知,

Ax b 0,Ax b

'

n

00

充分性: 由Ax

'

b,对于任何x R

n

,

(0) (Ax0 b,x) (Ax,x)| 0 0

即x是J(x)的驻点. §1-2

补充: 证明f(x)的不同的广义导数几乎处处相等.

证明:设f L(I),g,g L(I)为f(x)的广义导

2

2

1

2

偏微分方程数值习题解答

数,由广义导数的定义可知,对于任意 (x) C(I),有

0

b

ab

g1(x) (x)dx f(x) (x)dx

b

'

a

a

g2(x) (x)dx f(x) (x)dx

a

b

'

两式相减,得到

ba

(g1 g2) (x) 0

1

2

偏微分方程数值解试题(06B)

参考答案与评分标准

信息与计算科学专业

一（10分）、设矩阵A对称，定义J(x)?1(Ax,x)?(b,x)(x?Rn)，2?(?)?J(x0??x).若?'(0)?0,则称称x0是J(x)的驻点（或稳定点）.矩阵A对

称（不必正定），求证x0是J(x)的驻点的充要条件是：x0是方程组 Ax?b的解 解: 设x0?Rn是J(x)的驻点,对于任意的x?Rn,令

?(?)?J(x0??x)?J(x0)??(Ax0?b,x)??22(Ax,x), (3分)

?'(0)?0,即对于任意的x?Rn,(Ax0?b,x)?0,特别取x?Ax0?b,则有

(Ax0?b,Ax0?b)?||Ax0?b||2?0,得到Ax0?b. (3分) 反

之

,若

x0?Rn满足

Ax0?b,则对于任意的

1x,J(x0?x)??(1)??(0)?(Ax,x)?J(x0),因此x0是J(x)的最小值点. (4分)

2评分标准:?(?)的展开式3分, 每问3分,推理逻辑性1分

ddu??Lu??(p)?qu?f二（10分）、 对于两点边值问题：?dxdx'??u(a)?0,u(b)?0x?[a,b]x?(a,b)

其中p?C1([a,b]),