ltpda_compute_f
PURPOSE 
LTPDA_COMPUTE_F compute the frequency and resolution vectors for
SYNOPSIS 
function [f,r,m, L,rr,rrr] = ltpda_compute_f(fs, N, Kdes, Kmin, Jdes, olap)
DESCRIPTION 
CROSS-REFERENCE INFORMATION 
This function calls:
This function is called by:
- ltpda_lpsd LTPDA_LPSD implement LPSD algorithm for analysis objects.
- ltpda_ltfe LTPDA_LTFE implement weighted transfer-function estimated computed on a log frequency axis.
SOURCE CODE 
0001 function [f,r,m, L,rr,rrr] = ltpda_compute_f(fs, N, Kdes, Kmin, Jdes, olap)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024 fmin = fs/N;
0025 fmax = fs/2;
0026 j = 1;
0027 f(j) = fmin;
0028
0029
0030 rav = fs/N*(1+(1-olap)*(Kdes-1));
0031 rmin = fs/N*(1+(1-olap)*(Kmin-1));
0032 g = log(fmax)-log(fmin);
0033
0034
0035 while f < fmax
0036
0037 rr(j) = f(j)*g/(Jdes-1);
0038
0039 if rr(j) >= rav
0040 rrr(j) = rr(j);
0041 else
0042 rrr(j) = sqrt(rav*rr(j));
0043 end
0044
0045
0046 if rrr(j) <= rmin
0047 rrr(j) = rmin;
0048 end
0049
0050
0051 L(j) = round(fs/rrr(j));
0052
0053
0054 r(j) = fs/L(j);
0055 m(j) = f(j)/r(j);
0056
0057 tf = f(j) + r(j);
0058 if tf<fmax
0059 f(j+1) = tf;
0060 j = j + 1;
0061 else
0062 break
0063 end
0064 end
0065
0066
Generated on Wed 04-Jul-2007 19:03:10 by m2html © 2003