[4yp.git] / phaseNoise.m

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
Commit	Line	Data
1eeb62fb AIL	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