查看单个帖子
旧 2019-12-10, 20:48   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
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]'); 希望这可以帮助。



更多&回答...
poster 当前离线   回复时引用此帖