www.gusucode.com > 信道仿真的很好程序,可以实现不同参数情况下的信道仿真,比如超宽带系统的室内多径环境信道仿真 > 信道仿真/nakatest.m
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;