期刊文献+

不等距网格上求解ODE特征值问题若干高精度格式的计算与分析

CALCULATION AND ANALYSIS OF SEVERAL HIGH AUCCURACY SCHEMES FOR SOLVING ODE EIGENPROBLEM ON NON-UNIFORM GRID
原文传递
导出
摘要 本文针对不等距网格,从Raylei曲商(Raylei曲quotient)角度出发,构造了若干求解ODE特征值问题的高阶格式,并进行误差分析.文中高阶格式的构造是基于线性有限元及其对应的差分格式进行的.单纯的线性有限元及其对应的差分格式求解PDE特征值问题都只有二阶精度,我们利用质量集中和加权组合的思想通过将二者结合得到四阶精度的算法.本文从理论和实验的角度构造高阶格式并进行了相应的误差分析.通过在五种网格上计算四阶精度格式的误差阶系数,将四阶格式加权组合的新格式甚至可以达到六阶精度.最后用数值实验验证了构造的高阶格式的误差阶.同时,本文构造的两种四阶格式相对于传统的线性有限元方法,在同等量级误差的要求下,需要的网格数有量级的减少. From the perspective of Rayleigh quotient, we constructed several high-order accuracy schemes to solve the ODE eigenproblem on non-uniform grid, and analyze the error. Purely linear finite element scheme and finite difference scheme to solve corresponding ODE eigen- problem are only second order accuracy, Use a weighted combination of the two and error compensation or mass lumping method to construct the fourth-order accuracy schemes. We construct high-order accuracy schemes and analyse the corresponding error both in theo- retical and experimental levels. Through numerical experiments on five kinds of grid, we construct new schemes even up to sixth-order accuracy, by computing the coefficient of the specific error-order of the fourth-order accuracy schemes. At last, we use numerical ex- periment to verify the accuracy order of these new schemes. At the request of the same order of magnitude error, the number of grid that our two fourth-order schemes need have a magnitude of decrease compared to the traditional linear finite element method.
出处 《数值计算与计算机应用》 CSCD 2014年第2期131-152,共22页 Journal on Numerical Methods and Computer Applications
基金 国家自然科学基金(91230109)资助项目
关键词 特征值问题 不等距网格 高精度格式 质量集中 加权组合 eigenproblem non-uniform grid high-order accuracy schemes mass lump-ing weighted combination
  • 相关文献

参考文献4

二级参考文献36

  • 1LAI JunJiang,HUANG JianGuo,SHI ZhongCi.A lumped mass finite element method for vibration analysis of elastic plate-plate structures[J].Science China Mathematics,2010,53(6):1453-1474.
  • 2Li H, Sun J, Xu Y. Discrete Fourier analysis, cubature and interpolation on a hexagon and a triangle [J]. SIAM J. Numer. Anal. 2008, 46: 1653-1681.
  • 3Peter, Li and Shing-Tung Yau. On the SchrSdinger equation and the Eigenvalue problem[J]. Comm. Math. Phys., 1983, 88: 309-318.
  • 4Pinsky Mark A. The eigenvalues of an equilateral triangle[J]. SIAM J. Math. Anal., 1980, 11: 819-827.
  • 5Pinsky Mark A. Completeness of the eigenfunctions of the equilateral triangle[J]. SIAM J. Math. Anal., 1985, 16: 848-851.
  • 6Milan Prager. An investigation of eigenfunctions over the equilateral triangle and square[J]. Applications of mathematics, 1998, 43: 311-320.
  • 7McCartin B J. Eigenstructure of the equilateral triangle. I. The Dirichlet problem[J]. SIAM Rev., 2003, 45: 267~287.
  • 8Sun Jiachang and Li Huiyuan. Generalized Fourier transform on an arbitrary tirangular domain[J]. Advances in Computational Mathematics, 2005, 2: 223-248.
  • 9Laugesen R S and Siudeja B A. Dirichlet eigenvalue sums on triangles are minimal for equilaterals. preprint, 2011.
  • 10Laugesen R S, Pan Z C and Son S S. Neumann eigenvalue sums on triangles are (mostly) minimal for equilaterals. Mathematical Inequalities & Applications, preprint, 2011.

共引文献13

相关作者

内容加载中请稍等...

相关机构

内容加载中请稍等...

相关主题

内容加载中请稍等...

浏览历史

内容加载中请稍等...
;
使用帮助 返回顶部