用Matlab生成Ornstein-Uhlenbeck过程

先给出Ornstein-Uhlenbeck过程（简称OU过程）的SDE形式：

 是均值，代表均值回归速率，代表波动率（方差），为什么OU过程描述的是均值回归过程的直观理解可参见https://blog.csdn.net/qq_33254870/article/details/105137275

Matlab代码实现：

首先定义模拟生成函数（布朗运动增量用正态变量替代）

function [eta_t] = OU(mu, sigma, theta, dt, init_value)
eta_t = zeros(length(mu),1/dt);
for i=0:(1/dt)
    if i