小波分析在反应堆物理中的应用.docx

《核反应堆热工数值分析》课程报告 PAGE \* MERGEFORMAT8 小波分析简介及其在反应堆物理中的应用 小波分析是近二十年来迅速发展起来的全新数值分析方法，其基本思想是将函数用小波基函数来展开，计算其展开系数。由于小波基函数具有优良的紧支性，使得小波基函数可以聚集在研究的任意细节，被数学家和工程师们誉为“数学显微镜”[9]。即小波基函数可以很好的逼近各种剧烈变化的函数，尤其是局部震荡的函数，具有非常高的逼近精度。另外，小波分析具有多尺度、多分辨的特性[3]，能够提供多种基函数作为有限元插值函数，由此构造的小波基单元可以根据实际需要任意改变分析尺度，使在变化梯度小的求解域用大的分析尺度，而变化梯度大的求解域则采用小的分析尺度。这是一种优于传统单元网格加密和插值阶次升高的自适应有限元算法，这种变尺度算法数值稳定性好、运算速度快、求解精度高。本文简要介绍了小波基函数的基本性质、求解小波基函数的展开系数的数值方法、以及小波基函数在反应堆物理中的应用，包括用小波基函数离散中子通量密度中的角度变量、共振能量区间的能量变量，并结合有限元方法的思想，将所有的变量都用基函数展开，然后带入变分方程，得出一个代数方程组，从而可以求出展开系数。 理论基础 小波函数 小波函数基本概念 小波函数是一个在(-∞，+∞)区间积分为零的函数[9]，（容许条件）  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 称母小波；“小”指的是函数的局部性，积分为零说明含有波动性，而且正负两个方向的波动是均等的；小波函数经过平移、伸缩可以得到一族小波函数，即：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 其中a反应了特定函数的尺度（即伸缩情况），变量b指明了它沿着x轴的平移位置；当a和b取一系列离散值，如 a=2j,b=2jk(j,k∈Z)时，即可以得到一族小波函数，即：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 在一定条件下ψj,k|?j,k∈Z可以构成L2(R)的基，用它可以表示f(x)  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) Daubechies 小波 在数值计算中，通常使用Daubechies 小波[9]，Daubechies 小波具有紧支性和正交性；Daubechies 小波没有显式的数学表达式，其尺度函数和小波函数由以下两尺度方程给出：  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5)  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ MTEqn \c \* Arabic \* MERGEFORMAT 6)  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTSec \c \* Arabic \* MERGEFORMAT 0- SEQ