close
Warning:
Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.
- Timestamp:
-
Feb 2, 2011, 9:39:24 AM (14 years ago)
- Author:
-
Víctor de Buen Remiro
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v27
|
v28
|
|
74 | 74 | es que el método directamente fracasa y no es capaz de simular el siguiente punto. |
75 | 75 | |
76 | | == Diseños de generadores de candidatos factibles == |
77 | | |
78 | 76 | Se hace necesario por lo tanto disponer de generadores de candidatos que por construcción estén siempre |
79 | 77 | dentro del politopo, es decir generadores de candidatos factibles. Tal generador no podrá ser simétrico |
… |
… |
|
81 | 79 | eficiente. |
82 | 80 | |
83 | | === Definiciones === |
| 81 | == Definiciones previas == |
84 | 82 | Dado un punto [[LatexEquation(x)]] estrictamente interior al politopo |
85 | 83 | |
… |
… |
|
99 | 97 | será responsabilidad del analista no introducir demasiadas. |
100 | 98 | |
101 | | === Búsqueda de un punto inicial estrictamente interior === |
| 99 | == Búsqueda de un punto inicial estrictamente interior == |
102 | 100 | |
103 | 101 | Para encontrar un punto estrictamente interior hay que encontrar un punto que cumpla las restricciones |
… |
… |
|
134 | 132 | de parada razonables como pudiera ser un tiempo máximo y un valor mínimo para la log-densidad. |
135 | 133 | |
136 | | === Paseo aleatorio hiperesférico === |
| 134 | == Paseo aleatorio hiperesférico == |
137 | 135 | Se propone en primer lugar utilizar un generador con distribución uniforme en una hiperesfera centrada |
138 | 136 | en el último punto generado y que esté incluida estrictamente en el politopo lo cual será cierto si el |
… |
… |
|
157 | 155 | si se acerca a un vértice demasiado agudo. |
158 | 156 | |
159 | | ==== Función de densidad ==== |
| 157 | === Función de densidad === |
160 | 158 | La densidad del generador será porporcional al inverso del volumen de la hiperesfera y su logaritmo será, salvo una constante |
161 | 159 | |
… |
… |
|
170 | 168 | la cual habrá que calcular dos veces por cada candidato o precandidato generado. |
171 | 169 | |
172 | | ==== Función generatriz ==== |
| 170 | === Función generatriz === |
173 | 171 | |
174 | 172 | Para generar un candidato [[LatexEquation(y)]] a partir del actual [[LatexEquation(x)]] con esta |
… |
… |
|
204 | 202 | reservas enunciadas. |
205 | 203 | |
206 | | === Paseo aleatorio radial asimétrico === |
| 204 | == Paseo aleatorio radial asimétrico == |
207 | 205 | A continuación de define otra forma de generador factible similar a la anterior pero que podría |
208 | 206 | adaptarse mejor cuando el punto se acerca demasiado a la frontera, e incluso es válida si está en |
… |
… |
|
222 | 220 | [[Image(source:/tolp/OfficialTolArchiveNetwork/BysSampler/doc/image/RandWalk.InPolytope.chart.726488582_c.png)]] |
223 | 221 | |
224 | | ==== Función generatriz ==== |
| 222 | === Función generatriz === |
225 | 223 | Los pasos para muestrear un candidato serían los siguientes |
226 | 224 | 1. Primero se simulará un vector multinormal estandarizado [[BR]] [[BR]] |
… |
… |
|
252 | 250 | |
253 | 251 | |
254 | | ==== Función de densidad ==== |
| 252 | === Función de densidad === |
255 | 253 | |
256 | 254 | La función de densidad en este caso es simplemente proporcional al inverso de la |
… |
… |
|
271 | 269 | |
272 | 270 | |
273 | | === Paseo aleatorio mixto === |
| 271 | == Paseo aleatorio mixto == |
274 | 272 | El paseo aleatorio hiperesférico tiene la ventaja de ser muy rápido pues sólo se necesita la distancia |
275 | 273 | a la frontera, pero tiene el inconveniente que no es capaz de acercarse demasiado a la frontera y menos |