|
|
| (2 intermediate revisions by the same user not shown) |
| Line 1: |
Line 1: |
| − | <!-- metadata commented in wiki content
| + | Tabla.docx |
| − | | + | |
| − | | + | |
| − | ==HHA method to calculate the numerical integration and degree of fluctuation of fractally spaced data.==
| + | |
| − | | + | |
| − | '''Héctor Tabares-Ospina<sup>1</sup>*, Mauricio Osorio<sup>2</sup>'''
| + | |
| − | | + | |
| − | <span id='_Hlk5739322'><sup>1</sup>Universidad Nacional de Colombia, Sede Medellín, Facultad de Minas, [mailto:hatabare@unal.edu.co hatabare@unal.edu.co]. Institución Universitaria Pascual Bravo, Facultad de Ingeniería, Medellín, Colombia; [mailto:h.tabares@pascualbravo.edu.co h.tabares@pascualbravo.edu.co]</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;"><sup>2</sup>Universidad Nacional de Colombia, Sede Medellín, Escuela de matemáticas, maosorio[mailto:@unal.edu.co. @unal.edu.co.]</span></div>
| + | |
| − | -->
| + | |
| − | | + | |
| − | '''Abstract:''' This article proposes a new numerical method to calculate the integral of a function, using the trapezoid rule, but spacing the intervals fractally, with the purpose of discovering other complementary observations in relation to the degree of fluctuation of the function studied, a parameter that can be included in conventional numerical integration methods. The results are compared with the conventional method of generalized Trapezoid Rule, in which the intervals are equally spaced.
| + | |
| − | | + | |
| − | This is then a new contribution that extends universal knowledge on the subject.
| + | |
| − | | + | |
| − | '''AMS subject classifications: '''65D30, 65Yxx, 65M12.
| + | |
| − | | + | |
| − | '''Key words''' '':'' Numerical Integration; Computer aspects of numerical algorithms, Fractal dimension.
| + | |
| − | | + | |
| − | ==1. Introduction==
| + | |
| − | | + | |
| − | A semi-geometric element is “fractal” when it contains a repetitive structure in different scales with self-similarity characteristics. The fractal concept was developed by the mathematician Benoit Mandelbrot in the 1970s [1]. In fact, nature provides examples of fractals such as snowflakes, ferns, peacock feathers, and romanesco broccoli. Fractal theory has been used to analyze data and obtain relevant information in highly complex problems such as biology [2,3], health sciences [4–8], securities markets [9], network communications [10–12], and others. Other interesting applications are related with the electrical power. In [13], authors presented a new methodology to analyze the dynamics of the daily power demand by using the curve formed by attractors that move in the complex plane over the Mandelbrot set according to the law dictated by the load curve.
| + | |
| − | | + | |
| − | The fractal dimension is concept that indicate the degree of fluctuation of any function. In [14] authors explored fractal analysis in order to discriminate between pathological and healthy retinal texture.
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Other example is [15], where bentonite is studied when it is affected by a corrosion process in alkaline solution. The fractal dimensions of bentonite specimens at different stages of the corrosion process were obtained from </span> <math display="inline">{N}_{2}</math><span style="text-align: center; font-size: 75%;"> adsorption and swelling deformation tests, and the results of both methods suggested that the fractal dimension tends to increase with duration of the corrosion process. </span>
| + | |
| − | | + | |
| − | [16] used fractal dimension to describe the behavior of auriferous veins across the boreholes.
| + | |
| − | | + | |
| − | [17] used fractal dimension as an indicator of compactness. The fractal dimension of demolition waste was calculated by fitting the particle size demolition curves into the fractal model.
| + | |
| − | | + | |
| − | The exhaust stream of internal combustion engines contains soot agglomerates with fractal geometries. In [18] the fractal dimension of an agglomerate is estimated from the combination of its volume and a characteristic radius.
| + | |
| − | | + | |
| − | In [19] the fractal dimensions of aggregates and cement are discussed. A comprehensive understanding of the fractal properties of concrete meso-structure is important because it is associated with the complex and random mechanical behaviors of the concrete.
| + | |
| − | | + | |
| − | In [20], authors present a novel method to estimate pedestrian flow in uncontrolled environments by using the fractal dimension measured through the box-counting algorithm, which does not neither require the use of image pre-processing and nor intelligent algorithms.
| + | |
| − | | + | |
| − | In [21] fractal dimension is applied to study the growth of urban people of Beijing city, in China.
| + | |
| − | | + | |
| − | Finally, [22] present the spatial morphology of the soft agglomerate, such as maltodextrin, quantified by fractal dimension.
| + | |
| − | | + | |
| − | The previous Articles are the basis for the new contribution, considering the following hypothesis: The degree of fluctuation of a function is calculated using fractal dimension method. The fractal segments obtained also may be used to calculate the area under the curve, integrating the function using standard trapezoid rule, with fractally spaced data.
| + | |
| − | | + | |
| − | The most outstanding reference is [23], where a method to do numerical integration is studied when the original ‘‘signal’’ shows experimentally some kind of self-similarity.
| + | |
| − | | + | |
| − | Finally, this document is organized as follows. Section 2 provides a brief explanation of fractal dimension and the methods to do the numerical integration using trapezoid rule. Section 3 presents the method applied in this research and the algorithm used. Section 4 reports and discusses the most relevant results. Finally, the main conclusions of this research work are summarized in Sec. 5
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''2 MATHEMATICAL MODEL'''</span></div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''2.1 Fractal dimension'''</span></div>
| + | |
| − | | + | |
| − | <span id='_GoBack'></span><span style="text-align: center; font-size: 75%;">The concept of fractal dimension is introduced with respect to the fact that most objects of nature do not have an entire dimension, but they are in a fractional dimension and this dimension must be greater than its topological dimension. The fractal dimension is a dimensionless numerical measure that indicates how the length, area, or volume varies, depending on the size of the measurement element used. It is calculated with various box sizes, making a linear fit to a </span> <math display="inline">Log(N\left( r\right) )</math><span style="text-align: center; font-size: 75%;"> graph over </span> <math display="inline">Log(r)</math><span style="text-align: center; font-size: 75%;">. As the size of the boxes reduces their area, the points in the logarithmic graph align more and more and the fractal dimension can be calculated as the slope of the line that joins them, by the expression:</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>D=Log(N\left( r\right) )/Log(r)</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(1)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Where </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;"> is the size of the measure element used, and </span> <math display="inline">N\left( r\right)</math> <span style="text-align: center; font-size: 75%;"> is the length, area or volume related.</span>
| + | |
| − | | + | |
| − | In conclusion, D allows to determine the complexity of the set space under study. [24]
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''2.2 Numerical Integration'''</span></div>
| + | |
| − | | + | |
| − | A detailed study on the numerical methods of integration is beyond the scope of this unit. Full information is provided in (Mathews, J, Fink, K, 2000). A recount on the matter, are presented below:
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | :<span style="text-align: center; font-size: 75%;">''''''2.2.1''' Newton-Cótes methods.'''</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">It is sought to calculate the definite integral </span> <math display="inline">\int_{a}^{b}f\left( x\right) dx</math><span style="text-align: center; font-size: 75%;">, using the Newton-Cótes methods, which integrates an interpolation polynomial that approximates </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> in </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">. Therefore, these are methods, according to the interpolation polynomial considered. In the case of linear and quadratic interpolations, these methods are called Trapeze Method and Simpson Method, respectively.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | :<span style="text-align: center; font-size: 75%;">''''''2.2.2''' Trapezoid method'''</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The Trapezoid Method is a Newton-Cotes Method based on linear interpolation. With the purpose of integrating </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> from the point </span> <math display="inline">\left( a,f\left( a\right) \right)</math> <span style="text-align: center; font-size: 75%;">to </span> <math display="inline">\left( b,f\left( b\right) \right)</math> <span style="text-align: center; font-size: 75%;">, </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> is approximated by its linear interpolation polynomial in </span> <math display="inline">\left[ a,b\right]</math>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math display="inline">f\left( x\right) \approx {P}_{1}\left( x\right) =\frac{x-b}{a-b}f\left( a\right) +</math><math>\frac{x-a}{b-a}f(b)</math><span style="text-align: center; font-size: 75%;">, </span> <math display="inline">\forall x\in \left[ a,b\right]</math> </div>
| + | |
| − | | + | |
| − | And so
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>I=\int_{a}^{b}f\left( x\right) dx\approx \int_{a}^{b}{P}_{1}\left( x\right) dx=\frac{b-a}{2}\left( f\left( a\right) +f(b)\right)</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The value of the integral </span> <math display="inline">I</math><span style="text-align: center; font-size: 75%;"> approximates the area of the trapezoid determined by the lines </span> <math display="inline">x=</math><math>a,\, x=b</math><span style="text-align: center; font-size: 75%;">, the axis of abscissa and the line that joins the points </span> <math display="inline">\left( a,f\left( a\right) \right)</math> <span style="text-align: center; font-size: 75%;"> and </span> <math display="inline">\left( b,f\left( b\right) \right)</math> <span style="text-align: center; font-size: 75%;">. Another important aspect is also the knowledge of the degree of precision of the numerical solution.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The expression of the linear interpolation error, assuming that </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> is continuous and derivable twice in the interval </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>f\left( x\right) ={P}_{1}\left( x\right) +\, \epsilon \left( x\right)</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <math display="inline">\epsilon \left( x\right) =\frac{{f}^{''}\left( \delta \right) }{2}(x-</math><math>a)(x-b)</math><span style="text-align: center; font-size: 75%;">, </span> <math display="inline">a\leq \delta \leq b</math></div>
| + | |
| − | | + | |
| − | We will have then
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>I=\int_{a}^{b}f\left( x\right) dx=\frac{b-a}{2}\left( f\left( a\right) +f(b)\right) +E</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Where the error of numerical integration </span> <math display="inline">E</math><span style="text-align: center; font-size: 75%;"> is:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>E=\int_{a}^{b}\epsilon \left( x\right) dx=\frac{{f}^{''}(\delta )}{2}\int_{a}^{b}\left( x-a\right) \left( x-b\right) dx</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Integrating this last expression and naming </span> <math display="inline">h=</math><math>b-a</math><span style="text-align: center; font-size: 75%;">, it is concluded that:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>E=-\frac{{h}^{3}}{12}{f}^{''}\left( \delta \right) \Rightarrow \left| E\right| \leq \left| \frac{{h}^{3}}{12}{M}_{2}\right|</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <math display="inline">{M}_{2}</math><span style="text-align: center; font-size: 75%;">is the maximum value reached by the second derivative of the function in the given interval </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">. The demonstration of this general result can be found in advanced texts on numerical integration.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | :<span style="text-align: center; font-size: 75%;">'''2.2.3''' '''Composite Trapezoid method'''.</span></div>
| + | |
| − | | + | |
| − | If the interval in which the integral is carried out is large, the Simple Trapeze Method is usually very imprecise. To improve accuracy, it is possible to subdivide the interval into smaller ones and apply the Simple Method to each one.
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Thus, the compound or generalized Trapeze Method consists in taking a partition </span> <math display="inline">P=</math><math>\left\{ {x}_{0},{x}_{1},\ldots ,{x}_{n}\right\}</math> <span style="text-align: center; font-size: 75%;"> from </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">, </span> <math display="inline">\left( {x}_{0}=\right. </math><math>\left. a,\, \, {x}_{n}=b\right)</math> <span style="text-align: center; font-size: 75%;">, equiespaced, that is: </span> <math display="inline">{x}_{i+1}-</math><math>{x}_{i}=h</math><span style="text-align: center; font-size: 75%;">, </span> <math display="inline">\forall i=</math><math>1,\ldots ,n.</math><span style="text-align: center; font-size: 75%;"> We will have so:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>h=\frac{b-a}{n}</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Taking into account the basic properties of the definite integral:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{a}^{b}f\left( x\right) dx=\int_{{x}_{0}}^{{x}_{1}}f\left( x\right) dx+\int_{{x}_{1}}^{{x}_{2}}f\left( x\right) dx+\ldots +</math><math>\int_{{x}_{n-1}}^{{x}_{n}}f\left( x\right) dx</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">And applying to each integral the simple method:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{a}^{b}f\left( x\right) dx\approx \frac{h}{2}\left( f\left( {x}_{0}\right) +f\left( {x}_{1}\right) \right) +\frac{h}{2}\left( f\left( {x}_{1}\right) +f\left( {x}_{2}\right) \right) +\ldots +\frac{h}{2}\left( f\left( {x}_{n-1}\right) +f\left( {x}_{n}\right) \right)</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>=\frac{h}{2}\left( f\left( {x}_{0}\right) +2\left( f\left( {x}_{1}\right) +f\left( {x}_{2}\right) +\ldots +\right. \right. </math><math>\left. \left. f\left( {x}_{n-1}\right) \right) +f\left( {x}_{n}\right) \right)</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Therefore, we have the final expression for the Generalized Trapeze Method:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{a}^{b}f\left( x\right) dx\cong \frac{h}{2}\left( f\left( a\right) +2\sum _{i=1}^{n-1}f\left( {x}_{i}\right) +f\left( b\right) \right)</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">With regard to the integration error, it will be equal to the sum of the errors of each of the applications of the simple method:</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>E={E}_{1}+{E}_{2}+\ldots +{E}_{n}=-\frac{{h}^{3}}{12}{f}^{''}\left( {\epsilon }_{1}\right) -</math><math>\frac{{h}^{3}}{12}{f}^{''}\left( {\epsilon }_{2}\right) -\ldots -\frac{{h}^{3}}{12}{f}^{''}\left( {\epsilon }_{n}\right)</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">If we call </span> <math display="inline">{M}_{2}</math><span style="text-align: center; font-size: 75%;"> the maximum of the function </span> <math display="inline">{f}^{''}\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> in </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">, we will have:</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\left| E\right| \leq \left| \frac{{h}^{3}}{12}n{M}_{2}\right| =\left| \frac{\left( b-a\right) }{12}{h}^{2}{M}_{2}\right|</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(2)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">This means that the error is of order </span> <math display="inline">O\left( {h}^{2}\right)</math> <span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | :<span style="text-align: center; font-size: 75%;">''''''2.2.4''' Numerical integration for unequally spaced data.'''</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">In [25] is indicated that all numerical integration formulas have been based on equally spaced data. However, in practice there are many situations in which this assumption is not met and there are segments of unequal sizes. For example, experimentally obtained data, or fractally spaced data, are often of this type. In such cases, one method is to apply the trapezoid rule to each segment and add the results, as indicated in equation (2), where </span> <math display="inline">{h}_{i}=</math><span style="text-align: center; font-size: 75%;"> the width of segment </span> <math display="inline">i</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''Error of Trapezoid rule for unequal spacing'''</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Assume that </span> <math display="inline">y=</math><math>f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> be continuous, and possess continuous derivatives in </span> <math display="inline">\left[ {x}_{0},{x}_{n}\right] .</math><span style="text-align: center; font-size: 75%;"> Expanding </span> <math display="inline">y=</math><math>f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> in Taylor polynomials [26] around </span> <math display="inline">x=</math><math>{x}_{0}</math><span style="text-align: center; font-size: 75%;">, the result obtained is</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>y\left( x\right) ={y}_{0}+\frac{\left( x-{x}_{0}\right) }{1!}{y}_{0}^{'}+\frac{{\left( x-{x}_{0}\right) }^{2}}{2!}{y}_{0}^{''}+</math><math>\ldots</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(3)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Now, using (3)</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{{x}_{0}}^{{x}_{1}}ydx=\int_{{x}_{0}}^{{x}_{1}}\left[ {y}_{0}+\frac{\left( x-{x}_{0}\right) }{1!}{y}_{0}^{'}+\frac{{\left( x-{x}_{0}\right) }^{2}}{2!}{y}_{0}^{''}+\ldots \right] dx</math>
| + | |
| − | |}
| + | |
| − | </div>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>={\left[ {y}_{0}x+\frac{{\left( x-{x}_{0}\right) }^{2}}{2}{y}_{0}^{'}+\frac{{\left( x-{x}_{0}\right) }^{3}}{6}{y}_{0}^{''}+\ldots \right] }_{{x}_{0}}^{{x}_{1}}</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>={y}_{0}\left( {x}_{1}-{x}_{0}\right) +\frac{{\left( {x}_{1}-{x}_{0}\right) }^{2}}{2}{y}_{0}^{'}+</math><math>\frac{{\left( {x}_{1}-{x}_{0}\right) }^{3}}{6}{y}_{0}^{''}+</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(4)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">But</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{{x}_{0}}^{{x}_{1}}ydx=\frac{\left( {x}_{1}-{x}_{0}\right) }{2}\left( {y}_{0}+{y}_{1}\right)</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(5)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Putting </span> <math display="inline">{x=x}_{1}</math><span style="text-align: center; font-size: 75%;"> in (3)</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>y\left( {x}_{1}\right) \approx {y}_{1}={y}_{0}+\frac{{\left( {x}_{1}-{x}_{0}\right) }^{}}{1!}{y}_{0}^{'}+</math><math>\frac{{\left( {x}_{1}-{x}_{0}\right) }^{2}}{2!}{y}_{0}^{''}+\ldots</math>
| + | |
| − | |}
| + | |
| − | </div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Using the value of </span> <math display="inline">{y}_{1}</math><span style="text-align: center; font-size: 75%;"> in (5), the result obtained is</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{{x}_{0}}^{{x}_{1}}ydx=\frac{\left( {x}_{1}-{x}_{0}\right) }{2}\left( {y}_{0}+{y}_{1}\right) =\frac{\left( {x}_{1}-{x}_{0}\right) }{2}\left[ {y}_{0}+{y}_{0}+\frac{{\left( {x}_{1}-{x}_{0}\right) }^{}}{1!}{y}_{0}^{'}+\frac{{\left( {x}_{1}-{x}_{0}\right) }^{2}}{2!}{y}_{0}^{''}+\ldots \right]</math>
| + | |
| − | |}
| + | |
| − | </div>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{=y}_{0}\left( {x}_{1}-{x}_{0}\right) +\frac{{\left( {x}_{1}-{x}_{0}\right) }^{2}}{2}{y}_{0}^{'}+</math><math>\frac{{\left( {x}_{1}-{x}_{0}\right) }^{3}}{4}{y}_{0}^{''}+\ldots</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(6)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Now, for the error term subtracting (6) from (4)</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{{x}_{0}}^{{x}_{1}}ydx-\frac{{\left( {x}_{1}-{x}_{0}\right) }^{}}{2}\left( {y}_{0}+\right. </math><math>\left. {y}_{1}\right) =\frac{{\left( {x}_{1}-{x}_{0}\right) }^{3}}{6}{y}_{0}^{''}-</math><math>\frac{{\left( {x}_{1}-{x}_{0}\right) }^{3}}{4}{y}_{0}^{''}+\ldots</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>=\left( \frac{1}{6}-\frac{1}{4}\right) {\left( {x}_{1}-{x}_{0}\right) }^{3}{y}_{0}^{''}+\ldots =</math><math>-\frac{1}{12}{\left( {x}_{1}-{x}_{0}\right) }^{3}{y}_{0}^{''}+\ldots</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">This is the error for the interval </span> <math display="inline">\left[ {x}_{0},{x}_{1}\right]</math>
| + | |
| − | | + | |
| − | In the same way, the error for remaining intervals and the total error are
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>E=-\frac{1}{12}\left[ {\left( {x}_{1}-{x}_{0}\right) }^{3}{y}_{0}^{''}+{\left( {x}_{2}-{x}_{1}\right) }^{3}{y}_{1}^{''}+\ldots +\right. </math><math>\left. {\left( {x}_{n}-{x}_{n-1}\right) }^{3}{y}_{n-1}^{''}\right]</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>=\sum _{i=1}^{n}{\left( {x}_{i}-{x}_{i-1}\right) }^{3}{y}_{i-1}^{''}=O\left( {h}_{i}^{3}\right) \left[ {h}_{i}=\right]</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">This is the error of Trapezoid rule for unequal spacing</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The literature presents multiple articles on adaptations of the trapezoid rule to perform approximations to the integral of a function. However, there are few articles when the fractal dimension of the area is related. For example, in [28] were estimate the error in replacing an integral </span> <math display="inline">\int_{}^{}f.d\mu \,</math> <span style="text-align: center; font-size: 75%;">with respect to a fractal measure </span> <math display="inline">\mu</math> <span style="text-align: center; font-size: 75%;"> with a discrete sum </span> <math display="inline">\sum _{x\in E}^{}w\left( x\right) f(x)\,</math> <span style="text-align: center; font-size: 75%;">over a given sample set </span> <math display="inline">E</math><span style="text-align: center; font-size: 75%;"> with weights </span> <math display="inline">w</math><span style="text-align: center; font-size: 75%;">.The model is the classical Koksma-Hlawka theorem for integrals over rectangles, where the error is estimated by a product of a discrepancy that depends only on the geometry of the sample set and weights, and variance that depends only on the smoothness of </span> <math display="inline">f</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''3 METHODOLOGY'''</span></div>
| + | |
| − | | + | |
| − | In order to obtain universally valid results with respect to calculate the fractal dimension and the numerical integration, of a function divided with unequal spaced data, the procedure described below was followed. The numerical integration calculated with unequal segments will be compared with trapezoid rule equispaced and with the exact result. The algorithm proposed is:
| + | |
| − | | + | |
| − | '''Defining Input Variables'''
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Assume a function </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> defined in the interval </span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">, represented by the sequence </span> <math display="inline">{\left\{ {S}_{j}\left( x,y\right) \right\} }_{j=1}^{NTer}</math><span style="text-align: center; font-size: 75%;">, where </span> <math display="inline">NTer</math><span style="text-align: center; font-size: 75%;"> is the number of terms, and </span> <math display="inline">{\left\{ {P}_{j}\left( x,y\right) \right\} }_{j=1}^{NI}</math><span style="text-align: center; font-size: 75%;">, where </span> <math display="inline">NI</math><span style="text-align: center; font-size: 75%;"> is the number of intervals of the succession. It is interesting to calculate its fractal dimension, and area under the curve.</span>
| + | |
| − | | + | |
| − | <math display="inline">i:</math><span style="text-align: center; font-size: 75%;">Counter variable related with the iterations.</span>
| + | |
| − | | + | |
| − | <math display="inline">r</math><span style="text-align: center; font-size: 75%;">: Radius between two points.</span>
| + | |
| − | | + | |
| − | <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">: Interval.</span>
| + | |
| − | | + | |
| − | <math display="inline">inc</math><span style="text-align: center; font-size: 75%;">: Increment step to evaluate </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;"> from </span> <math display="inline">a</math><span style="text-align: center; font-size: 75%;"> to </span> <math display="inline">b</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <math display="inline">NTer:</math><span style="text-align: center; font-size: 75%;">Variable related with the Term Numbers.</span>
| + | |
| − | | + | |
| − | <math display="inline">j:</math><span style="text-align: center; font-size: 75%;"> Counter variable related with the </span> <math display="inline">NTer</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <math display="inline">NI:</math><span style="text-align: center; font-size: 75%;">Variable related with the Interval Numbers.</span>
| + | |
| − | | + | |
| − | <math display="inline">k:</math><span style="text-align: center; font-size: 75%;"> Counter variable related with the </span> <math display="inline">NI</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 1. Read </span> <math display="inline">f\left( x\right) ,\, a,\, b,\, inc</math><span style="text-align: center; font-size: 75%;">. </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 2. Calculate </span> <math display="inline">{\left\{ {S}_{j}.y=f({S}_{j}.x)\right\} }_{j=1}^{NTer}</math><span style="text-align: center; font-size: 75%;">, according to the following algorithm:</span>
| + | |
| − | | + | |
| − | j = 1 //Counter variable
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>Value=a</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>WHILE\left( Value\leq b\right) \, Do</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>{S}_{j}.x=Value</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>{S}_{j}.y=f({S}_{j}.x)</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>Value=Value+inc</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>j=j+1</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>EndWHILE</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>NTer=j-1</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 3. Making </span> <math display="inline">NI=1,\, i=</math><math>1</math><span style="text-align: center; font-size: 75%;"> as the number of subintervals in the first iteration. Counter variable </span> <math display="inline">k=</math><math>NI</math><span style="text-align: center; font-size: 75%;">. The first fractal partition of the curve is calculated with the points </span> <math display="inline">{P}_{(1)}=</math><math>{S}_{(1)}</math><span style="text-align: center; font-size: 75%;">, </span> <math display="inline">{P}_{(2)}=</math><math>{S}_{(NTer)}</math><span style="text-align: center; font-size: 75%;">. They join together and calculate the distance, radius, </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;"> between them.</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>r=\sqrt{{\left( {P}_{(k+1)}.x-{P}_{(k)}.x\right) }^{2}+{\left( {P}_{(k+1)}.y-{P}_{(k)}.y\right) }^{2}}</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(7)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 4. A first approach to the area under the curve formed by the trapezoid is calculated</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>{A}_{k}=\frac{h}{2}\left[ {P}_{(k)}.y+{P}_{(k+1)}.y\right]</math>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Where </span> <math display="inline">h=</math><math>{P}_{(k+1)}.x-{P}_{(k)}.x</math>
| + | |
| − | |}
| + | |
| − | | style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(8)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 5. Make </span> <math display="inline">i=i+</math><math>1</math><span style="text-align: center; font-size: 75%;">. A new radius </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;"> and equal to </span> <math display="inline">r=</math><math>r/2</math><span style="text-align: center; font-size: 75%;"> is calculated.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 6. Calculating the number of Intervals </span> <math display="inline">NI</math><span style="text-align: center; font-size: 75%;"> with the new radius </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;">. Searching intercepts from the starting point </span> <math display="inline">{P}_{(1)}=</math><math>{S}_{(1)}</math><span style="text-align: center; font-size: 75%;">, and using (7), a sequential search is made of the new points, which with </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;"> intercepts the function </span> <math display="inline">f\left( x\right)</math> <span style="text-align: center; font-size: 75%;">. The algorithm is the following</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>NI=1</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <math display="inline">L\_Index=0</math><span style="text-align: center; font-size: 75%;"> //Left_Index</span>
| + | |
| − | | + | |
| − | <math display="inline">R\_Index=1</math><span style="text-align: center; font-size: 75%;"> //Right_Index</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>WHILE\left( R\_Index\leq NTer\right) \, Do</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>Hypotenuse=\, \sqrt{{\left( {S}_{(R\_Index)}.x-{S}_{(L\_Index)}.x\right) }^{2}+{\left( {S}_{(R\_Index)}.y-{S}_{(L\_Index)}.y\right) }^{2}}</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>IF\left( Hypotenuse\geq \, r\right) Then</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">//</span> <span style="text-align: center; font-size: 75%;">The point of intersection was found</span>
| + | |
| − | | + | |
| − | <math display="inline">{P}_{(NI)}={S}_{(R\_Index)}</math>
| + | |
| − | | + | |
| − | <math display="inline">NI=NI+1</math><span style="text-align: center; font-size: 75%;"> // Increase variable about Interval Numbers.</span>
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>L\_Index=R\_Index</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>EndIF</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>R\_Index=R\_Index+1</math>
| + | |
| − | |}
| + | |
| − | {|class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>EndWHILE</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">This procedure replaces the standard “compass dimension” method [29] that manually determining the new point </span> <math display="inline">{P}_{(k+1)}</math><span style="text-align: center; font-size: 75%;">, with a compass centered on </span> <math display="inline">{P}_{(k)}</math><span style="text-align: center; font-size: 75%;"> and radius </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 7. The areas of each of the </span> <math display="inline">NI</math><span style="text-align: center; font-size: 75%;"> trapezoids made up of the </span> <math display="inline">{P}_{(N)}</math><span style="text-align: center; font-size: 75%;">unequally spaced points are calculated, using equation (8). The total area is calculated as </span> <math display="inline">A\_Total=</math><math>\sum _{k=1}^{NI}{A}_{(k)}</math>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 8. Save iteration </span> <math display="inline">i</math><span style="text-align: center; font-size: 75%;">, radius </span> <math display="inline">r</math><span style="text-align: center; font-size: 75%;">and number of subintervals </span> <math display="inline">NI.</math>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 9. Return to the step 5 </span> <math display="inline">m</math><span style="text-align: center; font-size: 75%;"> times, repeating the process recursively, reducing the size of the radius </span> <math display="inline">\frac{{r}_{i}}{2}\,</math> <span style="text-align: center; font-size: 75%;">in each new iteration.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Step 10. Once the </span> <math display="inline">m</math><span style="text-align: center; font-size: 75%;"> iterations are finished, for each one the points are calculated as </span> <math display="inline">\left( Log\left( 1/{r}_{(i)}\right) ,\, Log\left( {Longitude}_{(i)}\right) \right)</math> <span style="text-align: center; font-size: 75%;">, where </span> <math display="inline">{Longitude}_{(i)}=</math><math>{r}_{(i)}.{NTer}_{(i)}</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | Step 11. The regression line is calculated by the least squares method. The slope of the line obtained is an estimate of the fractal dimension of the curve studied.
| + | |
| − | | + | |
| − | ==4. NUMERICAL RESULTS AND DISCUSSION.==
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">For the following function, due to its degree of fluctuation or variation, it is interesting to calculate its integral and fractal dimension in the interval </span> <math display="inline">\left[ \mathit{\boldsymbol{0,\, \, 4.2}}\right]</math> <span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;"
| + | |
| − | |-
| + | |
| − | |
| + | |
| − | {| style="text-align: center; margin:auto;width: 100%;"
| + | |
| − | |-
| + | |
| − | | <math>\mathit{\boldsymbol{f}}\left( \mathit{\boldsymbol{x}}\right) =</math><math>\left( \mathit{\boldsymbol{2}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{5}}}-\right. </math><math>\left. \mathit{\boldsymbol{25}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{4}}}+\right. </math><math>\left. \mathit{\boldsymbol{112}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{3}}}-\right. </math><math>\left. \mathit{\boldsymbol{212}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{2}}}+\right. </math><math>\left. \mathit{\boldsymbol{144}}\mathit{\boldsymbol{x}}+\mathit{\boldsymbol{3}}\right)</math>
| + | |
| − | |}
| + | |
| − | | style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(9)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">In order not to saturate the tables of hundreds of results, the calculations will be shown to up to the fourth fractal partition.</span>
| + | |
| − | | + | |
| − | Integrating the function analytically, the area under the curve is
| + | |
| − | | + | |
| − | ==
| + | |
| − | {| class="formulaSCP" style="width: 100%; text-align: center;"
| + | |
| − | |-
| + | |
| − | | <math>\int_{\mathit{\boldsymbol{0}}}^{\mathit{\boldsymbol{4.2}}}\mathit{\boldsymbol{f}}\left( \mathit{\boldsymbol{x}}\right) \mathit{\boldsymbol{dx=}}{\left[ \frac{\mathit{\boldsymbol{2}}}{\mathit{\boldsymbol{6}}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{6}}}\mathit{\boldsymbol{-}}\frac{\mathit{\boldsymbol{25}}}{\mathit{\boldsymbol{5}}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{5}}}\mathit{\boldsymbol{+}}\frac{\mathit{\boldsymbol{112}}}{\mathit{\boldsymbol{4}}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{4}}}\mathit{\boldsymbol{-}}\frac{\mathit{\boldsymbol{212}}}{\mathit{\boldsymbol{3}}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{3}}}\mathit{\boldsymbol{+}}\frac{\mathit{\boldsymbol{144}}}{\mathit{\boldsymbol{2}}}{\mathit{\boldsymbol{x}}}^{\mathit{\boldsymbol{2}}}\mathit{\boldsymbol{+3}}\mathit{\boldsymbol{x}}\right] }_{\mathit{\boldsymbol{0}}}^{\mathit{\boldsymbol{4.2}}}\mathit{\boldsymbol{=53.212}}</math>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The function is numerically integrated in the interval </span> <math display="inline">\left[ 0,\, 4.2\right]</math> <span style="text-align: center; font-size: 75%;">, with the partitions indicated in </span> <math display="inline">NT</math><span style="text-align: center; font-size: 75%;"> subintervals, all with the same upper base length. The results are shown in figure 1, and are presented in table 1.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">In Figure 1(a), a grid with </span> <math display="inline">NT=</math><math>1</math><span style="text-align: center; font-size: 75%;">fractal segment with an upper base </span> <math display="inline">r=</math><math>4.21</math><span style="text-align: center; font-size: 75%;"> is superimposed on the interval </span> <math display="inline">\left[ 0,\, 4.2\right]</math> <span style="text-align: center; font-size: 75%;">. The resulting area of 7.464 is a first approximation to the area of the region (table 1, case a,). In order to improve it, in Figure 1 (b) a variation is made, partitioning the curve with 3 fractal segments with an upper base </span> <math display="inline">r=</math><math>2.105</math><span style="text-align: center; font-size: 75%;">, for a calculated total area of 38,195, with which a new coating of the ''R'' region (table 1, case b,) is formed. In Figure 1(c), a variation is presented, partitioning the curve with 9 fractal segments with an upper base </span> <math display="inline">r=</math><math>1.052</math><span style="text-align: center; font-size: 75%;"> for a calculated total area of 51.221 (table 1, case c). In Figure 1 (d) a last variation is made by partitioning the curve with 21 fractal segments with an upper base </span> <math display="inline">r=</math><math>0.526</math><span style="text-align: center; font-size: 75%;">. The integral of the curve is calculated by adding the area of the subintervals and equal to 54.306 (table 1, case d).</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image1.png|228px]]
| + | |
| − | | + | |
| − | (1) <span style="text-align: center; font-size: 75%;">Iter. 1, NT=1, r= 4.21</span>
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image2.png|234px]]
| + | |
| − | | + | |
| − | (2) <span style="text-align: center; font-size: 75%;">Iter. 2. NT=3, r= 2.105</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image3.png|228px]]
| + | |
| − | | + | |
| − | (3) <span style="text-align: center; font-size: 75%;">Iter. 3, NT=9, r= 1.052</span>
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image4.png|234px]]
| + | |
| − | | + | |
| − | (4) <span style="text-align: center; font-size: 75%;">Iter. 4, NT=21, r= 0.526</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 1. Area calculated by modified trapezoid rule.</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | colspan='7' style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Case a: </span> <math display="inline">\mathit{\boldsymbol{r}}=</math><math>\boldsymbol{4}.\boldsymbol{211}</math><span style="text-align: center; font-size: 75%;">''' '''</span>
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"><br/></span><span style="text-align: center; font-size: 75%;">'''Area'''</span>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{1}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{1}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}\mathit{\boldsymbol{=}}\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{1}}}-</math><math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>\mathit{\boldsymbol{A=}}\frac{\mathit{\boldsymbol{h}}\left( {\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{0}}}\mathit{\boldsymbol{+}}{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{1}}}\right) }{\mathit{\boldsymbol{2}}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A1'''</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4.2</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.553</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4.2</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7.464</span>
| + | |
| − | |-
| + | |
| − | | colspan='7' style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''<br/>'''</span><span style="text-align: center; font-size: 75%;">Case b: </span> <math display="inline">\mathit{\boldsymbol{r}}=</math><math>\boldsymbol{2}.\boldsymbol{105}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A1'''</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.151</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.282</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,151</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.757</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A2'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.557</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.151</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.373</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20.282</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.406</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.453</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A3'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.656</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.557</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.806</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.373</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.098</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.981</span>
| + | |
| − | |-
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Total'''</span>
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">38.195</span>
| + | |
| − | |-
| + | |
| − | | colspan='7' style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Case c: </span> <math display="inline">\mathit{\boldsymbol{r}}=</math><math>\boldsymbol{1},\boldsymbol{052}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A1'''</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.067</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">11.729</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.067</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.493</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A2'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.152</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.067</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20.370</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.729</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.085</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.364</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A3'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.28</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.152</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">29.007</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20.370</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.128</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.160</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A4'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.032</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.280</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">22.906</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">29.007</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.751</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.519</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A5'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.284</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.032</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.468</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">22.906</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.252</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.727</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A6'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.617</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.284</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.2416</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.468</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.332</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.448</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A7'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.636</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.617</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.4341</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.2416</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.018</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.477</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A8'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.689</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.636</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.2926</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.4341</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.052</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.806</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A9'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.186</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.689</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.6556</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.2926</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.497</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.223</span>
| + | |
| − | |-
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Total'''</span>
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | | style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">51.221</span>
| + | |
| − | |-
| + | |
| − | | colspan='7' style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Case d: </span> <math display="inline">\mathit{\boldsymbol{r}}=</math><math>\boldsymbol{0}.\boldsymbol{526}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A1'''</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.032</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.394</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.032</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.166</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A2'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.067</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.032</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.729</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.394</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.035</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.334</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A3'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.107</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.067</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.114</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.729</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.040</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.556</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A4'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.153</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.107</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20.456</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.114</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.046</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.841</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A5'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.209</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.153</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.811</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20.456</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.056</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.267</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A6'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.282</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.209</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">29.106</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.811</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.073</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.967</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A7'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.405</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.282</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">33.336</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">29.106</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.123</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.840</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A8'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.795</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.405</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">30.414</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">33.336</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.390</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.431</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A9'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.934</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.795</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">26.208</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">30.414</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.139</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.935</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A10'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.058</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.934</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">21.979</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">26.208</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.124</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.011</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A11'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.182</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.058</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.756</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">21.979</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.123</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.463</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A12'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.314</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.182</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">13.553</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.756</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.131</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.066</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A13'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.466</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.314</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.403</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">13.553</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.151</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.744</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A14'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.669</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.466</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.398</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.403</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.202</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.502</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A15'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.132</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.669</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.333</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.398</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.462</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.021</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A16'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.484</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.132</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.568</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.333</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.351</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.742</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A17'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.783</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.484</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10.141</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.568</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.298</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.498</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A18'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.219</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.783</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.569</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10.141</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.435</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.951</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A19'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.602</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.219</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.593</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.569</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.382</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.244</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A20'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.840</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.602</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.728</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.593</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.237</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.823</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A21'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.076</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.840</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.854</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.728</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.236</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.894</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Total '''</span>
| + | |
| − | |
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">54.306</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Table 1. Calculating area with fractal segments</span></div>
| + | |
| − | | + | |
| − | ==4.1 Calculating fractal dimension==
| + | |
| − | | + | |
| − | Next, the fractal dimension of the function (9), obtained with the trapezoid rule modified method, is calculated by dividing the curve with fractal segments, results recorded in Table 2, and represented in Figure 2.
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">In table 2, the slope </span> <math display="inline">m</math><span style="text-align: center; font-size: 75%;"> of the regression line that makes up the points </span> <math display="inline">\left( Log\left( 1/r\right) ,\, Log\left( Longitude\right) \right)</math> <span style="text-align: center; font-size: 75%;"> is calculated. </span> <math display="inline">m</math><span style="text-align: center; font-size: 75%;">is an estimate of the fractal dimension of the area under study.</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Iter.'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Fig'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''r'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''NI'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{Longitude=r.NI}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(1/r)}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(Longitude)}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1(a)</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4.211</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.211</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,624</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0,624</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1(b)</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2.105</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.317</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,323</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,800</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1(c)</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1.058</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">9</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.475</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0,022</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,976</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1(d)</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.526</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">21</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.055</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,278</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,043</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.2632</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">43</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.318</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,579</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,053</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.1312</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">86</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.318</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0,880</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,053</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">-</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0658</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">171</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.252</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,181</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1,051</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Table 2. Fractal dimension partitioning the curve using fractal segments.</span>
| + | |
| − | | + | |
| − | Note that the results are shown with three additional iterations.
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image5.png|450px]] </span></div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 2. Regression line.</span></div>
| + | |
| − | | + | |
| − | The fractal dimension is a parameter that characterizes the sinuous or irregular of an object. For the case under study. The fractal dimension determines the variability of the results in the different iterations. This means that when there is no variability, an error is entered into the calculations.
| + | |
| − | | + | |
| − | Note that the critical point of the Figure 2 was reached in iteration 4, from which no significant differences can be seen. For that reason, the fractal dimension calculated, using the least squares method, was delimited between 1 to 4 iterations, to eliminate the error due to invariability of the results in the interactions of 5 to 7. The result was rounded to 3 decimals, considering that the error obtained by the adjustment with least squares is in this range.
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The slope </span> <math display="inline">m</math><span style="text-align: center; font-size: 75%;"> of the regression line is an estimate of the fractal dimension of the curve under study. With respect to the topological dimension </span> <math display="inline">{(D}_{top})</math><span style="text-align: center; font-size: 75%;"> of a line, the number of dimensions needed to fix a point on the line is 1. And because the fractional dimension </span> <math display="inline">{(D}_{frac})</math><span style="text-align: center; font-size: 75%;"> must be larger than its dimension topological, this is calculated as </span> <math display="inline">{D}_{frac}=</math><math>{D}_{top}+m=1+m.</math><span style="text-align: center; font-size: 75%;">The calculation shows a measure of the fractal dimension of the curve with a value of 1.476, which means that it is a curve with a medium degree of fluctuation.</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''4.2 Numerical integration using trapezoid rule for equally spaced data'''</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The previous function is integrated in the same interval</span> <math display="inline">\left[ a,b\right]</math> <span style="text-align: center; font-size: 75%;">, with the partitions indicated in </span> <math display="inline">n</math><span style="text-align: center; font-size: 75%;"> subintervals, all of the same length </span> <math display="inline">h=</math><math>(b-a)/n</math><span style="text-align: center; font-size: 75%;">, results reported in table 3 and represented in figure 3.</span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">In Figure 3(a), a grid with 1 trapezium in the interval </span> <math display="inline">\left[ 0,\, 4.2\right]</math> <span style="text-align: center; font-size: 75%;"> is superimposed on the region. The resulting area of 7.477, is a first approximation to the area of the region. In order to improve it, in Figure 3(b), a variation is made with 3 trapezoids in the intervals for a total calculated area of 32.502, with which forms a new superimposed of the R region. In Figure 3(c), a variation with 9 trapezoids is presented for a total calculated area of 52.281. In Figure 3(d), a last variation is made with 21 trapezoids. The approximate area under the curve is calculated by adding the area of the trapezoids and equal to 54.490 (table 4).</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image6.png|246px]]
| + | |
| − | | + | |
| − | (1) <span style="text-align: center; font-size: 75%;">n=1, h=4.2</span>
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image7.png|246px]]
| + | |
| − | | + | |
| − | (2) <span style="text-align: center; font-size: 75%;">n=3, h=1.4</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image8.png|258px]]
| + | |
| − | | + | |
| − | (3) <span style="text-align: center; font-size: 75%;">n=9, h=0.466, </span>
| + | |
| − | | style="text-align: center;vertical-align: top;width: 50%;"|[[Image:Draft_Tabares-Ospina_730588725-image9.png|258px]]
| + | |
| − | | + | |
| − | (4) <span style="text-align: center; font-size: 75%;">n=21, h=0.2</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 3. Area calculating by trapezoid rule equiespaced.</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | colspan='7' style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Case a: '''</span> <math display="inline">\mathit{\boldsymbol{r}}\boldsymbol{=0.200}</math><span style="text-align: center; font-size: 75%;">''' '''</span>
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"><br/></span><span style="text-align: center; font-size: 75%;">'''Area'''</span>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{1}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{1}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>{\mathit{\boldsymbol{h}}\mathit{\boldsymbol{=}}\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{1}}}-</math><math>{\mathit{\boldsymbol{h}}}_{\mathit{\boldsymbol{0}}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>\mathit{\boldsymbol{A=}}\frac{\mathit{\boldsymbol{h}}\left( {\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{0}}}\mathit{\boldsymbol{+}}{\mathit{\boldsymbol{f}}}_{\mathit{\boldsymbol{1}}}\right) }{\mathit{\boldsymbol{2}}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A1'''</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.200</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.00</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24.176</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.000</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.717</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A2'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">33.228</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.176</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.740</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A3'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">34.187</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">33.228</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.741</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A4'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">30.279</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">34.187</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.446</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A5'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">30.279</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.427</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A6'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.00</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.192</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.119</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A7'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.124</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.192</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.831</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A8'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.563</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.124</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.768</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A9'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.855</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.563</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.041</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A10'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.855</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.685</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A11'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.728</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.672</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A12'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.580</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.728</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.930</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A13'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.979</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.580</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.356</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A14'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10.311</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">7.979</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.829</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A15'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10.311</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.231</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A16'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.584</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.458</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A17'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.796</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.584</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.438</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A18'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.400</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.635</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.796</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.143</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A19'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.600</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.447</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9.635</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.608</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A20'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.800</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.999</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.447</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.944</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''A21'''</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.200</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.00</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.560</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.999</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.356</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Total '''</span>
| + | |
| − | |
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">54.490</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Table 3. Calculating area with equispace intervals.</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The area of the </span> <math display="inline">R</math><span style="text-align: center; font-size: 75%;"> region bounded by the function (9) is compared by the standard simple trapezoid rule, partitioning the curve with fractal segments'''.'''</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Method'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Result'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">\mathit{\boldsymbol{E}}\mathit{\boldsymbol{RP=}}\frac{\left| {\mathit{\boldsymbol{A}}}^{\mathit{\boldsymbol{\ast }}}\mathit{\boldsymbol{-A}}\right| }{\mathit{\boldsymbol{A}}}</math><span style="text-align: center; font-size: 75%;">'''*100'''</span>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Exact Area</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">53.212</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Our trapezoid rule method for fractally (unequally) spaced data </span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">54.306</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2.014</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Trapezoid rule for equallly spacing data</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">54.490</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2.345</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Table 4. Comparing different methods for numeric integration.</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The value of the function integral (9) used as a reference is </span> <math display="inline">{\mathit{\boldsymbol{A}}}^{\mathit{\boldsymbol{\ast }}}=</math><math>53.212</math><span style="text-align: center; font-size: 75%;">. Table 4 presents the results obtained, comparing the values of the integral of the curve studied, the relative percentage error (ERP). It is concluded that both methods are equally comparable</span>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''4.3 Calculating the numerical integration and fractal dimension of other polynomial functions'''</span></div>
| + | |
| − | | + | |
| − | Next, the validation of the proposed method is presented, with the following three polynomial functions, all with different degrees of irregularity or fluctuation.
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Polynomial functions'''</span>
| + | |
| − | | style="border-top: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;vertical-align: top;"|
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;vertical-align: top;"|<math>y={x}^{4}-13{x}^{3}+49{x}^{2}-49x+45</math>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(10)</span>
| + | |
| − | |-
| + | |
| − | | style="vertical-align: top;"|<math>y={\left( x-3\right) }^{5}-8{\left( x-3\right) }^{3}+10\left( x-3\right) +60</math>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(11)</span>
| + | |
| − | |-
| + | |
| − | | style="vertical-align: top;"|<math>y=-3{x}^{4}+20{x}^{3}-41{x}^{2}+22x+88</math>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(12)</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Table 5. Other functions.</span></div>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">The polynomial functions (10), (11) and (12) will be integrated in the intervals </span> <math display="inline">\left[ 0,\, 7\right]</math> <span style="text-align: center; font-size: 75%;">, </span> <math display="inline">\left[ 0,\, 5.9\right]</math> <span style="text-align: center; font-size: 75%;"> and </span> <math display="inline">\left[ 0,\, 4\right]</math> <span style="text-align: center; font-size: 75%;">. The exact areas are 274.869, 349.230 and 318.976. Dividing the functions using fractal procedure, in the iteration number 5, there are 41, 24 and 19 segments respectively.</span>
| + | |
| − | | + | |
| − | Figure 4 shows the partitions of the areas of the curves under study, using equally spaced data in the first case, and unequally spaced in the second.
| + | |
| − | | + | |
| − | <br/>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Trapezoid rule equally spaced'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Trapezoid rule fractally spaced'''</span>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image10.png|234px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(1) Function (10), NI = 41 equispaced.</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image11.png|234px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(2) Function (10), NI=41 unequally spaced</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image12.png|228px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(3) Function (11), NI=24 equispaced</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image13.png|234px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(4) Function (11), NI= 24 unequally spaced</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image14.png|228px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(5) Function (12), NI=19 equispaced </span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image15.png|228px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">(6) Function (12), NI=19 unequally spaced</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 4. Intervals equispaced and not equispaced for other functions.</span>
| + | |
| − | | + | |
| − | Tables 6, 7 and 8 show the calculation of the integral and fractal dimension of the polynomial functions under study (10), (11) and (12) respectively.
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Iter'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{r}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{NI}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Area'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{Longitude=r.NI}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(1/r)}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(Longitude)}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.000</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">315.000</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.000</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.845</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.845</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.500</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">302.463</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.544</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.146</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.750</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">287.793</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.000</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.243</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.146</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.875</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">20</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">277.548</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.500</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.057</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.243</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.437</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">41</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">274.791</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.937</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.359</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.253</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.218</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">83</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">274.832</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">18.156</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.660</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.259</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.109</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">165</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">274.699</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">18.046</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.961</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.256</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Table 6. Fractal dimension, polynomial function </span> <math display="inline">y=</math><math>{x}^{4}-13{x}^{3}+49{x}^{2}-49x+45</math><span style="text-align: center; font-size: 75%;">.</span>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Iter'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{r}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{NT}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Area'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{Longitude=r.NT}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(1/r)}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(Longitude)}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.218</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">300.898</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.218</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.914</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.914</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.109</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">270.113</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.218</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.613</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.914</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.054</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">339.050</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">10.273</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.312</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.011</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.027</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">346.767</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">11.300</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.011</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.053</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.513</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">348.170</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.327</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.289</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.090</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.256</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">48</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">347.850</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.327</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.590</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.090</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.128</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">96</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">348.443</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.327</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.891</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.090</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Table 7. Fractal dimension, polynomial function </span> <math display="inline">y=</math><math>{\left( x-3\right) }^{5}-8{\left( x-3\right) }^{3}+10\left( x-3\right) +60</math>
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Iter'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{r}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{NT}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Area'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{Longitude=r.NT}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(1/r)}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{LOG(Longitude)}}</math>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.693</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">239.831</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.693</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.671</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.671</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.346</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">294.116</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.693</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.370</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.671</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">3</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.173</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">306.375</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.693</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">-0.069</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.671</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">4</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.586</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">9</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">315.353</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.279</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.231</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.722</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">5</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.293</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">318.351</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.573</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.532</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.746</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">6</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.146</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">38</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">318.650</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.573</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.833</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.746</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">7</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.073</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">76</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">318.942</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.573</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.134</span>
| + | |
| − | | style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.746</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Table 8. Fractal dimension, polynomial function </span> <math display="inline">y=</math><math>-3{x}^{4}+20{x}^{3}-41{x}^{2}+22x+88</math>
| + | |
| − | | + | |
| − | Next, the fractal dimension of the functions (10, 11 and 12), are calculated by dividing the curve with fractal segments, results are represented in Figures 5, 6 and 7.
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image16.png|282px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 5. Regression line, function (10).</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image17.png|282px]] </span>
| + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 6. Regression line, function (11).</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;"> [[Image:Draft_Tabares-Ospina_730588725-image18.png|318px]] </span></div>
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">Figure 7. Regression line, function (12).</span></div>
| + | |
| − | | + | |
| − | Note that in figures 5, 6 and 7 the critical points were reached in iterations 4, 5 and 5 respectively, from which no significant differences can be seen. For that reason, the fractal dimension calculated, using the least squares method, was delimited between 1 to 4, 1 to 5, 3 to 5 iterations respectively, to eliminate the error due to invariability of the results in the interactions of 5 to 7 (Figure 5), 6 to 7 (Figure 6), 1 to 2 and 6 to 7 (Figure 7), respectively.
| + | |
| − | | + | |
| − | Table 9 presents the results obtained, comparing the values of the integral of the curve studied, the relative percentage error (ERP), using the methods Trapezoid rule standard and Trapezoid rule modified, partitioning the curve with fractal segments.
| + | |
| − | | + | |
| − | As evidenced in Table 9, the method proposed in this article, dividing the trapezoids from the fractal fragmentation of its upper side, presents results comparable with those obtained with the traditional composite trapezoid method, in which the lower sides of the sets are equally spaced.
| + | |
| − | | + | |
| − | {| style="width: 100%;border-collapse: collapse;"
| + | |
| − | |-
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Function'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Exact Area'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{NI}}</math>
| + | |
| − | | + | |
| − | | + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Area equally spaced'''</span>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>\mathit{\boldsymbol{ERP}}</math>
| + | |
| − | | style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Area fractally spaced'''</span>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<math>\mathit{\boldsymbol{ERP}}</math>
| + | |
| − | | style="border: 1pt solid black;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">'''Fractal Dim.'''</span>
| + | |
| − | |-
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(10)</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">274.869</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">41</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">275.340</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.171</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">274.791</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.028</span>
| + | |
| − | | style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1.396</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(11)</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">349.230</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">24</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">349.040</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.054</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">348.170</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.304</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1.162</span>
| + | |
| − | |-
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">(12)</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">318.976</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">19</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">318.430</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.171</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">318.351</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.196</span>
| + | |
| − | | style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">1.124</span>
| + | |
| − | |}
| + | |
| − | | + | |
| − | | + | |
| − | <span style="text-align: center; font-size: 75%;">Table 9. Area comparison.''' '''</span>
| + | |
| − | | + | |
| − | With respect to the fractal dimension of the curves studied, their values are 1.396, 1.162 and 1.124 respectively, an index that effectively reveals their degree of fluctuation. The polynomial function (10) has a maximum point over two deep basins at different slopes. With respect to (10), the direction of the function (11) changes more smoothly between its two points of maximum and minimum (all of them almost at the same level), respectively. Finally (12) it is the function with smoother changes, with negative bias to the left of the study interval.
| + | |
| − | | + | |
| − | ==5. Conclusions==
| + | |
| − | | + | |
| − | In the present paper the working hypothesis was validated: The degree of fluctuation of a function is calculated using fractal dimension method. The fractal segments obtained also may be used to calculate the area under the curve, integrating the function using standard trapezoid rule, with fractally spaced data, from which other relevant observations may be obtained from its analysis, such as:
| + | |
| − | | + | |
| − | (1) The great advantage of the proposed method is that it unifies the calculation of the numerical integration of a function, together with its corresponding fractal dimension, a task not performed with conventional numerical integration methods.
| + | |
| − | | + | |
| − | (2) To calculate the fractal dimension, the standard method uses a compass that with a specific radius, an arc is drawn to find the point that intercepts the curve. In this article a compass was not used, because an algorithm was developed to find the point that intercepts the curve studied with a specific radius.
| + | |
| − | | + | |
| − | (3) The fractal dimension is a quantitative estimator of spatial complexity.
| + | |
| − | | + | |
| − | (4) It is confirmed that the fractal dimension satisfactorily describes the degree of irregularity of a function, an impossible measure to make using Euclidean (classical) geometry that is limited to regular shapes.
| + | |
| − | | + | |
| − | (5) The calculation of the fractal dimension of the studied surface cannot be obtained exactly but estimated, being very sensitive to numerical or experimental noise, and particularly to the limitations in the amount of data
| + | |
| − | | + | |
| − | <div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
| + | |
| − | <span style="text-align: center; font-size: 75%;">'''Acknowledgments'''</span></div>
| + | |
| − | | + | |
| − | This work was supported by the Agencia de Educación Superior de Medellín (SAPIENCIA) under the specific agreement celebrated with the Institución Universitaria Pascual Bravo. The project is part of the Energy System Doctorate Program of the Universidad Nacional de Colombia, Sede Medellín, Facultad de Minas.
| + | |
| − | | + | |
| − | '''Conflicts of Interest:''' The authors declare no conflict of interest.
| + | |
| − | | + | |
| − | ==References==
| + | |
| − | | + | |
| − | [1] M.F. Barnsley, R.L. Devaney, B.B. Mandelbrot, H.-O. Peitgen, D. Saupe, R.F. Voss, The Science of Fractal Images, Springer New York, New York, NY, 1988. doi:10.1007/978-1-4612-3784-6.
| + | |
| − | | + | |
| − | [2] S.H. Strogatz, Nonlinear Dynamics and Chaos, CRC Press, 2018. doi:10.1201/9780429492563.
| + | |
| − | | + | |
| − | [3] G.A. Losa, Fractals and Their Contribution to Biology and Medicine, Medicographia. 34 (2012) 365–374.
| + | |
| − | | + | |
| − | [4] T.A. Garcia, G.A. Tamura Ozaki, R.C. Castoldi, T.E. Koike, R.C. Trindade Camargo, J.C. Silva Camargo Filho, Fractal dimension in the evaluation of different treatments of muscular injury in rats, Tissue Cell. 54 (2018) 120–126. doi:10.1016/j.tice.2018.08.014.
| + | |
| − | | + | |
| − | [5] J.O. Rodríguez V, S.E. Prieto B, S.C. Correa H, M.Y. Soracipa M, L.R. Mendez P, H.J. Bernal C, N.C. Hoyos O, L.P. Valero, A. Velasco R, E. Bermudez, Nueva metodología de evaluación del Holter basada en los sistemas dinámicos y la geometría fractal: confirmación de su aplicabilidad a nivel clínico, Rev. La Univ. Ind. Santander. Salud. 48 (2016) 27–36. doi:10.18273/revsal.v48n1-2016003.
| + | |
| − | | + | |
| − | [6] N. Popovic, M. Radunovic, J. Badnjar, T. Popovic, Fractal dimension and lacunarity analysis of retinal microvascular morphology in hypertension and diabetes, Microvasc. Res. 118 (2018) 36–43. doi:10.1016/j.mvr.2018.02.006.
| + | |
| − | | + | |
| − | [7] J. de D. Hernández Velázquez, S. Mejía-Rosales, A. Gama Goicochea, Fractal properties of biophysical models of pericellular brushes can be used to differentiate between cancerous and normal cervical epithelial cells, Colloids Surfaces B Biointerfaces. 170 (2018) 572–577. doi:10.1016/j.colsurfb.2018.06.059.
| + | |
| − | | + | |
| − | [8] P. Moon, J. Muday, S. Raynor, J. Schirillo, C. Boydston, M.S. Fairbanks, R.P. Taylor, Fractal images induce fractal pupil dilations and constrictions, Int. J. Psychophysiol. 93 (2014) 316–321. doi:10.1016/j.ijpsycho.2014.06.013.
| + | |
| − | | + | |
| − | [9] B.B. Mandelbrot, R.L. Hudson, The (mis) Behaviour of Markets: A Fractal View of Risk, Ruin and Reward, Profile Books, London, 2004.
| + | |
| − | | + | |
| − | [10] R. Kumar, P.N. Chaubey, On the design of tree-type ultra wideband fractal Antenna for DS-CDMA system, J. Microwaves, Optoelectron. Electromagn. Appl. 11 (2012) 107–121. doi:10.1590/S2179-10742012000100009.
| + | |
| − | | + | |
| − | [11] Y.-J. Ma, M.-Y. Zhai, Fractal and multi-fractal features of the broadband power line communication signals, Comput. Electr. Eng. 72 (2018) 566–576. doi:10.1016/j.compeleceng.2018.01.025.
| + | |
| − | | + | |
| − | [12] D. Ye, M. Dai, Y. Sun, W. Su, Average weighted receiving time on the non-homogeneous double-weighted fractal networks, Phys. A Stat. Mech. Its Appl. 473 (2017) 390–402. doi:10.1016/j.physa.2017.01.013.
| + | |
| − | | + | |
| − | [13] H. Tabares-Ospina, F. Angulo, M. Osorio, A new methodology to analyze the dynamic of daily power demand with attractors into the Mandelbrot set, Fractals. 28 (2020) 2050003–1, 2050003–9. doi:10.1142/S0218348X20500036.
| + | |
| − | | + | |
| − | [14] J. Colomer, A. Naranjo, V. Janvier, T. Mossi, Evaluation of fractal dimension effectiveness for damage detection in retinal background, J. Comput. Appl. Math. 337 (2018) 341–353.
| + | |
| − | | + | |
| − | [15] F. Guosheng, X. Weimin, Y. Feng, Y. Yi, W. Yuan, Surface fractal dimension of bentonite affected by long-term corrosion in alkaline solution, Appl. Clay Sci. 175 (2019) 94–101.
| + | |
| − | | + | |
| − | [16] R. Swain, S. Roy, P. Mukherjee, B. Sawkar, Fractal dimension and its translation into a model of gold spatial proxy, Ore Geol. Rev. 110 (2019) 1–10.
| + | |
| − | | + | |
| − | [17] X. Yongfu, Fractal dimension of demolition waste fragmentation and its implication of compactness, Powder Technol. 339 (2018) 922–929.
| + | |
| − | | + | |
| − | [18] E. Martos, F. Lapuerta, M. Expósito, Juan. Sanmiguel-Rojas, Overestimation of the fractal dimension from projections of soot agglomerates, Powder Technol. 311 (2017) 528–536.
| + | |
| − | | + | |
| − | [19] F. Yang, X. Zhang, R. Ma, S. Yang, X, Wand, Fractal dimension of concrete meso-structure based on X-ray computed tomography, Powder Technol. 350 (2019) 91–99.
| + | |
| − | | + | |
| − | [20] J. Jiménez, A. Anzola, J. Jimenez-Triana, Pedestrian counting estimation based on fractal dimension, Heliyon. (2019). doi:10.1016/j.heliyon.2019. e01449.
| + | |
| − | | + | |
| − | [21] L. Chen, Y. Huang, Modeling growth curve of fractal dimension of urban form of Beijing, Physica A. 523 (2019) 1038–1056.
| + | |
| − | | + | |
| − | [22] E. Pashminehazar, R. Kharaghani, A. Tsotsas, Determination of fractal dimension and prefactor of agglomerates with irregular structure, Powder Technol. 343 (2019) 765–774.
| + | |
| − | | + | |
| − | [23] M.S.M. Navascués, Numerical integration of affine fractal functions, J. Comput. Appl. Math. 252 (2013) 169–176.
| + | |
| − | | + | |
| − | [24] D. Raghavendra, B. Dutt, Computing Fractal Dimension of Signals using Multiresolution Boxcounting Method, Int. J. Inf. Math. Sci. (2010) 50–65.
| + | |
| − | | + | |
| − | [25] R. Chapra, S. Canale, Numerical methods for engineers, Fifty, The McGraw-Hill Companies, Mexico, 2006.
| + | |
| − | | + | |
| − | [26] S. Mamun-Ur-Rashid, M. Hossain, M. Parvin, Numerical Integration Schemes for Unequal Data Spacing, Am. J. Appl. Math. 2 (2017) 48–56.
| + | |
| − | | + | |
| − | [27] K. Mathews, J. Fink, Numerical methods using MATLAB, third, Madrid, 2000.
| + | |
| − | | + | |
| − | [28] R. Malmquist, J. Stricharts, Numerical integration for fractal measures, J. Fractal Geom. 5 (2018) 165–226.
| + | |
| − | | + | |
| − | [29] G. Ranjan, S. Mishra, J, Palai, Analysing roughness of surface through fractal dimension: A review, Image Vis. Comput. 89 (2019) 21–34.
| + | |
