eeddabcd-CSDN博客

原创 matlab仿真双缝衍射强度分布曲线

%page 182%figure_5_33%双缝衍射强度分布曲线I_0=10;a=1;d=3;k=6;f=5;x=-10:0.05:10;alpha=0.5*k*a*x/f;delta=k*d*x/f;I=4*I_0*( ( sin(alpha) ).^2./alpha.^2 ).*(cos(delta/2)).^2;plot(x,I) %实线为衍射因子与干涉因子共同作用后的实际曲线hold onI_1=4*I_0*( sin(alpha) ).^2./alpha.^.

2021-11-07 14:20:38 1236

原创矩孔衍射在x轴上的强度分布曲线

%page 167%figure_5_13I_0=10;a=-4*pi:0.05*pi:4*pi;I=I_0*(sin(a)./a).^2;plot(a,I);title('矩孔衍射在x轴上的强度分布曲线')xlabel('a')ylabel('intensity')

2021-11-07 12:39:53 255

原创基于matlab的等光程差面

%page 77%等光程差面%Equal optical path difference surfacelamda=630*1e-9;d=0.5*1e-3;m=3;[x,y]=meshgrid([-10:0.15:10],[-10:0.15:10]);u=(x.^2+y.^2)./( (d/2)^2-(m*lamda/2).^2 )+1;z=(u.*(m*lamda/2).^2).^0.5;mesh(x,y,z);%colormap autumnhold onn=5;[x.

2021-11-07 11:50:10 371

原创 matlab仿真光拍合成

不同频率的两个单色光波叠加当两个光波的角速度w很接近时，二者叠加结果会出现“光学拍”现象利用matlab仿真光拍传播%light beatw1=10;w2=9.8;k1=10;k2=9.5;w_b=0.5*(w1+w2);w_m=0.5*(w1-w2);k_b=0.5*(k1+k2);k_m=0.5*(k1-k2);a=12;z=0:0.05:30;for t=0:200 E=2*a*cos(k_m*z-w_m*t).*cos(k_b*z-w_b*t); p

2021-11-07 11:47:25 630

原创薄膜光学中单层膜介质膜反射率随其光学厚度的变化（利用matlab仿真）

%Variation of reflectivity of dielectric film with its optical thicknesszeta_0=0;zeta=zeta_0;n_0=1;n_g=1.5;lamda=632.8*1e-9; %wavelength:632.8nmh=0:1e-9:lamda;%n=[1.0 1.2 1.4 1.5 1.7 2.0 3.0];n_1=1.0;n_2=1.2;n_3=1.4;n_4=1.5;n_5=1.7;n_6=2.0;n_7=3.0.

2021-10-17 12:50:03 1739

原创 matlab仿真——不同反射率下透射光条纹的强度分布曲线

对于多光束干涉情况下的透射光，当反射率R很大时，条纹对比度很好；而反射率R较小时，条纹对比度很差。根据精细度系数F以及透射光强与入射光强比值，绘制出不同反射率下的透射光条纹的强度分布曲线。phase=pi:0.01:2*(1+1)*pi+pi;R=[0.046 0.138 0.27 0.64 0.87];F=4*R./((1-R).^2);plot([2*pi,2*pi],[0,1],'k')hold onplot([4*pi,4*pi],[0,1],'k')hold onIt_Ii_1

2021-10-09 23:59:06 1579

原创 matlab拟合杨氏双缝干涉光强分布曲线

I_0为光强极值，d为双缝间距，lamda为波长，D1为双缝所在平面到光屏距离。%Intensiy distributionI_0=10;d=0.02;D1=50;lamda=5e-7;x=-67:67;k=pi*d/(lamda*D1);I=4*I_0*cosd(k*x).*cosd(k*x)plot([0,0],[0,45],'k')hold onplot(x,I,'b')hold offtext(40,42,'光强分布曲线I')text(2,43,'y轴')...

2021-10-06 23:11:39 1374

原创利用matlab实现驻波仿真

%stationary waveA=6;k=2;w=0.05;x=0:0.01:15;u_0=pi/8;for t=0:250 E_1=A*cos(k*x+w*t) E_2=A*cos(k*x-w*t+u_0) E=E_1+E_2 plot(x,E_1,'r') hold on plot(x,E_2,'c') hold on plot(x,E,'k') hold off ylim([-15,15]) titl.

2021-10-06 00:03:12 3711 1

原创利用matlab表示全反射下的s波与p波的位相跃变曲线

a_1=40:0.01:90;n=1/1.5;c_s=2*atand((-1)*(sind(a_1).^2.-n^2).^0.5./cosd(a_1))c_p=2*atand((-1)*(sind(a_1).^2.-n^2).^0.5./(n^2*cosd(a_1)))%axis ijplot(a_1,c_p)hold onplot(a_1,c_s)hold offgrid ontitle('s波与p波全反射下的位相跃变')xlabel('a_1角度变化')ylabel('s波.

2021-10-05 23:17:25 849

原创 matlab的rotate函数的使用

利用rotate函数可以实现定轴旋转rotate(vector,[axis],angle)函数括号内，第一项为矢量，就是要旋转的矢量；第二项为坐标轴，在空间选定某坐标轴后，矢量绕此坐标轴顺时针旋转；第三项为旋转的角度。...

2021-10-05 18:28:45 8190

原创利用matlab拟合光波的p波与s波反射系数与透射系数

入射光的电矢量总是可以分解为垂直于入射面的s波与平行于入射面的p波，求得s波与p波的反射系数及透射系数后，可以用matlab拟合出系数曲线。%figure 1.23 calculation of rs rp ts tpu_1=0:0.5:90;u_2=asind(sind(u_1)/1.5);r_s=sind(u_2-u_1)./sind(u_2+u_1);t_s=2*sind(u_2).*cosd(u_1)./sind(u_1+u_2);r_p=tand(u_1-u_2)./tand(u_

2021-10-05 00:49:42 2020

原创利用matlab实现光波（电磁波）的仿真

根据电磁波余弦曲线公式用matlab进行三维做图，做出电磁波仿真。%Electromagnetic wave simulation%假定电磁波沿x轴方向传播，电磁波在y轴方向和z轴方向有分量Ey_m=6; %电磁波y分量幅值Ez_m=8; %电磁波z分量幅值E=(Ey_m^2+Ez_m^2)^0.5; %求出电磁波的总幅值w=10; %电磁波的角速度为10a_1=pi/4; %y分量电磁波的初相位a_2=pi/2; %z分量电磁波的初相位%a_1和a_2的大小差值，决...

2021-10-05 00:25:38 7002 1

m0_51282055的博客