qpsk频谱图matlab仿真

QPSK (Quadrature Phase Shift Keying，四相位调制)是一种数字调制方式，其频谱图可以通过Matlab进行仿真。

首先，生成一个长度为N的随机比特序列，然后将其转换为QPSK信号。QPSK调制将每两个比特映射到一个复平面上的点，具体映射方式如下：

00 -> 1+j

01 -> -1+j

10 -> -1-j

11 -> 1-j

可以使用以下代码实现比特序列到QPSK信号的转换：

N = 1000; % 比特序列长度
bits = randi([0 1],1,N); % 生成随机比特序列
symbols = zeros(1,N/2); % 初始化符号序列
for i = 1:2:N
    if bits(i) == 0 && bits(i+1) == 0
        symbols((i+1)/2) = 1 + 1j;
    elseif bits(i) == 0 && bits(i+1) == 1
        symbols((i+1)/2) = -1 + 1j;
    elseif bits(i) == 1 && bits(i+1) == 0
        symbols((i+1)/2) = -1 - 1j;
    else
        symbols((i+1)/2) = 1 - 1j;
    end
end

接下来，对符号序列进行QPSK调制，得到QPSK信号：

Fs = 10000; % 采样率
fc = 1000; % 载波频率
t = 0:1/Fs:(N/2-1)/Fs; % 时间轴
I = real(symbols); % 实部作为I路信号
Q = imag(symbols); % 虚部作为Q路信号
s = sqrt(2)*I.*cos(2*pi*fc*t) - sqrt(2)*Q.*sin(2*pi*fc*t); % QPSK信号

最后，使用Matlab自带的fft函数计算QPSK信号的频谱，并画出频谱图：

L = length(s); % 信号长度
f = (-Fs/2:Fs/L:Fs/2-Fs/L); % 频率轴
S = fftshift(fft(s)); % 计算频谱
figure;
plot(f,abs(S)/L);
xlabel('频率(Hz)');
ylabel('幅值');
title('QPSK信号频谱图');

这样就可以得到QPSK信号的频谱图了。