?function?y=QPSK（signalFbFcFsCOD）

signal0=signal;
k0=length（signal0）;
signal=[signal0randint（1k0*（1/COD-1））];
Fr=Fs/Fb;
k=length（signal）;
Source1=signal（1:2:k）;
S1=kron（Source1ones（1Fr））;
Source2=signal（2:2:k）;
S2=kron（Source2ones（1Fr））;

L=length（S1）;
fai=zeros（1L）;

for?i=1:L
????if?（S1（i）==0）&&（S2（i）==0）
????????fai（i）=pi/4;
????elseif?（S1（i）==1）&&（S2（i）==0）
????????fai（i）=3*pi/4;
????elseif?（S1（i）==1）&&（S2（i）==1）
????????fai（i）=5*pi/4;
????elseif?（S1（i）==0）&&（S2（i）==1）
????????fai