chromaticDispersionTest.m

   1 Rsym = 2.5e10; % symbol/sec
   2 sps = 4; % samples/symbol
   3 fs = Rsym * sps; % sampling freq (Hz)
   4 Tsamp = 1 / fs;
   5
   6 x = 10 * normpdf(1:30, 10, 3);
   7
   8 xCD = chromaticDispersion(x, 20, 1550, 30, Tsamp);
   9
  10 figure(1);
  11 subplot(211);
  12 plot(abs(x) .^ 2);
  13 axis([0 30 -Inf Inf]);
  14 title('Before dispersion');
  15 ylabel('Power (arb. unit)');
  16 subplot(212);
  17 plot(abs(xCD) .^ 2)
  18 axis([20 50 -Inf Inf]);
  19 title('After dispersion');
  20 xlabel('Time (arb. unit)');
  21 ylabel('Power (arb. unit)');
  22
  23 %saveas(gcf, 'dispersionTest.png');