工光大作业1.docx

综合练习大作业一一、模拟两个频率相同、振动方向不同的单色光波的叠加1.基本分析 对于两个频率相同、振动方向不同的单色光波有：其合振动的大小和方向随时间变化，合成后为椭圆偏振光。2.MATLAB程序：下面采用0.5微米波长光波，振幅均为4，振动方向相互垂直的两束光波进行仿真分析，并列出不同相位差下的不同情况，得到二维和三维曲线。clear all; %清除变量wave=0.5e-6; %定义波长c=3e+8; %定义光速T=wave/c; %计算周期k=2*pi/wave; %计算波数w=2*pi/T; %计算角频率Ax=4;Ay=4; %定义振幅t=linspace(0,T,100); %分序列以绘制曲线z=linspace(0,4,100); i=1;for phase=0:pi/4:7*pi/4 %定义phase为x、y轴八种不同相位差Ex=Ax*cos(w*t-k*z); %光波波矢表达式Ey=Ay*cos(w*t-k*z+phase); subplot(4,4,i); %生成4*4个绘图窗口并在第i块绘制坐标系p=i/4;i=i+1;plot3(z,Ex,Ey); %以z Ex Ey为坐标轴绘制三维曲线zlabel(z); %标识z轴xlabel(x); %标识x轴ylabel(y); %标识y轴title([相位差= num2str(p) π]) %为三维坐标轴添加标题end n=9; for phase=0:pi/4:7*pi/4 %定义phase为x、y轴八种不同相位差Ex=Ax*cos(w*t-k*z); %光波波矢表达式Ey=Ay*cos(w*t-k*z+phase);subplot(4,4,n); %生成4*4个绘图窗口并在第n块绘制坐标系m=(n-9)/4;n=n+1; plot(Ex,Ey); %以Ex Ey为坐标轴绘制二维曲线xlabel(x); %标识x轴ylabel(y); %标识y轴title([相位差= num2str(m) π]) end3.程序运行结果二、模拟两个频率不同、振动方向相同的单色光波的叠加1.基本分析对于两个频率不同但接近、振动方向相同。振幅相同的单色光波的叠加会形成光学拍。沿用（一）中采用的光波公式，有：2.MATLAB程序下面采用频率为2000000000和2050000000，振幅均为1，振动方向相同的两束光波进行仿真分析，得到原始波、合成波和振幅、强度曲线。clear all; %清除变量c=3e+8; %定义光速f1=2000000000; %定义两个不同的频率f2=2050000000;w1=2*pi*f1; %计算角频率w2=2*pi*f2;k1=w1/c; %计算波数k2=w2/c;a=1; %定义振幅 ka=(k1+k2)/2;km=(k1-k2)/2;wa=(w1+w2)/2;wm=(w1-w2)/2;z=0:0.001:20; %分序列以绘制曲线t=0; %取出t=0时图像for i=1:1:20001 %循环 计算叠加公式结果E1(i)=a*cos(k1*z(i)-w1*t); %光波1E2(i)=a*cos(k2*z(i)-w2*t); %光波2E(i)=E1(i)+E2(i); %合成光波A(i)=2*a*cos(km*z(i)-wm*t); %合成光振幅I(i)=A(i)^2; %合成光光强endhold on;figure(1)plot(z,E1,r); %光波1图像，红色plot(z,E2,k); %光波2图像，黑色title(光波1（红色） 光波2（黑色）);figure(