aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorJay Belanger2007-08-04 04:03:50 +0000
committerJay Belanger2007-08-04 04:03:50 +0000
commitde6a5d7ba9b76912bbc4a3faddd59ecbe903a923 (patch)
treee3ed4f1ad93ddc5a7f54b97d1588ea660c4d3504
parent8ab437fd8ad9a1a50d5b7358ec5b67095228f497 (diff)
downloademacs-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/ChangeLog3
-rw-r--r--man/calc.texi32
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
72007-08-01 Alan Mackenzie <acm@muc.de> 102007-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
23962to be determined. For a typical set of measured data there will be 23962to be determined. For a typical set of measured data there will be
23963no single @expr{m} and @expr{b} that exactly fit the data; in this 23963no single @expr{m} and @expr{b} that exactly fit the data; in this
23964case, Calc chooses values of the parameters that provide the closest 23964case, Calc chooses values of the parameters that provide the closest
23965possible fit. 23965possible fit. The model formula can be entered in various ways after
23966the key sequence @kbd{a F} is pressed. If the letter @kbd{P}
23967is pressed after @kbd{a F} but before the model description is entered,
23968the data as well as the model formula will be plotted after the formula
23969is 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}.
24454Gaussian. 24458Gaussian.
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
24462Logistic @emph{s} curve.
24463@texline @math{a/(1+e^{b(x-c)})}.
24464@infoline @mathit{a/(1 + exp(b (x - c)))}.
24465@item b
24466Logistic 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
24470Hubbert 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
24459All of these models are used in the usual way; just press the appropriate 24475All 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
24462values of the parameters substituted. (You may find it easier to read 24478values of the parameters substituted. (You may find it easier to read
24463the parameter values from the vector that is placed in the trail.) 24479the parameter values from the vector that is placed in the trail.)
24464 24480
24465All models except Gaussian and polynomials can generalize as shown to any 24481All models except Gaussian, logistics, Hubbert and polynomials can
24466number of independent variables. Also, all the built-in models have an 24482generalize as shown to any number of independent variables. Also, all
24483the built-in models except for the logistic and Hubbert curves have an
24467additive or multiplicative parameter shown as @expr{a} in the above table 24484additive or multiplicative parameter shown as @expr{a} in the above table
24468which can be replaced by zero or one, as appropriate, by typing @kbd{h} 24485which can be replaced by zero or one, as appropriate, by typing @kbd{h}
24469before the model key. 24486before the model key.
@@ -24603,7 +24620,7 @@ to convert the model into this form. For example, if the model
24603is @expr{a + b x + c x^2}, then @expr{f(x) = 1}, @expr{g(x) = x}, 24620is @expr{a + b x + c x^2}, then @expr{f(x) = 1}, @expr{g(x) = x},
24604and @expr{h(x) = x^2} are suitable functions. 24621and @expr{h(x) = x^2} are suitable functions.
24605 24622
24606For other models, Calc uses a variety of algebraic manipulations 24623For most other models, Calc uses a variety of algebraic manipulations
24607to try to put the problem into the form 24624to 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
24662shows that it's actually similar to the quadratic model but with an 24679shows that it's actually similar to the quadratic model but with an
24663exponential that can be brought to the top and moved into @expr{Y}. 24680exponential that can be brought to the top and moved into @expr{Y}.
24664 24681
24665An example of a model that cannot be put into general linear 24682The logistic models cannot be put into general linear form. For these
24683models, and the Hubbert linearization, Calc computes a rough
24684approximation for the parameters, then uses the Levenberg-Marquardt
24685iterative method to refine the approximations.
24686
24687Another model that cannot be put into general linear
24666form is a Gaussian with a constant background added on, i.e., 24688form 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
24668this, your best bet is to replace enough of your parameters with 24690this, your best bet is to replace enough of your parameters with