攻读博士、硕士学位研究生试卷（作业）封面。
2012 至 2013 学年度第 1 学期）
题目双曲线方程的两个常用格式特点
科目偏微分方程的数值解法。
姓名张明发。
专业凝聚态物理。
入学年月2012-9-5
双曲线方程的两个常用格式特点。
摘要对迎风格式和二阶迎风格式的整理，列出了它们截断误差和稳定性的的特点。通过在实际计算，来比较二者的差异，找出最合理的计算方法。
关键词迎风格式二阶迎风格式。
引言迎风格式在实际计算中引起来普遍的重视，从而产生了很多好的方法和技巧。从简单的一阶线性双曲线方程开始，构造差分格式，分析其稳定性及其他性质。构造差分格式最简单最常用的方法是差商代替微商，迎风格式就是在双曲线方程中关于空间偏导数用在特征线方向一侧的单边差商来代替。
一、迎风格式。
截断误差： o(t+h)
稳定性： |a|l1 (l=t/h).
前两式可合写为。
二、顺风格式。
截断误差： o(t+h)
稳定性： 绝对不稳定。
若差分方程与微分方程特征线走向一致， 则条件稳定； 否则不稳定。
三、beam-warming格式。
3.1 a>0,用a、b、c三点的值进行抛物插值得u(q)
截断误差： o(t2+h2+th)
稳定性： al2 (l=t/h)
这里的讨论适合于a>0的情形。
3.2 a<0,用c、d、e三点进行抛物插值得u(q)
截断误差： o(t2+h2+th)
稳定性： |a|l2 (a<0, l=t/h).
b-w格式的缺点：
1. 须借助其它格式计算或。
2. 会引起解的非物理振荡。
例、考虑定解问题。
1)用迎风格式求解。
2)用beam-warming格式求解。
3)比较两种方法的异同。
解：取h=0.01, =计算至tn=0.5.
1）迎风格式求解。
2）beam-warming格式求解。
3)比较：迎风格式把解抹平了， 在两端点都没有振荡。数值解较精确解在波的上游偏小，下游偏大。
而beam-warming格式的数值解在波的上游较迎风格式更接近精确解，但在波的下游出现了振荡，很大距离后才与解析解重合。
四、程序。迎风格式。
program main
parameter (dj=150)
real a,h,lambda,x,x0,x1,dett,t0,t,t1
dimension u(0:dj),u1(0:dj)
integer i,j,k
x0=-0.5;x1=1.0;t0=0;t1=0.50;a=1.0;h=0.01;lambda=0.5;det=lambda*h;x=x0+h
do i=1,dj-1,1
if(x>0.0) then
u(i)=0.0
else u(i)=1.0
endifx=x+h
enddou(0)=1.0;u(dj)=0.0;t=t0+det
u1(0)=1.0
do while(t<=t1)
open(1,file=''
x=x0+h
do j=1,dj,1
u1(j)=u(j)-a*lambda*(u(j)-u(j-1
x=x+henddo
do k=1,dj,1
u(k)=u1(k)
enddot=t+det
enddox=x0
write(*,t-det
do j=0,dj,1
write(1,10)x,u1(j)
x=x+henddo
10 format(1x,f5.3,5x,f8.5)
close(1)
endbeam-warming格式。
program main
parameter (dj=150)
real a,h,lambda,x,x0,x1,dett,t0,t,t1
dimension u(0:dj),u1(0:dj)
integer i,j,k
x0=-0.5;x1=1.0;t0=0;t1=0.50;a=1.0;h=0.01;lambda=0.5;det=lambda*h;x=x0+h
do i=1,dj-1,1
if(x>0.0) then
u(i)=0.0
else u(i)=1.0
endifx=x+h
enddou(0)=1.0;u(dj)=0.0;t=t0+det
do while(t<=t1)
open(1,file=''
x=x0+h
u1(1)=u(1)-a*lambda*(u(2)-u(0))/2.0+a**2.0*lambda**2.0*(u(2)-2*u(1)+u(0))/2.0
do j=2,dj,1
u1(j)=u(j)-a*lambda*(u(j)-u(j-1))-a*lambda/2.0*(1.0-a*lambda)*(u(j)-2*u(j-1)+u(j-2))
x=x+henddo
do k=1,dj-1,1
u(k)=u1(k)
enddot=t+det
enddox=x0
write(*,t-det
do k=0,dj,1
write(1,10)x,u(k)
x=x+h
enddo10 format(1x,f5.3,5x,f8.5)
close(1)end
《微分几何》作业
一。填空题。1.曲面的第一基本形式为。2.空间曲线的基本公式是。3.曲面在任一点 u，v 的单位法向量公式为 4.空间曲线的切向量为。5.曲线的主法向量总是指向曲线 凹入方向。6.曲面上正常点满足的条件为。7.曲线的挠率表达式为。8.曲面上曲率线满足的微分方程为。9.在曲面s 上，线的微分方程是。1...
微分几何期中作业
证明 一个正则参数曲面是旋转面的充分必要条件是，它的所有法线都与一条固定的直线相交。证明过程 必要性 当我们假设旋转面的方程为。此时，如果对r分别对u,v求导，则可以得到。所以，通过对其做向量乘法可得 因此，我们可知单位法向量为 该曲线的法线的参数方程是。这里的 是法线上的参数，如果，则该处的法线与...
作业14微分中值定理
1 填空题。1 设，则在区间内方程有 2 个实根 在区间内有 2 个实根 2 函数在区间 0，1 上的有限增量公式中的等于 3 在区间上，函数满足罗尔微分中值定 4 在区间上，函数满足拉格朗日微分中值定理中的 5 在区间上，函数及满足柯西微分中值定理中的 2 设在上可导，证明 存在，使 证 设，则在...