数值积分方法的比较教学研究与实验刘鹏飞.docx

Comparative teaching research and experiment of numerical integration methodLiu Pengfei, XU NainanAbstract: Numerical integration is a very important teaching content in the course of calculation method or numerical analysis. The numerical integration method is also an important method to solve practical problems.In this paper, we discuss the basic ideas of numerical integration method and some in the actual teaching comparative research and experiment teaching. Keywords:Numerical integration Calculation methodNumerical analysis Numerical experimentsText: The invention of calculus is a great achievement in the history of human science. However, in practical problems, there are still many limitations in solving integral. If the product function is often very complex, can not find the corresponding original function, or the product function is not a function expression, but in the form of a table, its original function is not meaningful.Therefore, it is significant to investigate the numerical integration method of approximate calculation.1. Basic ideas of numerical integration method By the integral mean value theorem, it is known that (x), b (Qbaf) dx= (B-A) f (N) is set up in a (x), which is a continuous function f (N). Problem is difficult to accurately calculate the average height of F (n) values, the different algorithms of F (n), you can get different numerical quadrature method, such as left (right) formula of rectangular, trapezoidal, Simpson formula. Generally, a, b some node XK at the height of F (XK) weighted average to approximate the average height of F (n), f (XK) linear combination to approximate the definite integral, mechanical quadrature formula: Qbaf(x)dxUEnk=0Akf(xk) Directly using the function computed values in the discrete node function integral, avoids the difficulty in seeking the original function, only need to be determined with the product function and product Ak quadrature coefficients are independent between partitions. In order to ensure the accuracy of the mechanical quadrature formula, it is natural to want accurate simple function as much as possible. Namely, it is required that all the m polynomials are exact, and the m+1 time polynomials are not necessarily accurate. was about coefficient of n + 1 order linear equations.Because the coefficient determinant Vandermonde determinant, and it is not zero,solution Ak is unique. 2. Teaching analysis of several common numerical integration method2.1 The interpolation type quadrature formula In a, b, to node n Lagrange interpolation polynomial (x) as f(x) as the approximation function, can be obtained by interpolation quadrature formula:即 The interpolation type integral formula has the accuracy of N sub algebra, and the formula of Ak0 is stable.2.2 Newton - Cotes formulas The a b is divided into N parts, interpolation quadrature formula for: Called Newton- Cotes formula, the Cotes coefficient is independent of the product function and the integral interval, as long as the interval and other fractions n, can be obtained. Obviously, when n= 1 is a trapezoidal formula; when n= 2 is Simpson formula.2.3 Complex formula By the knowledge of definite integral, the definite integral is only related to the integral interval. The higher the number of polynomials is, the higher the degree of smoothness of the function is, and the Runge phenomenon occurs. Such as n 7, the Newton- Cotes formula is not stable. Therefore, people turn to the integral interval, similar to the piecewise interpolation, the integral interval is divided into a number of small areas, the use of the Newton- Cotes formula between each cell number, and then the results of each cell is added as a function of the whole range of the approximation, which is the basic idea of complex. Commonly used compound formula, complex trapezoid formula:Complex Simpson formula:2.4 Successive bisection technology and Romberg formula How to determine the appropriate n so that the difference between the approximate value and the true value of the allowable range, generally speaking, it is more difficult. And successive half technology is in quadrature process according to the requirements of precision and automatically determine n choose whether or not to meet the requirements of precision, the difference of two minutes before and after twice to estimate the error, such not only reduce the step size, and can keep the original calculating results, reduce the amount of computation. For be obtained: Called Romberg formula. The acceleration effect of Romberg formula is very significant, and the calculation amount is very small under the same precision.2.5 Gauss type quadrature formula To further improve the algebraic precision of quadrature formula, through the proper selection of interpolation nodes and quadrature coefficients, the algebraic accuracy up to a maximum of (2n + 1). The quadrature nodes and quadrature coefficients as parameters. The equations can be obtained by using the orthogonal polynomials. You can use orthogonal polynomials with different weighting functions of Gauss type quadrature formulas of different types are3. Comparison of several classes of numerical integration methods in Teaching Due to practical problems should be taken into account when the calculation precision and calculation, thus different quadrature formulas have different characteristics: Trapezoidal quadrature formula and Simpson quadrature formula although the calculation is simple, easy to use, but poor accuracy, but for less smooth integrand function sometimes ratio, high precision method is more effective. The effect of the periodic function is more prominent, especially the trapezoidal formula. Newton-Cotes, the n7 formula is unstable, while the complex trapezoid formula and Simpson formula in the retention of the advantages of low order formula, and can get higher accuracy, so the most widely used in the actual calculation. Romberg method using two points method is simple and easy to program. When the node encryption improve integral approximation degree, previous calculation results can be used for the back, very beneficial to reduce the calculation amount, and the estimation error is relatively simple, can get some integral sequence, if in the convergence of the control, and to check the rows and columns, to function in different state can with the fastest convergence sequence to approximate the integral. High precision Gauss type quadrature formulas, especially other methods cannot match the calculation of infinite integral and defect integral. But because of the node is irregular, when the node is increased, the function value of the previous calculation can not be used, and the N value of the node value and the coefficient of the different values are also stored in the computer.4.Numerical integration method in the teaching of numerical experiments In the teaching process of numerical integration method, the numerical experiments are carried out to improve the students understanding of the teaching contents and to master the basic method of numerical integration. The Matlab program is written in each method, and the advantages and disadvantages of various numerical integration methods are grasped and understood by comparison.4.1 Comparison of the accuracy of different algorithms with the same number of nodes In the course of teaching, the accuracy of various methods can be determined, but it is not conducive to students understanding. Using the powerful function of the popular software, the students can be used to operate the machine. It can be used to realize the accuracy of different numerical integration methods, and it can also improve the students experimental ability and practical ability.For example the precision requirement for.Comparison with the exact solutions of the trapezoidal formula, Simpson formula, Gauss formula. Use Matlab to write programs to run after the results showPrecision Trapezoidal SimpsonGauss0.7853981635 0.7849814972 0.7853981535 0.7853981635 Obviously, the trapezoidal formula with low accuracy, low Simpson formula and the Gauss quadrature formula is better, the highest precision in a given precision and exact solutions of (Precision) agreement. And then describes the method of various quadrature Matlab commands by calling and format, compared with their own programming, and make the students master the calculation method of numerical integration, in the calculation of actual problems skillfully.4.2 Comparison of the different points of the same algorithm. integrate Using trapezoid formula considering the same algorithm sync long more or different integral interval of scores in the comparison of calculation accuracy and computation time, let the students fully understand in numerical calculation of integral interval the importance of correct treatment. Matlab programming operation result shows that:Interval number:10204080The calculation results:0.63050.62280.62090.6204(精确解为: 0.6203)Computation time:0000.0100Graphical presentation: In short, along with the rise of the mathematics experiment, it has played a positive role in promoting the whole mathematics curriculum teaching reform, but also put forward a severe test to the reform of teaching methods,