乘性散斑噪声理论及其仿真实例分析-Matlab源代码

一、乘性散斑噪声理论分析与仿真模拟结果

1.1 乘性散斑噪声理论

噪声是影响激光相干场成像像质的重要因素之一，噪声来源主要有：与激光回波信号相关的乘性散斑噪声，与激光回波信号无关的加性背景噪声和探测器噪声等。

其乘性散斑噪声可用下式来表示：

$$I(x_0,y_0)=O(x_0,y_0)A(x_0,y_0)=o_i\exp \left(i{\phi }_i\right)\cdot a_k\exp \left(i{\phi }_k\right)$$

其中，$O(x_0,y_0)$为理想物光场复振幅，$A(x_0,y_0)$为散斑噪声复振幅。$o_i$与${\phi }_i$分别表示物光场的振幅与相位。$a_k$和${\phi }_k$分别为散斑噪声的振幅和相位。

1.2 仿真模拟结果

根据上述公式，可对乘性噪声进行仿真模拟，如图所示。为了进行对比说明，本文还同时给出了加性散斑噪声的模拟，也如下图所示。

二、Matlab源程序

2.1 加性散斑噪声模拟Matlab源程序

加性散斑噪声仿真模拟Matlab源程序如下：

clc
clear
close all 

%% 加性散斑噪声模拟
ph=peaks(1024)*5;                 %预设原始相位
U=exp(j*ph); 
figure,imshow(angle(U),[])    % 无噪声包裹相位

mu=0.8;                          %噪声系数为0；此时相当于不加入噪声
noise=rands(1024,1024).*mu.*pi;   %噪声
ph=ph+noise;                      %加性噪声

U=exp(j*ph); 
ph0=angle(U);                     %对应的包裹相位
figure,
imshow(ph0,[]);
title('加性散斑噪声')

运行后，可得到一下结果：

2.2 乘性散斑噪声模拟Matlab源程序

乘性散斑噪声仿真模拟Matlab源程序如下：

clc
clear
close all 

%% 乘性散斑噪声模拟
N=1024;
ph=peaks(N)*5;                 %预设原始相位
U=exp(j*ph); 
figure,imshow(angle(U),[])    % 无噪声包裹相位

mu=0.8;                          %噪声系数为0；此时相当于不加入噪声
fai=rand(N,N).*mu.*pi;    %随机相位值
A0=rand(N,N);          %表面振幅分布的随机值
Noise=A0.*exp(j*fai);
figure,imshow(angle(Noise),[])    % 散斑噪声


% 乘性散斑噪声
ph_noise=Noise.*U;
ph_noise__2=angle(ph_noise);                     
figure,
imshow(ph_noise__2,[]);
title('乘性散斑噪声')

---

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

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