matlab给出光学传输矩阵T，如何生成一张灰度图片对应的散斑图？

生成散斑图需要用到传输矩阵T和一个随机的相位场。下面是一份MATLAB代码示例：

% 首先定义传输矩阵T和图像大小
N = 256; % 图像大小
z = 0.1; % 传输距离
lambda = 0.6328e-6; % 波长
k = 2*pi/lambda; % 波数

% 生成传输矩阵T
L = 0.01; % 采样间隔
x = (-N/2:N/2-1)*L; % x坐标
[X,Y] = meshgrid(x,x); % 生成网格
r = sqrt(X.^2 + Y.^2); % 极径
theta = atan2(Y,X); % 极角
T = exp(1i*k*z) ./ (1i*lambda*z) .* exp(1i*k/2/z*r.^2) .* exp(-1i*k*z) .* exp(-1i*k/2/z*r.^2);

% 生成随机相位场
phi = rand(N);

% 计算散斑图
E_in = ones(N); % 入射光强度为常数
E = fftshift(fft2(ifftshift(sqrt(E_in).*exp(1i*phi)))); % 傅里叶变换
E_out = ifftshift(ifft2(fftshift(T.*E))); % 逆傅里叶变换
I = abs(E_out).^2; % 散斑图强度

% 显示散斑图
figure;
imshow(I, []);

这段代码首先定义了传输矩阵T和图像大小，然后根据传输距离、波长等参数生成传输矩阵T。接着生成一个随机的相位场，将其与入射光强度相乘后进行傅里叶变换，再与传输矩阵T相乘并进行逆傅里叶变换得到输出光强度分布，最后计算强度并显示散斑图。