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

### 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; ```