Code covered by the BSD License  

Highlights from
MATLAB for Digital Communication

image thumbnail

MATLAB for Digital Communication

by

 

The MATLAB programs in "MATLAB/Simulink for Digital Communication" authored by Won Y. Yang et. al

[J,JJ]=Jkb(K,beta)
function [J,JJ]=Jkb(K,beta)
% The 1st kind of kth-order Bessel function
tmpk= ones(size(beta));
for k=0:K
   tmp= tmpk;  JJ(k+1,:)= tmp;
   for m=1:100
      tmp=-tmp.*beta.*beta/4/m/(m+k);
      JJ(k+1,:)= JJ(k+1,:)+tmp; % Eq.(E1.17.4)
      if norm(tmp)<0.001, break; end
   end
   tmpk=tmpk.*beta/2/(k+1);  
end 
J=JJ(K+1,:);

Contact us