Вычитание каждого элемента из другого элемента в столбце В Matlab?

Question

Я пытаюсь вычесть каждый элемент из другого элемента в векторе в matlab. Я знаю, что diff (X) дает разность элементов. Пусть

x = {1 
    3 
    3 
    6 
    8 
    10 
    15 
    23} 

DiffX = diff (x); 

приведет

DiffX = {2 
     0 
     3 
     2 
     2 
     5 
     8} 

Однако то, что я ищу это различие каждого элемента из каждого элемента что-то вроде этого

DiffX = {0 
     2 
     2 
     5 
     7 
     9 
     14 
     22 
     -2 
     0 
     0 
     3 
     5 
     7 
     12 
     20 
     . 
     . 
     . 
     . 
     . 
     0} 

Который приходит от

DiffX = {x(1) - x(1), x(2) - x(1), x(3) - (1), x(4) - x(1), ... , x(8) - x(1), x(1) - x(2), x(2) - x(2), x(3) - x(2), ......., x(8) - x(8)} 

Я не знаю, как это сделать, чтобы получить ответ. Помощь очень ценится.

Большое спасибо

Answer 1

Это может быть сделано с bsxfun(@plus, x,-x'). Это добавляет поэтапную двоичную операцию + к массивам x и -x '(отрицательный x транспонирован).

---

В (было) указано в комментарии, то же самое может быть достигнуто за счет bsxfun(@minus, x, x').

---

Если результат должен быть вектор, вы должны изменить его на reshape(bsxfun(@minus, x', x),[],1). См. Также комментарий Роберта Р ниже.

Кстати, синтаксис в вашем вопросе немного запутан. Как было сказано ранее, { и } подразумевают, что результатом является массив ячеек, а применение diff указывает на нормальный массив. Кроме того, обозначение для присвоения значений x предполагает нормальный векторный результат, но отсутствие ; подразумевает, что результатом будет вектор-строка. Во всяком случае, играйте с символами транспонирования, и должен появиться желаемый результат.

Answer 2

bsxfun это лучший способ пойти, но это, как вы это сделать с помощью петель (не делать это таким образом, это только для иллюстрации концепции петель):

я предлагаю, что вы храните результат в виде матрицы 2D, и в этом случае решение петля проста:

x = {1 3 3 6 8 10 15 23} 

D = zeros(numel(x)); %// This preallocation step is essential for efficiency in Matlab 
for A = 1:numel(x) 
    for B = 1:numel(x) 
     D(A,B) = x{A} - x{B}; 
    end 
end 

Answer 3

Просто интересно, может быть, вы можете достичь что-то подобное с нормальными вектор разностных операций?

Начиная от чего-то вроде этого:

x = [1 3 3 6 8 10 15 23] 

X = repmat(x,size(x,2),1) 

делает разницу X и транспонированные вы должны получить разницу между всеми комбинациями элементов в x:

X - X' 

На этом этапе вы должны выбрать соответствующим элементам, которые вас интересуют (верхняя треугольная матрица?) и сглаживаются в вектор.