phi_nk.m

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

function y = Phi_nk(alpha, beta, s, xi_k) 
%%%%%%%  phi_nk(s)  

a=alpha(:);
b=beta(:);
I=ones(size(a));
[n1 , m1] =size (a) ;
[n2 , m2] =size(b) ;
if (m1~=1) | (m2~=1),
    error('x and nu must be vectors') ;
end ;
if(n1~=n2)
    error('x and nu must be of the same length');
end ;
if(s==0) 
    s = eps;
end
phi_nk_s = sum( (a*s)./(I-b*s) - log( I-b*s)  ) - log( s ) - s*xi_k;
y = phi_nk_s;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%