2DPSK的MATLAB代码实现

clear all;
close all;
fs=30000;
Time_Hold_On=0.1;
Num_Unit=fs*Time_Hold_On;
High_Level=ones（1Num_Unit）;
Low_Level=zeros（1Num_Unit）;
w=300;
A=1;
m=10;%随机二进制数个数
SNR=100;%信噪比


%[NWn]=buttord（2*pi*502*pi*150325‘s‘）;
%[ba]=butter（NWn‘s‘）;
%[bzaz]=impinvar（bafs）;
%dt=filter（bzazdt）;
%升余弦滤波器
beta=0.4;
span=6;
sps=40;

b=rcosdesign（betaspansps‘sqrt‘）;
%fvtool（b‘Analysis‘‘impulse‘）;
%z = upfirdn（dt b 40 sps）;%以40上采样，以sps下采样


%Sign_Set=[011001101110010101];
Sign_Set=randint（1m2）;
Lenth_Of_Sign=length（Sign_Set）;
c=zeros（1Lenth_Of_Sign）;
st=zeros（1Num_Unit*Lenth_Of_Sign）;
sign_orign=zeros（1Num_Unit*Lenth_Of_Sign）;
sign_orign_relative=zeros（1Num_Unit*Lenth_Of_Sign）;
sign_result=zeros（1Num_Unit*Lenth_Of_Sign）;
sign_result_2=zeros（1Num_Unit*Lenth_Of_Sign）;
sign_result_3=zeros（1Num_Unit*Lenth_Of_Sign）;
t=0:1/fs:Time_Hold_On*Lenth_Of_Sign-1/fs;
x=0:1/fs:（Time_Hold_On*Lenth_Of_Sign-1/fs）/Lenth_Of_Sign;

%将绝对码转换为相对