mxa.m

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

 
 function [A]=MXA(M,T0,Tb,detaPHI,Be,Bo,G,theta)
 A=zeros(2*M);
 for I=1:(2*M)
     f1=(I-M-1)/T0;
     for J=1:(2*M)
         f2=(J-M-1)/T0;
         H1f1=(exp(-j*2*pi*f1*Tb)+exp(-j*theta)*exp(j*detaPHI))/2;%H1(i)
         H1f2=(exp(-j*2*pi*f2*Tb)+exp(-j*theta)*exp(j*detaPHI))/2;%H1(j)
         Hof1=exp(-(f1/(0.65*Bo)).^4);%Ho(i)
         H2f1=(exp(-j*2*pi*f1*Tb)-exp(-j*theta)*exp(j*detaPHI))/2;%H2(i)
         H2f2=(exp(-j*2*pi*f2*Tb)-exp(-j*theta)*exp(j*detaPHI))/2;%H2(j)
         Hof2=exp(-(f2/(0.65*Bo)).^4);%Ho(j)
         fe=(J-I)/T0;
         ff=2.43*fe/Be;
         Hef=945/(j*ff.^5+15*ff.^4-105*j*ff.^3-420*ff.^2+945*j*ff+945);
         A(I,J)=(H1f2*conj(H1f1)-H2f2*conj(H2f1))*(G^2)*Hef*Hof2*conj(Hof1);
     end
 end