创建一个列数不断变化的矩阵

Question

我正在尝试制作一个自制版本的 peakfinder.m，使其与多个阵列一起工作，而不是一次只使用一个，以提高时间效率。 (http://www.mathworks.com/matlabcentral/fileexchange/25500-peakfinder)

我有一个 2D 矩阵，我需要在其中查找符号在 2 维中是否发生变化。

dx0 = diff(x0,1,2); % Find derivative
dx0(dx0 == 0) = -eps; % This is so we find the first of repeated values
ind = find(dx0(:,1:end-1).*(dx0(:,2:end)) < 0)+1; % Find where the derivative changes sign

现在我的问题是它确实找到了导数改变符号的位置，但它是一个大向量。因此，如果标志在同一行中发生两次变化（或没有连续变化），我无从得知。

因此，如果 x0 的大小为 1000x10，我希望 ind 的大小为 1000xY，其中 Y 是它在每行中更改符号的次数。我还需要知道 x0 的哪些值有符号变化。所以每一行的样式都是：

2 4 7
4 8
2 5 6 8

等等

这可能吗？或者我应该更改代码，如果它没有变化，它会放置一个 0，如果它发生变化，它会放置一个 1，考虑到我将使用它发生变化的值？

Answer 1

cellfun 方法-

b1 = padarray(sign(dx0(:,1:end-1))~=sign(dx0(:,2:end)),[0 1],'pre')
out = cellfun(@find,mat2cell(b1,ones(1,size(b1,1)),size(b1,2)),'uni',0)

以上代码假设您拥有padarray，这似乎是最近添加到 MATLAB 的图像处理工具箱中的。所以，如果你没有它，你可以像这样连接zeros（与false）-

b1 = sign(dx0(:,1:end-1))~=sign(dx0(:,2:end))
b1 = [false(size(b1,1),1) b1]
out = cellfun(@find,mat2cell(b1,ones(1,size(b1,1)),size(b1,2)),'uni',0)

使用cellfun with nonzeros function 的替代解决方案-

b1 = padarray(sign(dx0(:,1:end-1))~=sign(dx0(:,2:end)),[0 1],'pre')
out = cellfun(@nonzeros,mat2cell(bsxfun(@times,b1,1:size(b1,2)),ones(1,size(b1,1)),size(b1,2)),'uni',0)

out 包含sign 跨行变化的位置，可以使用celldisp(out) 显示。

sign 变化中的counts 可以使用 -

来计算

counts = cellfun(@numel,out)