% 生成OOK信号 fs = 1e6; % 采样率 T = 1/fs; % 采样时间间隔 f = 10e3; % 载波频率 duration = 1; % 信号持续时间 t = 0:T:duration-T; % 时间序列 data = randi([0 1], 1, length(t)); % 随机生成0和1的数据 signal = data.*sin(2*pi*f*t); % OOK信号 % 光纤传输 len = 10; % 光纤长度（km） lamda = 1550; % 中心波长（nm） c = 3e8; % 光速（m/s） D = 17; % 色散系数（ps/nm/km） beta2 = -D*(lamda*1e-9)^2/(2*pi*c); % 色散参数 L = len*1e3; % 光纤长度（m） wavelength = lamda*1e-9; % 光波长（m） span = 10; % 传输距离间隔（km） numSpans = L/span; % 总传输距离间隔数 spanLen = span*1e3; % 单个传输距离间隔长度（m） dispComp = exp(1j*0.5*beta2*wavelength^2*L*T^2); % 色散补偿系数 signal_out = zeros(size(signal)); % 接收信号 for i = 1:numSpans startIdx = (i-1)*spanLen/T+1; endIdx = i*spanLen/T; signal_span = signal(startIdx:endIdx); % 当前距离间隔内的信号 signal_span = ifft(fft(signal_span).*dispComp); % 色散补偿 signal_out(startIdx:endIdx) = signal_span; % 累加接收信号 end % 绘制信号波形 figure; subplot(2,1,1); plot(t, signal); xlabel('Time (s)'); ylabel('Amplitude'); title('Original Signal'); subplot(2,1,2); plot(t, abs(signal_out)); xlabel('Time (s)'); ylabel('Amplitude'); title('Received Signal after Dispersion Compensation');

这段代码的作用是生成一个OOK信号，并通过模拟光纤传输的方式进行色散补偿，最终绘制出接收信号的波形。其中，OOK信号的频率为10kHz，持续时间为1秒，随机生成0和1的数据，通过乘以正弦波的方式生成信号。接着，通过设定光纤长度、中心波长、色散系数等参数，计算出色散参数，并通过循环模拟光纤传输的过程，每隔10km对信号进行一次色散补偿，并累加接收信号。最后，绘制出原始信号和接收信号的波形图。