外区域上Dirichlet–Neumann算子的对角化和广义逆的正则逆表示.doc

外区域上Dirichlet-Neumann算子的对角化和广义逆的正则逆表示 【摘要】：本文是在法方导师FrancoisAlouges和中方导师陈果良教授的共同指导下完成的，全文共分成五个部分，第一部分是外区域上的Dirichlet-Neumann算子的对角化，这部分工作是在法国期间由法方教授FrancoisAlouges指导完成，其余四个部分关于互补广义逆的性质推导及应用则是中方导师陈果良教授悉心指导的结果．若Ω是R~d(d=2ou3)中的一个闭的凸集，在实际应用中经常需要求如下的能量积分(例如铁磁学中关于磁能的计算)其中φ满足在Ω上调和并且在无穷远处趋于零．利用Dirichlet-Neumann算子的定义，上式可以改写为然而众所周知，已有的求解Dirichlet-Neumann算子的近似算法的时问是与边界(?)Ω上的点的个数成平方关系，当点的个数足够多时，这将导致实际计算的困难，区域是三维的情形这一困难尤为突出．本文的想法是利用Dirichlet-Neumann算子DN的特征值和特征向量(λ_n,(?)_n)n≥1来逼近上述能量方程．若那么很自然的我们有此时能量方程转化为我们可以用有限项来逼近上式，通过本文第一部分的实例我们知道当区域为规则的圆或椭圆时，上式将很快收敛，即能量是由Dirichlet-Neumann算子DN的最小的若干个特征值和特征向量来决定．从而上述问题转化为求解Dirichlet-Neumann算子DN的最小的若干个特征值和相应的特征向量．在木文的第一章中我们首先分析了当区域为圆和球的情形，通过Fourier分析和调和分析，我们得出了如下结论：当区域为圆时，Dirichlet-Neumann算子DN的特征值为λ_n=|n|，其重数为2，相应的特征向量为(?)_n=e~(inθ)，即Fourier基．当区域为球时，Dirichlet-Neumann算子DN的特征值为λ_n=n，但其重数为2n+1，相应的特征向量为Y_l~m，为调和函数的基．利用无限元的思想，通过类比三角流我们给出了求解Dirichlet-Neumann算子DN最小若干个特征值和相应的特征向量的一个线性算法．通过大量的实例，与已有的方法相比我们的算法精度更高，且计算时间与区域上的点的个数成线性关系，这导致虽然当点数较少时，我们的算法需要较多的时间，但当点数足够多时，我们的方法将体现出其优越之处．矩阵理论在数值计算、线性规划、数据分析、网络优化等重大领域有着极其广泛的应用．随着科技进步，矩阵理论在实际中应用越来越广泛，矩阵理论的研究也显得越来越重要．在实际应用中，经常遇到如下的线性方程组：Ax=b，其中AC~(m×n)，xC~n，bC~m(1)若A为方阵且非奇异，则(1)一定有唯一解为x=A~(-1)b．若(1)为相容方程组(即bR(A))，则一定有解．但当其解有无数个时，通常我们关心一个范数极小的解．若(1)为矛盾方程组，通常我们求解使残差向量b-Ax范数最小的那个解，即最小二乘解．为了统一解决上述问题，广义逆理论应运而生，并最终解决了这些问题．1920年，E．H．Moore在美国数学会通报上给出了任意矩阵的广义逆的如下定义：设AC~(m×n)，则满足的矩阵GC~(n×m)称为A的广义逆矩阵，记做A~+，其中P_R(A)和P_R(G)分别是R(A)和R(G)上的正交投影矩阵．1955年，R．Penrose在美国剑桥哲学学会学报上以非常简单、直观的形式叙述了广义逆矩阵A~+满足的四个条件：设AC~(m×n)，则满足的矩阵GC~(n×m)称为A的广义逆矩阵，记做A~+．可知上述两个广义逆的定义是等价的，可参阅参考文献[81]．上述定义的广义逆，通常称为Moore-Perrose广义逆，简记为M-P逆，也称作伪逆．在此基础上，又衍生出了其它许多类型的广义逆．设(?)≠η(?){1，2，3，4}，则所有满足η中条件的矩阵集称为A的η-逆．根据广义逆的相关理论，矛盾线性方程组的最小二乘解与A的{1，3}-逆有关，相容线性方程组的极小范数解与A的{1，4}-逆有关，矛盾线性方程组的极小范数最小二乘解则与A的M-P广义逆有关．随着广义逆研究的深入，又产生了群逆，Drazin逆，Boot-Duffin逆，加权广义逆，α-β广义逆，A_(T,S)~(2)逆等其它广义逆．它们的部分定义罗列如下，相关性质可参阅参考文献[5，81，82，83，84，85]等．设AC~(n×n)，满足rank(A~(k+1))=rank(A~k)的最小正整数k称为A的指标，记为ind(A)=k．则满足：的矩阵GC~(n×n)称为A的Drazin逆，记做A~(d)．当ind(A)=1时，Drazin逆称为群逆，记做A~#．设AC~(n×n)，L(?)C~n，且AP_L+P_L非奇异，则称