Re: [問題] 提取矩陣中大於0的數值到另一個矩陣

作者: celestialgod (天)   2015-05-04 16:34:25
※ 引述《shadow5402 (shadow5402)》之銘言:
: 假設我目前有一個varible x
: x =
: 0 0 1 2
: 0 3 2 0
: 1 3 0 0
: 2 0 0 0
: 我希望將其中不等於0的部分移動到另外一varible y
: 變成
: y =
: 1 3 1 2
: 2 3 2 0
: 目前我的寫法是
: for k = 1:4
: for a = 1:2
: if x(k,1) > 0
: y(a,1) = x(k,1);
: end
: if x(k,2) > 0
: y(a,2) = x(k,2);
: end
: if x(k,3) > 0
: y(a,3) = x(k,3);
: end
: if x(k,4) > 0
: y(a,4) = x(k,4);
: end
: end
: end
: 但結果卻是
: y =
: 2 3 2 2
: 2 3 2 2
: 請問我是哪裡寫錯了?又該怎麼改呢?
: 謝謝。
你最主要的錯誤是不能用a那個迴圈...
會重複assign x中非零的值到y去
y = zeros(max(sum(x>0)), size(x, 1));
for k = 1:size(x, 2)
tmp = x(:, k) > 0;
if any(tmp)
y(1:sum(tmp), k) = x(tmp, k);
end
end
% 上面迴圈中,可以用另外一個迴圈代替,
% tmp = x(x(k,:) > 0, k);
% for a = 1:length(tmp)
% y(a, k) = tmp(a);
% end
tmp = arrayfun(@(y) y{:}(y{:}~=0), mat2cell(x, size(x, 1), ones(1,
size(x,2))), ...
'UniformOutput', false);
maxSize = max(cellfun(@length, tmp));
y2 = cell2mat(cellfun(@(x) [x; zeros(maxSize-length(x), 1)], tmp, ...
'UniformOutput', false));
isequal(y, y2) % 1
作者: shadow5402 (shadow5402)   2015-05-04 16:42:00
非常謝謝你的幫忙!不好意思,我剛剛發現結果跟預期的有點不同。會變成 y = [1 3 1 2; 2 2 3 0],而不是[1 3 1 2; 2 32 0]我想到了,改成y(1:sum(tmp), k) = x(tmp, k)就可以了感謝你的幫忙!或是把tmp = x(k, x(k,:) > 0);改為tmp = x(x(:,k)>0, k);不過我現在遇到一個問題,如果我要處理的矩陣不再是4x4,而是5x4,4x3甚至是120x4這類型的矩陣呢?我試了幾種方式,目前失敗中Orz呃..我成功了XD,用第二種方式就可以了,根本不用改第一種也行,改成tmp = x(:,k) > 0;

Links booklink

Contact Us: admin [ a t ] ucptt.com