diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-01-21 22:24:43 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2010-01-21 22:24:43 +0000 |
| commit | dbc850d4eb1969c3c1818faec8b0ac3838d7134c (patch) | |
| tree | 04b69a51043b6ef0f9f809938d8c8676b9071432 /useful_functions.m | |
| parent | 32c30704e9debda5604342d2b316c0ef948eb111 (diff) | |
| download | multi_mode_eit-dbc850d4eb1969c3c1818faec8b0ac3838d7134c.tar.gz multi_mode_eit-dbc850d4eb1969c3c1818faec8b0ac3838d7134c.zip | |
proper calculation of light transmission
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 22 |
1 files changed, 16 insertions, 6 deletions
diff --git a/useful_functions.m b/useful_functions.m index e675ba1..9127cab 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -287,7 +287,7 @@ function [dEdz_linear, dEdz_left, dEdz_right]=dEdz( atom_field_problem) endfor endfunction -function total_absorption=total_field_absorption(atom_field_problem) +function relative_transmission=total_relative_transmission(atom_field_problem) % summed across all frequencies field absorption coefficient global atom_properties; [dEdz_linear, dEdz_left, dEdz_right]=dEdz( atom_field_problem); @@ -300,11 +300,21 @@ function total_absorption=total_field_absorption(atom_field_problem) E_field.linear = atom_field_problem.E_field.linear; E_field.right = atom_field_problem.E_field.right; E_field.left = atom_field_problem.E_field.left; - total_absorption_v = ... - abs(E_field.linear).^2 - abs(E_field.linear - (1i)*dEdz_linear).^2 ... - +abs(E_field.right).^2 - abs(E_field.right - (1i)*dEdz_right).^2 ... - +abs(E_field.left).^2 - abs(E_field.left - (1i)*dEdz_left).^2; - total_absorption=sum(total_absorption_v); + + modulation_freq = atom_field_problem.modulation_freq ; + pos_freq_indexes=(modulation_freq>0); + transmited_intensities_vector = ... + abs(E_field.linear + (1i)*dEdz_linear).^2 ... + +abs(E_field.right + (1i)*dEdz_right).^2 ... + +abs(E_field.left + (1i)*dEdz_left).^2; + transmited_intensities_vector = transmited_intensities_vector(pos_freq_indexes); + input_intensities_vector = ... + abs(E_field.linear).^2 ... + +abs(E_field.right).^2 ... + +abs(E_field.left).^2; + input_intensities_vector = input_intensities_vector(pos_freq_indexes); + + relative_transmission = sum(transmited_intensities_vector) / sum(input_intensities_vector); endfunction % create full list of atom modulation frequencies from positive light frequency amplitudes |