Aquí tiene un listado de un programa que simula el comportamiento estadístico de un cultivo acuícola en el cual hay variaciones en supervivencia, ocurrencia de enfermedades y en la talla de los organismos

 

 

 

variable	Distribución	Parámetros
Numero de enfermedades	Poisson	Número promedio de ocurrencias en un ciclo
Numero de sobrevivientes	Binomial	Número de ensayos (organismos), Probabilidad de éxito (porcentaje que sobreviven)
Peso de un organismo	Normal	Peso promedio Desviación estándar

 

El programa calcula el peso del producto y lo multiplica por el precio por kg en cada ciclo.

 

 

 

%programa que simula sistema acuicola%

pv=80;             %precio por kilo%

pe=2.0 ;           %peso promedio%

de=0.2 ;           %desviacion en peso%

psu=0.8;            %tasa de superviviencia%

n=1000;             %numero de semillas%

pp=0.2;             %ocurrencia promedio de enfermedad por ciclo%

pse=0.2;             %probabilidad supervivencia en enfermedad%

nc= 10;                %numero de ciclos%

ganancia=1:nc;         %valores iniciales de ganancia%

ns=1:nc;

 

for i=1:nc

    enf=poissrnd(0.1);

    if (enf <1)

    ns(i)=binornd(n,psu) ;       %numero de sobrevivientes sin enfermedad en ciclo i%

else

    ns(i)=binornd(n,pse);       %numero de sobrevivientes con enfermedad%

end

    pesos=normrnd(pe,de,ns(i),1);  %matriz de pesos%

    speso=sum(pesos);           %suma de pesos%

    ganancia(i)=speso*pv;

end

subplot(1,2,1);

plot(ganancia)

title('ganancia')

subplot(1,2,2);

plot(ns)

title('sobrevivientes')