Bezdepbonusu7Программирование: статьи и ответы
Меню

Вопрос

Размещение значений (чисел) в многослойной матрице ячеек в MATLAB

Автор вопроса: ImanОтветов: 1

Предположим, что Q - матрица, которая имеет 3 ячейки, и в каждой ячейке есть 2 другие ячейки, что означает:

Q={ { [] [] } ; { [] [] }; { [] [] } }

Более того, если у нас есть "a" и "b", в каждом из которых по 3 члена, и мы хотели бы разместить

"a(1,1)" в "Q{1}{1}",

"b(1,1)" в "Q{1}{2}",

"a(2,1)" в "Q{2}{1}",

"b(2,1)" в "Q{2}{2}",

"a(3,1)" в "Q{3}{1}",

"b(3,1)" в "Q{3}{2}",

Например, если

a = [2; 3; 4];
b = [1; 5; 8]

Тогда Q должен быть как

Q={{2 1};
   {3 5};
   {4 8}}

Обратите внимание, что нам нужен векторизованный код, а не код цикла for, так как последний у меня уже есть, как показано ниже:

for i=1:size(Q,2)

     Q{i}{1}=a(i,:)
     Q{i}{2}=b(i,:)

end

Спасибо.

  • А в чем проблема с вашей петлей? Dan 16.07.2014
  • Я очистил ваш Q (сохранив его структуру идентичной тому, как он был у вас), и это противоречит вашему циклу, в котором есть ячейки в ячейках. Есть ли причина иметь внутренние матрицы в виде ячеек? то есть могут ли эти скаляры когда-либо не быть скалярами? А также почему бы просто не Q = [a,b] забыть о ячейках? Dan 16.07.2014
  • Как насчет этого - Q = mat2cell([a b].',2,ones(1,numel(a))).'? Divakar 16.07.2014
  • @Dan Есть причина поместить каждого члена в матрицу подъячейки в каждую суперячейку, поскольку я только что упростил матрицы Q, a, b. На самом деле Q и другие матрицы в моем коде содержат миллиарды чисел, и каждая подъячейка представляет собой другую характеристику. Более того, я хотел бы написать его в строке, а не в цикле, так как это может быть очень дорого и занимать много времени для таких огромных данных. Кстати, спасибо за внимание. Iman 16.07.2014
  • @Divakar Большое спасибо. Это выглядит очень хорошо. Мне было интересно, не могли бы вы сообщить мне, как я могу разделить каждый элемент a и b на отдельные ячейки. Например: Q={{2}{1};{3}{5};{4}{8}} Iman 16.07.2014
  • @Iman Может быть, это - mat2cell([a b].',[1 1],ones(1,numel(a))).' Divakar 16.07.2014

matlabvectorizationcellcell-arraymulti-layer

Ответы

Ответ 1 · рейтинг 2 · Divakar

Код

Q = mat2cell(num2cell([a b]),ones(1,numel(a)),2)

Пример

Код с отображением ввода и вывода

a = [2; 3; 4]; %// Inputs
b = [1; 5; 8];

Q = mat2cell(num2cell([a b]),ones(1,numel(a)),2); %// Output

celldisp(Q) %// Display results

Вывод при запуске кода

Q{1}{1} =
     2
Q{1}{2} =
     1
Q{2}{1} =
     3
Q{2}{2} =
     5
Q{3}{1} =
     4
Q{3}{2} =
     8

Сравнительный анализ

Функция для метода цикла

function out = loop1(a,b)

out = cell(size(a,1),1);
for i=1:size(out,1)
    out{i}{1}=a(i,:);
    out{i}{2}=b(i,:);
end

return;

Функция векторного метода

function out = vec1(a,b)

out = mat2cell(num2cell([a b]),ones(1,numel(a)),2);

return;

Код сравнительного анализа

N_arr = [50 100 200 500 1000 2000 5000 10000 50000]; %// array elements for N

timeall = zeros(2,numel(N_arr));
for k1 = 1:numel(N_arr)

    N = N_arr(k1);
    a = randi(9,N,1);
    b = randi(9,N,1);

    f = @() loop1(a,b);
    timeall(1,k1) = timeit(f);
    clear f

    f = @() vec1(a,b);
    timeall(2,k1) = timeit(f);
    clear f
end

%// Graphical display of benchmark results
figure,
hold on
plot(N_arr,timeall(1,:),'-ro')
plot(N_arr,timeall(2,:),'-kx')
legend('Loop Method','Vectorized Method')
xlabel('Datasize (N) ->'),ylabel('Time(sec) ->')

Результаты

Выводы

Похоже, векторизованный метод - это то, что нужно, поскольку он показывает почти удвоенную производительность (с точки зрения времени выполнения) по сравнению с циклическим подходом в широком диапазоне размеров данных.

Перевод вопроса и ответов сообщества Stack Overflow, лицензия CC BY-SA. Оригинал на Stack Overflow.

Похожие материалы