dapsktogray_frequency.m

通信系统中的最新的分块差分调制技术的源码

function [out]=dapsktogray_frequency(data,a,ncon)
lendata=length(data);
m=log2(ncon);
temp=1;
 for k=1:lendata
    theta=mod(angle(data(k)),2*pi);
    amplitude=abs(data(k));
    if ncon>16
        [bitp]=psk_16_gray(theta);
    else
        [bitp]=psk_8_gray(theta);
    end    
    
if amplitude>=(1+1/a)/2&amplitude<=(a+1)/2
    bita=0;
else
    bita=1;
end
out(temp:temp+m-2)=bitp;
out(temp+m-1)=bita;
temp=k*m+1;
end