多频外差

% 程序开始
clc;
close all;
clear;

% 图片的初始化
width = 1280;  
heigth = 800;

% 三频率
freq = [70 64 59];


C = cell(3,4);  
for i=1:3
    for j=1:4
        C{
   i,j} = zeros(heigth,width);
    end
end

% 利用余弦函数计算12张图的灰度值
% 图像的生成
% 三种频率，四组相位
for i = 1:3 % 对应三种不同的频率
    for  j = 0:3 % 对应四种相位
        for k = 1:width
            C{
   i,j+1}(:,k) = 128+127*sin(2*pi*k*freq(i)/width+j*pi/2);
        end
        
    end
end

% 对灰度值进行归一化处理
for i = 1:3
    for j = 1:4
        C{
   i,j} = mat2gray(C{
   i,j});
    end
end

% 显示12张图
for i = 1:3
    for j = 1:4
        n = 4*(i-1)+j;
        h = figure(n);
        imshow(C{
   i,j});
    end
end

% 初始化三组处理后的图片灰度矩阵
% phi也是分块矩阵
% 存储相位主值图像
phi = cell(3,1);
for i = 1:3
    phi{
   i