Engineering Statistics Handbook Episode 6 Part 6 ppt

4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (3 of 8) [5/1/2006 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (4 of 8) [5/1/2006 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (5 of 8) [5/1/2006 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (6 of 8) [5/1/2006 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (7 of 8) [5/1/2006 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd8113.htm (8 of 8) [5/1/2006 10:23:02 AM] Advantages Rational function models have the following advantages. Rational function models have a moderately simple form.1. Rational function models are a closed family. As with polynomial models, this means that rational function models are not dependent on the underlying metric. 2. Rational function models can take on an extremely wide range of shapes, accommodating a much wider range of shapes than does the polynomial family. 3. Rational function models have better interpolatory properties than polynomial models. Rational functions are typically smoother and less oscillatory than polynomial models. 4. Rational functions have excellent extrapolatory powers. Rational functions can typically be tailored to model the function not only within the domain of the data, but also so as to be in agreement with theoretical/asymptotic behavior outside the domain of interest. 5. Rational function models have excellent asymptotic properties. Rational functions can be either finite or infinite for finite values, or finite or infinite for infinite values. Thus, rational functions can easily be incorporated into a rational function model. 6. Rational function models can often be used to model complicated structure with a fairly low degree in both the numerator and denominator. This in turn means that fewer coefficients will be required compared to the polynomial model. 7. Rational function models are moderately easy to handle computationally. Although they are nonlinear models, rational function models are a particularly easy nonlinear models to fit. 8. Disadvantages Rational function models have the following disadvantages. The properties of the rational function family are not as well known to engineers and scientists as are those of the polynomial family. The literature on the rational function family is also more limited. Because the properties of the family are often not well understood, it can be difficult to answer the following modeling question: Given that data has a certain shape, what values should be chosen for the degree of the numerator and the degree on the denominator? 1. Unconstrained rational function fitting can, at times, result in undesired nusiance asymptotes (vertically) due to roots in the denominator polynomial. The range of values affected by the function "blowing up" may be quite narrow, but such asymptotes, when they occur, are a nuisance for local interpolation in the neighborhood of the asymptote point. These asymptotes are easy to 2. 4.8.1.2. Rational Functions http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd812.htm (2 of 4) [5/1/2006 10:23:03 AM] detect by a simple plot of the fitted function over the range of the data. Such asymptotes should not discourage you from considering rational function models as a choice for empirical modeling. These nuisance asymptotes occur occasionally and unpredictably, but the gain in flexibility of shapes is well worth the chance that they may occur. General Properties of Rational Functions The following are general properties of rational functions. If the numerator and denominator are of the same degree (n=m), then y = a n /b m is a horizontal asymptote of the function. ● If the degree of the denominator is greater than the degree of the numerator, then y = 0 is a horizontal asymptote. ● If the degree of the denominator is less than the degree of the numerator, then there are no horizontal asymptotes. ● When x is equal to a root of the denominator polynomial, the denominator is zero and there is a vertical asymptote. The exception is the case when the root of the denominator is also a root of the numerator. However, for this case we can cancel a factor from both the numerator and denominator (and we effectively have a lower-degree rational function). ● Starting Values for Rational Function Models One common difficulty in fitting nonlinear models is finding adequate starting values. A major advantage of rational function models is the ability to compute starting values using a linear least squares fit. To do this, choose p points from the data set, with p denoting the number of parameters in the rational model. For example, given the linear/quadratic model we need to select four representative points. We then perform a linear fit on the model Here, p n and p d are the degrees of the numerator and denominator, respectively, and the and Y contain the subset of points, not the full data set. The estimated coefficients from this fit made using the linear least squares algorithm are used as the starting values for fitting the nonlinear model to the full data set. Note: This type of fit, with the response variable appearing on both sides of the function, should only be used to obtain starting values for the nonlinear fit. The statistical properties of models like this are not well understood. 4.8.1.2. Rational Functions http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd812.htm (3 of 4) [5/1/2006 10:23:03 AM] The subset of points should be selected over the range of the data. It is not critical which points are selected, although you should avoid points that are obvious outliers. Example The thermal expansion of copper case study contains an example of fitting a rational function model. Specific Rational Functions Constant / Linear Rational Function1. Linear / Linear Rational Function2. Linear / Quadratic Rational Function3. Quadratic / Linear Rational Function4. Quadratic / Quadratic Rational Function5. Cubic / Linear Rational Function6. Cubic / Quadratic Rational Function7. Linear / Cubic Rational Function8. Quadratic / Cubic Rational Function9. Cubic / Cubic Rational Function10. Determining m and n for Rational Function Models11. 4.8.1.2. Rational Functions http://www.itl.nist.gov/div898/handbook/pmd/section8/pmd812.htm (4 of 4) [5/1/2006 10:23:03 AM] [...]... Examples: http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (2 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1 Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (3 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1 Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (4 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1 Constant... 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1 Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (5 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1 Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (6 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.2 Linear / Linear Rational Function Function Family: Rational Statistical Type: Nonlinear... Horizontal asymptote at: and vertical asymptote at: Additional Examples: http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8122.htm (2 of 6) [5/1/20 06 10:23:05 AM] 4.8.1.2.2 Linear / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8122.htm (3 of 6) [5/1/20 06 10:23:05 AM] . Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (5 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1. Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (6 of 6) [5/1/20 06. Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (2 of 6) [5/1/20 06 10:23:04 AM] 4.8.1.2.1. Constant / Linear Rational Function http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8121.htm (3 of 6) [5/1/20 06. Polynomial http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8113.htm (5 of 8) [5/1/20 06 10:23:02 AM] 4.8.1.1.3. Cubic Polynomial http://www.itl.nist.gov/div898 /handbook/ pmd/section8/pmd8113.htm (6 of 8) [5/1/20 06 10:23:02