gusucode.com > 信道仿真的很好程序,可以实现不同参数情况下的信道仿真,比如超宽带系统的室内多径环境信道仿真 > 信道仿真/normtest.m
function q=normtest(data,bins) %Gayatrii Prabhu and P. M. Shankar y=sort(data); y=(y-mean(y))./std(y); %mean 0, var 1 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=normcdf(intval,0,1); %cdf at each sub-interval sample 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;