一维波方程空间半离散化的边界可观性论文大学论文.docx
《一维波方程空间半离散化的边界可观性论文大学论文.docx》由会员分享,可在线阅读,更多相关《一维波方程空间半离散化的边界可观性论文大学论文.docx(11页珍藏版)》请在冰豆网上搜索。
![一维波方程空间半离散化的边界可观性论文大学论文.docx](https://file1.bdocx.com/fileroot1/2023-3/1/312998a5-399b-42f6-803c-fcd83673f45c/312998a5-399b-42f6-803c-fcd83673f45c1.gif)
一维波方程空间半离散化的边界可观性论文大学论文
一维波方程空间半离散化的边界可观性
摘要:
我们考虑在有界区间上具有齐次Dirichlet边界条件的一维波方程的空间半
离散化。
我们分析了边界可观性问题,也就是当网格尺寸hT0时,通过观测边
谱导出的伪模式,我们证明不存在一致界。
然而,在离散系统的低频谱生成的解子空间中,一致界是存在的。
当hT0时,这个有限维空间变大并最终覆盖整个
空间。
这样,当网格尺寸趋于零时,熟知的连续系统的可观性性质是离散观测估计的极限。
我们考虑了有限差分和有限元两种半离散化。
1.引言
考虑一维波方程:
(1.1)
系统(1.1)在能量空间H0(O,LML2(O,L)中是适定的。
确切地说,对任何
(Uo,U1)-h0(O,L)xL2(O,L),存在唯一解u。
(1.1)解的能量由
(1.2)
给出,并且能量是关于时间守恒的,即:
E(t)=E(0).
众所周知,当时间满足TA2L,解的总能量可以通过集中观测一个边界点的
能量一致估计,比如说在x=L处。
确切地说,对任何T>2L,存在C(T)>0,使得对(1.1)的每个有限能量解
(1.4)
成立。
当同时在两个端点进行观测时,上述不等式对所有的T>L成立。
在本文,我们主要研究不等式(1.4)。
形如(1.4)的不等式和波方程的边界可控性密切相关。
文[9]和[11]对波方程和板方程的相关问题进行了系统分析。
在本文,我们研究几个波方程空间半离散化的(1.4)相似形式。
为了说明上面提及的问题,让我们先考虑有限差分半离散化。
给定N^N,令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)的解,能量Eh关于时间是守恒的,即:
(1.8)
本文的主要目标就是研究如下(1.4)的离散形式
(1.9)
注1.1让我们讨论一下为何选择-UN(t)/h作为法向导数Ux(L,t)的近似。
都知道泰勒展开式给出Ux(L,t)的最简单逼近是
或者,用上述记号表示为
由于Dirichlet边界条件,考虑到山十=0,可有
另一方面,就像我们将在2.1节看到的那样,当频率固定和hT0时,对于每
个特征向量(特征函数),-UN(t)/hTUx(L,t)。
这也说明-UN(t)/h也是法向导
数Ux(L,t)的较好近似。
根据不等式(1.4),人们期望当时间T>2L时,存在一个与h无关的常数
C=C(T),使得不等式(1.9)对于(1.6)的每个解和每个0vhv1成立。
本文的第一个结果断定这是办不到的:
定理1.1对任何T>0,有
(1.10)
就像我们将要看到的那样,这是由数值格式产生的高频谱形成的伪模式造成的后果。
在文[3-5]中由R.Glowinski等人早已发现这种现象,文[3-5]与高维波方程的边界精确可控性和所谓实用数值HUM方法有密切关系。
在这些文献中,为了抑
制高频谱的病态效应,提出两种方法:
(a)当计算控制函数时,为了使二次泛函取得最小值,Tychonoff正则化过程被引入;(b)为了缩短离散系统解的分量的波长,滤波法被引入。
通过各种数值实验说明这两种方法的有效性。
为了证明定理1.1,我们分析了系统(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)
这些特征值可以明确地计算出来:
(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.1对于系统(1.11)的任何特征向量,如下等式成立:
(1.16)
这个等式给出了特征向量的总能量和由量|Un(t)/h|2表示的集中在端点x=L出观
测到的能量的关系。
另一方面,对所有hAO和(1.11)的所有特征值,容易验证
(1.17)
但是,(1.17)并不能排除(1.16)右侧的的常数爆破的可能。
事实上,容易验证
(1.18)
因而,爆破出现了。
这也立即说明了定理1.1的结果。
为了证明定理1.1的正面结果,我们不得不引出合适的(1.6)的解集。
给定任何O(1.19)
的特征向量生成的(1.6)的解集Ch(Y)。
更确切的就是
(1.20)
根据引理1.1,根据边界处的能量汇集,进入解集Ch(Y)的每个特征向量的能量
可以被一致估计。
只要时间T足够大,如下结果确保不等式(1.9)对于(1.6)的解集Ch(Y)中的每个解是一致成立的。
定理1.2设0vYc4。
则存在T(Y)二2L,对所有Tat仏)和在解集Ch(Y)中
的每个解,存在C=C(T,Y)使得当hT0时,(1.9)一致成立,此外
(a)
(b)
注1.2定理1.2断言只要T足够大,一致观测不等式(1.9)在解集Ch(Y)中成立。
事实上,当Yt4时,T仏)T4。
这是由于当特征值趋于临界值4时,相连特征
值的平方根的间隙趋于零。
然而,当Yt0时,观测常数T仏)收敛到2L,2L是
系统(1.1)的观测时间。
根据这个结果,并注意到对于TA2L和每个(1.6)如下形式的解
(1.21)
一致观测不等式(1.9)成立。
其中4(h)满足
(1.22)
这说明当hT0时,原始系统(1.1)的可观性以离散系统(1.6)的(1.21)-(1.22)的形式
的解的可观性的极限被恢复。
也可以看出一致观测不等式
(1.9)中的常数C仃,Y)收敛到L/2仃-2L),
L/2仃-2L)是利用离散乘子法得到的连续系统(1.1)的观测常数(见[11])。
注1.3容易看出系统(1.6)是可观的。
由于它是常微分系统,对所有T>0它
是可观的(见[10])。
因而,对所有h>0和T>0,存在C=C(T,h)使得(1.9)对于(1.6)
的每个解成立。
然而,为了得到一致观测不等式常数(当hT0时),我们不得
不滤掉高频部分(也就是考虑解集Ch(Y)中的解)并取足够大的时间T。
粗略地说,当hT0时,只要滤掉高频解,定理1.2可以保证半离散系统是
致可观的。
我们将给出定理1.2的两种证明方法。
第一种是采用经典乘子法,这种方法常用于证明波方程和板方程的可观性(见[9,11])。
第二种方法依赖于非调和
Fourier级数中的经典Ingham不等式。
你也许会认为这些结果的出现是由于我们选取了特殊的有限差分离散化导致的。
但事实不是这样的。
我们也考虑了有限元半离散化并得到了类似的结果,既有定理1.1的否定结果,也有定理1.2的正面结果。
值得提及的是,我们此处解释的关于波方程的离散化出现的现象早已从具有
剧烈震荡周期系数的一维波方程的可观性研究中被人们得知(见[1]和[2])。
在这两
种情形中,具有微局部结构或离散格式出生病态高频震荡的波的相互作用没有在连续模型中出现。
下文安排如下:
在第二节研究有限差分逼近。
特别地,我们提出并证明引言中叙述的结果。
在第二节研究有限元离散化。
在第四节简要比较有限差分和有限元半离散化。
2.有限差分半离散化
在本节,我们详细分析在引言中讨论的波方程(1.1)的有限差分空间半离散化的可观性问题。
首先,我们将进行详细的谱分析。
特别地,我们证明引理1.1以及它的直接
结果定理1.1。
然后利用乘子法详细证明定理1.2.我们也指出如何利用非调和
Fourier级数中的熟知结果证明同样的结论。
2.1谱分析
系统(1.6)的特征函数和特征值满足
(2.1)
这实际上是矩阵的特征值问题。
可以精确算出系统(2.1)的特征值和特征向量。
有(见[8],P.456):
(2.2)
(2.3)
特别地,可以看出离散系统的特征向量和连续系统的特征函数si"(兀‘)是相一致
的。
另一方面,对于固定的k,有
(2.4)。
这是连续系统的第k个特征值。
引理2.1对于(2.1)的任意特征值A相应的特征向量W,如下等式成立:
(2.5)。
如果Wk和Wl是特征值加对应的特征向量,有
证明:
用护j乘以(2.1)将上述等式对于1,2,…,N求和,则立即得到等式(2.5)(注意
到等式(2.5)可以解释为(A弟®)=)妙®))。
为了得到(2.6)式,需要指出的是由于A是对称矩阵,不同特征值对应的特
征向量是正交的和A正交的,即(俨,创)=0和(A护,3)=0。
这样有
(2.7)。
因而有
换句话说
(2.8)。
根据(2.7),上式和(2.6)等价。
对于较小的时间段,由于在x=L处所产生的几乎不可区分的时间谐波与特
征值密切相关,所以在分析边界可观性问题中连续特征值的间隙起着重要的作用。
对于连续模型,我们有
(2.9)。
这样,间隙兀/L和频域无关。
然而,就像我们将要看到的那样,在离散问题中,连续特征值的高频间隙逐渐增大并且当hT0时,它和h是同阶的。
我们有:
并且考虑到(N+1)/h=L,有
可以得到:
因而对某个jN,只要
(2.11)
成立,我们就有
(2.12).
这说明相应于指标k=N+1-j,N+2-j,…,N的第j个特征值间的距离是和h是
同阶的,级相差一个随着j增大而增大的数乘因子。
特别地,可以找到最大特征值间的距离的一个上界如下:
反过来,对于较小特征值,还可以找到它的下界。
事实上,
(2.13).
当k仍然有界或即使k无界但hk有上界6L(0v6<1),当hT0时,上述不等式
的右端收敛到兀/L。
注意到兀/L是连续模型连续特征值平方根的间隙。
2.2特征向量的边界可观性
本节的目的就是证明引理1.1的等式(1.16)。
根据(2.2)和(2.3)式中的特征向量和特征值的精确值,这个等式很容易直接证
明。
然而,我们将用乘子法来证明。
首先,我们将特征向量规范化使其满足
(2.14)
现在(2.5)变为:
(2.15)
这样有
(2.16)
用j(®j+-®jJ/2乘以等式(2.1)的两侧并将其从1加到N(注意这是乘子X护x的离
散版)。
根据(2.14)和(2.16)在左侧有:
在右侧我们有
因而有
(2.17)
换句话说就是
(2.18)
结合(2.18)和(2.15)有
这就完成了引理1.1的证明。
注2.1在注1.1我们已经说明选择-UN(t)/h作为Ux(L,t)逼近的合理性。
根据(2.3)
中特征向量的明显表达式,立即可以看出对任何固定的k,
对任何固定的
然而注意到当砧2
k,我们也可以在引理1.1的(1.16)式中取极限。
然后得到等式
T4(ashT0)时等式(1.16)退化。
这也是一个在波方程数值逼
近格式中由其导致的病态高频振荡的典型事例。
2.3非一致可观性的证明
本节用于证明定理1.1.就像在引言中指出的那样,它是引理1.1的直接结果。
事
实上,令u是相应于第N个特征向量的(1.6)的解,即:
(2.19)
根据引理1.1我们有
(2.20)
另一方面有
(2.21)
根据(2.20和(2.21)即可推出
(2.22)
然而根据(2.2):
结合(2.22)和(2.23),定理1得证。
2.4离散波方程的边界可观性:
乘子法
本节用乘子法来证明定理1.2。
就像我们将在后面2.5节将要看到的那样,通过Fourier级数技巧,可以更容易地得到本节的结果。
然而,由于在高维空间中相似问题分析时可能有潜在的应用价值,因而我们认为在本节提供的离散乘子法还是有特别意义的。
首先我们给出一些基本的等式。
引理(2.2)(能量守恒)