2DPSK调制解调simulink仿真及matlab程序

对DPSK调制解调过程进行仿真，分别使用了simulink和matlab程序

fs=60000;                             %采样频率
Rb=1200;                              %符号率
Time_Hold_On=1/Rb;                    %保持时间
Num_Unit=fs*Time_Hold_On;             %数据点数
High_Level=ones（1Num_Unit）;          %门限1
Low_Level=zeros（1Num_Unit）;          %门限0
f=1800;                               %载波频率
A=1;                                  %信号幅度
N=8;                                  %消息长度

%>>>>>>>>>>>码型变换>>>>>>>>>>>
sign_absolute=[1011010010];%产生绝对码
sign_ar=zeros（111）;
for I=1:10
    sign_ar（I+1）=xor（sign_ar（I）sign_absolute（I））;
end
sign_relative=sign_ar（[2:11]）;%相对码的产生
delay_relative=sign_ar（[1:10]）;%延时后的相对码的产生

%>>>>>>>>>>>initial_the_signal>>>>>>>>>>>>>>>>>>>>
Lenth_Of_Sign=length（sign_relative）;
st=zeros（1Num_Unit*Lenth_Of_Sign）;%调制后的信号
signal_absolute=zeros（1Num_Unit*Lenth_Of_Sign）;%基带信号（绝对码波形）
signal_relative=zeros（1Num_Unit*Lenth_Of_Sign）;%相对码波形
signal_delay_relative=zeros（1Num_Unit*Lenth_Of_Sign）;%延时相对码波形
signal_result=zeros（1Num_Unit*Lenth_Of_Sign）;%解调后的信号
t=0:1/fs:Lenth_Of_Sign*Time_Hold_On-1/fs;

%>>>>>>>>>>>>>>>基带信号（绝对码波形）>>>>>>>>>>>>>>>>>>>>
for I=1:Lenth_Of_Sign
    if sign_absolute（I）==1
        signal_absolute（（I-1）*Num_Unit+1:I*Num_Unit）=High_Level;
    else
        signal_absolute（（I-1）*Num_Unit+1:I*Num_Unit）=Low_Level;
    end
end


%>>>>>>>>>>>>>>>相对码波形>>>>>>>>>>>>>>>>>>>>
for I=1:Lenth_Of_Sign
    if sign_relative（I）==1
        signal_relative（（I-1）*Num_Unit+1:I*Num_Unit）=High_Level;
    else
        signal_relative（（I-1）*Num_Unit+1:I*Num_Unit）=Low_Level;
    end
end



%>>>>>>>>>>>>>>>>信号的调制>>>>>>>>>>>>>>>>>>>>>>
for I=1:Lenth_Of_Sign
    if sign_relative（I）==1
        st（（I-1）*Num_Unit+1:I*Num_Unit）=A*sin（2*pi*f*t（（I-1）*Num_Unit+1:...
            I*Num_Unit）+pi）;
    else
        st（（I-1）*Num_Unit+1:I*Num_Unit）=A*sin（2*pi*f*t（（I-1）*Num_Unit+1:...
            I*Num_Unit））;
    end
end
%>>>>>>>>>>>>>>>>延时后调制波形>>>>>>>>>>>>>>>>>>>>>>
for I=1:Lenth_Of_Sign
    if delay_relative（I）==1
        signal_delay_relative（（I-1）*Num_Unit+1:I*Num_Unit）=A*sin（2*pi*f*t...
            （（I-1）*Num_Unit+1:I*Num_Unit）+pi）;
            
    else
        signal_delay_relative（（I-1）*Num_Unit+1:I*Num_Unit）=A*sin（2*pi*f*t...
            （（I-1）*Num_Unit+1:I*Num_Unit））;
    end
end
%>>>>>>>>>>>>>>>>调制后信号通过高斯信道>>>>>>>>>>>>>>>>>>>>>>
st_awgn=awgn（st20）;
signal_delay_relative_awgn=awgn（signal_delay_relative20）;

%>>>>>>>>>>>>>相乘>>>>>>>>>>>>>>>>>>>>
xc=st_awgn.*signal_delay_relative_awgn;

%>>>>>>>>>>>>低通滤波部分>>>>>>>>>>>>>>
[NWn]=buttord（2*pi*10002*pi*5000325‘s‘）;%临界频率采用角频率表示
[ba]=butter（NWn‘s‘）;
[bzaz]=impinvar（bafs）;
dt=filter（bzazxc）;

%>>>>>>>>>>>>>>>>>>>>抽样判决>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
for I=1:Lenth_Of_Sign
    if dt（（Num_Unit/2）*（2*I-1））<0
        signal_result（（I-1）*Num_Unit+1:I*Num_Unit）=High_Level;
    else
        signal_result（（I-1）*Num_Unit+1:I*Num_Unit）=Low_Level;
    end
end


figure（1）
plot（tsignal_absolute‘b‘）;
a

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4140  2019-01-09 23:54  第十四题DPSK调制与解调电路仿真电路 高硕 李高歌 闫月宁\2dpsk.m

     文件      25988  2019-01-09 23:43  第十四题DPSK调制与解调电路仿真电路 高硕 李高歌 闫月宁\DPSK.slx

     文件    1789445  2019-01-14 09:47  第十四题DPSK调制与解调电路仿真电路 高硕 李高歌 闫月宁\第十四题DPSK调制与解调电路仿真电路.pdf

     目录          0  2019-01-14 09:47  第十四题DPSK调制与解调电路仿真电路 高硕 李高歌 闫月宁

----------- ---------  ---------- -----  ----

              1819573                    4