MSK调制的功率谱仿真MATLAB代码

MSK调制是一种常用的数字调制方式，其功率谱可以通过MATLAB进行仿真。下面是一个简单的MSK调制功率谱仿真MATLAB代码示例：

% MSK调制功率谱仿真MATLAB代码
% 定义参数
T = 1; % 符号周期
fs = 1000; % 采样率
f1 = 10; % 载波频率
fdev = 5; % 频率偏移量
N = T*fs; % 符号周期内的采样点数

% 产生随机二进制序列
data = round(rand(1,N));

% 生成MSK调制信号
s = zeros(1,N);
for n=1:N
    if data(n) == 1
        s(n) = cos(2*pi*f1*(n-1)/fs + pi/2);
    else
        s(n) = cos(2*pi*f1*(n-1)/fs - pi/2);
    end
end

% 计算功率谱
[P,f] = pwelch(s,[],[],[],fs,'twosided');
P = 10*log10(P/max(P)); % 将功率谱归一化为0-1，并转换为dB

% 绘制功率谱图
plot(f,P);
xlabel('Frequency (Hz)');
ylabel('Power (dB)');
title('MSK Modulation Spectrum');

代码中，首先定义了MSK调制的相关参数，包括符号周期、采样率、载波频率、频率偏移量等。然后生成随机二进制序列，并根据序列生成MSK调制信号。最后通过pwelch函数计算功率谱，并绘制功率谱图。注意，绘制的功率谱图是以dB为单位的，同时也进行了归一化处理，以便更好地观察功率谱的变化。