-2

我正在做我的研究。我无法编写代码来呈现以下矩阵。

A(:,:,1)=[ 1 2 3 4;
           3 2 1 4]
A(:,:,2)=[ 1 3 4 2;
           4 2 1 3]

C(:,:,1)=[ 0 0 1 2;
           1 1 1 0]
C(:,:,2)=[ 1 0 1 0;
           0 1 1 1]

矩阵 C 是显示矩阵 A 中值之前有多少个零的矩阵。

我有以下矩阵来表示矩阵 A 中的值。

value=1, 
B(1)=[1 1 1;
      1 1 1]
value=2,
B(2)=[2 2 2;
      2 2 2]
value=3,
B(3)=[3 3 3;
      4 4 4]
value=4,
B(4)=[4 4 4;
      3 3 3]

这样我将通过组合 A 和 C 得到以下矩阵。

Res(:,:,1)=[1 1 1 2 2 2 0 3 3 3 0 0 4 4 4;
            1 1 1 2 2 2 0 4 4 4 0 0 3 3 3;
            0 3 3 3 0 2 2 2 0 1 1 1 4 4 4;
            0 4 4 4 0 2 2 2 0 1 1 1 3 3 3]

Res(:,:,2)=[0 1 1 1 3 3 3 0 4 4 4 2 2 2 0;
            0 1 1 1 4 4 4 0 3 3 3 2 2 2 0;
            4 4 4 0 2 2 2 0 1 1 1 0 3 3 3;
            3 3 3 0 2 2 2 0 1 1 1 0 4 4 4]

我从其他论坛得到了这个。

CP = bsxfun(@minus, 3 , sum(C,2));
G = mat2cell(C,ones(size(C,1),1),ones(size(C,2),1),ones(size(C,3),1));
CP = mat2cell(CP,ones(size(CP,1),1),ones(size(CP,2),1),ones(size(CP,3),1));
D = B(A);
E = cellfun(@(x,n) [zeros(size(x,1),n) x], D, G,'UniformOutput',false);
for ii=1:depth;
    E(:,end,ii) = cellfun(@(x,n) [x,zeros(size(x,1),n)],E(:,end,ii),CP(:,:,ii),...
   'UniformOutput',false);
    NewMatrice(:,:,ii)=cell2mat(E(:,:,ii));
end

你有什么意见？

4

1 回答 1

2

对于初学者，我不明白您的矩阵如何同时A包含第一个A = [ 1 2 3 4; 3 2 1 4]和A(1)、A(2)和。A(3)A(4)

我想你有以下内容：

A = [1, 2, 3, 4; 3, 2, 1, 4];
B(:, :, 1) = [1, 1, 1; 1, 1, 1];
B(:, :, 2) = [2, 2, 2; 2, 2, 2];
B(:, :, 3) = [3, 3, 3; 4, 4, 4];
B(:, :, 4) = [4, 4, 4; 3, 3, 3];

假设这一点，您可以通过以下方式获得所需的结果：

Res = cell2mat(arrayfun(@(n)B(:, :, n), A, 'UniformOutput', false));

产生：

Res =

     1     1     1     2     2     2     3     3     3     4     4     4
     1     1     1     2     2     2     4     4     4     3     3     3
     3     3     3     2     2     2     1     1     1     4     4     4
     4     4     4     2     2     2     1     1     1     3     3     3

编辑＃1：

在更新您的问题之后，如果您有另一个矩阵C = [0, 0, 1, 2; 1, 1, 1, 0]表示要在B矩阵之前插入的前导零的数量，您可以通过以下方式获得所需的结果：

my_func = @(n)[zeros(size(B(:, :, A(n)), 1), C(n)), B(:, :, A(n))];
Res = cell2mat(arrayfun(my_func, reshape(1:numel(A), size(A)), 'Uniform', false))

这产生：

Res =

  Columns 1 through 15

     1    1    1    2    2    2    0    3    3    3    0    0    4    4    4
     1    1    1    2    2    2    0    4    4    4    0    0    3    3    3
     0    3    3    3    0    2    2    2    0    1    1    1    4    4    4
     0    4    4    4    0    2    2    2    0    1    1    1    3    3    3

编辑＃2：

在对您的问题进行第二次更新之后，如果您的A和C矩阵是 3 维的，这会增加代码的很多复杂性。我发现最容易使用for循环并将问题简化为二维，如下所示：

for k = 1:size(A, 3)
    A2 = A(:, :, k);
    C2 = C(:, :, k);
    my_func = @(n)[zeros(size(B(:, :, A2(n)), 1), C2(n)), B(:, :, A2(n))];
    Res(:,:,k) = cell2mat(arrayfun(my_func,reshape(1:numel(A2),size(A2)),'Un',0))
end

Res如果你想在每次迭代中使用不同的大小，我建议你制作Res一个单元格数组，即：

Res{k} = cell2mat(arrayfun(my_func,reshape(1:numel(A2),size(A2)),'Un',0))

PS
请注意，对于第二次迭代，C = [1, 0, 1, 0; 0, 1, 1, 1]您会遇到问题，因为第二行不会有与第一行相同数量的零。我相信你的意思是这样的：C = [1, 1, 1, 0; 0, 1, 1, 1]。

编辑#3：

这是生成随机矩阵的代码C，其中每一行总和为一个预定义的整数K（也可以任意选择）：

K = fix(size(A, 2) * (1 + rand));                         %# Pick a random K
C = rand(size(A));                                        %# Generate a random C
C = ceil(C * K ./ repmat(sum(C, 2), 1, size(A, 2))) - 1;  %# Normalize
C(:, end) = C(:, end) + K - sum(C, 2);                    %# Fix illegal rows
于 2012-09-06T12:57:58.400 回答