写一个matlab程序:MSK能够实现其数字基带波形，时域波形图及其功率谱图

以下是一个简单的 MATLAB 代码，可以生成 MSK 信号的数字基带波形、时域波形图和功率谱图。该代码使用 MATLAB 的 Communications Toolbox 中的函数。

% 设置参数
f_carrier = 1000; % 载波频率
f_symbol = 100; % 符号速率
T_symbol = 1/f_symbol; % 符号时间
T_sample = T_symbol/4; % 采样时间间隔
T = 10*T_symbol; % 信号持续时间
N = T/T_sample; % 采样点数
alpha = 0.5; % MSK 调制指数

% 生成数字基带信号
data = randi([0,1],N,1);
data(data==0) = -1;
symbols = data(1:2:end) + 1i*data(2:2:end);
baseband_signal = mskmod(symbols,f_symbol,alpha);

% 绘制数字基带波形
figure;
plot(real(baseband_signal),imag(baseband_signal),'o');
xlabel('I'); ylabel('Q');
title('MSK数字基带波形');
axis square;

% 生成载波信号
t = (0:N-1)*T_sample;
carrier = cos(2*pi*f_carrier*t);

% 生成 MSK 信号
msk_signal = real(baseband_signal.*exp(1i*2*pi*f_carrier*t));

% 绘制时域波形
figure;
plot(t,msk_signal);
xlabel('时间(s)'); ylabel('幅度');
title('MSK时域波形');

% 计算功率谱密度
[P,f] = periodogram(msk_signal,[],[],1/T_sample);
PdB = 10*log10(P);

% 绘制功率谱图
figure;
plot(f,PdB);
xlabel('频率(Hz)'); ylabel('功率谱密度(dB/Hz)');
title('MSK功率谱密度图');

代码中使用了 `mskmod` 函数生成 MSK 调制信号的数字基带波形，然后乘以一个载波信号得到 MSK 信号。最后，使用 `periodogram` 函数计算功率谱密度，并用 `plot` 函数绘制时域波形和功率谱图。