warped_cmfb_analysis.m

function [ X ] = warped_cmfb_analysis(wcmfb_parms, x)
%warped_cmfb_analysis - ðåàëèçóþò íåðàâíîïîëîñíûé áàíê ôèëüòðîâ àíàëèçà.
%   wcmfb_parms  --  ñòðóêòóðà, ñîäåðæàùàÿ ïàðàìåòðû áàíêà ôèëüòðîâ
%   x            --  àíàëèçèðóåìûé ñèãíàë.

% 1. Ïàðàìåòðû
M = wcmfb_parms.NumberOfChannel;    % ×èñëî êàíàëîâ áàíêà ôèëüòðîâ
N = wcmfb_parms.OrderOfFP;          % Ïîðÿäîê ôèëüòðà-ïðîòîòèïà (ÔÏ)
m = wcmfb_parms.OrderOfPPComponents;% Ïîðÿäîê ïîëèôàçíûõ êîìïîíåíò ÔÏ
alpha = wcmfb_parms.Alpha;
h = wcmfb_parms.CoeffOfFP;          % Êîýôôèöèåíòû ôèëüòðà-ïðîòîòèïà
Npt = length(x);
X = zeros(M,Npt);                   % ìàòðèöà êàíàëüíûé ñèãíàëîâ
allpass_delay_chain = zeros(N,1);

% 2. Êîýôôèöèåíòû îáîáùåííîé ïîëèôàçíîé ñòðóêòóðû
h = reshape(h,2*M,length(h)/(2*M));
h(:,2:2:size(h,2))=-h(:,2:2:size(h,2));     %G(z)=G(-z)
h = h(:);

% 3. Ïîäãîòîâêà ìîäóëÿöèîííîé ìàòðèöû
c = zeros(M,2*M);
for k =0:M-1,
    for l=0:2*M-1,
        c(k+1,l+1) = 2*cos( (2*k+1)*(pi/(2*M))*(l - (N-1)/2) + (-1)^(k)*(pi/4) );
    end;
end;
    
% 4. Ðåàëèçàöèÿ áàíêà ôèëüòðîâ àíàëèçà
I_2M = eye(2*M);
polyphase_sum = sparse(kron(ones(1,m),I_2M));
for n=1:Npt
    % 4.1 Öåïî÷êà ôàçîâûõ çâåíüåâ
    R1 = x(n);
    R2 = allpass_delay_chain(1);
    allpass_delay_chain(1) = R1;
    for i=2:N
        R3 = R2;
        R2 = allpass_delay_chain(i);
        R1 = (R2-R1)*alpha + R3;
        allpass_delay_chain(i) = R1; 
    end
    
    % 4.2 Óìíîæåíèå íà êîýôôèöèåíòû ôèëüòðà-ïðîòîòèïà
    PP_out = allpass_delay_chain.*h;
    
    Xk = polyphase_sum*PP_out;

    % 4.3 Êîñèíóñíàÿ ìîäóëÿöèÿ 
    X(:,n) = (c*Xk)';
end
end