diff options
| author | Jay Belanger | 2007-08-04 04:03:50 +0000 |
|---|---|---|
| committer | Jay Belanger | 2007-08-04 04:03:50 +0000 |
| commit | de6a5d7ba9b76912bbc4a3faddd59ecbe903a923 (patch) | |
| tree | e3ed4f1ad93ddc5a7f54b97d1588ea660c4d3504 | |
| parent | 8ab437fd8ad9a1a50d5b7358ec5b67095228f497 (diff) | |
| download | emacs-de6a5d7ba9b76912bbc4a3faddd59ecbe903a923.tar.gz emacs-de6a5d7ba9b76912bbc4a3faddd59ecbe903a923.zip | |
(Curve Fitting): Mention plotting the curves.
(Standard Nonlinear Models): Add additional models.
(Curve fitting details): Mention the Levenberg-Marquardt method.
| -rw-r--r-- | man/ChangeLog | 3 | ||||
| -rw-r--r-- | man/calc.texi | 32 |
2 files changed, 30 insertions, 5 deletions
diff --git a/man/ChangeLog b/man/ChangeLog index 18cd4b72bc6..662250740b9 100644 --- a/man/ChangeLog +++ b/man/ChangeLog | |||
| @@ -3,6 +3,9 @@ | |||
| 3 | * calc.texi (Basic Graphics): Mention the graphing of error | 3 | * calc.texi (Basic Graphics): Mention the graphing of error |
| 4 | forms. | 4 | forms. |
| 5 | (Graphics Options): Mention how `g s' handles error forms. | 5 | (Graphics Options): Mention how `g s' handles error forms. |
| 6 | (Curve Fitting): Mention plotting the curves. | ||
| 7 | (Standard Nonlinear Models): Add additional models. | ||
| 8 | (Curve fitting details): Mention the Levenberg-Marquardt method. | ||
| 6 | 9 | ||
| 7 | 2007-08-01 Alan Mackenzie <acm@muc.de> | 10 | 2007-08-01 Alan Mackenzie <acm@muc.de> |
| 8 | 11 | ||
diff --git a/man/calc.texi b/man/calc.texi index 9e50629a3b2..95737d6106a 100644 --- a/man/calc.texi +++ b/man/calc.texi | |||
| @@ -23962,7 +23962,11 @@ such as @expr{y = m x + b} where @expr{m} and @expr{b} are parameters | |||
| 23962 | to be determined. For a typical set of measured data there will be | 23962 | to be determined. For a typical set of measured data there will be |
| 23963 | no single @expr{m} and @expr{b} that exactly fit the data; in this | 23963 | no single @expr{m} and @expr{b} that exactly fit the data; in this |
| 23964 | case, Calc chooses values of the parameters that provide the closest | 23964 | case, Calc chooses values of the parameters that provide the closest |
| 23965 | possible fit. | 23965 | possible fit. The model formula can be entered in various ways after |
| 23966 | the key sequence @kbd{a F} is pressed. If the letter @kbd{P} | ||
| 23967 | is pressed after @kbd{a F} but before the model description is entered, | ||
| 23968 | the data as well as the model formula will be plotted after the formula | ||
| 23969 | is determined. | ||
| 23966 | 23970 | ||
| 23967 | @menu | 23971 | @menu |
| 23968 | * Linear Fits:: | 23972 | * Linear Fits:: |
| @@ -24454,6 +24458,18 @@ Quadratic. @mathit{a + b (x-c)^2 + d (x-e)^2}. | |||
| 24454 | Gaussian. | 24458 | Gaussian. |
| 24455 | @texline @math{{a \over b \sqrt{2 \pi}} \exp\left( -{1 \over 2} \left( x - c \over b \right)^2 \right)}. | 24459 | @texline @math{{a \over b \sqrt{2 \pi}} \exp\left( -{1 \over 2} \left( x - c \over b \right)^2 \right)}. |
| 24456 | @infoline @mathit{(a / b sqrt(2 pi)) exp(-0.5*((x-c)/b)^2)}. | 24460 | @infoline @mathit{(a / b sqrt(2 pi)) exp(-0.5*((x-c)/b)^2)}. |
| 24461 | @item s | ||
| 24462 | Logistic @emph{s} curve. | ||
| 24463 | @texline @math{a/(1+e^{b(x-c)})}. | ||
| 24464 | @infoline @mathit{a/(1 + exp(b (x - c)))}. | ||
| 24465 | @item b | ||
| 24466 | Logistic bell curve. | ||
| 24467 | @texline @math{ae^{b(x-c)}/(1+e^{b(x-c)})^2}. | ||
| 24468 | @infoline @mathit{a exp(b (x - c))/(1 + exp(b (x - c)))^2}. | ||
| 24469 | @item o | ||
| 24470 | Hubbert linearization. | ||
| 24471 | @texline @math{{y \over x} = a(1-x/b)}. | ||
| 24472 | @infoline @mathit{(y/x) = a (1 - x/b)}. | ||
| 24457 | @end table | 24473 | @end table |
| 24458 | 24474 | ||
| 24459 | All of these models are used in the usual way; just press the appropriate | 24475 | All of these models are used in the usual way; just press the appropriate |
| @@ -24462,8 +24478,9 @@ result will be a formula as shown in the above table, with the best-fit | |||
| 24462 | values of the parameters substituted. (You may find it easier to read | 24478 | values of the parameters substituted. (You may find it easier to read |
| 24463 | the parameter values from the vector that is placed in the trail.) | 24479 | the parameter values from the vector that is placed in the trail.) |
| 24464 | 24480 | ||
| 24465 | All models except Gaussian and polynomials can generalize as shown to any | 24481 | All models except Gaussian, logistics, Hubbert and polynomials can |
| 24466 | number of independent variables. Also, all the built-in models have an | 24482 | generalize as shown to any number of independent variables. Also, all |
| 24483 | the built-in models except for the logistic and Hubbert curves have an | ||
| 24467 | additive or multiplicative parameter shown as @expr{a} in the above table | 24484 | additive or multiplicative parameter shown as @expr{a} in the above table |
| 24468 | which can be replaced by zero or one, as appropriate, by typing @kbd{h} | 24485 | which can be replaced by zero or one, as appropriate, by typing @kbd{h} |
| 24469 | before the model key. | 24486 | before the model key. |
| @@ -24603,7 +24620,7 @@ to convert the model into this form. For example, if the model | |||
| 24603 | is @expr{a + b x + c x^2}, then @expr{f(x) = 1}, @expr{g(x) = x}, | 24620 | is @expr{a + b x + c x^2}, then @expr{f(x) = 1}, @expr{g(x) = x}, |
| 24604 | and @expr{h(x) = x^2} are suitable functions. | 24621 | and @expr{h(x) = x^2} are suitable functions. |
| 24605 | 24622 | ||
| 24606 | For other models, Calc uses a variety of algebraic manipulations | 24623 | For most other models, Calc uses a variety of algebraic manipulations |
| 24607 | to try to put the problem into the form | 24624 | to try to put the problem into the form |
| 24608 | 24625 | ||
| 24609 | @smallexample | 24626 | @smallexample |
| @@ -24662,7 +24679,12 @@ The Gaussian model looks quite complicated, but a closer examination | |||
| 24662 | shows that it's actually similar to the quadratic model but with an | 24679 | shows that it's actually similar to the quadratic model but with an |
| 24663 | exponential that can be brought to the top and moved into @expr{Y}. | 24680 | exponential that can be brought to the top and moved into @expr{Y}. |
| 24664 | 24681 | ||
| 24665 | An example of a model that cannot be put into general linear | 24682 | The logistic models cannot be put into general linear form. For these |
| 24683 | models, and the Hubbert linearization, Calc computes a rough | ||
| 24684 | approximation for the parameters, then uses the Levenberg-Marquardt | ||
| 24685 | iterative method to refine the approximations. | ||
| 24686 | |||
| 24687 | Another model that cannot be put into general linear | ||
| 24666 | form is a Gaussian with a constant background added on, i.e., | 24688 | form is a Gaussian with a constant background added on, i.e., |
| 24667 | @expr{d} + the regular Gaussian formula. If you have a model like | 24689 | @expr{d} + the regular Gaussian formula. If you have a model like |
| 24668 | this, your best bet is to replace enough of your parameters with | 24690 | this, your best bet is to replace enough of your parameters with |