|
|
[交流]
把单缝改多缝，缝之间距离相等
clear all;clc;ticdz=0.5e-4; %%define the distance step, unit [m]maxz=0.5e4; %%%%%%%%%%%%%%%%%%%%define the space % grid parameters, the unit is in meterLx=0.6e-2;N=1024; x=-Lx/2:Lx/N:Lx/2-Lx/N; dx=Lx/N; kx=[0:N/2-1 -N/2:-1]*2*pi/Lx;lambda=632.8e-9;k=2*pi/lambda;i=sqrt(-1);w0=0.5e-3;A0=1;E=zeros(1,N);phi=zeros(1,N);N1=fix((-w0/2-(-Lx/2))/dx);N2=fix((+w0/2-(-Lx/2))/dx);E(N1:N2)=A0.*ones(1,N2-N1+1);D=exp(i/(2*k).*(-kx.^2).*dz/2);%plot(x,E)temp=maxz./200; gg=1;for ii=1:maxz %%%%%%%%%%%fft processing UD_h_2_U=fft(E).*D; uD_h_2_U=ifft(UD_h_2_U); UDD_h_U=fft(uD_h_2_U).*D; E=ifft(UDD_h_U); % Ec(ii)=E(N/2);% z(ii)=ii*dz; if mod(ii,temp)==0 IE=sum(abs(E).^2)*dx waveE(gg, =E; zz(gg)=ii*dz; gg=gg+1 endendfigure;mesh(x,zz,abs(waveE).^2);view(0,90);axis tight;
发自小木虫Android客户端 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
|
回复此楼