对比gfsk和fsk信号的功率谱的matlab代码实现

要对比GFSK和FSK信号的功率谱，需要首先生成相应的信号。然后使用Matlab提供的函数进行功率谱的计算和绘图。

下面是一个简单的示例代码：

% 生成GFSK信号
fs = 1000;              % 采样频率
t = 0:(1/fs):1;         % 时间序列
fc = 100;               % 载波频率
dev = 20;               % 频偏
bt = 1;                 % 位时产品
alpha = 0.5;            % 频率调制指数
modulationIndex = (1 + alpha*dev/bt) * pi;
gfsk_signal = cos(2*pi*fc*t + cumsum(modulationIndex*sin(2*pi*dev*t)));

% 生成FSK信号
fsk_signal_1 = cos(2*pi*100*t);                        % 频率为100Hz的载波
fsk_signal_2 = cos(2*pi*(100+dev)*t);                  % 频率为100+dev Hz的载波
fsk_signal = (fsk_signal_1 + fsk_signal_2) / 2;

% 计算功率谱
window = hamming(length(gfsk_signal));         % 使用汉明窗口
nfft = 1024;                                  % 傅里叶变换点数
[P_gfsk,f_gfsk] = periodogram(gfsk_signal,window,nfft,fs);
[P_fsk,f_fsk] = periodogram(fsk_signal,window,nfft,fs);

% 绘制功率谱
figure;
subplot(2,1,1);
plot(f_gfsk, 10*log10(P_gfsk));
title('GFSK功率谱');
xlabel('频率(Hz)');
ylabel('功率谱(dB)');
subplot(2,1,2);
plot(f_fsk, 10*log10(P_fsk));
title('FSK功率谱');
xlabel('频率(Hz)');
ylabel('功率谱(dB)');

以上代码中，首先生成了GFSK和FSK信号，然后使用periodogram函数计算功率谱，并使用plot函数进行绘制。最后使用subplot函数将两个功率谱图像分别绘制在同一幅图上。

注意：生成的信号参数和窗口长度等可根据实际需求进行调整。