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。然后可以独立修改数组。

Julia在循环中更新变量

1 个答案: