1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
1;
load '/tmp/xi_vs_power.mat' ;
Er=(1+I*xi_right)*E_field_pos_freq.right;
El=(1+I*xi_left) *E_field_pos_freq.left;
Ex=(Er+El)/sqrt(2);
Ey=I*(Er-El)/sqrt(2);
%extra rotation to compensate rotation due to ellipticity
% actually no need for it since x-polarization shifts by positive phase
% and y-pol by negative phase
%el_rot=0*psi_el;
%Ex=cos(el_rot)*Ex-sin(el_rot)*Ey;
%Ey=sin(el_rot)*Ex+cos(el_rot)*Ey;
Ipos=(abs(Ey).^2)/2;
Ineg=(abs(Ex).^2)/2;
figure(1);
hold off;
plot(Ep.^2, real(xi_left-xi_right), '-');
title("differential real xi");
xlabel("two photon detuning");
figure(2);
hold off;
plot(Ep.^2, imag(xi_left-xi_right), '-');
title("differential imag xi");
xlabel("two photon detuning");
figure(3);
hold off;
plot(Ep.^2, imag(xi_left), '-', Ep, imag(xi_right), '-');
title("imag xi");
xlabel("two photon detuning");
figure(4);
hold off;
plot(Ep.^2, real(xi_left), '-', Ep.^2, real(xi_right), '-');
title("real xi");
xlabel("two photon detuning");
figure(5);
hold off;
%plot(Ep.^2, (Ipos-Ineg), '-');
semilogx(Ep.^2, (Ipos-Ineg), '-');
%semilogx(Ep.^2, (Ipos-Ineg)./(Ep.^2), '-');
title("BPD signal xi");
xlabel("two photon detuning");
%figure(1);
%hold off;
%plot(detuning_freq, imag(xi_linear), '-1;linear;');
%hold on;
%plot(detuning_freq, imag(xi_left), '-2;left;');
%plot(detuning_freq, imag(xi_right), '-3;right;');
%title("probe absorption");
%hold off;
%figure(2);
%hold off;
%plot(detuning_freq, real(xi_linear), '-1;linear;');
%hold on;
%plot(detuning_freq, real(xi_left), '-2;left;');
%plot(detuning_freq, real(xi_right), '-3;right;');
%title("probe dispersion");
%hold off;
|
