git.gag.com Git - debian/gnuradio/blob - gnuradio-examples/python/ofdm/plot_ofdm.m

   1 function plot_ofdm(fft_size, occ_tones)
   2
   3 ofdm = read_complex_binary('ofdm_frame_acq_c.dat');
   4 ofdm_split = split_vect(ofdm, occ_tones);
   5
   6 ofdm_derot = read_complex_binary('ofdm_frame_sink_c.dat');
   7 ofdm_derot_split = split_vect(ofdm_derot, occ_tones);
   8
   9 fftc = read_complex_binary('fft_out_c.dat');
  10 fftc_split = split_vect(fftc, fft_size);
  11
  12 size(ofdm_split)
  13 size(ofdm_derot_split)
  14 disp "DEROTATED SPLIT"
  15 ofdm_derot(1:100)
  16
  17
  18 figure(1)
  19 #set(gcf, 'Position', [50 50 1000 600]);
  20
  21 a = size(ofdm_split);
  22 if nargin == 3
  23     maxcount = count;
  24     if maxcount > a(1)
  25         error('plot_ofdm:tolong', 'Requested count size exceeds size of vectors');
  26     end
  27 else
  28     maxcount = a(1);
  29 end
  30
  31 for i = 1:size(ofdm_split)[0]
  32     x = ofdm_split(i,:);
  33     y = fftc_split(i+1,:);
  34
  35     subplot(2,2,1);
  36     plot(real(x), imag(x), 'bo');
  37     #set(gca, 'FontSize', 30, 'FontWeight', 'Bold');
  38     axis([-1.5, 1.5, -1.5, 1.5])
  39     #title('I&Q Constellation', 'FontSize', 36);
  40     #xlabel('Inphase', 'FontSize', 32);
  41     #ylabel('Quadrature', 'FontSize', 32);
  42
  43     subplot(2,2,3);
  44     plot(angle(x*j), 'bo');
  45     #set(gca, 'FontSize', 30, 'FontWeight', 'Bold');
  46     axis([0, occ_tones, -3.5, 3.5])
  47     #title('Equalized Symbol Angle', 'FontSize', 36);
  48     #xlabel('Bin Number (Occ. Tones)', 'FontSize', 32);
  49     #ylabel('Symbol Angle', 'FontSize', 32);
  50
  51     subplot(2,2,2);
  52     plot(angle(y*j), 'bo');
  53     #set(gca, 'FontSize', 30, 'FontWeight', 'Bold');
  54     axis([0, fft_size, -3.5, 3.5])
  55     #title('Unequalized Symbol Angle', 'FontSize', 36);
  56     #xlabel('Bin Number (FFT Size)', 'FontSize', 32);
  57     #ylabel('Symbol Angle', 'FontSize', 32);
  58
  59     Y = 20*log10(abs(y) ./ max(abs(y)));
  60     subplot(2,2,4);
  61     plot(Y, 'b-');
  62     #set(gca, 'FontSize', 30, 'FontWeight', 'Bold');
  63     axis([0, fft_size, -50, 1]);
  64     #title('Frequency Domain of Unequalized Rx', 'FontSize', 36);
  65     #xlabel('Bin Number (FFT Size)', 'FontSize', 32);
  66     #ylabel('Power (dB)', 'FontSize', 32);
  67
  68     #N = 20*log10(var(abs(x)-1));
  69
  70     disp(sprintf('Symbol Number: %d\n', i))
  71     #disp(sprintf('\tFreq Error: %f\n', anglesh_pn(1+(i-1)*fft_size)))
  72     pause
  73
  74 end