close
Warning:
Can't synchronize with repository "(default)" (/var/svn/mms does not appear to be a Subversion repository.). Look in the Trac log for more information.
- Timestamp:
-
Sep 14, 2010, 11:04:50 AM (14 years ago)
- Author:
-
Pedro Gea
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v1
|
v2
|
|
43 | 43 | B·Z == C |
44 | 44 | }}} |
45 | | viene dad por: |
| 45 | viene dada por: |
46 | 46 | {{{ |
47 | 47 | Z = Sigma2·B' · (B·Sigma2·B')**-1 · (C-B·Mu) + Mu |
… |
… |
|
112 | 112 | y también como valores más probables (modas) |
113 | 113 | la solución obtenida minimizando la distancia de Mahalanobis. |
| 114 | |
| 115 | === Nota matemática === |
| 116 | |
| 117 | En las ecuaciones y cálculos anteriores se han utilizado las siguientes relaciones: |
| 118 | {{{ |
| 119 | fNormal(mu, sigma, x+b) == fNormal(mu-b, sigma, x) |
| 120 | fNormal(mu, sigma, -x) == fNormal(-mu, sigma, x) |
| 121 | fNormal(mu, sigma, a*x) == fNormal(mu/a, sigma/a, x)/a ; con a>0 |
| 122 | }}} |
| 123 | donde: |
| 124 | {{{ |
| 125 | fNormal(mu, sigma, x) := 1/(Sqrt(2*Pi)*Sigma) * Exp(-(1/2)*((x-mu)/sigma)**2) |
| 126 | }}} |
| 127 | y también el siguiente resultado que nos permite multiplicar dos funciones de densidad normales: |
| 128 | {{{ |
| 129 | fNormal(muA, sigmaA, x) * fNormal(muB, sigmaB, x) == |
| 130 | fNormal((muA*sigmaB**2 + muB*sigmaA**2)/sigmaAB**2, sigmaA*sigmaB/sigmaAB, x) * fNormal(muA-muB, sigmaAB, 0) |
| 131 | }}} |
| 132 | donde: |
| 133 | {{{ |
| 134 | sigmaAB = Sqrt(sigmaA**2 + sigmaB**2) |
| 135 | }}} |