function deta=freq_pingheng(A,q,N)
%A is the frequences array,q is the fanwei of frequence ,N is the lenth of
%frequences
aa=A;
%計算跳頻序列平衡性參數
length_aa=length(aa);
bb=zeros(1,q);%儲存每個頻率出現次數的數組
%統計每個頻率出現的次數
cc=[0:(q-1)];
for i=1:q
    for j=1:length_aa
        if cc(i)==aa(j)
           aa(j)=-1;
            bb(i)= bb(i)+1;
        end        
    end
end
bb;         %輸出序列中每個頻率的出現次數
%計算參數...............
dd=zeros(q,0);%過渡變量
for k=1:q
    dd(k)=(bb(k)-N/q)^2;
end
dd;
ff=1/q*sum(dd);
deta=q/N*sqrt(ff);%平衡參數