Home

Introduction

Contents

Scripts

FeedBack

 Date 06/09/2010
Script S3_2_13A.m
Download Script S3_2_13A.m
%===============================================
%A quarter-wave plate 1A
%transmission axis of the first polarizer
%parallel to the x axis
%===============================================
%
%
%beam emerging from the first polarizer
v1=[1;0]
%matrix of the quarter-wave whose optic axis
%is parallel to x axis
m1=[1 0;0 -i]
%polarization form of the beam emerging
%from quarter-wave plate
vs=m1*v1
%
%intensity Is of the beam emerging
%from the quarter-wave
vst=vs'
Is=vst*vs
%
%matrix of the second rotating polarizer
max=25;
tetag=linspace(0,360,max)
teta=tetag*pi/180;
for i=1:max
%elements of the matrix
c1=cos(teta(i));
s1=sin(teta(i));
a11=c1^2;
a12=c1*s1;
a21=c1*s1;
a22=s1^2;
m2=[a11 a12;a21 a22];
%beam emerging from P2
vf=m2*m1*v1;
%the complex conjugate
vft=vf';
%intensity
Int(i)=vft*vf;
end
%array of the intensities
Int
%the corresponding plot
plot(tetag,Int,'r-o'),grid on
title('Intensity of the beam emerging from P2 varying tetag')
%===============================================
%
Top
Script S3_2_13B.m
Download Script S3_2_13B.m
%===============================================
%A quarter-wave plate 1B
%transmission axis of the first polarizer
%subtends an angle of 45° with x axis
%===============================================
%
%A beam of natural light is incident
%on the first polarizer
%
%vector of the first polarizer whose
%transmission axis subtends
%an angle of 45° with the x axis
tetg=45
tet=tetg*pi/180;
cc1=cos(tet);
ss1=sin(tet);
v1=[cc1;ss1]
%matrix of the quarter-wave whose optic axis
%is parallel to x axis
m1=[1 0;0 -i]
%polarizzation form of the beam
%emerging from quarter-wave plate
vs=m1*v1
%intensity of the beam
%emerging from the quarter_wave
vst=vs'
Is=vst*vs
%
%matrix of the rotating polarizer P2
max=25;
tetag=linspace(0,360,max)
teta=tetag*pi/180;
for i=1:max
%elements of the matrix
c1=cos(teta(i));
s1=sin(teta(i));
a11=c1^2;
a12=c1*s1;
a21=c1*s1;
a22=s1^2;
m2=[a11 a12;a21 a22];
%beam emerging from P2
vf=m2*m1*v1;
%complex conjugate
vft=vf';
%intensity for a given angle
Int(i)=vft*vf;
end
%array of the intensities
Int
%plot of the intensity emerging from P2 varying the angle of rotation
plot(tetag,Int,'r-o'),grid on
title('Intensity of the beam emerging from P2 varying tetag')
% axis([0,360,0.4,0.58])
%===============================================
%
Top