在能夠計算
選擇權評價-BS MODEL (with Excel) 及
隱含波動率 (with EXCEL VBA、C#) 之後，再來剩GREEKS了，要
計算四個影響價格的變數是Delta、Gamma、 Vega、Theta，完成這部份就可以得出目前選擇權部位組合的各項數值。關於這些希臘字母的定義和作用，都需要深入瞭解後才能對選擇權策略有所掌握，進而設計出選擇權策略。
此篇要紀錄的僅是關於系統開發，雖說大部份的交易軟體都會提供GREEKS資料，但總要是自己寫的才能確定正確性，因為券商軟體算的不一定對。當然也有許多人不在意些微差距，或根本沒在觀察GREEKS，但那樣操作選擇權是無法進步的，再說若要自己建立策略模組，這些也只是第一步而已。
大部份設計選擇權監視程式的軟體是EXCEL，因為最簡便，直接利用公式和一些函數設計就可以完成，若要即時串接部位的話，大部份是寫VBA程式和券商API串接，而行情的部份可以用DDE或API來接收。
但是EXCEL會有效能問題，如果接受行情多，再加上運算邏輯複雜，EXCEL很可能會有延遲的狀況，因此比較好的方式還是自己寫程式，可以
先由券商API所提供的範例程式來測試，可能是VB、C#、C++或Delphi等等，選一個好寫的來測試收行情、回報、部位。此篇L要紀錄的是以C#來實作選擇權GREEKS。
整段範例程式碼貼在最後面，需要先說明一下前因後果，先有一個DataTable是用來紀錄未平倉的選擇權部位，這個副程式所要處理的是計算目前所有未平倉選擇權部位的GREEKS並寫進DataTable，這個DataTable再與DataGridView串接，就可以顯示在FORM上
這個DataTable，取名dtOI，至少有以下的欄位:
dtOI.Columns.Add(new DataColumn("SettleM", typeof(string))); //月份
dtOI.Columns.Add(new DataColumn("Strike", typeof(double))); //履約價
dtOI.Columns.Add(new DataColumn("CP", typeof(string))); //Call or Put
dtOI.Columns.Add(new DataColumn("BS", typeof(string))); //Buy or Short
dtOI.Columns.Add(new DataColumn("Qty", typeof(string))); //口數
dtOI.Columns.Add(new DataColumn("Match", typeof(double))); //成交價
dtOI.Columns.Add(new DataColumn("Market", typeof(double))); //市價
dtOI.Columns.Add(new DataColumn("VOL(%)", typeof(double))); //隱含波動率
dtOI.Columns.Add(new DataColumn("DELTA", typeof(double))); //DELTA
dtOI.Columns.Add(new DataColumn("GAMMA", typeof(double)));//GAMMA
dtOI.Columns.Add(new DataColumn("THETA", typeof(double))); //THETA
dtOI.Columns.Add(new DataColumn("VEGA", typeof(double))); //VEGA
在執行副程式之前，除了最後五欄之外都是準備好的，然後有一些變數要準備
AssetPrice 期貨價格
Strike 履約價
InterestRate 利率
Expiry 到期因子
Target 選擇權的市價
還有之前紀錄隱含波動率的類別，一開始要NEW一下。
Option op = new Option();
然後就開始了，看片段來說明一下
以下PUT的部份就依此類推。程式碼如後所貼，有一些使用到的參數及關聯FORM的變數沒有仔細說明，不過整個架構已經是可以瞭解的了。以前 L 設計時找不到完整的參考資料，一步步的自己設計，現在
若有同樣需求的朋友看到這邊可以自由取用，歡迎參考及指正。
private void Greeks()
{
int i = 0;
AssetPrice = double.Parse(label16.Text);
while (i < dtOI.Rows.Count)
{
Strike = double.Parse(dtOI.Rows[i]["Strike"].ToString());
Target = double.Parse(dtOI.Rows[i]["Market"].ToString());
Qty = double.Parse(dtOI.Rows[i]["Qty"].ToString());
//期貨
if (dtOI.Rows[i]["CP"].ToString() == "")
{
if (dtOI.Rows[i]["BS"].ToString().Trim() == "B")
{
dtOI.Rows[i].SetField("DELTA", 200 * Qty );
}
if (dtOI.Rows[i]["BS"].ToString().Trim() == "S")
{
dtOI.Rows[i].SetField("DELTA", -200 * Qty );
}
}
//CALL civ(double AssetPrice, double Strike, double InterestRate, double Expiry, double Target)
if ((dtOI.Rows[i]["CP"].ToString().Trim()) == "C")
{
if (dtOI.Rows[i]["PId"].ToString().Trim() == "TXO")
{
if (dtOI.Rows[i]["SettleM"].ToString().Trim() == label10.Text.Substring(0, 6)) //近月
{
AssetPrice = double.Parse(label16.Text);
Expiry = Expirynear;
}
if (dtOI.Rows[i]["SettleM"].ToString().Trim() == label11.Text.Substring(0, 6)) //遠月
{
if (label22.Text != "")
{
AssetPrice = double.Parse(label22.Text);
}
else { AssetPrice = double.Parse(label16.Text); }
Expiry = Expiryfar;
}
}
else //週選
{
if (wnFid != "MX")
{
if (dtOI.Rows[i]["PId"].ToString().Trim().Substring(1, 2) == wnFid.Substring(1, 2)) //近週
{
AssetPrice = double.Parse(label23.Text);
Expiry = Expirywn;
}
}
if (wfFid != "MX")
{
if (dtOI.Rows[i]["PId"].ToString().Trim().Substring(1, 2) == wfFid.Substring(1, 2)) //遠週
{
AssetPrice = double.Parse(label22.Text);
Expiry = Expirywf;
}
}
}
civ = op.civ(AssetPrice, Strike, InterestRate, Expiry, Target);
d1 = op.d1(AssetPrice, Strike, InterestRate, Expiry, civ);
nd1 = Math.Exp(-0.5 * d1 * d1) / Math.Sqrt(2 * Math.PI);
nd2 = op.NormsDist(d1 - civ * Math.Sqrt(Expiry));
//DELTA NormsDist(double x) d1(double AssetPrice, double Strike, double InterestRate, double Expiry, double Volatility)
DELTA = op.NormsDist(d1);
GAMMA = nd1 / (AssetPrice * civ * Math.Sqrt(Expiry));
THETA = (-(AssetPrice * nd1 * civ / (2 * Math.Sqrt(Expiry))) - InterestRate * Strike * Math.Exp(-InterestRate * Expiry) * nd2) / 365;
VEGA = (AssetPrice * (Math.Sqrt(Expiry)) * nd1) / 100;
dtOI.Rows[i].SetField("VOL(%)", Math.Round(civ * 100, 4));
if (dtOI.Rows[i]["BS"].ToString().Trim() == "B")
{
dtOI.Rows[i].SetField("DELTA", Math.Round(DELTA * 50 * Qty , 0));
dtOI.Rows[i].SetField("GAMMA", Math.Round(GAMMA * 50 * Qty * AssetPrice / 100, 0));
dtOI.Rows[i].SetField("THETA", Math.Round(THETA * 50 * Qty, 0));
dtOI.Rows[i].SetField("VEGA", Math.Round(VEGA * 50 * Qty, 0));
}
else if (dtOI.Rows[i]["BS"].ToString().Trim() == "S")
{
dtOI.Rows[i].SetField("DELTA", Math.Round(-DELTA * 50 * Qty , 0));
dtOI.Rows[i].SetField("GAMMA", Math.Round(-GAMMA * 50 * Qty * AssetPrice / 100, 0));
dtOI.Rows[i].SetField("THETA", Math.Round(-THETA * 50 * Qty, 0));
dtOI.Rows[i].SetField("VEGA", Math.Round(-VEGA * 50 * Qty, 0));
}
}
//PUT piv(double AssetPrice, double Strike, double InterestRate, double Expiry, double Target)
if ((dtOI.Rows[i]["CP"].ToString().Trim()) == "P")
{
if (dtOI.Rows[i]["PId"].ToString().Trim() == "TXO")
{
if (dtOI.Rows[i]["SettleM"].ToString().Trim() == label10.Text.Substring(0, 6)) //近月
{
AssetPrice = double.Parse(label16.Text);
Expiry = Expirynear;
}
if (dtOI.Rows[i]["SettleM"].ToString().Trim() == label11.Text.Substring(0, 6)) //遠月
{
if (label22.Text != "")
{
AssetPrice = double.Parse(label22.Text);
}
else { AssetPrice = double.Parse(label16.Text); }
Expiry = Expiryfar;
}
}
else //週選
{
if (wnFid != "MX")
{
if (dtOI.Rows[i]["PId"].ToString().Trim().Substring(1, 2) == wnFid.Substring(1, 2)) //近週
{
AssetPrice = double.Parse(label23.Text);
Expiry = Expirywn;
}
}
if (wfFid != "MX")
{
if (dtOI.Rows[i]["PId"].ToString().Trim().Substring(1, 2) == wfFid.Substring(1, 2)) //遠週
{
AssetPrice = double.Parse(label22.Text);
Expiry = Expirywf;
}
}
}
piv = op.piv(AssetPrice, Strike, InterestRate, Expiry, Target);
d1 = op.d1(AssetPrice, Strike, InterestRate, Expiry, piv);
nd1 = Math.Exp(-0.5 * d1 * d1) / Math.Sqrt(2 * Math.PI);
nd2 = op.NormsDist(d1 - piv * Math.Sqrt(Expiry));
//DELTA NormsDist(double x) d1(double AssetPrice, double Strike, double InterestRate, double Expiry, double Volatility)
DELTA = op.NormsDist(d1) - 1;
GAMMA = nd1 / (AssetPrice * piv * Math.Sqrt(Expiry));
THETA = (-(AssetPrice * nd1 * piv / (2 * Math.Sqrt(Expiry))) + InterestRate * Strike * Math.Exp(-InterestRate * Expiry) * nd2) / 365;
VEGA = (AssetPrice * (Math.Sqrt(Expiry)) * nd1) / 100;
dtOI.Rows[i].SetField("VOL(%)", Math.Round(piv * 100, 4));
if (dtOI.Rows[i]["BS"].ToString().Trim() == "B")
{
dtOI.Rows[i].SetField("DELTA", Math.Round(DELTA * 50 * Qty , 0));
dtOI.Rows[i].SetField("GAMMA", Math.Round(GAMMA * 50 * Qty * AssetPrice / 100, 0));
dtOI.Rows[i].SetField("THETA", Math.Round(THETA * 50 * Qty, 0));
dtOI.Rows[i].SetField("VEGA", Math.Round(VEGA * 50 * Qty, 0));
}
else if (dtOI.Rows[i]["BS"].ToString().Trim() == "S")
{
dtOI.Rows[i].SetField("DELTA", Math.Round(-DELTA * 50 * Qty , 0));
dtOI.Rows[i].SetField("GAMMA", Math.Round(-GAMMA * 50 * Qty * AssetPrice / 100, 0));
dtOI.Rows[i].SetField("THETA", Math.Round(-THETA * 50 * Qty, 0));
dtOI.Rows[i].SetField("VEGA", Math.Round(-VEGA * 50 * Qty, 0));
}
}
//dtOI.Rows[i].SetField("VOL(%)", i);
i = i + 1;
}
}
