intefun.m

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

function y = InteFun(vtr_alpha, vtr_beta, u0, kai, v, ksi_k) 
%%%%%%%%  \phi_nk(s) 的一阶导数

a= vtr_alpha(:);
b=vtr_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
s = u0 + 0.5*kai*v*v + i*v;
%s = u0 + i*v;
phi_nk_s = sum( (a*s)./(I-b*s) - log( I-b*s )  ) - log(s) - s*ksi_k;
r = real( exp( phi_nk_s ) *( 1-i*kai*v ) );
if  isinf(r) | isnan(r)
    r=0;
end
y =r;

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