一维波方程空间半离散化的边界可观性论文_大学论文

1、界处的能量，解的总能量能否被一致估计的问题。由于存在数值格式产生的高频一维波方程空间半离散化的边界可观性摘要：我们考虑在有界区间上具有齐次Dirichlet边界条件的一维波方程的空间半离散化。我们分析了边界可观性问题，也就是当网格尺寸hT 0时，通过观测边谱导出的伪模式，我们证明不存在一致界。然而，在离散系统的低频谱生成的解 子空间中，一致界是存在的。当hT 0时，这个有限维空间变大并最终覆盖整个空间。这样，当网格尺寸趋于零时，熟知的连续系统的可观性性质是离散观测估 计的极限。我们考虑了有限差分和有限元两种半离散化。1.1.引言考虑一维波方程:(1.1)系统(1.1)在能量空间H0(O,LML

2、2(O,L)中是适定的。确切地说，对任何(Uo,U1)- H0(O,L)XL2(O,L)，存在唯一解 u。(1.1)解的能量由(1.2)给出，并且能量是关于时间守恒的，即：E(t) = E(0).众所周知，当时间满足T A2L ,解的总能量可以通过集中观测一个边界点的能量一致估计，比如说在x=L处。确切地说，对任何T2L，存在C(T)0， 使得对(1.1)的每个有限能量解(1.4)成立。当同时在两个端点进行观测时，上述不等式对所有的T L成立。在本文，我们主要研究不等式(1.4)。形如(1.4)的不等式和波方程的边界可控性密切相关。文9和11对波方程和 板方程的相关问题进行了系统分析。在本文，

3、我们研究几个波方程空间半离散化的(1.4)相似形式。为了说明上面提及的问题，让我们先考虑有限差分半离散化。给定NN，令h =L/(N +1)，引入网格(1.5)其中Xj = jh。然后，引入如下(1.1)的有限差分半离散化(1.6)在(1.6 )中，撇代表关于时间的导数。由于根据边界条件有U0=O,UN十=0,系统(1.6)是关于未知量U1,UN的N元线性微分方程组。显然，u是(1.1)的解，只要初始数据(Uj0,Uj1)是(1.1 )的初值的一个近似值,Uj(t)就是U(Xj,t)的一个近似值。系统(1.6)的能量由(1.7)给出，这是连续能量E的离散化形式。容易看出，对于(1.1)的解，能

4、量Eh关于时间是守恒的，即:(1.8)本文的主要目标就是研究如下(1.4)的离散形式(1.9)注1.11.1让我们讨论一下为何选择-UN(t)/h作为法向导数Ux(L,t)的近似。都知道 泰勒展开式给出Ux(L,t)的最简单逼近是或者，用上述记号表示为由于Dirichlet边界条件，考虑到 山十=0，可有另一方面，就像我们将在2.1节看到的那样，当频率固定和hT 0时，对于每个特征向量(特征函数)，-UN(t)/hT Ux(L,t)。这也说明-UN(t)/h也是法向导数Ux(L,t)的较好近似。根据不等式(1.4)，人们期望当时间T 2L时，存在一个与 h无关的常数C=C(T)，使得不等式(1

5、.9)对于(1.6)的每个解和每个0vhv1成立。本文的第一个结果断定这是办不到的：定理1.11.1对任何T0,有(1.10)就像我们将要看到的那样，这是由数值格式产生的高频谱形成的伪模式造成的后 果。在文3-5中由R. Glowinski等人早已发现这种现象，文3-5与高维波方程的 边界精确可控性和所谓实用数值 HUM方法有密切关系。在这些文献中，为了抑制高频谱的病态效应，提出两种方法：(a)当计算控制函数时，为了使二次泛函 取得最小值，Tychonoff正则化过程被引入；(b)为了缩短离散系统解的分量的波 长，滤波法被引入。通过各种数值实验说明这两种方法的有效性。为了证明定理1.1,我们分

6、析了系统(1.6)的谱，并且对于相应于(1.6)的特征值的特征向量利用离散乘子法得到精确的观测不等式。为了证明定理1.1的正面 部分，也就是在hT0时，形如(1.9)的不等式是一致成立的，我们利用离散乘子 法。如上所述，为了使这些不等式是一致成立的， 必须剔除由数值格式产生的高 频伪模式。利用(1.6)的低频谱生成的合适的(1.6)的解子空间，或者换句话说，(1.6) 的解的Fourier展开的合适截断，这个问题可以解决。因而，我们的方法同上面 提及的滤波法很相似。更详细地说，让我们考虑相应于(1.6)的特征值问题:(1.11)用71(h),-N(h)代表(1.11)的N个特征值:(1.12)

7、这些特征值可以明确地计算出来：(1.13)相应于特征值 (h)的特征函数 护也可以明确的算出来:(1.14)根据系统(1.11)的特征向量，(1.6)的解有如下的Fourier展开。确切地说，对于合 适的系数ak,bk-R , (1.6 )的每个解u=(ui,，UN)能写成:(1.15)其中ak,bk R可由初值明确地给出。在详细研究(1.6)的解的观测不等式前，分析特征向量的边界观测是十分有意 义的。如下引理给我们提供了解答:引理1.11.1对于系统(1.11)的任何特征向量，如下等式成立:(1.16)这个等式给出了特征向量的总能量和由量|UN (t)/h |2表示的集中在端点x=L出观测到

8、的能量的关系。另一方面，对所有hAO和(1.11)的所有特征值，容易验证(1.17)但是，(1.17)并不能排除(1.16)右侧的的常数爆破的可能。事实上，容易验证(1.18)因而，爆破出现了。这也立即说明了定理1.1的结果。为了证明定理1.1的正面结果，我们不得不引出合适的(1.6)的解集。给定任 何OY0它是可观的(见10)。因而，对所有h0和T0,存在C =C(T,h)使得(1.9)对于(1.6)的每个解成立。然而，为了得到一致观测不等式常数(当hT 0时)，我们不得不滤掉高频部分(也就是考虑解集 Ch(Y)中的解)并取足够大的时间T。粗略地说，当hT 0时，只要滤掉高频解，定理1.2可

9、以保证半离散系统是致可观的。我们将给出定理1.2的两种证明方法。第一种是采用经典乘子法，这种方法 常用于证明波方程和板方程的可观性(见9,11)。第二种方法依赖于非调和Fourier级数中的经典Ingham不等式。你也许会认为这些结果的出现是由于我们选取了特殊的有限差分离散化导 致的。但事实不是这样的。我们也考虑了有限元半离散化并得到了类似的结果, 既有定理1.1的否定结果，也有定理1.2的正面结果。值得提及的是，我们此处解释的关于波方程的离散化出现的现象早已从具有剧烈震荡周期系数的一维波方程的可观性研究中被人们得知（见1和2）。在这两种情形中，具有微局部结构或离散格式出生病态高频震荡的波的相

10、互作用没有在 连续模型中出现。下文安排如下：在第二节研究有限差分逼近。特别地，我们提出并证明引言 中叙述的结果。在第二节研究有限元离散化。在第四节简要比较有限差分和有限 元半离散化。2.2.有限差分半离散化在本节，我们详细分析在引言中讨论的波方程（1.1）的有限差分空间半离散化 的可观性问题。首先，我们将进行详细的谱分析。特别地，我们证明引理1.1以及它的直接(2.6)。结果定理1.1。然后利用乘子法详细证明定理1.2.我们也指出如何利用非调和Fourier级数中的熟知结果证明同样的结论。2.12.1谱分析系统（1.6）的特征函数和特征值满足(2.1)这实际上是矩阵的特征值问题。可以精确算出系

11、统（2.1）的特征值和特征向量。有（见8,P.456）:(2.2)(2.3)特别地，可以看出离散系统的特征向量和连续系统的特征函数si （兀）是相一致的。另一方面，对于固定的k，有(2.4)。这是连续系统的第k个特征值。引理2.1对于（2.1）的任意特征值A相应的特征向量W ,如下等式成立:(2.5)。如果W k和Wl是特征值加对应的特征向量，有(2.10).证明：用护j乘以(2.1)将上述等式对于1,2,，N求和，则立即得到等式(2.5)(注意到等式(2.5)可以解释为(A弟)=)妙)。为了得到(2.6)式，需要指出的是由于 A是对称矩阵，不同特征值对应的特征向量是正交的和A正交的，即(俨,

12、创)=0和(A护,3)=0。这样有(2.7)。因而有换句话说(2.8)。根据(2.7)，上式和(2.6)等价。对于较小的时间段，由于在 x=L处所产生的几乎不可区分的时间谐波与特征值密切相关，所以在分析边界可观性问题中连续特征值的间隙起着重要的作 用。对于连续模型，我们有(2.9)。这样，间隙兀/L和频域无关。然而，就像我们将要看到的那样，在离散问题中, 连续特征值的高频间隙逐渐增大并且当 hT 0时，它和h是同阶的。我们有:并且考虑到(N +1)/h=L，有可以得到:因而对某个j N，只要(2.11)成立，我们就有(2.12).这说明相应于指标k = N +1 - j, N + 2- j,N

13、的第j个特征值间的距离是和h是同阶的，级相差一个随着j增大而增大的数乘因子。特别地，可以找到最大特征值间的距离的一个上界如下: 反过来，对于较小特征值，还可以找到它的下界。事实上,(2.13).当k仍然有界或即使k无界但hk有上界6L(0v61)，当hT 0时，上述不等式的右端收敛到兀/ L。注意到兀/ L是连续模型连续特征值平方根的间隙。2.22.2特征向量的边界可观性本节的目的就是证明引理1.1的等式(1.16)。根据(2.2)和 (2.3 )式中的特征向量和特征值的精确值，这个等式很容易直接证明。然而，我们将用乘子法来证明。首先，我们将特征向量规范化使其满足(2.14)现在(2.5)变为

14、:(2.15)这样有(2.16)用j(j+-jJ/2乘以等式(2.1 )的两侧并将其从1加到N(注意这是乘子X护x的离散版)。根据(2.14)和(2.16)在左侧有:在右侧我们有(2.23)因而有(2.17)换句话说就是(2.18)这是连续波方程第k个特征函数Wk(x)的法相导数。对任何固定的然而注意到当砧2结合(2.18)和(2.15)有这就完成了引理1.1的证明。注2.12.1在注1.1我们已经说明选择-UN(t)/h作为Ux(L,t)逼近的合理性。根据(2.3)中特征向量的明显表达式，立即可以看出对任何固定的 k，k，我们也可以在引理1.1的(1.16)式中取极限。然后得到等式T 4(as hT 0)时等式(1.16)退化。这也是一个在波方程数值逼近格式中由其导致的病态高频振荡的典型事例。2.32.3非一致可观性的证明本节用于证明定理1.1.就像在引言中指出的那样，它是引理1.1的直接结果。事实上，令u是相应于第N个特征向量的(1.6)的解，即:(2.19)根据引理1.1我们有(2.20)另一方面有(2.21)根据(2.