2

我有以下问题：

max CEQ(w) st w in (0,1) 我对 CEQ(w) 一无所知，除了由 CEQ(w) = F(CEQ(w)) 形式的定点方程给出。如果我修正 aw，我可以使用 fzero 函数求解定点方程并获得 CEQ 的值。如果我选择不同的 w，我会得到另一个 CEQ 值。因此，我可以遍历 w 的所有可能值，然后选择给出最高 CEQ 的值。虽然这看起来很糟糕，我想知道我是否可以在 MATLAB 中更有效地做到这一点：我想将我的定点方程的解建模为 w 的函数，但我不知道如何实现它。

更准确地说，这是一个示例代码：

clear all
clc

NoDraws = 1000000;
T_hat = 12;
mu = 0.0058;
variance = 0.0017;
rf = 0.0036;
sim_returns(:,T_hat/12) = T_hat*mu + sqrt(T_hat*variance)*randn(NoDraws,1);

A = 5;
kappa=1;
l=0;
theta = 1 - l*(kappa^(1-A) - 1) *(kappa>1);

CEQ_DA_0 = 1.1; 
CEQ_opt = -1000;
w_opt = 0;

W_T = @(w) (1-w)*exp(rf*T_hat) + w*exp(rf*T_hat + sim_returns(:,T_hat/12));
for w=0.01:0.01:0.99
    W=W_T(w);
    fp = @(CEQ) theta*CEQ^(1-A)/(1-A) - mean( W.^(1-A)/(1-A)) + l*mean( ((kappa*CEQ)^(1-A)/(1-A) - W.^(1-A)/(1-A)) .* (W < kappa*CEQ)); 
    CEQ_DA = fzero(fp,CEQ_DA_0);
    if CEQ_DA > CEQ_opt
        CEQ_opt = CEQ_DA;
        w_opt = w;
    end
end

也就是说，在循环中，我修复了 aw，求解了不动点方程并存储了 CEQ 的值。如果某个其他 w 为 CEQ 提供了更大的值，则当前最优 w 将被新的 w 替换。我想要的（而不是循环部分）是这样的：

fp = @(CEQ,w) theta*CEQ^(1-A)/(1-A) - mean( W_T(w).^(1-A)/(1-A)) + l*mean( ((kappa*CEQ)^(1-A)/(1-A) - W_T(w).^(1-A)/(1-A)) .* (W_T(w) < kappa*CEQ)); 
CEQ_DA = @(w) fzero(fp,CEQ_DA_0);
[w_opt, fval]=fminbnd(CEQ_DA,0,1);
4

1 回答 1

1

您提出的解决方案非常接近。换句话说，您将定义fp为两个参数的函数，并希望 CEQ_DA 是 的函数w，它解决fp了 CEQ 的问题，给定的w. 唯一的问题是fzero不知道要解决 fp 的哪个参数，因为它无法通过名称匹配匿名函数参数和 fp 参数。

答案是在 fzero 中又多了一个匿名函数，fp(CEP,w)变成fp_w(CEP)，这对于 CEQ 是可解的

CEQ_DA = @(w) fzero(@(CEQ) fp(CEQ, w),CEQ_DA_0);
于 2015-04-30T15:17:28.363 回答