Scallop 損失式を使用して Rectangular、Hamming、および Blackman ウィンドウの Scallop 損失の値を計算する関数を作成したいと考えています。関数を作成しましたが、その答えが 0 しか返されません。エラーでしょうか?
function s_l = scallop loss(len)
window = rectwin(len);
num_total = 0;
den_total = 0;
for n = 0:len-1
A1 = exp(-((1i*(n)*pi)/len));
A2 = window(n+1)*A1;
num = abs(A2);
den = win(n+1);
num_total = num_total + num;
den_total = den_total + den:
end
result = 20*log(num_total/den_total);
s_l = result;