Учитывая M
и k
, вот как вы можете выполнить суммирование:
A = rand(M+1,1); %# Create M+1 random values
Y = sin(2*pi*(k-(0:M)))*A; %# Use a matrix multiply to perform the summation
РЕДАКТИРОВАТЬ: Вы даже можете создать функцию для Y
, которая принимает k
и A
в качестве аргументов:
Y = @(k,A) sin(2*pi*(k+1-(1:numel(A))))*A; %# An anonymous function
result = Y(k,A); %# Call the function Y