oddhkz_cem.m

来自「基于复指数调制多相滤波器的设计」· M 代码 · 共 38 行

function Hk_z = oddHkz_CEM(Theta,M_ch,m)

% M_ch: 
% m:
% Theta:

Num_lattice = floor(M_ch/2);                                           % ����������
Theta = Theta(1:m*Num_lattice);
Theta = reshape(Theta,Num_lattice,[]);

[iRow,iCol] = size(Theta);
if((iRow ~= Num_lattice) || (iCol ~= m))
    disp('Error Theta');
    return;
end

syms w z;
Theta = sym(Theta);
pi = 4*atan(1);
P_puredelay = [pi/4 pi/2*ones(1,m-1)];    % ����ʱ��Ԫ�IJ���Ϊ��ʼ���IJ���
P_puredelay = sym(P_puredelay);

% ��1������Theta�õ�G_k(z)��G_Mk(z)
G_kpz = [cos(Theta(:,1)) zeros(Num_lattice,m-1)];
G_Mkpz = [sin(Theta(:,1)) zeros(Num_lattice,m-1)];
G_oddk = [cos(P_puredelay(1)) zeros(1,m-1)];
G_oddMk = [sin(P_puredelay(1)) zeros(1,m-1)];

j = sqrt(-1);
z = exp(2*M_ch*j*w);

for k=1:Num_lattice
    for p=1:(m-1)
        G_kpz(k,p+1) = G_kpz(k,p)*cos(Theta(k,p+1)) + z^(-1)*G_Mkpz(k,p)*sin(Theta(k,p+1));
        G_Mkpz(k,p+1) = G_kpz(k,p)*sin(Theta(k,p+1)) - z^(-1)*G_Mkpz(k,p)*cos(Theta(k,p+1));
    end
end
% ���Ϲ���Num_lattice����