情节不作为/行动状态空间

1视图(30天)

显示旧的评论

卡洛斯Fellmann 2023年3月7日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1924680-plot-inaction-action-state-space

回答: Sugandhi 2023年3月20日

我想情节不作为/行动状态空间(A、K)。的问题是,公司选择投资或不投资每段时期,基于预期贴现利润。我用蛮力网格搜索解决问题。如下的代码。

这个想法是情节(K),基于公司投资的决定。

任何帮助将是感谢!

                          % % p 4
                         
                          %管家
                         
                          抽搐
                         
                          清晰的
                         
                          关闭所有
                         
                          clc
                         
                          全球% A0θβp伽马F_cost k0 valueFunction概率j kgrid % valueFunction_a valueFunction_i
                         
                          % % 1校准
                         
                          抽搐
                         
                          θ= 0.592;
                         
                          p = 1;
                         
                          β= 0.90;
                         
                          γ= 0.049;
                         
                          F = 0.039;
                         
                          概率= (0.997 - 0.003 0;
                         
                          0.0003 0.9994 0.0003;
                         
                          0 0.003 - 0.997);
                         
                          = (0.8134 - 1 1.2293)';
                         
                          expectedprod =“*概率;
                         
                          expectedprod = expectedprod ';
                         
                          % k1 = 0.01;
                         
                          % k2 = 100;
                         
                          nGridCapital = 500;
                         
                          A_ss =和(A) /(长度(A));
                         
                          托尔= 1 * 10 ^ (-10);
                         
                          % %稳态
                         
                          首都= (F / A_ss) ^ (1 / (1-theta));
                         
                          kmin = 0.01;% 0.1 *首都;
                         
                          kmax = 100;% 1.1 *首都;
                         
                          % %电网建设
                         
                          vGridCapital = linspace (kmin kmax nGridCapital);
                         
                          % vGridCapital = vGridCapital ';
                         
                          % %矩阵
                         
                          nGridCapital mValueFunction = 0(3日);
                         
                          nGridCapital mValueFunctionnew = 0(3日);
                         
                          nGridCapital mValueFunctionnew_i = 0(3日);
                         
                          nGridCapital mValueFunctionnew_a = 0(3日);
                         
                          nGridCapital mPolicyFunction = 0(3日);
                         
                          nGridCapital mPolicyFunction_i = 0(3日);
                         
                          nGridCapital mPolicyFunction_a = 0(3日);
                         
                          nGridCapital valueProvisional_i = 0(3日);
                         
                          nGridCapital valueProvisional_a = 0(3日);
                         
                          nGridCapital adjcost = 0(3日);
                         
                          nGridCapital投资= 0(3日);
                         
                          nGridCapital收入= 0(3日);
                         
                          范植伟= mValueFunction;
                         
                          % %的收入
                         
                          π= * (vGridCapital ^θ);%预构建收入每一点在首都网格和不同的生产力水平
                         
                          % %的迭代
                         
                          diff_i = 10;
                         
                          diff_a = 10;
                         
                          迭代= 0;
                         
                          diff = max (diff_i diff_a);
                         
                          而(diff > tol)
                         
                          %为资本价值循环
                         
                          为nCapital = 1: nGridCapital
                         
                          %值循环
                         
                          为刺激= 1:尺寸(1)
                         
                          %循环资本在接下来的时期
                         
                          为nCapitalNext = 1: nGridCapital
                         
                          %循环下一时期的生产力
                         
                          为prodnext = 1:尺寸(expectedprod, 1)
                         
                          adjcost (nCapital nCapitalNext) =(γ/ 2)* (vGridCapital (nCapitalNext) -vGridCapital (nCapital) / vGridCapital (nCapital)) ^ 2 * vGridCapital (nCapital);
                         
                          投资(nCapital nCapitalNext) = p * (vGridCapital (nCapitalNext) -vGridCapital (nCapital));
                         
                          收入(nCapital) =刺激* vGridCapital (nCapital)。^θ;
                         
                          valueProvisional_i (prodnext nCapitalNext) =(1测试版)*(收入(nCapital)) +β* (mValueFunction (prodnext nCapitalNext));
                         
                          valueProvisional_a (prodnext nCapitalNext) =(1测试版)*(收入(nCapital) - F * vGridCapital (nCapital) -adjcost (nCapital nCapitalNext)投资(nCapital nCapitalNext)) +β* (mValueFunction (prodnext nCapitalNext));
                         
                          结束
                         
                          结束
                         
                          %发现资本'和'最大限度地为每一个函数,
                         
                          %调整和不活动,考虑到这一时期的生产力
                         
                          mValueFunctionnew_i (prod nCapital) = max (valueProvisional_i (prod nCapital));
                         
                          mPolicyFunction_i (prod nCapital) = max (valueProvisional_i (prod nCapital));
                         
                          mValueFunctionnew_a (prod nCapital) = max (valueProvisional_a (prod nCapital));
                         
                          mPolicyFunction_a (prod nCapital) = max (valueProvisional_a (prod nCapital));
                         
                          %选择两者之间的最大新值函数和政策
                         
                          %的功能
                         
                          mValueFunctionnew = max (mValueFunctionnew_a mValueFunctionnew_i);
                         
                          mPolicyFunction = max (mPolicyFunction_i mPolicyFunction_a);
                         
                          结束
                         
                          结束
                         
                          % mPolicyFunction = vGridCapital (1、mPolicyFunction);%的政策功能转换回单位
                         
                          %的差异
                         
                          diff =规范(mValueFunctionnew-mValueFunction);
                         
                          %设置新值函数
                         
                          mValueFunction = mValueFunctionnew;
                         
                          迭代=迭代+ 1
                         
                          流('迭代= % d,一口= % 2.8 f \ n '、迭代diff)
                         
                          流(' \ n ')
                         
                          结束
                         
                          % %的情节
                         
                          图(2)
                         
                          情节(vGridCapital mPolicyFunction (1:)“红色”)
                         
                          持有在
                         
                          :情节(vGridCapital mPolicyFunction (2),“蓝”)
                         
                          持有在
                         
                          情节(vGridCapital mPolicyFunction (3:)“绿色”)
                         
                          传奇(“低”,”中“,“高”)
                         
                          toc