相位重建GS算法

function [ I ] = F_GS（ I1I2k ）
%UNtitleD4 Summary of this function goes here
%   Detailed explanation goes here
%   I1为输入光场的振幅
%   I2为输出光场的振幅
%   k为算法程序的迭代次数
%   I为输出的光场，为复数，包括相位和振幅

[mn]=size（I1）;
I1=double（I1）;
I2=double（I2）;

I1re=rand（mn）;
I1lm=rand（mn）;
I1x=I1re+i*I1lm;
I1x=I1x./abs（I1x）;%随机生成输入面的初始相位

I2q=I2.^2;
sumG=sum（sum（I2q））;
sse=zeros（1k）;
l=0;

for l=1:k