作者celestialgod (攸藍)
看板MATLAB
標題Re: [問題] 提取矩陣中大於0的數值到另一個矩陣
時間Mon May 4 16:34:25 2015
※ 引述《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
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.205.27.107
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1430728467.A.BBF.html
※ 編輯: celestialgod (123.205.27.107), 05/04/2015 16:41:38
1F：推 shadow5402: 非常謝謝你的幫忙！ 05/04 16:42
2F：推 shadow5402: 不好意思，我剛剛發現結果跟預期的有點不同。 05/04 17:02
3F：→ shadow5402: 會變成 y = [1 3 1 2; 2 2 3 0]，而不是[1 3 1 2; 2 3 05/04 17:03
4F：→ shadow5402: 2 0] 05/04 17:03
5F：推 shadow5402: 我想到了，改成y(1:sum(tmp), k) = x(tmp, k)就可以了 05/04 17:24
6F：→ shadow5402: 感謝你的幫忙！ 05/04 17:24
7F：推 shadow5402: 或是把tmp = x(k, x(k,:) > 0);改為 05/04 18:12
8F：→ shadow5402: tmp = x(x(:,k)>0, k); 05/04 18:12
這部分，很抱歉，我那時候沒很仔細看，沒發現你是要每個column...
我寫成了每個row，我稍作修改之後，就會是正確的結果了
9F：推 shadow5402: 不過我現在遇到一個問題，如果我要處理的矩陣不再是 05/04 18:16
10F：→ shadow5402: 4x4，而是5x4,4x3甚至是120x4這類型的矩陣呢？ 05/04 18:17
11F：→ shadow5402: 我試了幾種方式，目前失敗中Orz 05/04 18:17
12F：推 shadow5402: 呃..我成功了XD，用第二種方式就可以了，根本不用改 05/04 18:26
13F：推 shadow5402: 第一種也行，改成tmp = x(:,k) > 0; 05/04 18:32
※ 編輯: celestialgod (123.205.27.107), 05/04/2015 18:51:43