не могу вычислить выражение

[Faton_11]

Добрый вечер!
как вычислить выражение
z=(1-a).*x+a.*y;
если a меняется от 0 до 1 с шагом 0.0001?
т.е. нужно получить первый столбец со значениями z при а=0, второй столбец при а=0.0001, третий - при а=0.0002 и т.д.
x и y -случайные последовательности.

Была попытка создать цикл, но как правильно организовать массив - непонятно.

Код: Выделить всё

MU=0;
SIGMA=1;
x=normrnd(MU,SIGMA,[10001 1]);
y=normrnd(MU,SIGMA,[10001 1]);
%a=(0:0.0001:1)';
%z=(1-a).*x+a.*y;
 
n=1;
for a=0:0.0001:1
z=(1-a).*x+a.*y;
h[z;n]; % как организовать эту строчку?
n=n+1;
end;

А может есть способ попроще? Если знаете, подскажите, пожалуйста

[abobroff]

Код: Выделить всё

MU = 0;
SIGMA = 1;
x = normrnd(MU,SIGMA,[10001 1]);
y = normrnd(MU,SIGMA,[10001 1]);
a = 0:0.0001:1;

z = bsxfun(@times,1 - a,x) + bsxfun(@times,a,y);

[Faton_11]

Мне долго пришлось с этим мучаться, а вы нашли такое оптимальное решение! Спасибо!

[sandy]

А зачем в этом случае bsxfun? достаточно простого матричного умножения:

Код: Выделить всё

z = x*(1 - a) + y*a;