基于边界节点法MATLAB工具箱开发.docVIP

基于边界节点法MATLAB工具箱开发 　　0 引 言 　　Helmholtz方程作为简化的波动方程，长期以来备受关注.求解Helmholtz方程的数值方法主要有有限差分法、有限元法和边界元法等.这些方法主要基于网格近似或者需要背景网格积分，在处理高波数或高频率Helmholtz方程及高维问题时存在网格划分困难和内存消耗大等缺点.基于点的近似思想，无网格方法在一定程度上减少对网格的依赖，甚至不需要网格.目前，主要有2类无网格法：区域型无网格法和边界型无网格法.区域型无网格法主要有Kansa方法和MLPG法等；边界型无网格法由于只需要边界节点的信息，近年来备受关注.基本解法作为一个典型的边界型无网格法已被成功应用于工程领域，如数值求解Helmholtz方程等，基本解法已经取得很大的成功，但为克服核函数的奇异性，虚拟边界的选取至今仍无定论.CHEN等［1-3］提出的边界节点法不仅保持基本解法的优点，而且克服基本解法最大的缺点――虚拟边界难以选取.边界节点法的主要思想是将控制方程的通解作为核函数，避免基本解法中核函数的奇异性，将源点和配点取在物理边界上.边界节点法已被成功应用于Helmholtz?方程求解. 　　为使研究人员更好地了解和运用边界节点法，本文开发基于边界节点法的MATLAB工具箱.自MATLAB出现以来，各领域的专家为研究高技术和应用方面的问题，已经开发出符号计算工具箱和图形处理工具箱等各领域的工具箱，而边界节点法工具箱仍未见报道.本文介绍基于边界节点法开发的MATLAB工具箱，用于求解简化的Helmholtz方程，主要考虑齐次Helmholtz方程、非齐次Helmholtz方程以及Helmholtz方程控制的边界型反问题.运用本文开发的MATLAB工具箱，可非常简单地运用边界节点法拟合上述问题，将正则化方法加入到本文开发的MATLAB工具箱中，便于研究人员更好地了解正则化方法对边界节点法的影响. 　　1 Helmholtz方程和边界节点法 　　1.1 Helmholtz方程 　　 二维Helmholtz控制方程及其边界条件为式中： 　　Δ为Laplace算子； 　　?λ?为控制方程的波数； 　　?Ω?为二维单连通区域； 　　?Γ?1和Γ??2分别为Dirichlet边界和Neumann边界； 　　?B?1和B??2为边界条件算子； 　　?p?为二维区域节点坐标?(x,y)?； 　???f(p)?为方程的源项； 　　?g(p?1)和?h(p?2)?分别为给定的Dirichlet和Neumann边界条件值.当f(p)=0时，式（1）为齐次方程，反之为非齐次方程.当Γ?1∩Γ?2=?且Γ?1∪Γ?2=Ω时，只需要给出区域内部的数值解；当同时给出边界Γ?1上的?Dirichlet条件和Neumann?条件且边界Γ?2不可测时，式（1）～（3）构成一类典型的反问题，需给出不可测边界Γ?2上的数值解. 　　1.2 边界节点法 　　 边界节点法主要通过控制方程的非奇异通解的线性组合来近似方程的解，可表示为当N增加到一定程度时，式（6）的条件数会剧烈变化，所得解的收敛曲线会出现振荡现象.有效条件数能更好地表现矩阵的病态程度，故引入有效条件数.当求解的问题为反问题时，矩阵方程的解并不连续依赖于边界条件的值，且若边界条件含有扰动时，会对数值解带来致命影响.正则化方法已经被用于解决上述问题.［4-5］本文介绍的MATLAB工具箱主要采用以下几种正则化方法. 　　2 正则化方法和正则化参数 　　在引入正则化方法前简单介绍奇异值分解（Singular Value Decomposition,SVD）.对于式（6）中的矩阵式中：第一项为残差项；第二项为正则化项，包含解的先验信息，如光滑性和有界性等；参数?μ?≥0为正则化参数控制正则化项和残差项的比重，当?μ?=0时，式(16)转化为最小二乘问题.基于SVD的Tikhonov正则化方法可描述为 　　3 MATLAB工具箱简介 　　基于边界节点法求解Helmholtz方程的MATLAB工具箱主要由3类功能型的函数组成：(1)产生边界数据及内部检测点数据函数;(2)公用数据计算函数;(3)求解函数. 　　3.1 产生边界数据和内部检测点数据函数 　　3.1.1 产生边界数据函数 　　 bkm_2d_circle 　　功能：产生二维圆形边界数据. 　　用法：?［x，y］?=bkm_2d_circle(?O?x，O?y，r，n?). 　　描述：用于产生以?（O?x，O?y）?为圆心、?r?为半径的均匀分布的?n?个点. 　　bkm _2d_rect 　　功能：用于产生二维矩形边界节点. 　　用法：?［x，y］?=bkm_2d_rect(?O?x，