Trapezoidal RuleThe code approximates the definite integral of a function. The integral is calculated using the trapezoidal rule. Parameters of the function are the limits of integration and the number of intervals within the limits.The function to be integrated is another parameter and must be defined before running this program. For example, if we want to integrate the function f(x) = y = x3, we can define it, in Matlab, like this: function y = fun_for_integration(x) y = x^3; And then, we can create the integration function, like this: function y = trap2(lower_lim, upper_lim, interv, fun) % initialize the result y = 0; % 'step' is related to the size of each interval step = (upper_lim - lower_lim) / interv; % add up the area of each trapezoid for j = lower_lim : step : upper_lim y = y + feval(fun,j); end % compute integral y = (y - (feval(fun, lower_lim) + feval(fun, upper_lim))/2) * step; Now, we can call our integration function from another Matlab script or from the command window, for example z1 = trap2(0, 2, 10, 'fun_for_integration') z2 = trap2(0, 2, 100, 'fun_for_integration') z3 = trap2(0, 2, 1000, 'fun_for_integration') and we get the Matlab results z1 = 4.0400 z2 = 4.0004 z3 = 4.0000 Note that this is a numerical integration, and so we have to be very aware of the possible inaccuracies of the method. The first two parameters given to the function are the lower and upper limits, respectively. The third parameter is related to the size of each interval: the higher the number the better accuracy in the evaluation we can reach. The fourth parameter is the name of the function to be integrated (the instruction 'feval' is in charge of evaluating that function in the main body of the code). Now, the best part is that Matlab has its own function to do the integration using the trapezoidal rule ('trapz'), so we can save all of our programming thinking for other things... However, the built-in 'trapz' function, works a little different. We first define our interval and desired step in a variable vector x, and define the value of the corresponding function in a variable vector y. For example, we want to evaluate the same function as above, within the interval [0 2], in 0.1 increments, so we use x = 0 : .1 : 2. If we want a finer tunning of the integration function, we can set the step in our independent variable. Then we can call the Matlab function like this (three different cases) x = 0 : .1 : 2; y = x.^3; z4 = trapz(x, y) x = 0 : .01 : 2; y = x.^3; z5 = trapz(x, y) x = 0 : .001 : 2; y = x.^3; z6 = trapz(x, y) and we get the Matlab response, with different approximations z4 = 4.0100 z5 = 4.0001 z6 = 4.0000 From 'Trapezoidal Rule' to home From 'Trapezoidal Rule' to 'Matlab Cookbook'
|
