-
Notifications
You must be signed in to change notification settings - Fork 2
/
warped_cmfb_synthesis.m
57 lines (47 loc) · 1.52 KB
/
warped_cmfb_synthesis.m
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
51
52
53
54
55
56
function [ x ] = warped_cmfb_synthesis(wcmfb_params, X)
%warped_cmfb_synthesis_by_sum -- Âûïîëíÿåò ðåêîíñòðóêöèþ ñèãíàëà ïî åãî
% ñóáïîëîñíûì êîìïîíåíòàì ïóòåì ñóììèðîâàíèÿ.
% 1. Ïàðàìåòðû
eq_filt = wcmfb_params.PhaseCorr; % Ôèëüòð êîððåêöèè Ô×Õ
h = wcmfb_params.CoeffOfFP; % Êîýôôèöèåíòû ôèëüòðà-ïðîòîòèïà
M = wcmfb_params.NumberOfChannel; % ×èñëî êàíàëîâ áàíêà ôèëüòðîâ
N = wcmfb_params.OrderOfFP; % Ïîðÿäîê ôèëüòðà-ïðîòîòèïà (ÔÏ)
m = wcmfb_params.OrderOfPPComponents;% Ïîðÿäîê ïîëèôàçíûõ êîìïîíåíò ÔÏ
alpha = wcmfb_params.Alpha;
Sig_length = size(X,2);
APC_delay_1 = zeros(N,1);
APC_delay_2 = zeros(N,1);
x = zeros(Sig_length,1);
c = zeros(M,2*M);
% 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. Ïîäãîòîâêà ìîäóëÿöèîííîé ìàòðèöû
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. Ðåàëèçàöèÿ áàíêà ôèëüòðîâ ñèíòåçà
for n=1:Sig_length
% 4.1 Êîñèíóñíàÿ äåìîäóëÿöèÿ
Curr_X = c'*X(:,n);
% 4.2 Ïîëèôàçíàÿ ôèëüòðàöèÿ
Xk = repmat(Curr_X, m, 1);
Xk = Xk.*h;
% 4.3 Öåïî÷êà ôàçîâûõ çâåíüåâ
R1 = Xk(1);
for i=1:N-1
R3 = APC_delay_1(i);
APC_delay_1(i) = R1;
R2 = APC_delay_2(i);
R1 = (R2-R1)*alpha + R3;
APC_delay_2(i) = R1;
R1 = R1 + Xk(i+1);
end
x(n) = R1;
end
% 5. Âûðàâíèâàíèå Ô×Õ
x = filter(eq_filt,1,x);
end