chapter03_matlab格林函数_格林函数_格林函数MATLAB_格林函数的分析_bargrz_

% File: dss004.m % function [ux]=dss004(xl,xu,n,u) % % Function dss004 computes the first derivative, u , of a % x % variable u over the spatial domain xl le x le xu from classical % five-point, fourth-order finite difference approximations % % Argument list % % xl Lower boundary value of x (input) % % xu Upper boundary value of x (input) % % n Number of grid points in the x domain including the % boundary points (input) % % u One-dimensional array containing the values of u at % the n grid point points for which the derivative is % to be computed (input) % % ux One-dimensional array containing the numerical % values of the derivatives of u at the n grid points % (output) % % The mathematical details of the following Taylor series (or % polynomials) are given in routine dss002. % % Five-point formulas % % (1) Left end, point i = 1 % % 2 3 4 % a(u2 = u1 + u1 ( dx) + u1 ( dx) + u1 ( dx) + u1 ( dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u1 ( dx) + u1 ( dx) + u1 ( dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % b(u3 = u1 + u1 (2dx) + u1 (2dx) + u1 (2dx) + u1 (2dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u1 (2dx) + u1 (2dx) + u1 (2dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % c(u4 = u1 + u1 (3dx) + u1 (3dx) + u1 (3dx) + u1 (3dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u1 (3dx) + u1 (3dx) + u1 (3dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % d(u5 = u1 + u1 (4dx) + u1 (4dx) + u1 (4dx) + u1 (4dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u1 (4dx) + u1 (4dx) + u1 (4dx) + ...) % 5x 5f 6x 6f 7x 7f % % Constants a, b, c and d are selected so that the coefficients % of the u1 terms sum to one and the coefficients of the u1 , % x 2x % u1 and u1 terms sum to zero % 3x 4x % % a + 2b + 3c + 4d = 1 % % a + 4b + 9c + 16d = 0 % % a + 8b + 27c + 64d = 0 % % a + 16b + 81c + 256d = 0 % % Simultaneous solution for a, b, c and d followed by the solu- % tion of the preceding Taylor series, truncated after the u % 4x % terms, for u1 gives the following five-point approximation % x % 4 % u1 = (1/12dx)(-25u1 + 48u2 - 36u3 + 16u4 - 3u5) + O(dx ) (1) % x % % (2) Interior point, i = 2 % % 2 3 4 % a(u1 = u2 + u2 (-dx) + u2 (-dx) + u2 (-dx) + u2 (-dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u2 (-dx) + u2 (-dx) + u2 (-dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % b(u3 = u2 + u2 ( dx) + u2 ( dx) + u2 ( dx) + u2 ( dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u2 ( dx) + u2 ( dx) + u2 ( dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % c(u4 = u2 + u2 (2dx) + u2 (2dx) + u2 (2dx) + u2 (2dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u2 (2dx) + u2 (2dx) + u2 (2dx) + ...) % 5x 5f 6x 6f 7x 7f % % 2 3 4 % d(u5 = u2 + u2 (3dx) + u2 (3dx) + u2 (3dx) + u2 (3dx) % x 1f 2x 2f 3x 3f 4x 4f % % 5 6 7 % + u2 (3dx) + u2 (3dx) + u2 (3dx) + ...) % 5x 5f 6x 6f 7x 7f % % -a + b + 2c + 3d = 1 % % a + b + 4c + 9d = 0 % % -a + b + 8c + 27d = 0 % % a + b + 16c + 81d = 0 % % Simultaneous solution for a, b, c and d followed by the solu- % tion of the preceding Taylor series, truncated after the u % 4x % terms, for u1 gives the following five-point approximation % x % 4 % u2 = (1/12dx)(-3u1 - 10u2 + 18u3 - 6u4 + u5) + O(dx ) (2) % x % % (3) Interior point i, i ne 2, n-1 % % 2 3 % a(ui-2 = ui + ui (-2dx) + ui (-2dx) + ui (-2dx) % x 1f 2x 2f 3x 3f % % 4 5 6 % + ui (-2dx) + ui (-2dx) + ui (-2dx) + ...) % 4x 4f 5x 5f 6x 6f % % 2 3 % b(ui-1 = ui + ui ( -dx) + ui ( -dx) + ui ( -dx) % x 1f 2x 2f 3x 3f % % 4 5 6 % + ui ( -dx) + ui ( -dx) + ui ( -dx) + ...) % 4x 4f 5x 5f 6x 6f % % 2 3 % c(ui+1 = ui + ui ( dx) + ui ( dx) + ui ( dx) % x 1f 2x 2f 3x 3f % % 4 5 6 % + ui ( dx) + ui ( dx) + ui ( dx) + ...) % 4x 4f 5x 5f 6x 6f % % 2 3 % d(ui+2 = ui + ui ( 2dx) + ui ( 2dx) + ui ( 2dx) % x 1f 2x 2f 3x 3f % % 4 5 6 % + ui ( 2dx) + ui ( 2dx) + ui ( 2dx) + ...) % 4x 4f 5x 5f 6x 6f % % -2a - b + c + 2d = 1 % % 4a + b + c + 4d = 0 % % -8a - b + c + 8d = 0 % % 16a + b + c + 16d = 0 % % Simultaneous solution for a, b, c and d followed by the solu- % tion of the preceding Taylor series, truncated after the u % 4x % terms, for u1 gives the following five-point approximation % x % 4 % ui = (1/12dx)(ui-2 - 8ui-1 + 0ui + 8ui+1 - ui+2) + O(dx ) (3) % x % % (4) Interior point, i = n-1 % % 2 3 % a(un-4 = un-1 + un-1 (-3dx) + un-1 (-3dx) + un-1 (-3dx) % x 1f 2x 2f 3x 3f % % 4 5 6 % + un-1 (-3dx) + un-1 (-3dx) + un-1 (-3dx) + ... % 4x 4f 5x 5f 6x 6f % % 2 3 % b(un-3 = un-1 + un-1 (-2dx) + un-1 (-2dx) + un-1 (-2dx) % x 1f 2x 2f 3x 3f % % 4 5