请写一个MATLAB程序，采用傅立叶变换空间载波法从强度分布得到相位分布，即从空间上呈正弦分布的光强信息，恢复出波面的相位信息

以下是一个简单的MATLAB程序，用傅立叶变换空间载波法从强度分布得到相位分布：

% 读入干涉强度分布图
I = imread('intensity.bmp');
I = double(I);

% 对干涉图进行傅里叶变换
F = fft2(I);

% 计算干涉图的大小
[M, N] = size(I);

% 设定载波频率
fx = 0.1;
fy = 0.1;

% 生成傅里叶空间载波
[X, Y] = meshgrid(1:N, 1:M);
carrier = exp(-1i*2*pi*(fx*X+fy*Y));

% 对频域图像进行滤波
F_filtered = F.*carrier;

% 对滤波后的频域图像进行反傅里叶变换
I_filtered = real(ifft2(F_filtered));

% 相位解算
phase = atan2(imag(I_filtered), real(I_filtered));

% 相位校正
phase = phase - median(phase(:));

% 显示结果
figure;imshow(phase,[]);colormap(jet);colorbar;

需要注意的是，以上程序中的载波频率需要根据具体的实验条件进行调整，以保证解相位的精度和可靠性。此外，在使用傅立叶变换空间载波法进行相位解算时，还需要注意干涉图的稳定性和信噪比，以避免解相位误差的产生。