vect_ek.m

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

canshu;
ck=zeros(1,N);
dk=zeros(1,N);
fbconnection=[1 0 0 1 0];
a=m_sequence(fbconnection);
format long g;
[s1l,s2l]=sl(a,N,Tb,L,Eb,Be,Bo,b2,detaPHI,Long,G);
A=MXA(M,T0,Tb,detaPHI,Be,Bo,G);
[V,D]=eig(A);
[U,value]=MXU(V,D,A,M);
for k=1:N
    tk=detat+k*Tb/2;%%%%
    ck(k)=vect_c(s1l,s2l,L,N,Tb,Be,tk);
    v=vect_v(s1l,s2l,T0,tk,N,Tb,M,L,Be,Bo,detaPHI,G);
    b=zeros(2*M+1,1);
    b=U'*v;
    sum=0;
    for I=1:2*M+1
        sum=sum+(abs(b(I))^2)/value(I);
    end
    dk(k)=ck(k)-sum;
end
ek=r_th-dk;