1 % Contour plots
   2 
   3 ! perl mkat lcgt
   4 
   5 load lcgt_run1.out
   6 load lcgt_run2.out
   7 load lcgt_run3.out
   8 load lcgt_run4.out
   9 load lcgt_run5.out
  10 load lcgt_run6.out
  11 load lcgt_run7.out
  12 load lcgt_run8.out
  13 load lcgt_run9.out
  14 load lcgt_run10.out
  15 load lcgt_run11.out
  16 load lcgt_run12.out
  17 load lcgt_run13.out
  18 load lcgt_run14.out
  19 
  20 Nxy = sqrt( size(lcgt_run1, 1 ) );
  21 phase1 = zeros( Nxy, Nxy );
  22 phase2 = zeros( Nxy, Nxy );
  23 SPdc = zeros( Nxy, Nxy );
  24 SPlp= zeros( Nxy, Nxy );
  25 SPlm= zeros( Nxy, Nxy );
  26 SPls= zeros( Nxy, Nxy );
  27 APdc = zeros( Nxy, Nxy );
  28 APlp= zeros( Nxy, Nxy );
  29 APlm= zeros( Nxy, Nxy );
  30 APls= zeros( Nxy, Nxy );
  31 POdc = zeros( Nxy, Nxy );
  32 POlp= zeros( Nxy, Nxy );
  33 POlm= zeros( Nxy, Nxy );
  34 POls= zeros( Nxy, Nxy );
  35 N = size(lcgt_run13,1);
  36 SPLp= zeros( N, 2 );
  37 APLm= zeros( N, 2 );
  38 
  39 for ind = 1 : Nxy
  40 phase1(ind,1:Nxy) = lcgt_run1( 1 : Nxy : 1+Nxy*(Nxy-1), 1 )';
  41 phase2(1:Nxy,ind) = lcgt_run1( 1 : Nxy, 2 );
  42 SPdc(ind,1:Nxy) = lcgt_run1( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  43 SPlp(ind,1:Nxy) = lcgt_run2( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  44 SPlm(ind,1:Nxy) = lcgt_run3( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  45 SPls(ind,1:Nxy) = lcgt_run4( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  46 APdc(ind,1:Nxy) = lcgt_run5( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  47 APlp(ind,1:Nxy) = lcgt_run6( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  48 APlm(ind,1:Nxy) = lcgt_run7( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  49 APls(ind,1:Nxy) = lcgt_run8( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  50 POdc(ind,1:Nxy) = lcgt_run9( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  51 POlp(ind,1:Nxy) = lcgt_run10( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  52 POlm(ind,1:Nxy) = lcgt_run11( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';%
  53 POls(ind,1:Nxy) = lcgt_run12( ind : Nxy : ind+Nxy*(Nxy-1), 3 )';
  54 
  55 end
  56 SPLp=lcgt_run13;
  57 APLm=lcgt_run14;
  58 
  59 figure(1)
  60 axis([0 180 0 180])
  61 axis square
  62 box on
  63 %plot (0,0,'k :', 0,0, 'b', 0,0, 'r', 0,0, 'g')
  64 hold on
  65 plot ( 0,0, 'b', 0,0, 'r', 0,0, 'g')
  66 %legend ('dc=0', 'l+', 'l-', 'ls')
  67 legend ( 'l+', 'l-', 'ls')
  68 %contour(phase1, phase2, SPdc,1, 'k:');%1 means "conour level 1" ??
  69 [C,h] = contour(phase1, phase2, SPlp, 'b');
  70 clabel(C,h, 'FontSize', 10)
  71 [C,h] = contour(phase1, phase2, SPlm, 'r');
  72 clabel(C,h, 'FontSize', 10)
  73 [C,h] = contour(phase1, phase2, SPls, 'g');
  74 clabel(C,h, 'FontSize', 10)
  75 
  76 set(gca,'xtick',[0 45 90 135 180])
  77 set(gca,'ytick',[0 45 90 135 180])
  78 
  79 
  80 h=findobj('LineStyle',':');
  81 set (h, 'LineWidth', 3);
  82 h=findobj('LineWidth', 0.5);;
  83 set (h, 'LineWidth', 2);
  84 h=findobj('FontSize', 9);
  85 set (h, 'FontSize', 20);
  86 
  87 %xlabel('Demodulation Phase of f1','FontSize',16)
  88 %ylabel('Demodulation Phase of f2','FontSize',16)
  89 title('Double Demodulation at SP','FontSize',12)
  90 h=findobj('LineStyle',':');
  91 set (h, 'LineWidth', 3);
  92 h=findobj('LineWidth', 0.5);
  93 set (h, 'LineWidth', 1);
  94 h=findobj('FontSize', 9);
  95 set (h, 'FontSize', 12);
  96 h = gcf;
  97 saveas(h,'SP.fig');
  98 saveas(h,'SP.emf');
  99 %clf
 100 
 101 
 102 figure(2)
 103 axis([0 180 0 180])
 104 axis square
 105 box on
 106 %plot (0,0,'k :', 0,0, 'b', 0,0, 'r', 0,0, 'g')
 107 hold on
 108 plot ( 0,0, 'b', 0,0, 'r', 0,0, 'g')
 109 %legend ('dc=0', 'l+', 'l-', 'ls')
 110 legend ( 'l+', 'l-', 'ls')
 111 %contour(phase1, phase2, APdc,1, 'k:');
 112 [C,h] = contour(phase1, phase2, APlp, 'b');
 113 clabel(C,h, 'FontSize', 10)
 114 [C,h] = contour(phase1, phase2, APlm, 'r');
 115 clabel(C,h, 'FontSize', 10)
 116 [C,h] = contour(phase1, phase2, APls, 'g');
 117 clabel(C,h, 'FontSize', 10)
 118 
 119 set(gca,'xtick',[0 45 90 135 180])
 120 set(gca,'ytick',[0 45 90 135 180])
 121 
 122 
 123 h=findobj('LineStyle',':');
 124 set (h, 'LineWidth', 1);
 125 h=findobj('LineWidth', 0.5);
 126 set (h, 'LineWidth', 1);
 127 h=findobj('FontSize', 3);
 128 set (h, 'FontSize', 5);
 129 
 130 %xlabel('Demodulation Phase of f1','FontSize',16)
 131 %ylabel('Demodulation Phase of f2','FontSize',16)
 132 title('Double Demodulation at AP','FontSize',12)
 133 h=findobj('LineStyle',':');
 134 set (h, 'LineWidth', 3);
 135 h=findobj('LineWidth', 0.5);
 136 set (h, 'LineWidth', 1);
 137 h=findobj('FontSize', 9);
 138 set (h, 'FontSize', 12);
 139 h = gcf;
 140 saveas(h,'AP.fig');
 141 saveas(h,'AP.emf');
 142 %clf
 143 
 144 figure(3)
 145 axis([0 180 0 180])
 146 axis square
 147 box on
 148 %plot (0,0,'k :', 0,0, 'b', 0,0, 'r', 0,0, 'g')
 149 hold on
 150 plot ( 0,0, 'b', 0,0, 'r', 0,0, 'g')
 151 %legend ('dc=0', 'l+', 'l-', 'ls')
 152 legend ( 'l+', 'l-', 'ls')
 153 %contour(phase1, phase2, POdc,1, 'k:');
 154 [C,h] = contour(phase1, phase2, POlp, 'b');
 155 clabel(C,h, 'FontSize', 10)
 156 [C,h] = contour(phase1, phase2, POlm, 'r');
 157 clabel(C,h, 'FontSize', 10)
 158 [C,h] = contour(phase1, phase2, POls, 'g');
 159 clabel(C,h, 'FontSize', 10)
 160 %xlabel('Demodulation Phase of f1','FontSize',16)
 161 %ylabel('Demodulation Phase of f2','FontSize',16)
 162 title('Double Demodulation at PO','FontSize',12)
 163 
 164 set(gca,'xtick',[0 45 90 135 180])
 165 set(gca,'ytick',[0 45 90 135 180])
 166 
 167 
 168 h=findobj('LineStyle',':');
 169 set (h, 'LineWidth', 1);
 170 h=findobj('LineWidth', 0.5);
 171 %set (h, 'LineWidth', 1);
 172 set (h, 'LineWidth', 1);
 173 h=findobj('FontSize', 9);
 174 %h=findobj('FontSize', 3);
 175 set (h, 'FontSize', 12);
 176 %set (h, 'FontSize', 5);
 177 h = gcf;
 178 saveas(h,'PO.fig');
 179 saveas(h,'PO.emf');