Comment implémenter cette fonction dans MATLAB?

Question

J'utilise un livre avec une fonction que j'aimerais utiliser. Cependant, je ne pense pas que je reçois les bonnes valeurs de ma fonction.

Voici les instructions du livre:

est la fonction ici comme je l'ai créé dans Matlab:

function [ shortProbability ] = pShort(zkt, zktStar, short) 

if zkt > zktStar 
    shortProbability = 0; 
else 
    normalizer = 1/(1-exp(-short*zktStar)); 
    shortProbability = normalizer * (short * exp(-short*zkt)); 
end 

end 

Les valeurs que je suis BRANCHER sont:

zkt = 0:1:100 
zktStar = 50; 
short = 0.01; 

Cependant, mon graphique ne se comporte pas comme le un que je suis censé finir avec, ce qui est la suivante:

Je reçois ce à partir du graphique, ce qui semble correct, mais je ne pense pas qu'il est normalisé correctement:

Quelqu'un peut-il m'aider à corriger cette fonction?

Answer 1

Ceci est le exponential distribution. Vous pouvez utiliser EXPPDF et EXPCDF de la Boîte à outils statistiques:

normalizer = 1 ./ (expcdf(zktStar,1/short) - expcdf(0,1/short)); 
shortProbability = exppdf(zkt, 1/short) * normalizer; 

et il devrait être équivalent à ce que vous avez ...

---

C'est ce que je reçois quand je le combine avec le graphique de votre previous question:

Pour confirmer, nous calculons la surface sous chacun des courbes (assez proche de 1):

>> trapz(zkt,hitProbabilty) 
ans = 
    1 

>> trapz(zkt,shortProbability) 
ans = 
     1.0077