Matlab-MSK调制与解调 MSK信号产生的信息、I路、Q路信号对比图、接收端解调的信息恢复的对比图、I路、Q路载波信号与滤波后的信号对比图

本博客环境为Matlab2018 ，软件版本不同可能会有些出入，需要稍作修改。

MSK调制与解调
仿真要求：
画出MSK信号产生的信息、I路、Q路信号对比图以及接收端解调的信息恢复的对比图。
I路、Q路载波信号与滤波后的信号对比图，直接对比滤波前后的波形图。

代码如下：

clear all;clc;
N=20; % 输入的比特
T=1; % 比特周期
fc=2; % 载波频率
Fs=100; % 采样频率
bitstream=randi([0,1],1,N); % 随机产生的比特流
bitstream1= 2*bitstream-1; % 0 to -1; 1 to 1 单极性变成双极性

% 差分编码
b0=1;
for i=1:N
    encode_output(i)= b0*bitstream1(i);%差分编码的输出，对应原理图的bk
    b0=encode_output(i);
end

I=[];Q=[];
for i=1:N %奇数进I路，偶数进Q路
    if mod(i,2)~=0
        I=[I,encode_output(i)];
    else
        Q=[Q,encode_output(i)];
    end
end
% using plot to compare bitsteam, I, Q
bit_data=[];
for i=1:N
    bit_data=[bit_data,encode_output(i)*ones(1,T*Fs)];
end

I_data=[];Q_data=[];base_wave=-T:1/Fs:T-1/Fs;
for i=1:N/2
    I_data=[I_data,I(i)*cos(pi*base_wave/(2*T))];%分别与cos（πt/2T）相乘
    Q_data=[Q_data,Q(i)*cos(pi*base_wave/(2*T))];
end
% Q路的延时一个比特周期
number_delay=length(base_wave)/2; 
Q_data1=[zeros(1,number_delay),Q_data(1:length(Q_data)-number_delay)];
% plot
figure();
t=0:1/Fs:N*T-1/Fs;
subplot(3,1,1)
plot(t,bit_data);legend('Bitstream')
subplot(3,1,2)
plot(t,I_data);legend('I Bitstream')
subplot(3,1,3)
plot(t,Q_data1);legend('Q Bitstream')

% carrier signal
bit_t=0:1/Fs:N*T-1/Fs;%定义时间戳
I_carrier=cos(2*pi*fc*bit_t);%定义载波信号2*pi*fc*t
Q_carrier=cos(2*pi*fc*bit_t+pi/2);%定义载波信号2*pi*fc*t+pi/2
% transmit signal
MSK_signal=I_data.*I_carrier+Q_data1.*Q_carrier;%发送的MSK信号

snr=1; % 信噪比
% Receive singal
MSK_receive= awgn(MSK_signal,snr);%加入噪声
% demodulate
I_output=MSK_receive.*I_carrier;
Hd=myfilter;%调用低通滤波器前的系数
I_filter_ouput=filter(Hd,I_output);%滤波后的信号
Q_output= MSK_receive.*Q_carrier;
Q_filter_ouput= filter(Hd,Q_output);
for i=1:N/2 %抽样判别
    if I_filter_ouput((2*i-1)*number_delay)>0 %I路的奇数倍时；例如i=2时，第三个周期进行采样
        I_recover(i)= 1;
    else
        I_recover(i)=-1;
    end
    if Q_filter_ouput(2*i*number_delay)>0 %Q路是偶数倍周期时采样
        Q_recover(i)= 1;
    else
        Q_recover(i)=-1;
    end
end
bit_recover=[];
for i=1:N
    if mod(i,2)~=0
        bit_recover=[bit_recover,I_recover((i-1)/2+1)];
    else
        bit_recover=[bit_recover,Q_recover(i/2)];
    end
end

figure();
t=0:1/Fs:N*T-1/Fs;
subplot(4,1,1)
plot(t,I_carrier);legend('I carrier signal')
subplot(4,1,2)
plot(t,Q_carrier);legend('Q carrier signal')
subplot(4,1,3)
plot(t,I_filter_ouput);legend('I filter')
subplot(4,1,4)
plot(t,Q_filter_ouput);legend('Q filter')


% 差分译码
for i=1:N 
    if i==1
        bit_recover1(i)=bit_recover(i);
    else
        bit_recover1(i)= bit_recover(i)*bit_recover(i-1);
    end
end
% using plot to compare bitsteam, I, Q
recover_data=[];
for i=1:N 
    recover_data=[recover_data,bit_recover1(i)*ones(1,T*Fs)];
end

% plot
figure();
t=0:1/Fs:N*T-1/Fs;
bit_stream=[];
for i=1:N 
    bit_stream=[bit_stream,bitstream1(i)*ones(1,T*Fs)];
end
subplot(2,1,2)
plot(t,bit_stream);legend('original Bitstream')
subplot(2,1,1)
plot(t,recover_data);legend('recover Bitstream')

低通滤波器：
在命令行输入filterDesigner。

修改为你所需要的参数。Fpass我改为1，改为2或3也可以的。
修改完后点击最下端的Design Filter按钮来设计这个滤波器。

依次点击 File - Generate MATLAB Code - Filter Design Function 将其保存为函数文件。我将其文件命名为myfilter.m

在MATLAB里也可以打开这个文件。

运行结果如下：


载波信号与滤波后的信号对比图：