来个色散补偿的matlab代码例子

以下是一个简单的色散补偿的Matlab代码例子：

% 生成一段信号
f_s = 10e9; % 采样率
t = 0:1/f_s:1e-6;
f_c = 2.5e9; % 中心频率
s = cos(2*pi*f_c*t);

% 添加色散
beta2 = 10e-27; % 色散参数
D = beta2 * (2*pi*f_c)^2; % 色散参数
s_disp = s .* exp(1i * 0.5 * D * t.^2);

% 绘制时域和频域图像
figure;
subplot(2,1,1);
plot(t, real(s), 'b', t, real(s_disp), 'r');
xlabel('Time (s)');
ylabel('Amplitude');
legend('Original signal', 'Dispersion compensated signal');
subplot(2,1,2);
f = linspace(-f_s/2, f_s/2, length(s));
s_fft = fftshift(fft(s));
s_disp_fft = fftshift(fft(s_disp));
plot(f, abs(s_fft), 'b', f, abs(s_disp_fft), 'r');
xlabel('Frequency (Hz)');
ylabel('Magnitude');
legend('Original signal', 'Dispersion compensated signal');

解释：该代码生成一个中心频率为2.5 GHz的正弦波信号，并添加了一个色散参数beta2=10e-27。然后绘制了原始信号和色散补偿后的信号的时域和频域图像，以便比较它们的形状和频谱。