一、期权的概念

期权是一种金融衍生工具，赋予持有人在未来某一特定日期或之前以预定价格买入或卖出标的资产的权利，但不承担必须买卖的义务。期权主要分为看涨期权（Call Option）和看跌期权（Put Option）。

这个视频讲的不错：

1.1 期权的基本要素

1. 标的资产：期权合约所涉及的资产，可以是股票、指数、商品、货币等。
2. 行权价格（Strike Price）：期权持有人可以在期权到期日或之前买入或卖出标的资产的价格。
3. 到期日（Expiration Date）：期权合约的有效期限，超过此日期期权将失效。
4. 期权费（Premium）：购买期权时支付的费用，是期权的价格。

简单来说，（看涨期权）你可以拿钱去卖方买期权（相当于一张合约、一张合同），你买要花钱，这个就是期权的期权费。现在拥有了一个期权合约，这个合约的作用就是，约定了一个价格（成为行权价格、执行价格），在合约到期（到期日）之前，你有权利以这个价格向卖家买期权合约对应的标的物（标的物可以是股票、指数、商品等等东西），当然了你也可以不买；但是卖家（卖你期权的人）必须（有义务）满足你的权利（以约定的价格即执行价格购买标的物）。

你怎么赚钱呢？

• 买期权的时候，期权合约不是约定了一个执行价格嘛，就是你可以以这个价格买标的物，但是你买了期权后过了几天，标的物的价格上涨（市场价格浮动嘛）了，那你就可以以低价买高价的东西了（卖家必须卖给你，他交了押金的）；
• 或者你不买标的物，你发现期权的价格上涨了，你就可以把期权卖给别人。
• 总之就是赚差价，但要注意，期权合约是有时间限制的，超过了约定的到期日，这个合约就没用用了。

当然了，要是标的物价格下跌了、期权价格下跌了，你就亏啦，但最多也就亏个本金嘛，就是你买期权合约花的钱；理论上赚钱是无上限的。

1.2 看涨期权和看跌期权

• 看涨期权（Call Option）：赋予持有人在未来某一日期或之前以特定价格买入标的资产的权利。买入看涨期权的投资者预期标的资产价格将上涨。
• 看跌期权（Put Option）：赋予持有人在未来某一日期或之前以特定价格卖出标的资产的权利。买入看跌期权的投资者预期标的资产价格将下跌。

这两种期权又叫认购期权和认沽期权。

前面我详细讲解的是看涨期权。

认沽期权也是类似的，标的物价格下跌，你就可以以高价（执行价格）卖给期权卖家了；或者期权价格上涨，你也可以转差价。

欧式期权，只能在到期日行使权利（以执行价格买卖标的物），美式期权则不必。但大多数人，不选择行使权利，而是直接买卖期权（平仓），这样利润更大（不展开解释了）。

期权还有时间价值，简单来说就是离到期日越远，期权的价格、标的物的价格有更多的波动可能性，潜在的价值更高。

此外，你也可以做期权卖家，去卖这2种期权。但这时，你的亏损理论上无限大。

1.3 期权的操作策略

1. 买入看涨期权（Long Call）：预期标的资产价格会上涨，投资者买入看涨期权。
2. 买入看跌期权（Long Put）：预期标的资产价格会下跌，投资者买入看跌期权。
3. 卖出看涨期权（Short Call）：预期标的资产价格不会上涨，投资者卖出看涨期权。
4. 卖出看跌期权（Short Put）：预期标的资产价格不会下跌，投资者卖出看跌期权。

1.4 期权的定价模型

期权定价的主要模型有：

1. Black-Scholes模型：这是最早也是最著名的期权定价模型，基于无套利理论和对数正态分布假设。
2. 二叉树模型（Binomial Model）：使用树状结构模拟标的资产价格变化，适用于美式期权。
3. 蒙特卡罗模拟（Monte Carlo Simulation）：通过大量随机模拟路径计算期权价格，适用于复杂期权。

1.5 期权的风险和收益

• 风险：买入期权的最大损失是期权费，而卖出期权的潜在损失是无限的。
• 收益：买入看涨期权的收益是标的资产价格减去行权价格和期权费；买入看跌期权的收益是行权价格减去标的资产价格和期权费。

二、期权定价模型和隐含波动率计算

常见的期权定价模型有 Black-Scholes 模型，其公式如下：

对于欧式看涨期权（Call Option）的定价公式是：

$$C=S_0\cdot N(d_1)-K\cdot e^{-rT}\cdot N(d_2)$$

对于欧式看跌期权（Put Option）的定价公式是：

$$P=K\cdot e^{-rT}\cdot N(-d_2)-S_0\cdot N(-d_1)$$

其中，$d_1$ 和 $d_2$ 的计算公式为：

$$d_1=\frac{\ln \left(\frac{S_0}{K}\right)+\left(r+\frac{{\sigma }^2}{2}\right)T}{\sigma \sqrt{T}}$$

$$d_2=d_1-\sigma \sqrt{T}$$

符号解释：

• $C$: 看涨期权价格
• $P$: 看跌期权价格
• $S_0$: 标的资产当前价格
• $K$: 行权价格
• $r$: 无风险利率
• $\sigma$: 标的资产的波动率（隐含波动率）
• $T$: 到期时间
• $N(\cdot )$: 标准正态分布函数的累积分布函数

这些公式是基于以下假设的：

1. 标的资产价格遵循几何布朗运动。
2. 无风险利率和波动率在期权有效期内保持不变。
3. 不考虑红利支付。
4. 市场是无摩擦的（即没有交易成本、税收等）。

---

根据前面的公式就可以求解隐含波动率，常用的方法有牛顿法、二分法。

三、代码示例

以纳斯达克100指数期权为例，简写NDX。

假设今天是2024-6-15，标的物价格为：19659.8，无风险利率为5.53%，期权价格为$\frac{Ask+Bid}{2}$

到期日为2024-6-21的期权数据：（Skrike：执行价格，Ticker：合约代码，Bid：买方愿出最高价，Ask：买方可就是最低价，Last：最近成交价，Volm：交易量）

Calls						Puts
Strike	Ticker	Bid	Ask	Last	Volm	Strike	Ticker	Bid	Ask	Last	Volm
19530	NDX 6/21/24 C19530	214.1000061	229.8000031	165.3999939	0	19530	NDX 6/21/24 P19530	63	65.40000153	67.30000305	2
19540	NDX 6/21/24 C19540	206.6000061	222	160.6499939	0	19540	NDX 6/21/24 P19540	59.79999924	68.19999695	70.09999847	3
19550	NDX 6/21/24 C19550	199.3000031	214.3999939	174.4799957	29	19550	NDX 6/21/24 P19550	62.5	70.90000153	72.69999695	2
19560	NDX 6/21/24 C19560	192.1999969	207.6000061	140	2	19560	NDX 6/21/24 P19560	65.69999695	73.80000305	75.90000153	2
19570	NDX 6/21/24 C19570	184.8999939	200.8999939	132.3500061	0	19570	NDX 6/21/24 P19570	68.59999847	76.80000305	79.19999695	8
19575	NDX 6/21/24 C19575	181.5	197.5	164	25	19575	NDX 6/21/24 P19575	69.69999695	78.40000153	80.30000305	4
19580	NDX 6/21/24 C19580	178.1999969	194.3000031	159.6199951	6	19580	NDX 6/21/24 P19580	71.30000305	80	86.84999847	10
19590	NDX 6/21/24 C19590	171.1999969	187.5	140.1999969	1	19590	NDX 6/21/24 P19590	74.80000305	83.40000153	85.40000153	6
19600	NDX 6/21/24 C19600	164.6000061	180.5	158.7200012	18	19600	NDX 6/21/24 P19600	77.80000305	86.59999847	95	54
19610	NDX 6/21/24 C19610	158	173.8999939	160	6	19610	NDX 6/21/24 P19610	81.30000305	90.19999695	87.62999725	20
19620	NDX 6/21/24 C19620	152	161.5	125.5	41	19620	NDX 6/21/24 P19620	85	93.80000305	103.1100006	43
19625	NDX 6/21/24 C19625	148.8999939	158.3000031	148	18	19625	NDX 6/21/24 P19625	86.90000153	95.80000305	105.1500015	22
19630	NDX 6/21/24 C19630	145.6999969	155.1999969	84.19999695	0	19630	NDX 6/21/24 P19630	88.69999695	97.59999847	0	0
19640	NDX 6/21/24 C19640	139.6000061	149	135.2700043	16	19640	NDX 6/21/24 P19640	92.59999847	104.5	106.5	10
19650	NDX 6/21/24 C19650	133.5	142.8999939	130.8999939	21	19650	NDX 6/21/24 P19650	96.5	105.8000031	109	12
19660	NDX 6/21/24 C19660	128	137.1000061	126.5500031	5	19660	NDX 6/21/24 P19660	100.6999969	109.8000031	130.5	1
19670	NDX 6/21/24 C19670	121.9000015	131.1999969	111.8799973	18	19670	NDX 6/21/24 P19670	104.9000015	114.0999985	145.3999939	0
19675	NDX 6/21/24 C19675	119.5999985	128.3000031	109.1999969	6	19675	NDX 6/21/24 P19675	107.1999969	116.1999969	157.3999939	0
19680	NDX 6/21/24 C19680	116.0999985	125.5	90.41999817	0	19680	NDX 6/21/24 P19680	109.6999969	118.5999985	0	0
19690	NDX 6/21/24 C19690	110.6999969	120	81.30000305	0	19690	NDX 6/21/24 P19690	114.4000015	123.0999985	0	0
19700	NDX 6/21/24 C19700	105.5	114.5999985	95.47000122	10	19700	NDX 6/21/24 P19700	119.0999985	127.9000015	160.5500031	6
19710	NDX 6/21/24 C19710	100.3000031	109.4000015	91.80000305	4	19710	NDX 6/21/24 P19710	123.9000015	132.8999939	0	0
19720	NDX 6/21/24 C19720	96	104.4000015	59.31999969	0	19720	NDX 6/21/24 P19720	128.8999939	137.8999939	161.3999939	3
19725	NDX 6/21/24 C19725	93.19999695	101.9000015	84.05000305	1	19725	NDX 6/21/24 P19725	131.3999939	140.3000031	0	0
19730	NDX 6/21/24 C19730	90.69999695	99.5	80.09999847	6	19730	NDX 6/21/24 P19730	134.1000061	143.1999969	0	0
19740	NDX 6/21/24 C19740	86.30000305	94.69999695	0	0	19740	NDX 6/21/24 P19740	139.1999969	148.5	0	0
19750	NDX 6/21/24 C19750	81.30000305	90.09999847	80.80000305	12	19750	NDX 6/21/24 P19750	144.6999969	154.1000061	198.3000031	1
19760	NDX 6/21/24 C19760	76.90000153	85.69999695	64.5	1	19760	NDX 6/21/24 P19760	150.1000061	165.8000031	0	0
19770	NDX 6/21/24 C19770	73.19999695	81.30000305	63.40000153	4	19770	NDX 6/21/24 P19770	155.8000031	171.3999939	0	0
19775	NDX 6/21/24 C19775	71.09999847	79.30000305	68.40000153	1	19775	NDX 6/21/24 P19775	158.6999969	174.1999969	0	0
19780	NDX 6/21/24 C19780	69	77.19999695	60	1	19780	NDX 6/21/24 P19780	161.6000061	177	0	0

---

import pandas as pd
import numpy as np
from scipy.stats import norm
from scipy.optimize import newton
from datetime import datetime
import matplotlib.pyplot as plt


class OptionAnalysis:
    def __init__(self, file_path):
        self.file_path = file_path
        self.S = 19659.80
        self.r = 0.0553
        self.current_date = datetime(2024, 6, 15)
        self.data = pd.read_excel(file_path, sheet_name=0, header=[0, 1])

    def bs_price(self, S, K, T, r, sigma, option_type='call'):
        d1 = (np.log(S / K) + (r + 0.5 * sigma ** 2) * T) / (sigma * np.sqrt(T))
        d2 = d1 - sigma * np.sqrt(T)
        if option_type == 'call':
            return S * norm.cdf(d1) - K * np.exp(-r * T) * norm.cdf(d2)
        return K * np.exp(-r * T) * norm.cdf(-d2) - S * norm.cdf(-d1)

    def implied_volatility(self, option_price, S, K, T, r, option_type='call'):
        def difference(sigma):
            return self.bs_price(S, K, T, r, sigma, option_type) - option_price
        return newton(difference, 0.2, tol=1e-5)

    def process_data(self):
        strike_calls = self.data[('Calls', 'Strike')].dropna().values
        bid_calls = self.data[('Calls', 'Bid')].dropna().values
        ask_calls = self.data[('Calls', 'Ask')].dropna().values
        price_calls = (bid_calls + ask_calls) / 2

        strike_puts = self.data[('Puts', 'Strike')].dropna().values
        bid_puts = self.data[('Puts', 'Bid')].dropna().values
        ask_puts = self.data[('Puts', 'Ask')].dropna().values
        price_puts = (bid_puts + ask_puts) / 2

        expiry_date = datetime.strptime(self.data.index.name, '%m-%d-%y')
        T = (expiry_date - self.current_date).days / 365

        implied_vols_calls = [self.implied_volatility(price, self.S, K, T, self.r, 'call') for price, K in zip(price_calls, strike_calls)]
        implied_vols_puts = [self.implied_volatility(price, self.S, K, T, self.r, 'put') for price, K in zip(price_puts, strike_puts)]

        self.data[('Calls', 'Implied Volatility')] = pd.Series(implied_vols_calls, index=self.data.index[:len(implied_vols_calls)])
        self.data[('Puts', 'Implied Volatility')] = pd.Series(implied_vols_puts, index=self.data.index[:len(implied_vols_puts)])

        return strike_calls, implied_vols_calls, strike_puts, implied_vols_puts

    def plot_data(self):
        strike_calls, implied_vols_calls, strike_puts, implied_vols_puts = self.process_data()

        plt.figure(figsize=(14, 6))

        # Calls
        plt.subplot(1, 2, 1)
        plt.plot(strike_calls, implied_vols_calls, 'go-', label='Call Implied Volatility')
        plt.xlabel('Strike Price')
        plt.ylabel('Implied Volatility')
        plt.title('Call Implied Volatility vs Strike Prices')
        plt.grid(True)
        plt.legend()

        # Puts
        plt.subplot(1, 2, 2)
        plt.plot(strike_puts, implied_vols_puts, 'mo-', label='Put Implied Volatility')
        plt.xlabel('Strike Price')
        plt.ylabel('Implied Volatility')
        plt.title('Put Implied Volatility vs Strike Prices')
        plt.grid(True)
        plt.legend()

        plt.suptitle(f'NDX: Options Implied Volatility', color='blue')
        plt.tight_layout()
        plt.savefig('ndx_img/Options_volatility.png')


if __name__ == '__main__':
    file_path = 'NDX.xlsx'
    analysis = OptionAnalysis(file_path)
    analysis.plot_data()