Matlab - Grafica Estimada

   
Vista:

Grafica Estimada

Publicado por Luis (9 intervenciones) el 10/09/2017 23:00:37
Buenas tardes,

Me encuentro graficando una funcion beta generando numeros aleatorios por acep-rechzo y tengo la grafica teorica de dicha funcion

Función por aceptación rechazo

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
function [b,tp] = beta34v2(n)
% Escribe un codigo que genere numeros aleatorios con distribucion 
% B(2, 2) ~ (f(x) = 6x(1 - x), 0 >= x >= 1) mediante el metodo de 
% aceptacion-rechazo simple
% Input:
% n :   número de muestras
% Output:
% b :   muestras aleatorias generadas y aceptadas
% tp :  total de puntos generados, aceptados y rechazados 
a = 0.5;
b = []; tp = [];
for i = 1:n,                            % Hasta que tenga todas las muestras
    u1 = rand;                          % Generar u1
    u2 = a * rand;                      % Generar u2
    while (u2 > 6 * u1 * (1-u1) )  % Mientras (U * a * g(x))> f(x) 
        % Bloque para generar x
        tp = [tp ; u1];                 % Recoger todas las muestras
        u1 = rand;                      % Volver a generar u1
        u2 = a * rand;                  % Volver a generar u2
    end
    x = u1;                             % x generado
    b = [b ; x];                        % Recoge las muestras aceptadas
    tp = [tp ; u1];                     % Recoge las muestras aceptadas
end


Grafica teorica

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear all
figure(gcf)
x = 0:0.01:1;
syms X
a=2;b=2;
y = zeros(size(x));
 
N=int(X^(a+b)*exp(-X), 0, inf);
D1=int(X^(a)*exp(-X), 0, inf);
D2=int(X^(b)*exp(-X), 0, inf);
for i = 1:size(x,2)
% Aplicacion de la Funcion de Densidad de Probabilidad
 
y(i)=N/(D1*D2)*(x(i)^(a-1)*(1-x(i))^(b-1));
end
 plot(x,y)
 set(gca,'xtick',0:0.1:1),
 set(gca,'ytick',0:0.2:2.6),
title('Densidad de Probabilidad Teorica')

Desearia de ser posible que me ayuden a realizar la grafica estimada en el mismo plot de la teorica para ver la diferencia.

Agradezco de antemano cualquier ayuda

varianza
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Revisar política de publicidad