[291105]: comp / comp_gabmixdual_fac.m  Maximize  Restore  History

Download this file

53 lines (41 with data), 974 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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
function gammaf=comp_gabmixdual_fac(gf1,gf2,L,a,M)
%COMP_GABMIXDUAL_FAC Computes factorization of mix-dual.
% Usage: gammaf=comp_gabmixdual_fac(gf1,gf2,a,M)
%
% Input parameters:
% gf1 : Factorization of first window
% gf2 : Factorization of second window
% L : Length of window.
% a : Length of time shift.
% M : Number of channels.
%
% Output parameters:
% gammaf : Factorization of mix-dual
%
% GAMMAF is a factorization of a dual window of gf1
%
% This function does not verify input parameters, call
% GABMIXDUAL instead
%
% See also: gabmixdual, comp_fac, compute_ifac
% AUTHOR : Peter L. Søndergaard.
% TESTING: OK
% REFERENCE: OK
LR=prod(size(gf1));
R=LR/L;
b=L/M;
N=L/a;
c=gcd(a,M);
d=gcd(b,N);
p=b/d;
q=N/d;
gammaf=zeros(p*q*R,c*d);
G1=zeros(p,q*R);
G2=zeros(p,q*R);
for ii=1:c*d
G1(:)=gf1(:,ii);
G2(:)=gf2(:,ii);
S=G2*G1';
Gpinv=M*S\G2;
gammaf(:,ii)=Gpinv(:);
end;

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:





No, thanks