Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.

Close

[e2c226]: comp / comp_dgtreal_long.m Maximize Restore History

Download this file

comp_dgtreal_long.m    39 lines (31 with data), 982 Bytes

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
function cout=comp_dgtreal_long(f,g,a,M)
%COMP_DGTREAL_LONG Full-window factorization of a Gabor matrix.
% Usage: c=comp_dgtreal_long(f,g,a,M);
%
% Input parameters:
% f : Factored input data
% g : Window
% a : Length of time shift.
% M : Number of channels.
% Output parameters:
% c : M x N*W*R array of coefficients, where N=L/a
%
% Do not call this function directly, use DGT instead.
% This function does not check input parameters!
%
% References: so07-2 st98-8
% AUTHOR : Peter L. Søndergaard.
% TESTING : TEST_DGT
% REFERENCE : OK
L=size(f,1);
W=size(f,2);
N=L/a;
M2=floor(M/2)+1;
gf=comp_wfac(g,a,M);
% Compute the window application
% We know the output is real, but comp_dgt_walnut cannot detect this, so
% we force the output to be real.
cout=real(comp_dgt_walnut(f,gf,a,M));
% FFT with only positive frequencies
cout=fftreal(cout)/sqrt(M);
cout=reshape(cout,M2,N,W);