我知道要得到10个0,一个人可以做到。
A = zeros(10, 1);要得到10个1,一个人可以做到。
A = ones(10, 1);任意的数字呢?比方说,我想要10 3。我想出了一种方法。
A = linspace(3, 3, 10);这样满意吗?有什么更优雅的方法吗?
发布于 2014-11-03 12:18:09
一些替代办法:
repmat:
A= repmat(3,5 7);填充值3的%// 5x7矩阵下面是所有建议的方法之间的时序比较。正如您所看到的,@Dennis' solutions 是最快的(至少在我的系统上)。
功能:
function y = f1(x,m,n) %// Dennis' "ones"
y = x*ones(m,n);
function y = f2(x,m,n) %// Dennis' "zeros"
y = x+zeros(m,n);
function y = f3(x,m,n) %// Luis' "repmat"
y = repmat(x,[m n]);
function y = f4(x,m,n) %// Luis' "dirty"
y(m,n) = 0;
y(:) = x;
function y = f5(x,m,n) %// Luis' "indexing"
y(1:m,1:n) = x;
function y = f6(x,m,n) %// Divakar's "matrix porod"
y = x*ones(m,1)*ones(1,n);正方形数组的基准脚本:
x = 3;
sizes = round(logspace(1,3.7,10)); %// max limited by computer memory
for s = 1:numel(sizes)
n = sizes(s);
m = sizes(s);
time_f1(s) = timeit(@() f1(x,m,n));
time_f2(s) = timeit(@() f2(x,m,n));
time_f3(s) = timeit(@() f3(x,m,n));
time_f4(s) = timeit(@() f4(x,m,n));
time_f5(s) = timeit(@() f5(x,m,n));
time_f6(s) = timeit(@() f6(x,m,n));
end
loglog(sizes, time_f1, 'r.-');
hold on
loglog(sizes, time_f2, 'g.:');
loglog(sizes, time_f3, 'b.-');
loglog(sizes, time_f4, 'm.-');
loglog(sizes, time_f5, 'c.:');
loglog(sizes, time_f6, 'k.:');
xlabel('Array size')
ylabel('Time')
legend('ones', 'zeros', 'repmat', 'dirty', 'indexing', 'matrix prod')对于列数组:只需更改以下行:
sizes = round(logspace(1,3.7,10)).^2; %// max limited by computer memory
n = 1;
m = sizes(s);对于行数组:
sizes = round(logspace(1,3.7,10)).^2; %// max limited by computer memory
n = sizes(s);
m = 1;结果为双核CPU,2-GB内存,Windows,Matlab R2010b:



发布于 2014-11-03 12:17:14
有两种基本方法可以做到这一点:
A = ones(10,1)*3
B = zeros(10,1)+3第一种是最常用的例子,但如果我没有弄错,第二种则略好一些。总之,这只是品味的问题。
当然,如果您有一个现有的矩阵,还有一种简单的方法:
C = zeros(10,1)
C(:) = 3;为了完整起见,@Luis建议的爬行动物解决方案也很好。
发布于 2014-11-03 13:52:05
作为另一种选择,基于matrix multiplication的方法(在MATLAB上被认为是相当快的)也可以用于2D或multidimensional array的分配工作。
因此,假设m为rows,n为columns,x为所有元素的赋值,代码为-
y = x*ones(m,1)*ones(1,n);https://stackoverflow.com/questions/26714067
复制相似问题