dpsk的matlab实现

%%E7_1_DPSKModem.m程序清单
%仿真DPSK信号的产生及相干解调过程
%将成形滤波器系数量化为10比特后，写入tra_lpf.coe文件中
%将接收滤波器系数量化为10比特后，写入rec_lpf.coe文件中
%将DPSK已调数据量化为8比特后，写入Dpsk.txt文件中
ps=1*10^6;   %码速率为1MHz
a=0.8;       %成形滤波器系数为0.8
B=（1+a）*ps;  %中频信号处理带宽
Fs=8*10^6;  %采样速率为8MHz
fc=2*10^6;   %载波频率为2MHz
N=20000;      %仿真数据的长度

t=0:1/Fs:（N*Fs/ps-1）/Fs;%产生长度为N频率为fs的时间序列
s=randint（N12）;       %产生随机数据作为原始数据
%并将绝对码变换为相对码
ds=ones（1N）;
for i=2:N
    if s（i）==1 
        ds（i）=-ds（i-1）;
    else
        ds（i）=ds（i-1）;
    end
end
%对相对码数据以Fs频率采样
Ads=upsample（dsFs/ps）;

%设计平方升余弦滤波器
n_T=[-2 2];
rate=Fs/ps;
T=1;
Shape_b = rcosfir（an_TrateT）;%figure（4）;freqz（Shape_b）
%对采样后的数据进行升余弦滤波;
rcos_Ads=filter（Shape_b1Ads）;

%产生载频信号
f0=sin（2*pi*fc*t）;       
%产生DPSK已调信号
dpsk=rcos_Ads.*f0;       

%与相干载波相乘