nakatest.m

来自「信道仿真的很好程序」· M 代码 · 共 41 行

function [q,m]= nakatest(data,bins,m)
%Gayatrii Prabhu and P. M. Shankar
y=sort(data); %received data
% m=mean(y.^2)^2/mean((y.^2-mean(y.^2)).^2); %Nakagami m parameter
omega=mean(y.^2); %Nakagami parameter

inter = max(y) - min(y); %interval for chi-square test 
intstep = inter/bins; %size of sub-interval
intval = [min(y):intstep:max(y)]; %samples at the sub-intervals
k = [zeros(1,bins)]; %number of samples in each sub-interval

%counting the samples in each sub-interval
for i = 1:length(y)
   for j = 1:bins
      if y(i) <= intval(j+1)
         k(j) = k(j) + 1;
         break;
      end;
   end;
end;

Fy=gammainc(((intval.^2).*(m/omega)),m); %Nakagami cdf
%%%%%%%%%%%%%%%add plot%%%%%%%%%%%%%%%%
plot(Fy);

p = zeros(1,bins);
for i = 1:bins
   p(i) = Fy(i+1) - Fy(i); %probabilities at the sub-interval samples
end;

np = p.*length(y);

q = 0; 
for i = 1:bins
   q = q + ((k(i)-np(i))^2)/np(i); %chi-square test statistic
end;