Question

我感觉我在我的解决方案中缺少一些直观的东西，用于生成部分变化的块对角网格。在任何情况下，我想摆脱我的功能循环（为了挑战...）

给定参数元组，间隔数和百分比变化：

params = [100 0.5   1
           24   1 0.9];
nint   = 1; 
perc   = 0.1;

所需的输出应为：

pspacegrid(params,perc,nint)
ans =
   90.0000    0.5000    1.0000
  100.0000    0.5000    1.0000
  110.0000    0.5000    1.0000
  100.0000    0.4500    1.0000
  100.0000    0.5000    1.0000
  100.0000    0.5500    1.0000
  100.0000    0.5000    0.9000
  100.0000    0.5000    1.0000
  100.0000    0.5000    1.1000
   21.6000    1.0000    0.9000
   24.0000    1.0000    0.9000
   26.4000    1.0000    0.9000
   24.0000    0.9000    0.9000
   24.0000    1.0000    0.9000
   24.0000    1.1000    0.9000
   24.0000    1.0000    0.8100
   24.0000    1.0000    0.9000
   24.0000    1.0000    0.9900

您可以在其中看到此掩码所表示的值发生变化：

mask =
     1     0     0
     1     0     0
     1     0     0
     0     1     0
     0     1     0
     0     1     0
     0     0     1
     0     0     1
     0     0     1
     1     0     0
     1     0     0
     1     0     0
     0     1     0
     0     1     0
     0     1     0
     0     0     1
     0     0     1
     0     0     1

函数pspacegrid()是：

function out = pspacegrid(params, perc, nint)
% PSPACEGRID Generates a parameter space grid for sensitivity analysis

% Size and number of variation steps
sz     = size(params);
nsteps = nint*2+1;

% Preallocate output
out    = reshape(permute(repmat(params,[1,1,nsteps*sz(2)]),[3,1,2]),[],sz(2));

% Mask to index positions where to place interpolated 
[tmp{1:sz(2)}] = deal(true(nsteps,1));
mask           = repmat(logical(blkdiag(tmp{:})),sz(1),1);

zi = cell(sz(1),1);
% LOOP per each parameter tuple
for r = 1:sz(1)
    % Columns, rows, rows to interpolate and lower/upper parameter values
    x     = 1:sz(2);
    y     = [1; nint*2+1];
    yi    = (1:nint*2+1)';
    z     = [params(r,:)*(1-perc); params(r,:)*(1+perc)];
    % Interpolated parameters 
    zi{r} = interp2(x,y,z, x, yi);
end
out(mask) = cat(1,zi{:});

Answer 1

我想我得到了它，建立了你的预循环代码：

params = [100 0.5   1
           24   1 0.9];
nint   = 1; 
perc   = 0.1;

sz     = size(params);
nsteps = nint*2+1;

% Preallocate output
out    = reshape(permute(repmat(params,[1,1,nsteps*sz(2)]),[3,1,2]),[],sz(2));

%Map of the percentage moves
[tmp{1:sz(2)}] = deal(linspace(-perc,perc,nint*2+1)');
mask = repmat(blkdiag(tmp{:}),sz(1),1) + 1; %Add one so we can just multiply at the end

mask.*out

因此，不要让你的面具复制我制作的面具，而是复制每个元素的百分比移动，这是一个重复的模式，基本元素是这样的：

linspace(-perc,perc,nint*2+1)'

然后就像将1添加到整个事物并乘以out矩阵一样简单

我按如下方式测试了它：

me = mask.*out;
you = pspacegrid(params, perc, nint);

check = me - you < 0.0001;

mean(check(:))

当我摆弄输入时似乎工作。但是我的功能确实出错了，我不得不将true(...)更改为ones(...)。这可能是因为我在网上运行它可能使用Octave而不是Matlab。

块对角线变化网格

1 个答案: