1

学校の演習では、要素ごとの関数を使用して 10x10 の乗算表を作成し、できるだけ短くする必要があります。これは私が書いたコードです (機能していますが長すぎます)。このコードにいくつかの工夫を提案してください。前もって感謝します (:

base=zeros(10);
oneten=[1:1:10];
base(1,:)=1.*oneten;
base(2,:)=2.*oneten;
base(3,:)=3.*oneten;
base(4,:)=4.*oneten;
base(5,:)=5.*oneten;
base(6,:)=6.*oneten;
base(7,:)=7.*oneten;
base(8,:)=8.*oneten;
base(9,:)=9.*oneten;
base(10,:)=10.*oneten
4

5 に答える 5

4

次のようにします。

(1:10)' * (1:10)

編集 --->

N が大きい場合、Daniel、私、Luis Mendo、および David によって提案されたソリューションの速度をテストしました。

N = 100;    % number of iterations
runtime_a = zeros(N, 1);    % runtime of Daniel's solution
runtime_b = zeros(N, 1);    % runtime of the obvious solution
runtime_c = zeros(N, 1);    % runtime of Luis Mendo's solution
runtime_d = zeros(N, 1);    % runtime of Luis Mendo's solution
runtime_e = zeros(N, 1);    % runtime of David's solution
n = 5000;    % number of elements
one_to_n = 1:n;
for hh = 1:N
  % Solution by Daniel R.
  tic, a = bsxfun(@times, one_to_n, one_to_n'); runtime_a(hh) = toc;
  clear a
  tic, b = one_to_n' * one_to_n; runtime_b(hh) = toc;
  clear b
  % Solution by Luis Mendo
  tic, c = cell2mat(arrayfun(@(x) (x:x:n*x).', one_to_n, 'uni', false)); runtime_c(hh) = toc;
  clear c
  % Solution by Luis Mendo.
  tic, d = cumsum(repmat(one_to_n, [n 1])); runtime_d(hh) = toc;
  clear d
  % Solution by David
  tic, [A, B] = meshgrid(one_to_n); e = A.*B; runtime_e(hh) = toc;
  clear e
end

% Check mean and standard deviation:
mean([runtime_a, runtime_b, runtime_c, runtime_d, runtime_e])
std([runtime_a, runtime_b, runtime_c, runtime_d, runtime_e])

結果は次のとおりです。

% mean:
0.105048900691251   0.188570704057917   0.491929288458701   0.787045026437718   0.979624197407329
% standard deviation:
0.034274873626281   0.077388368324427   0.163983982925543   0.285395301735065   0.372693172505310

したがって、明らかに、N が大きい場合、ダニエルの解が最速です。

于 2013-11-03T20:33:55.830 に答える
2

楽しみのために: 他の可能性は

  • cumsum(repmat(1:10,[10 1]))
    
  • cell2mat(arrayfun(@(n) (n:n:10*n).',1:10,'uni',false))
    
于 2013-11-03T21:21:17.610 に答える
1
oneten=[(1:10)]
base = bsxfun(@times,oneten,oneten')

この場合、事前割り当て ( base=zeros(10);) は不要です。

そして、理解しやすい別のソリューション:

base=zeros(10);
oneten=[(1:10)];
for k=oneten
    base(k,:)=k.*oneten;
end
于 2013-11-03T20:38:25.990 に答える
0
[A,B]=meshgrid(1:10);
A.*B

要素ごとの乗算を使用

于 2013-11-03T21:41:31.470 に答える