Matlab插值模拟在光的偏振中的应用
- 格式：doc
- 大小：90.50 KB
- 文档页数：5
Matlab 插值模拟在光的偏振中的应用
姓名：刘胜 学号：2013302540209
摘要：在光的偏振试验中，由于对数据处理比较多，所以，在查阅资料和在请教同学的基础上，自行编写了一组程序，来模拟光的偏振中的实验图谱。
关键词：matlab ，插值法,马吕斯定律，光的椭圆偏振。
正文：
实验内容及其步骤如下： 验证马吕斯定律
如果光源中的任一波列（用振动平面E 表示）投射在起偏器P 上（如下图）， 只有相当于它的成份之一的Ey （平行于光轴方向的矢量）能够通过，另一成份
Ex （=E cos θ）则被吸收。
与此类似，若投射在检偏器A 上的线偏振光的振幅为E0，则透过A 的振幅为E0 cos θ（这里θ是P 与A 偏振化方向之间的夹角）。
由于光强与振幅的平方成正比，可知透射光强I 随θ而变化的关系为
这就是马吕斯定律。
这是利用matlab 写的程序，运行结果如下： >> clear all;
>> x=[0:pi/9:2*pi];
>> y=[134 135 104 50 16 9 23 68 118 138 124 99 53 16 9 24 67 115 140]; >> xi=[0:pi/720:2*pi];
>> yi=interp1(x,y,xi,'spline'); >> plot(x,y,'*',xi,yi);
并且总结了几点心得： 1. 关于插值的应用。
clear all;
>> x=[0:pi/9:2*pi];
>> y=[134 135 104 50 16 9 23 68 118 138 124 99 53 16 9 24 67 115 140]; >> xi=[0:pi/180:2*pi];
>> yi=interp1(x,y,xi,'spline'); >> plot(x,y,'*',xi,yi);
注意:1.矩阵坐标只能用弧度表示；
2. 对应矩阵的数目一定要对称；interp 函数的调用格式yi=interp1(x,Y ,xi);yi=interp1(Y,xi);
3. yi=interp1(x,y,xi,method);
4. 利用星号在plot 中做数据的标记，spline 是三次条样插值。
结果如下：
4. 波片
(1) 预备知识1：两个互相垂直的、同频率的简谐振动的合成 两简谐振动的位移方程分别为
)cos()
cos(2211ϕωϕω+=+=t A y t A x (3)
合成的运动是两个分运动之和，消去参数t ，得到轨迹方程
)(sin )cos(212212212
2
2212ϕϕϕϕ-=--+A A xy A y A x (4)
上式表明，一般情况下，合振动矢量末端运动轨迹是椭圆，该椭圆在2122A A ⨯的矩形范围内。
如果(3)式是两线偏振光，则叠加后一般成为椭圆偏振光。
下面讨论相位差12ϕϕϕ-=∆为几种特殊值的价况。
① πϕk 2=∆ （k =0, ±1, ±2, …），(4)式变为
21A y
A x =
(5)
合振动矢量末端运动轨迹是上述矩形的对角线，与y 轴的夹角为α（请想一想：α＝？）。
这就是说，相互垂直的、同相位的两线偏振光合成后仍为线偏振光，但E 振动有新的取向。
② πϕ)12(+=∆k （k =0, ±1, ±2, …），(4)式变为
21A y
A x -= (6)
(6)式表明，相互垂直的、相位相反的两线偏振光合成后也是线偏振光，E 振动方向与y 轴的夹角也是α，但在y 轴的另一侧。
③
2)
12(π
ϕ+=∆k （k =0, ±1, ±2, …），(4)式变为
122
2
2
12=+A y A x (7)
此时，合振动矢量末端运动轨迹是正立在上述矩形内的椭圆。
即相互垂直的、相位差是2π的两线偏振光合成为椭圆偏振光，椭圆的长半轴和短半轴分别等于两线偏振光的振幅。
④ 根据情况③且A A A ==21，(4)式变为
222A y x =+ (8)
即振幅相等的、相互垂直的、相位差是2π的两线偏振光合成为圆偏振光。
下面是实验后的数据处理：
>> clear all;
>> a=[0:pi/9:2*pi];
>> m=[59 75 84 74 54 37 19 21 35 52 68 82 75 55 33 21 21 35 58]; >> b=ones(size(m)); >> ix=87; >> iy=17.5; >> for t=1:19; k=1;
E(k)=sqrt((ix*iy)/(ix+iy-m(t))); b(t)=E(k); k=k+1; end
>> polar(a,b); >>
270
0刘胜 2013302540209
2700
刘胜 2013302540209
>> clear all;
>> x=[0:pi/9:2*pi];
>> m=[59 75 84 74 54 37 19 21 35 52 68 82 75 55 33 21 21 35 58]; >> y=ones(size(m));
>> ix=87;
>> iy=17.5;
>> for t=1:19;
k=1;
E(k)=sqrt((ix*iy)/(ix+iy-m(t)));
y(t)=E(k);
k=k+1;
end
>> ai=[0:pi/180:2*pi];
>> bi=interp1(x,y,ai,'spline');
>> polar(x,y,'*');
>> hold on;
>> polar(ai,bi);
>> hold off;
90
270
注意：
（1）hold on可以将两个图像相结合。
（2）极坐标条件下最好给出范围，不能像plot的一次性描点，for循环蛮好用。
（3）数组的替换注意。
总结：经过这次试验让我更加深刻认识物理实验的乐趣和提高了我的动手能力，以后的实验中我也会好好应用学过的知识来学习。
参考文献：《matlab仿真技术与应用》《matlab图形处理》《matlab数学计算范例教程》。
