清华大学有限元分析课件3-4 Gauss Quadrature_607804617.pptxVIP

第3章 一维问题3.4 高斯积分Gauss Quadrature[Fish]-4.63.4 Gauss QuadratureGauss quadrature is one of the most efficient techniques for functions that are polynomials or nearly polynomialsApproximate the integral by?Wi — weightsξi — integration points3.4 Gauss QuadratureGauss Quadrature: Choose the weights and integration points that yield an exact integral of a polynomial of a given order3.4 Gauss QuadratureG. Q.A system of nonlinear algebraic equations for the unknown matrices M and WM — n parametersW — n parametersf(ξ) — polynomial of order p defined by p+1 parametersp＋1 ＝ 2n : An n-point Gauss formula can integrate a (2n?1)-order polynomial exactlyThe number of integration points need to integrate a polynomial of order p exactly is given byIn one dimension the Gauss quadrature formulas are optimal, but in multi-dimensional cases, the Gauss quadrature rules are no longer necessarily optimal3.4 Gauss QuadraturePosition of Gauss points and corresponding weights3.4 Gauss Quadrature% get gauss points in the parent element domain [-1, 1]% and the corresponding weightsfunction [w,gp] = gauss(ngp) if ngp == 1 gp = 0; w = 2;elseif ngp == 2 gp = [-0 0; w = [1, 1];elseif ngp == 3 gp = [-0.7745966692, 0.7745966692, 0.0]; w = [0.5555555556, 0.5555555556, 0.8888888889];end