Matlab优雅地添加行和列

Question

假设我们有以下随机矩阵：

1 2 3 4
5 6 7 8
9 8 7 6
5 4 3 2

我想将其转换为以下内容：

1 0 2 0 3 0 4 0
0 0 0 0 0 0 0 0
5 0 6 0 7 0 8 0
0 0 0 0 0 0 0 0
9 0 8 0 7 0 6 0
0 0 0 0 0 0 0 0
5 0 4 0 3 0 2 0
0 0 0 0 0 0 0 0

出于某种原因，我不能使用mathjax格式，所以它看起来有点可怕，对不起。点，是，我想在当前行和列之间添加零行和列，以便我将其大小增加2倍。

我提出了以下代码，但它只适用于非常小的矩阵，如果我在一个大图像上使用它，由于内存限制问题它无法完成。

clear all

I=imread('image.png');
I=rgb2gray(I);

B=zeros(2*size(I));

[x, y]=find(-inf<I<inf);

xy=[x,y];
nxy=xy;

%coord change
nxy=2*xy-1;

B(nxy(:,1),nxy(:,2))=I(xy(:,1),xy(:,2)); 

我预计会很快，因为它是用maltlab函数完全矢量化的，但它很难失败。有没有其他优雅的方法来做到这一点？

Answer 1

如果您查看索引向量，对于2x2矩阵，这类似于I([1 1 2 2] ,[1 2 1 2] );，这意味着您将每行和每列索引两次。正确的解决方案是B(1:2:end,1:2:end)=I;，它为每隔一行和每隔一列编制索引。

Answer 2

这也可以通过一个班轮来完成，比如原来的矩阵称为A，然后

kron(A,[1,0;0,0])