正压大气Helmholtz方程小波数值解法研究.pdfVIP

正压大气Helmholtz 方程小波数值解法* 汪萍 戴新刚 马柱国 中国科学院东亚区域气候-环境重点实验室 大气物理研究所东亚中心 北京 100029 摘 要 用Wavelet-Galerkin方法求解了准地转正压大气的Helmholtz方程 检验了FFTJacobian和QR三 种矩阵分解方法的正交性和重构误差 发现无论是正交性还是矩阵重构 FFT分解方法的误差都最大 其 余两种方法的误差较小且很接近 Jacobian方法略好于QR对Helmholtz方程和Poisson方程的小波数值 求解检验表明 FFT方法会造成很大的误差 Jacobian和QR的精度较高 此外 用Daubechies 6小波基 -3 底求解的精度要比迭代算法高一个量级,达到 10 随着小波阶数的提高数值解的精度也会有进一步的提 高 关键词 正压大气 Helmholtz 方程 Wavelet -Galerkin 方法 矩阵分解 1引言 Helmholtz方程是求解大气运动数值预报模式最常遇到的问题[1,2]已经有几种成熟的 [3] 数值解法 如迭代法 Fourier变换法 变分法和有限元法等 每种解法各有其特点 可 以根据边界条件或计算要求进行选择 利用小波基求解微分方程的方法近些年来发展很快 [4] 小波基是局域基底 在波数空间和物理空间都具有良好的局域性[57] 用小波表示函数 [8,9,10] 或信号可以大大压缩信号的自由度 同时可以保持很高的重构精度 这些优良性质是 吸引人们用小波求解微分方程的主要原因 目前基于小波基求解微分方程主要有三种方法 [11,12] [13 15] Wavelet-Galerkin方法 小波配置法 和 Wavelet-Capacitance matrix [11,16,17] method 第一种方法的主要困难是边界条件较难处理 第二和第三种方法解可以在边 界附近获得较高的精度 但不宜用于大规模数值计算 Qian和 Weiss首次用 Wavelet-Galerkin方法将一个 Poisson方程的解表示成矩阵形式并求出了数值解 但这种 解法需要进一步完善 本文主要结合周期边界条件下Helmholtz方程的Wavelet-Galerkin 数值解法对有关的代数和数值求解精度等问题做一详细研究 下面的第2节介绍周期小波基 第3节用小波变换法导出小波谱形式的Helmholtz方程 第 4节研究矩阵的正交分解与误差 第 5节在周期边界条件下计算不同的矩阵分解方法下 Helmholtz方程的数值解 第6节检验小波阶数对方程数值解精度的影响 最后是结论 2周期小波基和小波变换 利用小波函数求解微分方程的主要困难之一是复杂边界的处理 传统的小波定义在实数 轴上 有限区间上的小波基构造仍然是一个正在发展的方向[18] 为简单起见 我们选择定 义在[0,1]区间上的Daubechies 周期小波作为展开基底[19] 周期小波可以通过对传统小波做 Poisson 求和得到[19] 如果 (x)和 (x)是传统小波函数和尺度函数 则周期尺度函数和周 j,k j,k 期小波函数可以表示为 * 国家重点基础研究发展规划项目 G1999043400 和国家自然科学基金 No 联合资助 1 ~ jj , k (x) jj ,k (x l) 1 l ZZ ~ y