sl.m

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

function [s1l,s2l]=sl(a,N,Tb,L,Eb,Be,Bo,b2,detaPHI,Long,G);
x=zeros(1,2*L+1);
s1l=zeros(1,2*L+1);
s2l=zeros(1,2*L+1);
for l=-L:1:L
    f=l/(N*Tb);
    pl=sqrt(Eb/Tb)*Tb*sinc(f*Tb);
    Hf=1;%exp(-j*2*pi^2*b2*Long*f.^2);
    Hof=1;%exp(-(f/(0.65*Bo)).^4);
    H1f=(exp(-j*2*pi*f*Tb)+exp(j*detaPHI))/2;
    H2f=(exp(-j*2*pi*f*Tb)-exp(j*detaPHI))/2;
    sum=0;
    for k=0:N-1
        sum=sum+a(k+1)*exp(-j*2*pi*l*k/N);
    end
    x(l+L+1)=1/(N*Tb)*pl*sum;
    s1l(l+L+1)=x(l+L+1)*Hf*G*Hof*H1f;
    s2l(l+L+1)=x(l+L+1)*Hf*G*Hof*H2f;
end