Working Kerr effect; PDM; speedups; removed unused files

[4yp.git] / passband.m

diff --git a/passband.m b/passband.m
deleted file mode 100644 (file)
index 9e88af6..0000000
--- a/passband.m
+++ /dev/null
@@ -1,79 +0,0 @@
-function passband(rolloff, M, numSymbs)
-  %% Set defaults for inputs
-  if nargin < 3
-    numSymbs = 1000;
-  end
-  if nargin < 2
-    M = 2;
-  end
-  if nargin < 1
-    rolloff = 0.5;
-  end
-
-  %% https://www.mathworks.com/help/comm/examples/passband-modulation-with-adjacent-channel-interference.html
-  Rsym = 2.5e10; % symbol rate (sym/sec)
-
-  span = 6; % filter span
-  sps = 4; % samples per symbol
-
-  fs = Rsym * sps; % sampling freq (Hz)
-
-  t = (0 : 1 / fs : numSymbs / Rsym + (1.5 * span * sps - 1) / fs)';
-
-
-  EbN0_db = 0:0.2:10;
-  EbN0 = 10 .^ (EbN0_db ./ 10);
-
-  Es = 1;
-  Eb = Es / log2(M);
-  N0 = Eb ./ EbN0;
-
-  EsN0 = EbN0 .* log2(M);
-  EsN0_db = 10 .* log10(EsN0);
-
-  plotlen = length(EbN0);
-
-  ber = zeros(1, plotlen);
-
-
-  data = randi([0 M - 1], numSymbs, 1);
-  modData = pskmod(data, M, 0, 'gray');
-
-  xBaseband = txFilter(modData, rolloff, span, sps);
-
-  fc = 1e12; % Carrier freq (Hz)
-  carrier = sqrt(2) * exp(j * 2 * pi * fc * t);
-
-  xPassband = xBaseband .* carrier(1:length(xBaseband));
-
-  for i = 1:plotlen
-    snr = EbN0_db(i) + 10 * log10(log2(M)) - 10 * log10(sps);
-    noiseEnergy = 10 ^ (-snr / 10);
-
-    yPassband = awgn(xPassband, snr, 'measured');
-    yLO = yPassband .* conj(carrier(1:length(yPassband)));
-    rBaseband = rxFilter(yLO, rolloff, span, sps);
-
-    rSampled = rBaseband(sps*span/2+1:sps:(numSymbs + span/2) * sps);
-
-    demodData = pskdemod(rSampled, M, 0, 'gray');
-    [bitErrors, ber(i)] = biterr(data, demodData);
-  end
-
-  figure(1);
-  clf;
-
-  %% Plot simulated results
-  semilogy(EbN0_db, ber, 'r', 'LineWidth', 2);
-  hold on;
-
-  theoreticalPSK(EbN0_db, M, 'b', 'LineWidth', 1);
-  legend('Simulated RRC', 'Discrete');
-
-  title(strcat(num2str(M), '-PSK RRC with Gray code'));
-  grid on;
-  xlabel('$E_b/N_0$ (dB)');
-  ylabel('BER');
-
-  formatFigure;
-end