问题描述
我目前正在尝试为我的固定收益证券课程模拟一个随机过程(具体来说,Vasicek 模型)。我需要对相同的参数多次运行模拟,以显示利率可以采取的不同路径。我目前使用的代码是:
import numpy as np
import matplotlib.pyplot as plt
def vasicek(r0,K,theta,sigma,T,N,seed=777):
np.random.seed(seed)
dt = T/float(N)
rates = [r0]
for i in range(N):
dr = K*(theta-rates[-1])*dt + sigma*np.random.normal()
rates.append(rates[-1] + dr)
return range(N+1),rates
我的问题是该函数在绘制范围和比率时返回相同的路径。每次运行此函数时, np.random.normal() 不应该返回不同的值吗?如何更改此代码以使其执行此操作?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)