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
0052
0053
0054
0055
0056
0057
0058
0059 L(j) = round(fs/rrr(j));
0060
0061
0062 r(j) = fs/L(j);
0063 m(j) = f(j)/r(j);
0064
0065 tf = f(j) + r(j);
0066 if tf<fmax
0067 f(j+1) = tf;
0068 j = j + 1;
0069 else
0070 break
0071 end
0072 end
0073
0074
Generated on Fri 08-Jun-2007 16:09:11 by m2html © 2003