0

2つの入力変数があります。

N=10;
M=4;

Matlabを使用して、次の行列を生成したいと思います-forループを使用せずに、任意のMとNに対して：

    %N = 1   2   3   4   5   6   7   8  9
Mat1 = [2,  3,  4,  5,  6,  7,  8,  9, 10;      %M=1 -> Mat1(1, i) = N(i)+1
        3,  6,  9, 12, 15, 18, 21, 27, 30;      %M=2 -> Mat1(2, i) = N(i)+N(i)*2
        4,  8, 12, 16, 20, 24, 28, 32, 36;      %M=3 -> Mat1(2, i) = N(i)+N(i)*3
        5, 10, 15, 20, 25, 30, 35, 40, 45]      %M=4 -> Mat1(2, i) = N(i)+N(i)*4

   %N = 1   2   3   4   5   6   7   8  9
Mat2 = [1, 2,  3,  4,  5,  6,  7,  8,  9;       %M=1 -> Mat2(2, i) = N(i)
        1, 4,  7, 10, 13, 16, 19, 25, 28;       %M=2 -> Mat2(2, i) = N(i)+N(i)*2-2
        1, 5,  9, 13, 17, 21, 25, 29, 33;       %M=3 -> Mat2(2, i) = N(i)+N(i)*3-3
        1, 6, 11, 16, 21, 26, 31, 36, 41]       %M=4 -> Mat2(2, i) = N(i)+N(i)*4-4

概要：

% i=1:N

%Mat1(M, i) = N(i)+N(i)*M    , if M>1
%             N(i)+1         , if M=1

%Mat2(M, i) = N(i)+N(i)*M-M  , if M>1
%             N(i)           , if M=1

2つのforループを使用してコードを作成しましたが、ループを使用せずにこれらの結果を取得する方法があるかどうかを確認したいと思います。bsxfunとを試してみましarrayfunたが、正しい結果を得ることができませんでした。

clear Mat1 Mat2
N=10;
M=9;

Mat1 = ones(M, N);
Mat2 = ones(M, N);
for i=1:M
    for j=1:N
        if  i==1
            Mat1(i, j) = j+1;
            Mat2(i, j) = j;
        elseif j ==1  
            Mat1(i, j) = j+j*i;
        elseif i~=1 && j~=1
            Mat1(i, j) = j+j*i;
            Mat2(i, j) = j+j*i-i;
        end
        
    end
end

ありがとうございました

4

4 に答える 4

3

私の試み：

N=10;
M=4;
mat1 = bsxfun(@(x,y) x*(y.*(y>1)+1)+(y==1), 1:N,(1:M)')
mat2 = bsxfun(@(x,y) x*(y.*(y>1)+1)-(y>1).*y, 1:N,(1:M)')
于 2013-02-21T10:34:13.987 に答える
1

私はこれがあなたが探していることをするかもしれないと思います：

最初のマトリックス：

Mat1 = (2:N+1)'*(1:N);
Mat1 = Mat1(1:M,:);
Mat1(1,:) = 2:N+1;

これにより、（N=10、M=4）が得られます。

Mat1 =

 2     3     4     5     6     7     8     9    10    11
 3     6     9    12    15    18    21    24    27    30
 4     8    12    16    20    24    28    32    36    40
 5    10    15    20    25    30    35    40    45    50

2番目のマトリックス：

tmp = (0:N-2)'*(1:N-1);
Mat2 = repmat((4:3:3*N),M-1,1) + tmp(1:M-1,:);
Mat2 = [ones(M,1), [2:N; Mat2]];

これは次のようになります。

Mat2 =

 1     2     3     4     5     6     7     8     9    10
 1     4     7    10    13    16    19    22    25    28
 1     5     9    13    17    21    25    29    33    37
 1     6    11    16    21    26    31    36    41    46

これらのマトリックスを構築するためのルールが正確に何であるかがあまり明確ではないため、これは少しパズルでした。しかし、必要なのがこれらの行列だけである場合、これはそれらを取得するための高速な方法であるはずです。

これらのコード行はすべてN、ループを使用してコーディングする場合と同じ行列を生成します。Mfor

于 2013-02-20T11:26:55.790 に答える
1

その間、私は使い方を学びましbsxfunたrepmat。Mat1ジェネレーターとMat2ジェネレーターの両方のコンパクトなコードを入手しました。

N=10;
M=4;

Mat1 = bsxfun(@times, (1:M)', (1:N)) +...
      [ones(1, N)  ;  horzcat(  ones(M-1,1) , repmat(2:N,M-1,1)   )   ];

Mat2 = [ones(M,1), [(2:N); bsxfun(@plus, (2:M)',(2:N)) +...
        horzcat(  zeros(M-1,1),  bsxfun(@times, (2:M)',(1:N-2))  )   ] ];

Mat1 =

 2     3     4     5     6     7     8     9    10    11
 3     6     9    12    15    18    21    24    27    30
 4     8    12    16    20    24    28    32    36    40
 5    10    15    20    25    30    35    40    45    50

Mat2 =

 1     2     3     4     5     6     7     8     9    10
 1     4     7    10    13    16    19    22    25    28
 1     5     9    13    17    21    25    29    33    37
 1     6    11    16    21    26    31    36    41    46
于 2013-02-21T09:36:31.497 に答える
0

私はすべてのテストケースでこれらを試したわけではないので、それが完全に一般的であるかどうかはわかりません。私がこれを行っていたものを調査し、最初のマトリックスの解にも到達できるはずだと思います。上記のサンプルマトリックスにタイプミスがあるかもしれませんが、これで正しい方向に進むことができれば幸いです。

 M=4;
 N=10;

 Mat2=repmat([1:N],M,1);
 Mat2(2:M,:)=Mat2(2:M,:)+bsxfun(@times, Mat2(2:M,:), [2:M]');
 Mat2(2:M,:)=bsxfun(@minus, Mat2(2:M,:), [2:M]');

私たちに与える：

Mat2 =

 1     2     3     4     5     6     7     8     9    10
 1     4     7    10    13    16    19    22    25    28
 1     5     9    13    17    21    25    29    33    37
 1     6    11    16    21    26    31    36    41    46

M=5およびN=10の場合

Mat2 =

 1     2     3     4     5     6     7     8     9    10
 1     4     7    10    13    16    19    22    25    28
 1     5     9    13    17    21    25    29    33    37
 1     6    11    16    21    26    31    36    41    46
 1     7    13    19    25    31    37    43    49    55

これはMat1だと思います。

 Mat1=repmat([1:N],M,1);
 Mat1(2:M,:)=Mat1(2:M,:)+bsxfun(@times, Mat1(2:M,:), [2:M]');
 Mat1(1,:)=Mat1(1,:)+1;

Mat1 =

 2     3     4     5     6     7     8     9    10    11
 3     6     9    12    15    18    21    24    27    30
 4     8    12    16    20    24    28    32    36    40
 5    10    15    20    25    30    35    40    45    50
于 2013-02-20T11:59:32.637 に答える