Phase Noise

Oscillator Phase Noise Model

Description

function Sout = add_phase_noise( Sin, Fs, phase_noise_freq, phase_noise_power )
 
  Oscillator Phase Noise Model
  
   INPUT:
      Sin - input COMPLEX signal
      Fs - sampling frequency ( in Hz ) of Sin
      phase_noise_freq - frequencies at which SSB Phase Noise is defined (offset from carrier in Hz)
      phase_noise_power - SSB Phase Noise power ( in dBc/Hz )
 
   OUTPUT:
      Sout - output COMPLEX phase noised signal
 
   NOTE:
      Input signal should be complex
 
   EXAMPLE ( How to use add_phase_noise ):
          Assume SSB Phase Noise is specified as follows:
       -------------------------------------------------------
       | Offset From Carrier | Phase Noise |
       -------------------------------------------------------
       | 1 kHz | -84 dBc/Hz |
       | 10 kHz | -100 dBc/Hz |
       | 100 kHz | -96 dBc/Hz |
       | 1 MHz | -109 dBc/Hz |
       | 10 MHz | -122 dBc/Hz |
       -------------------------------------------------------
 
       Assume that we have 10000 samples of complex sinusoid of frequency 3 KHz
       sampled at frequency 40MHz:
        
        Fc = 3e3; % carrier frequency
        Fs = 40e6; % sampling frequency
        t = 0:9999;
        S = exp(j*2*pi*Fc/Fs*t); % complex sinusoid
 
       Then, to produse phase noised signal S1 from the original signal S run follows:
 
        Fs = 40e6;
        phase_noise_freq = [ 1e3, 10e3, 100e3, 1e6, 10e6 ]; % Offset From Carrier
        phase_noise_power = [ -84, -100, -96, -109, -122 ]; % Phase Noise power
        S1 = add_phase_noise( S, Fs, phase_noise_freq, phase_noise_power );

• MATLAB release: MATLAB 6.5.1 (R13SP1)
• Other requirements: -

[Ming-Yu Hsieh]

This model is very adequate for ofdm system.

[Madhukar Ramamurthy]

 

[thomas höhne]

Great work. Easy to adapt to own needs. The plot of the generated PSD could be done better.

[john kedziora]

Works very well!

[chuntel Lin]

thank you

[CT Lin]

thank you

[Agus Suhendar]

i have problem, hope u can help me..
after qam mapping and ifft the signal is in time domain and complex form, with BW=7 MHz the problem is how to sampling the complex time domain signal with sampling factor=8/7 or Fs=8MHz
thanxs b4..

[Nick A]

This approach is very elegant. However, it limits the amount of low frequency phase noise extracted from the PSD by the original sampling frequency. Any ideas about circumventing this limitation?
Thanks
Nick

[wendy]

how come i run this program but it says Not enough input arguments. ?

[xu zhenhua]

xie xie

[sia yousef]

Hi I wonder how this model can satisfy the real shape of the oscillator that firs we have a 1/f^3 and then 1/f^2 and then noise floor.
and it doesn't have close to carrier phase noise.
I don't know how based on which model you generate phase noise?
thanks