diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2011-11-16 16:33:49 -0500 |
|---|---|---|
| committer | Eugeniy E. Mikhailov <evgmik@gmail.com> | 2020-09-21 16:29:52 -0400 |
| commit | f25ef4ac36b03b1c4e26ef0aa14905734c2db9fe (patch) | |
| tree | 8fd3a754993193f81488b72136bdf2f086584ceb /faraday/basis_transformation.m | |
| parent | b06d36c8981fb5a3859099b84e1ce442b5364b00 (diff) | |
| download | multi_mode_eit-f25ef4ac36b03b1c4e26ef0aa14905734c2db9fe.tar.gz multi_mode_eit-f25ef4ac36b03b1c4e26ef0aa14905734c2db9fe.zip | |
clean up of unused files and link to the master copy
Diffstat (limited to 'faraday/basis_transformation.m')
| l---------[-rw-r--r--] | faraday/basis_transformation.m | 58 |
1 files changed, 1 insertions, 57 deletions
diff --git a/faraday/basis_transformation.m b/faraday/basis_transformation.m index 8b8eb04..28d8beb 100644..120000 --- a/faraday/basis_transformation.m +++ b/faraday/basis_transformation.m @@ -1,57 +1 @@ -1; - -% matrix of circular to linear transformation -% [x, y, z]' = lin2circ * [r, l, z]' -function transformation_matrix = circ2lin() - transformation_matrix = ... - [ ... - [ 1/sqrt(2), 1/sqrt(2), 0]; ... - [ 1i/sqrt(2),-1i/sqrt(2), 0]; ... - [ 0, 0, 1] ... - ]; -endfunction - - -% matrix of linear to circular transformation -% [r, l, z]' = lin2circ * [x, y, z]' -function transformation_matrix = lin2circ() - transformation_matrix = ... - [ ... - [ 1/sqrt(2), -1i/sqrt(2), 0]; ... - [ 1/sqrt(2), 1i/sqrt(2), 0]; ... - [ 0, 0, 1] ... - ]; -endfunction -% linear basis rotation -% x axis untouched -% z and y rotated by angle theta around 'x' axis -% [x_new, y_new, z_new]' = oldlin2newlin * [x_old, y_old, z_old]' -function oldlin2newlin_m = oldlin2newlin(theta) - oldlin2newlin_m = [ ... - [ 1, 0, 0]; ... - [ 0, cos(theta), -sin(theta)]; ... - [ 0, sin(theta), cos(theta)]... - ]; -endfunction - -% rotate x polarized light by angle phi around -% light propagation axis (Z) -function [E_field_x, E_field_y] = rotXpolarization(phi, E_field_linear) - % important negative frequency behave as they rotate in opposite direction - E_field_x=cos(phi)*E_field_linear; - E_field_y=sin(phi)*E_field_linear; -endfunction - -% transform x,y,z linearly polarized light in the lab/light system coordinate -% to left, right, linear along z atom system of coordinate -% atom magnetic field is along new axis Z wich is at angle theta with respect to -% light propagation direction -function E_field_pos_freq=xyz_lin2atomic_axis_polarization(theta, E_field_lab_pos_freq) - coord_transf_m = lin2circ() * oldlin2newlin(theta); - E_field_pos_freq.right = coord_transf_m(1,1)*E_field_lab_pos_freq.x + coord_transf_m(1,2)*E_field_lab_pos_freq.y + coord_transf_m(1,3)*E_field_lab_pos_freq.z; - E_field_pos_freq.left = coord_transf_m(2,1)*E_field_lab_pos_freq.x + coord_transf_m(2,2)*E_field_lab_pos_freq.y + coord_transf_m(2,3)*E_field_lab_pos_freq.z; - E_field_pos_freq.linear = coord_transf_m(3,1)*E_field_lab_pos_freq.x + coord_transf_m(3,2)*E_field_lab_pos_freq.y + coord_transf_m(3,3)*E_field_lab_pos_freq.z; -endfunction - - -% vim: ts=2:sw=2:fdm=indent +../basis_transformation.m
\ No newline at end of file |