vect_v.m

来自「一种新的光纤通信调制技术」· M 代码 · 共 19 行

function [v]=vect_v(s1l,s2l,T0,tk,N,Tb,M,L,Be,Bo,detaPHI,G);
v=zeros(2*M+1,1);
for I=1:2*M+1   
    sum=0;
    for l=-L:L
        f=(I-M-1)/T0;
        H1f=(exp(-j*2*pi*f*Tb)+exp(j*detaPHI))/2;
        H2f=(exp(-j*2*pi*f*Tb)-exp(j*detaPHI))/2;
        Hof=1;%exp(-(f/(0.65*Bo)).^4);
        H1f=H1f*Hof*G;
        H2f=H2f*Hof*G;
        fe=l/(N*Tb)-f;
        ff=2.43*fe/Be;
        Hef=945/(j*ff.^5+15*ff.^4-105*j*ff.^3-420*ff.^2+945*j*ff+945);
        sum=sum+(s1l(l+L+1)*conj(H1f)-s2l(l+L+1)*conj(H2f))*Hef*exp(j*2*pi*l*tk/(N*Tb));
    end
    v(I)=sum;
end