向量化连续幂矩阵的创建

[poster]

令x=1:100和N=1:10 。我想创建矩阵x^N使得i 列包含条目[1 ii^2 ... i^N]

我可以使用for循环轻松完成此操作。但是有办法使用矢量化代码吗？



回答：

我会去：

x = 1:100; N = 1:10; Solution = repmat(x,[length(N)+1 1]).^repmat(([0 N])',[1 length(x)]); 另一个解决方案（可能效率更高）：

Solution = [ones(size(x)); cumprod(repmat(x,[length(N) 1]),1)]; 甚至：

Solution = bsxfun(@power,x,[0 N]'); 希望这可以帮助。



更多&回答...