2

假设我有一个包含 3 个元素的单元格数组;第一个和第三个单元格是矩阵,而第二个是空的:

>> P
P = 
  [3x8 double]    []    [3x8 double]

我可以通过对矩阵中的每个单独元素进行插值来填充第二个矩阵,但肯定有办法一次性做到这一点吗?

编辑:平均值肯定适用于P。但假设我有类似的东西

>> Q
Q = [3x8 double]    [3x8 double]     [3x8 double]    []    []    [3x8 double]

因此,鉴于我的四个已知矩阵,我可以使用更高级的插值(如样条)来填充缺失的矩阵。

4

3 回答 3

4

您可以在对 INTERP1 的单个函数调用中执行所有插值,但它需要将数据设置为所需的形状,然后在完成后返回:

%# cell array of matrices (some are missing)
sz = [3 8];         %# all of same size
Q = {ones(sz), [], 3*ones(sz), [], [], 6*ones(sz)};

%# indices of missing matrices
idx = ~cellfun(@isempty,Q);
x = 1:numel(Q);

%# merge cells into a multidimensional matrix, call INTERP1, then divide back
QQ = Q(idx);
QQ = permute(cat(3,QQ{:}), [3 1 2]);
QQ = interp1(x(idx), QQ, x);            %# one call to interpolation function
QQ = reshape(num2cell(permute(QQ, [2 3 1]), [1 2]), 1,[]);

请记住,矢量化有时会以牺牲可读性为代价......

于 2012-05-31T23:23:51.357 回答
2

像这样的东西适用于您的玩具箱:

l = length(P{1}(:)); m = size(P{1}, 1); n = size(P{1}, 2);
myfun = @(x) interp1([1 3], [P{1}(x) P{3}(x)], 2)
P{2} = reshape(arrayfun(myfun, 1:l), m, n);

更通用的代码(已验证!!)通常可以工作。您可以在一个 interp1 调用中指定要使用的插值类型:

[m, n] = size(P{1}); l = m*n;
avbl = find(not(cellfun('isempty', P)));
missing = find(cellfun('isempty', P));

extractor = @(k) cellfun(@(x, k) x(k), {P{avbl}}, cellmat(1, length(avbl), 1, 1, k))
myfun = @(x) interp1(avbl, extractor(x), missing)
intermediate = reshape(arrayfun(myfun, 1:l, 'UniformOutput', false), m, n);
ext2 = @(k) cellfun(@(p) p(k), intermediate);
P(missing) = arrayfun(ext2, 1:length(missing), 'UniformOutput', false);

如果你想让我解释一些部分,请告诉我。它基本上只是调用 cellfun 几次,然后重塑输出。

于 2012-05-31T22:37:25.870 回答
0

你的意思是:

     p{2} = (p{1} + p{3}) / 2;  % note the curly braces

?

于 2012-05-31T22:39:41.040 回答