朱莉娅在循环中更新变量

Question

我有一个代码来编程值函数迭代。代码是从 Matlab 翻译而来的。第一次迭代很好，但在第二次迭代中，旧的值函数在某个时候被错误地更新了，尽管第二次迭代仍然从 matlab 产生等效的结果。

首先，我在这里声明变量：

sigma = 1.5;
delta = 0.1;
beta  = 0.95;
alpha = 0.30;
ks    = 1.0;
csy   = 0.2;
kmin  = 0.1*ks;
kmax  = 1.1*ks;
nbk   = 20;
devk  = (kmax-kmin)/(nbk-1);
k     = linspace(kmin,kmax,nbk)';
v0    = zeros(nbk,1);
v0    = (csy*k.^alpha).^(1-sigma);
v0    = v0';
v     = zeros(nbk,1);
ik1   = zeros(nbk,1);
iter  = 1;
crit  = 1;
tol   = 1e-6;

其次，我执行这个(这是一个while循环(，但我更喜欢一个接一个地执行它(：

 for i=1:nbk
      imin =  max(ceil(((1-delta)*k[i]-kmin)/devk)+1.0,1);
      imax =  min(floor((k[i]^alpha+(1-delta)*k[i]-kmin)/devk)+1.0,nbk);
      imin=trunc(Int, imin);
      imax=trunc(Int, imax);
      c = k[i]^alpha+(1-delta)*k[i]-k[imin:imax];
      u = (c.^(1-sigma)-1)/(1-sigma);
      (v[i],itmp)= findmax(u+beta*v0[imin:imax]);
      ik1[i] = imin-1+itmp;
   end;
  v0  = v;

第一次运行第二部分会正确打印 v 和错误：

[-9.41322; -6.19362; -4.35414; -3.01282; -1.95198; -1.07545; -0.330322; 0.316098; 0.885574; 1.39337; 1.85067; 2.29198; 2.69438; 3.06364; 3.40432; 3.72013; 4.0141; 4.28875; 4.54618; 4.79853]
error:5.124021319754114

第二次运行它会导致错误，因为现在 v0 变成了我不知道为什么等于 v(尽管 v 仍然正确生成，就像在 Matlab 中一样(：

[-7.06807; -4.7177; -3.24724; -2.16729; -1.30218; -0.578569; 0.0429673; 0.586699; 1.09022; 1.54641; 1.98317; 2.37955; 2.7416; 3.07412; 3.38098; 3.66532; 3.92974; 4.18466; 4.42981; 4.66936]
error:0.0

所以误差为零，两次迭代后不应该是这种情况。有什么想法吗？

Answer 1

当你设置v0 = v时，两个变量都引用内存中的同一区域，即当你尝试修改一个变量时，你实际上同时修改了两个变量。

使用copy!(v0, v)将v的内容复制到v0中。然后可以独立修改数组。