算法：CIR模型的参数估计量

Step 1：取初始参数初始值K = 0.21459，b = 0.08579，θ = 0.07830，n = 100，X₀ = 0.1，T = 0.5 (代表半年)，其n为[0, T]的划分份数， $dt=\frac{1}{n}$ ；

Step 2：将方程(1)式离散化为：

$X_t=X_{t-1}+0.21459\left(0.08579-X_{t-1}\right)dt+0.0783X_{t-1}\epsilon \sqrt{dt}$

生成X_t的100条样本轨道；其中ε~N(0, 1)；

Step 3：取出t₁、t₂时刻与轨道随机相交的一个点以及t₂、t₃时刻与轨道相交的一组数据点；

Step 4：将离散化数据代入 $\widehat{K},\widehat{b},\widehat{\theta }$ 估计量的表达式，得到参数的估计值。