用MATLAB写一段完整的QPSK调制解调的代码，要求解调中用到锁相环，需要画出信号的星座图以及眼图。

以下是一个简单的QPSK调制解调的MATLAB代码，包括锁相环和信号的星座图和眼图。

%% QPSK调制解调的MATLAB代码

clc;
clear all;
close all;

% 信号参数设置
N = 1000;                       % 信号长度
fc = 10e3;                      % 载波频率
fs = 100e3;                     % 采样率
Ts = 1/fs;                      % 采样时间
t = (0:N-1)*Ts;                 % 时间向量
M = 4;                          % 调制阶数
k = log2(M);                    % 每个符号的比特数
numBits = k*N;                  % 总比特数
dataIn = randi([0 1], numBits, 1);   % 生成随机二进制数据流
dataInMatrix = reshape(dataIn, k, length(dataIn)/k)';    % 将二进制数据流分成k位一组

% QPSK调制
dataMod = qammod(dataInMatrix, M);   % QPSK调制
txSig = reshape(dataMod, 1, length(dataMod));  % 转换为行向量

% 产生QPSK调制信号
carrier = cos(2*pi*fc*t) + 1i*sin(2*pi*fc*t);    % 产生载波信号
txSig = real(txSig.*carrier);                   % QPSK调制信号

% 信道噪声
EbNo = 10;                      % 信噪比
snrdB = EbNo + 10*log10(k) - 10*log10(fs);   % SNR
rxSig = awgn(txSig, snrdB, 'measured');   % 加入高斯白噪声

% QPSK解调
dataOut = qamdemod(rxSig, M);    % QPSK解调
dataOutMatrix = reshape(dataOut, length(dataOut)/k, k);    % 将解调后的数据分成k位一组
dataOut = dataOutMatrix(:);     % 将数据转换为列向量

% 锁相环
NCO = exp(-1i*2*pi*fc*t);       % 构造本振
loopGain = 0.1;                 % 锁相环增益
phaseError = angle(rxSig.*NCO); % 相位误差
filteredError = filter([loopGain], [1 -1], phaseError);   % 滤波
vco = exp(1i*cumsum(filteredError));    % 伪相锁定VCO
rxSig = rxSig./vco;             % 除去本振

% 信号星座图
scatterplot(dataMod);
title('QPSK调制信号星座图');

% 眼图
delay = 8;                      % 延迟
span = 4;                       % 窗口长度
eyediagram(rxSig(delay:end), span);
title('QPSK调制信号眼图');

注意，这里的锁相环只是一个简单的实现，实际应用中可能需要更复杂的算法来保证解调的精度和稳定性。