这个程序如何理解

miki321

clc
clear

h(1:128)=cos((0:127)*2*pi*6/256);
h(31:71)=h(31:71)+0.5*cos((-10:30)*2*pi*52/256);
h(129:256)=cos((0:127)*2*pi*25/256);
N = 256;
H = [fft(h) fft(h) fft(h)];
al = [0:(N/2)-1, -(N/2):-1];
for n=1:N
W=exp(-2*pi^2*al.^2./((n-(N/2)-1).^2));
S(n,:) = ifft(H(n+(N/2):n+(N/2)+N-1) .* W);
end

S(N/2+1,:) = ifft(H(N+1:2*N) .* [1 zeros(1,N-1)]);
S_TRANS = abs(S);
figure;
pcolor(S_TRANS)

程序是时间序列的S变换，其公式为
S(t,f)=∫_(-∞)^∞▒〖H(a+f)exp⁡((-2π^2 a^2)/f^2 )  exp⁡(+2πiat)da〗
上面程序中
H = [fft(h) fft(h) fft(h)]; 这句应该怎么理解，为什么要延拓信号的傅里叶变换？
al = [0:(N/2)-1, -(N/2):-1];这句为什么要设置成这种形式？

miki321

真的木有来帮忙看看啊？

ChaChing

个人水平有限,感觉这不是编程问题!?
或许应移至信号版块!

firecat_2

确实是属于信号处理问题，去信号版问问应该能有结果

happy

楼主应该先说明这个程序是干嘛的？