一种混合多高斯相位的理论与模拟仿真-含Matlab源代码

一、混合多高斯相位模拟结果

写在最前面：仿真结果如图1与图2所示:

图1 混合多高斯相位模拟演示一

图2 混合多高斯相位模拟演示二

二、混合多高斯相位理论公式

$$P\text{hase}=A_1{e}^{-\frac{1}{2{{\sigma }_1}^2}\left[{\left(x-u_1\right)}^2+{\left(y-v_1\right)}^2\right]}+A_2{e}^{-\frac{1}{2{{\sigma }_2}^2}\left[{\left(x-u_2\right)}^2+{\left(y-v_2\right)}^2\right]}+A_3{e}^{-\frac{1}{2{{\sigma }_3}^2}\left[{\left(x-u_3\right)}^2+{\left(y-v_3\right)}^2\right]}$$

其中，A控制高斯相位的深度，即大小；σ控制高斯相位的束腰大小，u与v分别控制高斯相位在图像矩阵中的位置。

三、Matlab源代码

clear;
close all;
clc;
%% 混合多高斯相位---图1

N=1024;     %取样数
h=0.532e-3; %波长(mm)
L0=26;      %模拟图像宽度(mm)

n=1:N;                    
x=-L0/2+L0/N*(n-1);	    		
y=x;
[xx,yy] = meshgrid(x,y);  % 物平面坐标网格

dh=20*exp(-90*((xx-5).^2+(yy+2).^2)/L0/L0)+...
    28*exp(-90*((xx+8).^2+(yy-5).^2)/L0/L0)+...
    30*exp(-40*((xx+2).^2+(yy+2).^2)/L0/L0)+...
    30*exp(-120*((xx-8).^2+(yy-8).^2)/L0/L0);
figure,mesh(dh)
figure,imshow(dh,[])

%% 混合多高斯相位---图2

N=1024;     %取样数
h=0.532e-3; %波长(mm)
L0=26;      %模拟图像宽度(mm)

n=1:N;                    
x=-L0/2+L0/N*(n-1);	    		
y=x;
[xx,yy] = meshgrid(x,y);  % 物平面坐标网格

dh=20*exp(-90*((xx-5).^2+(yy+2).^2)/L0/L0)+...
    28*exp(-90*((xx+8).^2+(yy-5).^2)/L0/L0)+...
    30*exp(-40*((xx+2).^2+(yy+2).^2)/L0/L0)+...
    30*exp(-120*((xx-8).^2+(yy-8).^2)/L0/L0)+...
    20*exp(-60*((xx+8).^2+(yy+8).^2)/L0/L0);
figure,mesh(dh)
figure,imshow(dh,[])

---

⭐️◎⭐️◎⭐️◎⭐️ · · · **博 主 简 介** · · · ⭐️◎⭐️◎⭐️◎⭐️ ♪

▁▂▃▅▆▇ 博士研究生生 ，研究方向主要涉及定量相位成像领域，具体包括干涉相位成像技术(如**全息干涉☑**、散斑干涉☑等)、非干涉法相位成像技术(如波前传感技术☑，相位恢复技术☑)、此外，还对各种相位解包裹算法☑，相干噪声去除算法☑等开展过深入的研究。

---