MSK的matlab仿真

关于MSK的matlab仿真，包括调制，使用非相干解调进行解调！！！ 关于MSK的matlab仿真，包括调制，使用非相干解调进行解调！！！ 关于MSK的matlab仿真，包括调制，使用非相干解调进行解调！！！

close all;
clear all;
fs = 50000;%采样率  
Time_Hold_On = 0.001;   %Rb = 1000
Num_Unit = fs * Time_Hold_On;  %T中采样点数
fc =3000;                                             %初始化载波频率 f1 = 3250 ；f2 = 2750
f1 = fc + 1/Time_Hold_On/4;
f2 = fc - 1/Time_Hold_On/4;
Sign_Set = [-1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1];
%Sign_Set = [-1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1 -1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1 -1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1 -1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1 -1 -1 1 1 -1 1 -1 -1 -1 -1 -1 -1 1 1 -1 -1];
Lenth_Of_Sign_Set = length （ Sign_Set ）; %计算信号长度 

Sign_SetI = zeros （ 1  Lenth_Of_Sign_Set ）;
Sign_SetQ = zeros （ 1  Lenth_Of_Sign_Set ）;
Sign_SetI1 = zeros （ 1 Num_Unit * Lenth_Of_Sign_Set ）;
Sign_SetQ1 = zeros （ 1 Num_Unit * Lenth_Of_Sign_Set ）;
Sign_SetI2 = zeros （ 1 Num_Unit * Lenth_Of_Sign_Set ）;
Sign_SetQ2 = zeros （ 1 Num_Unit * Lenth_Of_Sign_Set ）;
Sign_Setdiff =  ones （ 1  Lenth_Of_Sign_Set ）;
j=1;

%差分编码
Sign_Setdiff（1） = -1;
for （i = 2:Lenth_Of_Sign_Set）
    Sign_Setdiff（i） =（Sign_Set（i） * Sign_Set（i-1））;
end
for I=1:2:Lenth_Of_Sign_Set                            %信号分离成两路奇数偶数信号
    Sign_SetI（j）= Sign_Set（I+1）;Sign_SetI（j+1）= Sign_Set（I+1）;
    Sign_SetQ（j）=Sign_Set（I）;Sign_SetQ（j+1）= Sign_Set（I）;
    j=j+2;
end

for I = 1 : Lenth_Of_Sign_Set
    Sign_Set1（ （（I-1）*Num_Unit + 1） : I*Num_Unit）  = Sign_Set（I）;%生成m1（t）
    Sign_Setdiff1（ （（I-1）*Num_Unit + 1） : I*Num_Unit）  = Sign_Setdiff（I）;%生成mq（t）
end 

for I = 2 : Lenth_Of_Sign_Set
    Sign_SetI1（ （（I-1）*Num_Unit + 1） : I*Num_Unit）  = Sign_SetI（I-1）;%生成m1（t）
    Sign_SetQ1（ （（I-1）*Num_Unit + 1） : I*Num_Unit）  = Sign_SetQ（I）;%生成mq（t）
end 

t = 0 : 1/fs : Time_Hold_On * Lenth_Of_Sign_Set- 1/fs; %Rb * 16 = 0.01*16 = 0.16
Sign_SetI2 = Sign_SetI1 .* cos （pi * t / 2 / Time_Hold_On）;
Sign_SetI2 = Sign_SetI2 .* cos （2 * pi * fc * t ）;
Sign_SetQ2 = Sign_SetQ1 .* sin （pi * t