Çok önceden hazırlamış olduğum bu konuyu biraz da isteksizlik ve üşengeçlik nedeniyle sizlerle geç paylaşıyorum. Neyse geç olsun da güç olmasın diyerekten bu duyarlılık göstergelerinin ne olduğunu ve matlabda nasıl hesaplandığını sizlerle paylaşayım.
Delta : Bir opsiyonun deltası demek, onun dayanak varlığının fiyatındaki değişimin, opsiyonun fiyatını nasıl etkilediğini gösteren bir parametredir. Yani, opsiyonun fiyatı opsiyona konu olan varlığın fiyatına karşı ne kadar duyarlıdır, bunu ölçen bir göstergedir.
Gamma : Bu ise, delta nın duyarlılığını ölçer. Yani, opsiyona konu olan varlığın fiyatındaki değişmeler, opsiyonun deltasını nasıl etkiler, bunu gösteren bir parametredir.
Vega : Opsiyona konu olan varlığın fiyatındaki volatilitenin, opsiyon fiyatına olan etkisini ölçen bir parametredir.
Theta : Vadeye kalan gün sayısının, opsiyonun fiyatı üzerinde nasıl bir etkisi olduğunu gösteren bir parametredir.
Rho : Faiz oranlarının opsiyonun fiyatında ne kadar etkisi olduğunu gösteren bir parametredir. Yani faiz oranlarındaki 1 birimlik değişim, opsiyonun fiyatında ne kadarlık bir artış veya azalışa neden olduğunu gösteren bir parametredir.
Bu risk göstergelerinin nasıl hesaplandığına gelince. Bunlar, opsiyon fiyatlama modeli olan Black&Scholes formülünün türevleri alınarak hesaplanır. Bu denklemin türevlerinin nasıl hesaplandığını wikipediadan bulabilirsiniz ( http://en.wikipedia.org/wiki/Greeks_(finance) ). Opsiyon fiyatlama modeli Black&Scholes ile ilgili yazıma da buradan ulaşabilirsiniz.
Gelelim bu parametrelerin wikipedia linkinde gösterilen türev alma hesaplamaları ışığında matlab kodlarına :
Matlab kodu için tıklayınız.
prezzi=xlsread(‘unicredit.xlsx’)
k=1.9
r=0.04
t=0.25
sigma=std(prezzi)
for i=1:length(prezzi)
d1(i)=(log(prezzi(i)/k)+(r+(sigma^2)/2)*t)/(sigma*sqrt(t));
d2(i)=d1(i)-sigma*sqrt(t);
Nd1(i)=normcdf(d1(i),0,1);
dellta(i)=Nd1(i);%valore di delta
Nd2(i)=normcdf(d2(i),0,1);
C(i)=prezzi(i)*Nd1(i)-k*Nd2(i)*exp(-r*t);%prezzo dell’opzione
theta(i)=-(prezzi(i)*((exp(-0.5*d1(i)^2))/sqrt(2*pi))*sigma)/(2*sqrt(t))-r*k*Nd2(i)*exp(-r*t);%valore di theta
gamma(i)=((exp(-0.5*d1(i)^2))/sqrt(2*pi))/(prezzi(i)*sigma*sqrt(t));%valore di gamma
vega(i)=prezzi(i)*sqrt(t)*((exp(-0.5*d1(i)^2))/sqrt(2*pi));%valore di vega
rho(i)=k*t*normcdf(d2(i),0,1)*exp(-r*t);
end
plot(prezzi,dellta)
figure
plot(prezzi,theta)
figure
plot(prezzi,gamma)
figure
plot(prezzi,vega)
figure
plot(prezzi,rho)