function [xPN, phiTx_phiLO]  = phaseNoise(x, linewidthTx, linewidthLO, Tsamp)
  %% fIF = 15e3; % fIF = fc - fLO
  %% carrier = exp(-j * (2 * pi * fIF * t - phi);
  %% Neglecting fIF for now...

  dphiTx = normrnd(0, sqrt(2 * pi * linewidthTx * Tsamp), length(x), 1);
  dphiLO = normrnd(0, sqrt(2 * pi * linewidthLO * Tsamp), length(x), 1);
  phiTx = cumsum(dphiTx);
  phiLO = cumsum(dphiLO);

  phiTx_phiLO = phiTx - phiLO;

  pn = exp(-j * phiTx_phiLO);

  xPN = x .* pn;
end