用matlab2018b写一个msk调制示例代码，调制信号进行频谱图绘制

好的，我可以为您提供一个MSK调制的matlab示例代码，并绘制出调制信号的频谱图。下面是代码：

% MSK调制
fs = 10000; % 采样率
dt = 1/fs; % 采样时间间隔
t = 0:dt:1; % 时间序列
fc = 500; % 载波频率
kf = 2000; % 调频灵敏度
fdev = kf*0.5; % 频偏
m = 0.5; % 调制指数

% 生成调制信号
bitrate = 100; % 比特率
data = randi([0 1],1,length(t)*bitrate); % 生成随机数据
data_p = 2*data - 1; % 极性编码
phi = zeros(1,length(t)); % 相位序列初始化
for i = 1:length(data)
    phi(i+1) = phi(i) + 2*pi*bitrate*(1+m*data_p(i))*dt; % 相位连续变化
end
phi = phi(1:end-1);
s_t = sqrt(2/T)*cos(2*pi*fc*t + phi); % MSK调制信号
s_f = fftshift(fft(s_t)); % 频谱图

% 绘制调制信号的频谱图
df = fs/length(s_t);
f = -fs/2:df:fs/2-df;
figure(1)
plot(f,abs(s_f),'LineWidth',1.5)
xlabel('频率(Hz)')
ylabel('|S_t(f)|')
title('MSK调制信号的频谱图')
grid on

希望这可以帮到您。如果您有更多问题，请随时问我。