phaseNoise.m

   1 function [xPN, phiTx_phiLO]  = phaseNoise(x, linewidthTx, linewidthLO, Tsamp)
   2   %% fIF = 15e3; % fIF = fc - fLO
   3   %% carrier = exp(-j * (2 * pi * fIF * t - phi);
   4   %% Neglecting fIF for now...
   5
   6   dphiTx = normrnd(0, sqrt(2 * pi * linewidthTx * Tsamp), length(x), 1);
   7   dphiLO = normrnd(0, sqrt(2 * pi * linewidthLO * Tsamp), length(x), 1);
   8   phiTx = cumsum(dphiTx);
   9   phiLO = cumsum(dphiLO);
  10
  11   phiTx_phiLO = phiTx - phiLO;
  12
  13   pn = exp(-j * phiTx_phiLO);
  14
  15   xPN = x .* pn;
  16 end