使用 scipy.minimze 最大化夏普比率

Question

我正在尝试使用 scipy.minimize

最大化夏普比率

$1$

我这样做是为了找到 CAPM 的证券市场线

$2$

所以我有一个等式：

$3$

$4$

可选（如果不允许空头头寸）：

$5$

所以我正在尝试解决这个问题：

def target_func(x, cov_matix, mean_vector, virtual_mean):
    f = float(-(x.dot(mean_vector) - virtual_mean) / np.sqrt(x.dot(cov_matix).dot(x.T)))
    return f


def optimal_portfolio_with_virtual_mean(profits, virtual_mean, allow_short=False):
    x = np.zeros(len(profits))
    mean_vector = np.mean(profits, axis=1)
    cov_matrix = np.cov(profits)
    cons = ({'type': 'eq',
             'fun': lambda x: np.sum(x) - 1})
    if not allow_short:
        bounds = [(0, None,) for i in range(len(x))]
    else:
        bounds = None
    minimize = optimize.minimize(target_func, x, args=(cov_matrix, mean_vector, virtual_mean,), bounds=bounds,
                                 constraints=cons)
    return minimize

但我总是得到成功：错误（超出迭代限制）。我尝试设置 maxiter = 10000 选项，但没有帮助。

如有任何帮助，我将不胜感激

P.S。我用 python 2.7

Answer 1

我不知道为什么，但是当我替换

时它工作得很好

x = np.zeros(len(profits))

有

x = np.ones(len(profits))

使用 scipy.minimze 最大化夏普比率

Maximize Sharpe's ratio using scipy.minimze

python

mathematical-optimization

scipy

nonlinear-optimization