基于魏格纳函数和焦散线方法的自加速光束matlab模拟与仿真

1.程序功能描述

2.测试软件版本以及运行结果展示

3.核心程序

4.本算法原理

5.完整程序

1.程序功能描述

基于魏格纳函数和焦散线方法的自加速光束matlab模拟与仿真。通过魏格纳函数法，来产生多种自加速的光束，设计自加速光束方法，模拟出相应的图结果。

2.测试软件版本以及运行结果展示

MATLAB2022A版本运行

3.核心程序

...................................................................
figure;
subplot(121);
imshow(Phase2,[])
colorbar 
title('Phase');
subplot(122);
imshow(Amp2,[])
colorbar 
title('amplititude');
%保存图片
Phase2_save = imresize(Phase2,[1080,1080],'nearest');
Amp2_save   = imresize(Amp2,[1080,1080],'nearest');
imwrite(Phase2_save,'Phase2_save.bmp','bmp'); 
imwrite(Amp2_save,'Amp2_save.bmp','bmp'); 
%%
%平面光场图,魏格纳函数
%W = Int(0..zmax){E(x+y)E(x-y)exp[2ixy]}
Wxkx    = func_wigner(X');
%光场
amp     = abs(sum(Wxkx,2));%根据魏格纳函数得到幅度，公式7
fai     = (sum([t].*Wxkx,2))./sum(Wxkx,2);%根据魏格纳函数得到相位，公式7
E       = amp'.*exp(sqrt(-1)*fai);
E       = abs(E);%根据相位和幅度得到E
E       = Zmax*[E/max(max(E))]';%E归一化

%计算光场
%根据公式11进行叠加
[s2,t2] = meshgrid(s,2*t);
%其实这里之前也没大问题，之前是用了这个sin代入了，其实只要把计算得到E代入就可以了，所以之前你才觉得去掉E还是不变这个问题
I       = abs(airy(0, E-s2,1)+airy(0,-E+s2,1));%非凸轨迹

figure;
subplot(211);
surface(t2,s2,I)  
shading interp
axis([min(t),max(t),-75,75]);
colormap(hot)
%%
%fit and xmax
%获得xmax值
[RR,CC] = size(I);
for i = 1:RR
    [VV,II]  = max(I(i,:));
    ysamp(i) = (II-RR/2)/RR*max(s);
end
ysamp=ysamp-mean(ysamp);
%进行拟合
xdata = 2*Z;
ydata = ysamp;
x     = lsqcurvefit(@myfun1,[3;1;1],xdata,ydata);%sin类型拟合
x
subplot(212);
plot(Z,x(1)*sin(x(2)*t+x(3)),'b');
hold on
plot(xdata(1:10:end),ysamp(1:10:end),'r.');
legend('拟合','主瓣位置');
axis([0,20,-8,8]);
grid on
xlabel('z/mm');
ylabel('x/um');
24_008m