[4yp.git] / chromaticDispersionTest.m

Rsym = 2.5e10; % symbol/sec
sps = 4; % samples/symbol
fs = Rsym * sps; % sampling freq (Hz)
Tsamp = 1 / fs;

x = 10 * normpdf(1:30, 10, 3);

xCD = chromaticDispersion(x, 20, 1550, 30, Tsamp);

figure(10001);
subplot(121);
plot(abs(x) .^ 2);
axis([0 30 -Inf Inf]);
title('Before dispersion');
ylabel('Power (arb. unit)');
xlabel('Time (arb. unit)');
subplot(122);
plot(abs(xCD) .^ 2)
axis([20 50 -Inf Inf]);
title('After dispersion');
xlabel('Time (arb.~unit)');
ylabel('Power (arb.~unit)');

subplot(121);
formatFigure;
subplot(122)
formatFigure;
set(figure(10001), 'Units', 'centimeters', ...
    'Position', [0 0 24 8], 'PaperPositionMode', 'auto');
subplot(121);
set(gca, 'XTick', 0:5:30);
subplot(122);
set(gca, 'XTick', 20:5:50);

%saveas(gcf, 'dispersionTest.png');
Commit	Line	Data
1eeb62fb AIL	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
f9a73e9e AIL	10	figure(10001);
f9a73e9e AIL	11	subplot(121);
1eeb62fb AIL	12	plot(abs(x) .^ 2);
	13	axis([0 30 -Inf Inf]);
	14	title('Before dispersion');
	15	ylabel('Power (arb. unit)');
f9a73e9e AIL	16	xlabel('Time (arb. unit)');
f9a73e9e AIL	17	subplot(122);
1eeb62fb AIL	18	plot(abs(xCD) .^ 2)
	19	axis([20 50 -Inf Inf]);
	20	title('After dispersion');
f9a73e9e AIL	21	xlabel('Time (arb.~unit)');
	22	ylabel('Power (arb.~unit)');
	23
	24	subplot(121);
	25	formatFigure;
	26	subplot(122)
	27	formatFigure;
	28	set(figure(10001), 'Units', 'centimeters', ...
	29	'Position', [0 0 24 8], 'PaperPositionMode', 'auto');
	30	subplot(121);
	31	set(gca, 'XTick', 0:5:30);
	32	subplot(122);
	33	set(gca, 'XTick', 20:5:50);
1eeb62fb AIL	34
1eeb62fb AIL	35	%saveas(gcf, 'dispersionTest.png');