形状変換

現在19×19の行列が231250個あります。

この中身は全て1と０の数字しかなく対角部分は数字が同じで対角線上は必ず1が入るようにしています。

そして対角になる部分だけを取り出して行列の形状変換をお願いしたいです。

行列が多いので３×３で下記で説明するので、それを１９×１９を２３１２５０個でどのように書けばよいかスクリプトを教えていただきたいです。

3×3行列が３個ある場合で説明

１　１　０　　　　　１　０　１　　　　１　１　０

１　１　０　　　　　０　１　１　　　　１　１　１

０　０　１　　　　　１　１　１　　　　０　１　１

この１行２列目と１行３列目と２行３列目の数字だけ取り出して１行にまとめる

（１，０，０）

（０，１，１）

（１，０，１）

このように三次元のものを２次元に変更したいです

１９×19の231250個分の場合171×231250行列になるはずです。

よろしくお願いいたします

rng(1,'philox')

X = randi([0 1], 19, 19, 231250);

乱数使用で１と0の行列は作りました。

Answer 1

[Shunichi Kusano]

reshapeによる配列変形で楽に処理可能です。

clear
wh = 19;
num = 10000;
X = randi([0 1], wh, wh, num); % ダミーデータ(どうせ使わないので対角性は考慮していない)

% 下対角成分だけを抽出するためのマスク
mask = double(~triu(ones(wh)));
mask(mask==0) = nan;

% マスキング
X_ = X .* mask;

% 変形
X_ = reshape(X_,wh^2,num);
X_(isnan(X_(:,1)),:) = []; % マスクした要素を行ごと消去
X_ = X_'; % 必要に応じて