用matlab做最小频移键控调制（msk）

最小频移键控调制 (MSK) 是一种连续相位调制技术，通常用于数字通信系统中。在 MSK 中，相邻两个符号的相位差为 $\pm 90^\circ$，也就是说，MSK 信号的频谱中只有两个相邻峰值，它们之间的频率差为符号速率的一半。

下面是用 MATLAB 实现 MSK 调制的代码：

% 定义参数
fs = 10000;         % 采样率
fc = 500;           % 载波频率
T = 1/fs;           % 采样周期
Ts = 1/10;          % 符号周期
N = Ts/T;           % 每个符号采样点数
fc1 = fc - 1/(4*Ts); % 频率偏移量
fc2 = fc + 1/(4*Ts); % 频率偏移量

% 生成随机比特流
data = randi([0 1], 1, 100);

% 生成 MSK 信号
t = 0:T:Ts*N*length(data)-T;
msk = zeros(1, length(t));
for i = 1:length(data)
    if data(i) == 0
        msk((i-1)*N+1:i*N) = cos(2*pi*fc1*t((i-1)*N+1:i*N));
    else
        msk((i-1)*N+1:i*N) = cos(2*pi*fc2*t((i-1)*N+1:i*N));
    end
end

% 绘制 MSK 信号时域图和频谱图
figure(1)
subplot(2,1,1)
plot(t, msk)
xlabel('Time (s)')
ylabel('Amplitude')
title('MSK Signal in Time Domain')
subplot(2,1,2)
f = linspace(-fs/2, fs/2, length(msk));
msk_fft = fftshift(abs(fft(msk)));
plot(f, msk_fft)
xlabel('Frequency (Hz)')
ylabel('Amplitude')
title('MSK Signal in Frequency Domain')

上述代码中，首先定义了一些参数，包括采样率、载波频率、符号周期等。然后生成了一个随机的比特流，用于产生 MSK 信号。在生成 MSK 信号时，根据比特流的值选择不同的频率偏移量，从而得到 MSK 信号。最后，绘制了 MSK 信号在时域和频域的图像。

需要注意的是，上述代码中只是实现了 MSK 调制，如果需要解调，还需要进行相应的解调处理。