diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2012-10-17 22:30:58 -0400 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2012-10-17 22:30:58 -0400 |
| commit | bc5f1eea9b5718c3e32d4cc0ea8b384983be6d93 (patch) | |
| tree | 641cc56ec1a0b119c637e5ccb5dd3917be5522bb | |
| parent | d30e6365e55c7b7cea59cd0b52c032b500e61029 (diff) | |
| download | wgmr-bc5f1eea9b5718c3e32d4cc0ea8b384983be6d93.tar.gz wgmr-bc5f1eea9b5718c3e32d4cc0ea8b384983be6d93.zip | |
split coupling function and its demo to separate scripts
| -rw-r--r-- | demo_prism_disk_coupling.m | 29 | ||||
| -rw-r--r-- | prism_disk_coupling.m (renamed from coupling_angles.m) | 29 |
2 files changed, 36 insertions, 22 deletions
diff --git a/demo_prism_disk_coupling.m b/demo_prism_disk_coupling.m new file mode 100644 index 0000000..e809968 --- /dev/null +++ b/demo_prism_disk_coupling.m @@ -0,0 +1,29 @@ +%% Calculates incident angle for proper coupling into the disc via prism + +% angle of the prism faces in degrees +prism_angle_in_degrees = 45; + +prism_angle = prism_angle_in_degrees*pi/180; + +%% prism index of refraction +% Rutile (TiO2) see http://refractiveindex.info/?group=CRYSTALS&material=TiO2 +n_rutile_o = 2.4885; % p - polarization +n_rutile_e = 2.75324; % s - polarization + +n_prism=n_rutile_o ; % for horizontal or p polarization + +%% disk material index of refraction +% Magnesium Fluoride (MgF2) see http://refractiveindex.info/?group=CRYSTALS&material=MgF2 +n_MgF2_o = 1.3751; +n_MgF2_e = 1.38679; + +% Measured +n_MgF2_p = 1.465; % p-polarization + + +n_disk=n_MgF2_p; +coupling_description='Rutile prism, MgF_{2} disk, p-polarization'; + +% calculate coupling angle and draw coupling solution with annotations +prism_disk_coupling(prism_angle_in_degrees, n_disk, n_prism, coupling_description); + diff --git a/coupling_angles.m b/prism_disk_coupling.m index c660a64..40e0980 100644 --- a/coupling_angles.m +++ b/prism_disk_coupling.m @@ -1,29 +1,12 @@ +function prism_disk_coupling(prism_angle_in_degrees, n_disk, n_prism, coupling_description) %% Calculates incident angle for proper coupling into the disc via prism - -% angle of the prism faces in degrees -prism_angle_in_degrees = 45; +% prism_angle_in_degrees - angle of the prism faces in degrees +% coupling_description - short annotation of the situation +% for example: +% coupling_description='Rutile prism, MgF_{2} disk, p-polarization'; prism_angle = prism_angle_in_degrees*pi/180; -%% prism index of refraction -% Rutile (TiO2) see http://refractiveindex.info/?group=CRYSTALS&material=TiO2 -n_rutile_o = 2.4885; % p - polarization -n_rutile_e = 2.75324; % s - polarization - -n_prism=n_rutile_o ; % for horizontal or p polarization - -%% disk material index of refraction -% Magnesium Fluoride (MgF2) see http://refractiveindex.info/?group=CRYSTALS&material=MgF2 -n_MgF2_o = 1.3751; -n_MgF2_e = 1.38679; - -% Measured -n_MgF2_p = 1.465; % p-polarization - - -coupling_description='Rutile prism, MgF_{2} disk, p-polarization'; -n_disk=n_MgF2_p; - %% critical angle for beam from prism to disk % recall n_d*sin(theta_disk)=n_prism*sin(theta_prism) where angles are counted from normal to the face % and we want theta_disk to be 90 degrees for the total internal reflection @@ -165,3 +148,5 @@ axis([-1,1,-0.5,1.5],'equal'); %% output of the plot to the file print('prism_disk_coupling.eps','-depsc2'); print('prism_disk_coupling.png'); + +end |