nk_mean_var.m

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

E_nk=zeros(1,N);
var_nk=zeros(1,N);
for k=1:N
    tk=detat+k*Tb/4;%%%%
    v=vect_v(s1l,s2l,T0,tk,N,Tb,M,L,Be,Bo,detaPHI,G);
    b=zeros(2*M+1,1);
    b=U'*v;
    alpha=zeros(2*M+1,1);
    beta=zeros(2*M+1,1);
    for I=1:(2*M+1)
        alpha(I)=((abs(b(I)))^2)/value(I);
        beta(I)=2*value(I)*sgma^2;
    end
    
    for I=1:(2*M+1);
        E_nk(k)=E_nk(k)+alpha(I)+beta(I);
        var_nk(k)=var_nk(k)+beta(I)*(2*alpha(I)+beta(I));
    end
end