基于基尔霍夫衍射理论编写matlab程序计算波带片聚焦,但为什么我感觉程序是对的,结果算出来是错的呢?请问有人知道为啥吗
%% I. 清理环境变量
%%
clc
clear all;
%% II. 绘制zp
lambda=0.0028; %入射波长,单位微米, 632.8nm
k0=2*pi/lambda; %波数
focallength=50000; %波带片焦长,单位微米
deltar=15.0; %最外环宽度,最小特征尺寸,微米
Number=floor(focallength*lambda/(4.0*deltar*deltar)); %波带片的总环数 floor函数即向下取整函数
Radius=sqrt(Number*lambda*focallength); %波带片半径,单位微米 4.2093e+03
xp1=linspace(-100,100,400)