5 Rsym = 2.5e10; % symbol rate (sym/sec)
7 span = 6; % filter span
8 sps = 4; % samples per symbol
10 fs = Rsym * sps; % sampling freq (Hz)
13 t = (0 : 1 / fs : numSymbs / Rsym + (1.5 * span * sps - 1) / fs).';
17 EbN0 = 10 .^ (EbN0_db ./ 10);
23 EsN0 = EbN0 .* log2(M);
24 EsN0_db = 10 .* log10(EsN0);
26 plotlen = length(EbN0);
28 ber = zeros(1, plotlen);
30 data = randi([0 M - 1], numSymbs, 1);
31 %%modData = dpskmod(data, M, pi / M, 'gray');
32 modData = dpskmod(data, M, 0, 'gray');
34 x = txFilter(modData, rolloff, span, sps);
37 %%linewidthLO = 1e6; % Hz
38 linewidthLO = Rsym * 1e-3;
44 [xPN, pTxLO] = phaseNoise(x, linewidthTx, linewidthLO, Tsamp);
47 snr = EbN0_db(i) + 10 * log10(log2(M)) - 10 * log10(sps);
48 noiseEnergy = 10 ^ (-snr / 10);
50 y = awgn(xPN, snr, 'measured');
52 r = rxFilter(y, rolloff, span, sps);
54 r = normalizeEnergy(r, numSymbs, 1 + noiseEnergy);
56 rSampled = r(sps*span/2+1:sps:(numSymbs+span/2)*sps);
57 [rSaPhEq, phiests] = phaseNoiseCorr(rSampled, M, avgSa);
59 adaptiveFilterOut = adaptiveCMA(rSaPhEq.');
61 demodData = dpskdemod(rSaPhEq, M, 0, 'gray').';
62 demodAdapt = dpskdemod(adaptiveFilterOut, M, 0, 'gray');
64 [~, ber(i)] = biterr(data, demodData);
68 plot(repelem(-phiests, sps));
71 legend('estimate', 'actual');
84 %% Plot simulated results
85 semilogy(EbN0_db, ber, 'r', 'LineWidth', 2);
88 theoreticalPSK(EbN0_db, M, 'b', 'LineWidth', 1);
89 legend({'Simulated phase noise + correction', 'Theoretical AWGN'}, ...
90 'Location', 'southwest');
93 title({'QPSK with phase nosie and correction', ...
94 strcat(num2str(numSymbs * log2(M) / 1e3), '~kbit, LO~', ...
95 num2str(linewidthLO / 1e6), '~MHz, Av~', num2str(avgSa), ...
98 xlabel('$E_b/N_0$ (dB)');