基于matlab的仿真MSK信号的代码

以下是基于MATLAB的MSK信号仿真的代码示例：

% 设置参数
f0 = 10^6; % 载波频率
Tb = 1/f0; % 每个比特的时间
fs = 10*f0; % 采样率
Ts = 1/fs; % 采样时间间隔
N = 100; % 信号长度为100个比特
A = 1; % 调制深度为1

% 生成随机比特序列
bits = randi([0 1], [1 N]);

% 产生MSK信号
t = (0:N-1)*Tb;
s = zeros(1, length(t));
for i=2:N
    if bits(i)==1
        s((i-1)*Tb/Ts+1:i*Tb/Ts) = A*cos(2*pi*f0*t((i-1)*Tb/Ts+1:i*Tb/Ts)-pi/2);
    else
        s((i-1)*Tb/Ts+1:i*Tb/Ts) = A*cos(2*pi*f0*t((i-1)*Tb/Ts+1:i*Tb/Ts)+pi/2);
    end
end

% 绘制MSK信号时域波形
figure;
plot(t, s);
xlabel('Time (s)');
ylabel('Amplitude');
title('MSK Signal in Time Domain');

% 绘制MSK信号频域波形
f = (-fs/2:fs/N:fs/2-fs/N);
S = fftshift(fft(s));
figure;
plot(f, abs(S));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('MSK Signal in Frequency Domain');

这段代码中，我们首先设置了一些参数，比如载波频率、每个比特的时间、采样率等等。然后，我们生成了一个随机比特序列，长度为100个比特。接着，根据MSK调制的原理，我们根据比特序列产生了MSK信号。最后，我们分别绘制了MSK信号的时域波形和频域波形。

需要注意的是，这段代码中产生的MSK信号是理想的，没有加入任何噪声，所以只能用于学习和演示用途。在实际应用中，需要考虑到噪声的影响，可以使用AWGN信道模型来模拟信道噪声。