利用时域或频域均衡的方法，从数字域完成色散补偿matlab代码

以下是一个简单的基于频域均衡的色散补偿的 Matlab 代码示例：

% 设置参数
fs = 10e9; % 采样率
N = 1024; % FFT点数
f = linspace(-fs/2,fs/2,N); % 频率向量
T = 1/fs; % 采样周期
L = 80; % 信道长度，单位为km
D = 17; % 色散参数，单位为ps/(nm*km)
beta2 = 21.7e-27; % 二阶色散参数，单位为s^2/km
lambda = 1550e-9; % 波长

% 生成信号
t = (0:N-1)*T; % 时间向量
f0 = 193.1e12; % 中心频率
s = exp(1j*2*pi*f0*t); % 原始信号

% 传输信号，加入色散
h = exp(1j*pi/2*D*(lambda*f).^2*L); % 信道冲激响应
s_ch = ifft(fft(s).*fft(h)); % 传输信号

% 频域均衡
H = exp(-1j*pi/2*D*(lambda*f).^2*L); % 频域均衡器
s_eq = ifft(fft(s_ch)./fft(H)); % 均衡后的信号

% 绘图
figure;
subplot(2,1,1);
plot(t*1e12, abs(s_ch).^2);
xlabel('时间 (ps)');
ylabel('功率 (W)');
title('传输信号');

subplot(2,1,2);
plot(t*1e12, abs(s_eq).^2);
xlabel('时间 (ps)');
ylabel('功率 (W)');
title('均衡后的信号');

该代码在时域生成一个中心频率为 193.1 THz 的正弦信号，然后通过频域均衡器进行色散补偿，最后绘制传输信号和均衡后的信号的功率图。