我正在研究 Matlab 中的一种算法,该算法需要定期更新矩阵的某些元素,并寻找最好的方法。这是我要实现的目标的描述:
- 我有一个 MxN 数组
A
和 1xN 向量B
。 - 基本上,向量
B
是一个逻辑索引,它描述了我需要选择哪一列,A
即C = A(:,B)
。 - 不幸的是,逻辑向量
B
因某些进程而异。这意味着其中的列数C
不是固定的。 - 其他一些处理将用作输入并生成另一个与ie大小相同的
C
数组D
C
size(D) == size(C)
- 现在,我需要“重塑”
D
以使其具有与A
. 棘手的部分是上面#2中没有选择的那些列A
需要替换为NaN
s。当然,我可以使用循环的粗略方式来做到这一点。但我希望以 Matlab 的方式执行此操作,即线性或逻辑索引、矢量化等。这是我目前陷入困境的地方。
一些让事情更清楚的例子:
让我们说
A = [1 2 3; 4 5 6; 7 8 9]
B = [1 0 1]
C = A(:,B) = [1 3; 4 6; 7 9]
经过一些处理,我会得到D = [2 5; 6 7; 3 3]
. 现在,我需要通过ie填充将其“重塑”D
为相同的大小。A
NaN
D = [2 NaN 5; 6 NaN 7; 3 NaN 3]
到目前为止我尝试过的,
Atmp = NaN(size(A));
Btmp = find(repmat(B,[size(B,1),1]));
Atmp(Btmp) = D(Btmp); %-> error because D is smaller than A.