октаву функция векторизовать приложения

Question

Для заданного (м по п) матрица poly Я хотел бы генерировать свои строки со следующим выражением:

poly(i, :) = [X(i) X(i)^2 X(i)^3 ... X(i)^p]

Я дан (м на 1) X вектор и значение p. Мое текущее решение:

poly(:,1) = X; 
for i = 2:p 
    poly(:,i) = X.^i; 
end; 

Мой вопрос: есть ли способ их для дальнейшего векторизации это? Я также сгенерировал массив ячеек функций, которые можно было применить к матрице подряд за строкой, но мне все еще пришлось петли.

ТИА

Answer 1

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

Шаг через раствор:

poly = ones(size(X,1),p ); 
poly = poly .* X; 
powers = 1:p; 
poly = poly .^ powers; 

Однострочник:

poly = (ones(size(X,1),p) .* X) .^[1:p];