STFT

花瓣雨

tfr= zeros (N,tcol) ;  
for icol=1:tcol,
    ti= t(icol);
    tau=-min([round(N/2)-1,Lh,ti-1]):min([round(N/2)-1,Lh,xrow-ti]);
    indices= rem(N+tau,N)+1;
    tfr(indices,icol)=sig(ti+tau,1).*conj(h(Lh+1+tau));
end;
这个是STFT的被分析信号和窗口函数的乘积过程的实现过程。
其中N为数据点数，tcol为时间点数；tfr被分析信号和窗口函数共扼的乘积；
我想问这一步tau=-min([round(N/2)-1,Lh,ti-1]):min([round(N/2)-1,Lh,xrow-ti]);和 indices= rem(N+tau,N)+1; 怎么解释？

tjucruiser

应该是信号边界的判断吧
你这里没有对信号进行边界延拓，所以信号的最开始部分有二分之窗宽，以及结束部分的最后二分之窗宽不能进行STFT