Procesos estocásticos

Notas procesos estocásticos
Ejemplo seleccionar la muestra Agosto 12 k15 http://www.mathworks.com/help/stats/examples/selecting-a-sample-size.html

Ejemplo Agosto 12 k15 Fitting Custom Univariate Distributions  http://www.mathworks.com/help/stats/examples/fitting-custom-univariate-distributions.html

Ejemplo Agosto 12 k15 Simulating Dependent Random Variables Using Copulas http://www.mathworks.com/help/stats/examples/simulating-dependent-random-variables-using-copulas.html

mACHINE LEARNIG CANAL
MatLab en línea
SCIlab ayuda
Equivalentes MATLAB-SCILAB
Generar numeros aleatorios
Distribuciones discretas    http://www.mathworks.com/help/stats/discrete-distributions.html
Cumultive distribution funtions    http://www.mathworks.com/help/stats/cdf.html
Julio 20 k15 Monte Carlo Simulation    http://www.mathworks.com/discovery/monte-carlo-simulation.html
Simulación Montecarlo VIDEOS http://www.mathworks.com/videos/search.html?q=%22monte+carlo+simulation%22
Explicación Metodo de Montecarlo vía cadenas de Markov  https://www.youtube.com/watch?v=LI38DSFBcKk


clear
tini=clock;
p= normcdf (4.5, 4, 0.64)
clc
p
clc
1-p
A= [0 1 0 0 0 0 0 0 0; 0 0 1 0 0 0 0 0 0; 0 0 0 1 1 0 0 0 0; 1 0 0 0 0 0 0 0 0; 0 0 0 0 0 1 1 0 0; 0 0 0 0 1 0 0 0 0; 0 0 0 0 0 0 0 1 0; 0 0 0 0 0 0 0 0 1; 0 0 0 0 0 0 1 0 0]
b = 0.3
[m, n] = size(A)
D = logical(inv(eye(m) - b*A))
C = D & D'

A=rand(6), B=inv(A), B*A

n=100; A=rand(n); B=rand(n); C=zeros(n);

 C=B*A; tend=clock; mflops=(2*n^3)/1000000;

mflops/etime(tend,tini)
x=-4:.01:4; y=sin(x); plot(x,y), grid, title('Funtion sin(x)')
whos
date
clock
pol=[1 -4 5 -2]
roots(pol)

# junio 21 2015
# ver
# Difiero en algunos resultados que el video presenta pero dan las mismos lamda 0
A= [.6 .5 .3; .2 .3 .2; .2 .2 .5]
[V,E] = eig(A) #Calcula los autovectores o vectores propios

ahora aplicando el metodo QR Shifted
b= A(3,3)
# eye (3,3 es la idéntica de 3x3
A= A-b*eye(3,3)
[Q,R]=qr(A)
A= R*Q
A=A+b*eye(3,3)

b= A(3,3)
# eye (3,3 es la idéntica de 3x3
A= A-b*eye(3,3)
[Q,R]=qr(A)
A= R*Q
A=A+b*eye(3,3)

b= A(3,3)
# eye (3,3 es la idéntica de 3x3
A= A-b*eye(3,3)
[Q,R]=qr(A)
A= R*Q
A=A+b*eye(3,3)

b= A(3,3)
# eye (3,3 es la idéntica de 3x3
A= A-b*eye(3,3)
[Q,R]=qr(A)
A= R*Q
A=A+b*eye(3,3)