对QPSK调制格式的光信号加入色散，并进行色散补偿，输出补偿前后的星座图、眼图、时域波形、频域波形的matlab代码

以下是一个简单的QPSK调制信号加入色散并进行色散补偿的示例matlab代码，同时输出补偿前后的星座图、眼图、时域波形和频域波形。

% QPSK调制信号加入色散并进行色散补偿的示例matlab代码

% 设置QPSK调制参数
M = 4;                          % 调制阶数
k = log2(M);                    % 每个符号的比特数
n = 30000;                      % 生成的符号数
data = randi([0 1],n*k,1);      % 生成数据比特流
data_matrix = reshape(data,k,length(data)/k)';
data_symbols = bi2de(data_matrix);   % 将比特流转换成符号

% 进行QPSK调制
qpskMod = comm.QPSKModulator;
qpskDemod = comm.QPSKDemodulator;
txSig = qpskMod(data_symbols);

% 加入色散
span = 20;                      % 矩形滤波器的长度（符号数）
rolloff = 0.2;                  % 瑞利因数
sps = 4;                        % 每个符号的样本数
Fs = 100e9;                     % 采样率
h = rcosdesign(rolloff,span,sps);
txSig = upfirdn(txSig,h,sps);

% 色散补偿
D = 17e-12;                     % 色散时间延迟
lambda = 1550e-9;               % 波长
Dp = 17e-6;                     % 色散参数
txSig = dspec(txSig,Fs,D,Dp,lambda);

% 输出补偿前后的星座图
subplot(2,2,1)
plot(real(txSig),imag(txSig),'r.')
axis([-2 2 -2 2])
title('星座图-补偿前')
xlabel('I')
ylabel('Q')

rxSig = txSig;

subplot(2,2,2)
plot(real(rxSig),imag(rxSig),'b.')
axis([-2 2 -2 2])
title('星座图-补偿后')
xlabel('I')
ylabel('Q')

% 输出补偿前后的眼图
eyediagram(txSig,2*sps)

rxSig = txSig;

eyediagram(rxSig,2*sps)

% 输出补偿前后的时域波形
subplot(2,2,3)
plot(0:1/Fs:(length(txSig)-1)/Fs,real(txSig),'r')
title('时域波形-补偿前')
xlabel('时间 (s)')
ylabel('幅度')

subplot(2,2,4)
plot(0:1/Fs:(length(rxSig)-1)/Fs,real(rxSig),'b')
title('时域波形-补偿后')
xlabel('时间 (s)')
ylabel('幅度')

% 输出补偿前后的频域波形
subplot(2,2,3)
plot(linspace(-Fs/2,Fs/2,length(txSig)),10*log10(abs(fftshift(fft(txSig)))))
title('频域波形-补偿前')
xlabel('频率 (Hz)')
ylabel('功率谱密度 (dBm/Hz)')

subplot(2,2,4)
plot(linspace(-Fs/2,Fs/2,length(rxSig)),10*log10(abs(fftshift(fft(rxSig)))))
title('频域波形-补偿后')
xlabel('频率 (Hz)')
ylabel('功率谱密度 (dBm/Hz)')

这段代码将生成一个QPSK调制信号并加入色散，然后进行色散补偿，并输出补偿前后的星座图、眼图、时域波形和频域波形。