min(abs(fp-ft));
0125             f(n)=fp(k);
0126             <span class="comment">% Determine nominal smoothing window length</span>
0127             Nwt=round(min(Lsmin*exp(s*(n-1)),2*(fs/2-f(n))/df));
0128             <span class="keyword">if</span> mod(Nwt,2) == 0
0129                 Nwt=Nwt+1;
0130             <span class="keyword">end</span>
0131             Nw(n)=Nwt;
0132             win=feval(winstr,Nw(n));
0133             win = win/sum(win);
0134             <span class="keyword">if</span> k+(