特定の行列の置き換えに方法ついて

Aという10行50列のデータのうち、Bという3行50列に格納されている行番号のデータのみを０に置きかえる方法はありますでしょうか。

Bには、0<10が入っています。

イメージとしては、B = [1;3;5, 2;4;6, 8;7;4,....];

の場合、Aの１列目の1,3,5行が０に、2列目の2,4,6行が0に...

と出来れば思います。

Answer 1

[Atsushi Ueno]

A = rand(10,50); % Aという10行50列のデータ
B = randi(size(A,1),3,50); % Bという3行50列のデータ(0<10)
B(1:3,1:3) = [1,2,8;3,4,7;5,6,4]; % イメージとしては、B = [1;3;5, 2;4;6, 8;7;4,....];
for j = 1:size(B,2)
for k = 1:size(B,1)
    A(B(k,j),j) = 0; % Bという3行50列に格納されている行番号のデータ(0<10)のみを０に置きかえる
end
end
A % Aの１列目の1,3,5行が０に、2列目の2,4,6行が0に...

A = 10×50

0 0.7718 0.7895 0.1420 0 0.5693 0.4373 0.4954 0 0 0.7592 0.8207 0.3226 0.3680 0.1719 0.5081 0.7613 0 0 0.4406 0.8574 0.9051 0 0.0606 0.3038 0.6038 0.7309 0.4929 0.1167 0.1499 0.0597 0 0.8033 0 0.2491 0.2651 0.1170 0.2995 0.6122 0.3208 0.9350 0 0.0262 0.5147 0.1819 0.9510 0.2624 0.9878 0 0.7489 0.9119 0.9811 0.4622 0.6916 0.8696 0.6368 0 0.2488 0.6778 0.3306 0 0.0289 0.3459 0.4727 0.9457 0.3140 0 0.1217 0.7499 0.5453 0 0 0.9879 0.5238 0.2765 0.5372 0 0 0.2808 0.1379 0 0 0.6291 0 0.5940 0 0.8516 0.0977 0 0.7673 0.8414 0 0 0.6882 0 0.2365 0.6656 0.9974 0.6479 0.7240 0.1648 0 0 0 0 0.2182 0 0.5141 0.6374 0.4489 0.7500 0 0.1799 0.4598 0 0.5162 0.1979 0 0.5856 0 0 0.5431 0.8756 0 0.6322 0 0 0.4227 0.3913 0.7080 0.7437 0.0055 0.3072 0.1676 0.1763 0.6099 0.9582 0.4855 0.2154 0 0.7996 0 0 0.0206 0.4523 0.5729 0.5933 0.2590 0.3823 0.4316 0.3775 0 0.0405 0 0.0574 0.0673 0.6523 0 0.3517 0 0.4356 0.2352 0 0.9575 0.0908 0.0005 0.0950 0.0451 0.9584 0.3804 0.8418 0.5794 0.5128 0.2612 0.5995 0.5951 0.2944 0.4499 0.2436 0 0.0717 0.6202 0 0.9999 0.4440 0.7030 0.0393 0.7488 0.3140 0.5701 0.6889 0.8160 0.5022 0.1222 0.8667 0.5064 0.3532 0 0.9160 0.6394 0 0.8477 0.9690 0.3574 0.6484 0.7077 0 0.9140 0 0.1547 0.5234 0.5710 0 0.7823 0.0430 0 0.0299 0.5471 0 0 0.0358 0.8269 0.1276 0.4183 0.9646 0 0 0.8720 0.2193 0 0.6064 0.3184 0.4654 0 0 0 0.0294 0.4956 0.0299 0.4631 0.6826 0.9556 0.5362 0.8360 0.6736 0.4998 0.9703 0 0.8535 0.6667 0.7108 0.4505 0.7293 0 0 0.5153 0.0007 0.5709 0 0.9282 0.6700 0.4001 0.6387 0.1857 0.4951 0.7932 0 0 0 0.8176 0.9923 0.4814 0.7679 0.6291 0 0.3932 0 0 0 0.7987 0 0.8770 0 0.2043 0 0 0.2565 0.5743 0.0042 0.5482 0 0.7701 0 0.5745 0.0743 0.8791 0.8381 0 0.4316 0.5902

[Atsushi Ueno]

for 文を使ったら負けシリーズ第一弾「for 文を使うな!」

A = rand(10,50); % Aという10行50列のデータ
B = randi(size(A,1),3,50); % Bという3行50列のデータ(0<10)
B(1:3,1:3) = [1,2,8;3,4,7;5,6,4]; % イメージとしては、B = [1;3;5, 2;4;6, 8;7;4,....];

B = [B(:), repelem(1:size(A,2),3)']; % Bを縦に伸ばして列番号を付ける
for k = 1:size(B,1) % Bという3行50列に格納されている行番号のデータ(0<10)のみを０に置きかえる
    A(B(k,1),B(k,2)) = 0;
end
A % Aの１列目の1,3,5行が０に、2列目の2,4,6行が0に...

A = 10×50

0 0.7497 0.1425 0.8855 0.9225 0.4005 0 0 0.7609 0.9561 0.3972 0.6424 0 0.7142 0 0.4531 0.3648 0.1243 0.6602 0 0 0 0.8697 0.4928 0.5289 0.9496 0 0 0.3475 0.8097 0.3496 0 0.7837 0 0.5754 0.9276 0.8328 0 0.4465 0.2761 0.4433 0.9242 0.8425 0 0 0.7860 0.9364 0.4157 0 0.5302 0.7360 0.2727 0 0 0.0275 0 0 0.0525 0.4596 0.8992 0 0.3540 0.1962 0.8181 0.5055 0.4052 0.2151 0.1939 0.8956 0.4176 0.6184 0 0.0377 0.2790 0.7206 0.8905 0 0 0 0.5874 0 0.2425 0 0.2361 0.9852 0.4720 0.4507 0.7465 0.6323 0 0.1978 0 0 0.6620 0.7684 0 0.8433 0.3068 0.9981 0 0.0916 0.3044 0 0 0.3423 0.3669 0 0.3928 0.4333 0.6058 0.2946 0 0.9750 0 0 0 0.0598 0.2841 0.0097 0.5426 0 0.2919 0.9235 0.3797 0.5070 0.5660 0.7406 0.0758 0 0.2263 0 0 0.9833 0.0835 0.7185 0 0 0.2856 0.2963 0.7383 0.9326 0.2200 0.2816 0.3468 0.9429 0.5471 0.0082 0.8278 0 0 0.7733 0 0.0236 0.1386 0.2934 0.0711 0.5768 0.1185 0.5326 0.4534 0 0.3542 0.4656 0.4532 0 0.5724 0.7682 0.8891 0.3198 0 0.1935 0.4248 0.7839 0.8530 0.5289 0.7125 0.6174 0.1189 0.7328 0.1458 0.8384 0.6291 0 0.0466 0 0.2859 0.8237 0.3882 0.8707 0 0.0331 0.6833 0.0712 0.8605 0.0686 0 0.3896 0.2588 0.6668 0 0 0.5430 0.2066 0 0.2986 0 0.6372 0.0532 0 0.2095 0.4579 0.7434 0 0.7851 0.8136 0.0469 0.5070 0.2482 0.1310 0 0.1963 0.4046 0 0.5282 0.4422 0.1093 0.7941 0 0.3302 0.9649 0.9885 0.9928 0.5302 0.4871 0 0.9072 0.7696 0.3190 0 0 0.1859 0.5809 0.7748 0.2628 0 0 0.2630 0 0 0.5350 0 0.2762 0.8124 0.5039 0.7834 0 0.3574 0.1599 0.6846 0.6716 0.3821 0 0 0.7196 0.2743 0.9117 0.5957 0 0.7799 0.0261 0.7194 0.6962 0.6026 0 0.2470 0.4075 0 0.1083 0 0.0391 0.1544 0.8159 0.3287 0.8662 0.8170 0.4367 0.0146 0 0 0.5257 0.2836 0.5829 0.1786 0.0537 0.0696 0.5383 0 0.7752 0.2460 0.8630

[Atsushi Ueno]

for 文を使ったら負けシリーズ第二弾「for 文を使うな! 2」

A = rand(10,50); % Aという10行50列のデータ
B = randi(size(A,1),3,50); % Bという3行50列のデータ(0<10)
B(1:3,1:3) = [1,2,8;3,4,7;5,6,4]; % イメージとしては、B = [1;3;5, 2;4;6, 8;7;4,....];

B = sub2ind(size(A), B(:), repelem(1:size(A,2),3)'); % Bを添字から線形インデックスに変換する
A(B) = 0 % Aを線形インデックスで参照して０にする

A = 10×50

0 0.1357 0.8594 0.0817 0.9012 0.9353 0.5676 0.0249 0.8765 0.9944 0 0 0.9033 0.0198 0.4167 0.7282 0.0071 0 0.1230 0 0 0 0.6772 0.0765 0.2088 0.2291 0.7593 0.1517 0.3637 0 0.7943 0 0.8825 0.7574 0.1642 0 0.0121 0 0.8586 0.7842 0 0.2693 0.6650 0.6643 0 0 0.6630 0.0152 0.9455 0.1156 0.9668 0.7847 0.1166 0 0.8238 0.0791 0.1323 0.5051 0.7100 0 0 0.1594 0.3565 0.0363 0.5042 0.2520 0 0 0 0.5180 0.7266 0.8572 0.1699 0.9795 0.8283 0.2424 0.3086 0 0.7127 0 0.8867 0.2969 0.5249 0.6126 0.2115 0.8572 0.5446 0.4756 0 0.7542 0.3852 0 0 0.0422 0 0.9528 0.6489 0.3605 0.8136 0.1312 0.3351 0.1636 0.0892 0.2708 0 0.6772 0 0.8242 0 0.5022 0.8300 0.1373 0 0 0 0 0.2354 0.2821 0.5271 0.5323 0 0.4808 0.0168 0.2407 0.5110 0.2110 0.5067 0.7097 0 0.5665 0.6760 0.4823 0 0 0.8300 0.8586 0.5059 0.2325 0.9588 0.9664 0.8051 0 0.4566 0 0 0 0.2838 0.5499 0.2581 0.5155 0.1718 0 0.6915 0 0.2333 0.4086 0.8220 0.0670 0.3443 0 0.2418 0.4943 0 0.6556 0 0.5261 0.1519 0.5556 0.8818 0.3320 0 0.8688 0.5160 0.8026 0.9965 0.4707 0 0.4135 0.0141 0.9135 0.6842 0.4798 0 0.7064 0 0 0.4340 0.8516 0.2051 0.1118 0.3980 0 0.1107 0 0.6546 0.0834 0.6001 0.6534 0 0 0.4385 0 0.1284 0.1555 0.9127 0.8374 0.9816 0.1734 0 0.2541 0.2558 0.5228 0 0.9257 0 0 0 0.9809 0.3137 0 0.5758 0.1306 0.3039 0.5355 0.1078 0.1498 0 0 0.4261 0.7349 0.5894 0.6804 0.1740 0.6219 0.9862 0.7325 0 0.5245 0.7625 0.0855 0.0919 0.6812 0.6527 0 0.7233 0.6858 0.1016 0 0.4064 0.6110 0 0 0.4780 0.6361 0.6517 0 0.7445 0.6611 0 0.6320 0.6082 0.4250 0 0.4658 0.8303 0 0.8014 0 0 0 0.0750 0.4096 0.6157 0.1658 0.4854 0.6870 0.3486 0.1738 0 0 0.6959 0.4733 0 0 0.0842 0.7070 0.9346 0.5912 0.6954 0.7525 0 0.1295 0.9545 0.4123 0 0.7158 0.8816 0 0.9037 0.7394

[flower]

やりたいことが実現できました。バリエーション豊富でどれも思いつかないものでした...すごいです。

ありがとうございました。