STFT
tfr= zeros (N,tcol) ;for icol=1:tcol,
ti= t(icol);
tau=-min():min();
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():min();和 indices= rem(N+tau,N)+1; 怎么解释？ 应该是信号边界的判断吧
你这里没有对信号进行边界延拓，所以信号的最开始部分有二分之窗宽，以及结束部分的最后二分之窗宽不能进行STFT
页:
[1]