%?*******************************************
clear?all
close?all
clc
%%?*************初始相位**************
N?=?512;
G?=?2;
phi0?=?peaks（N）*G;?%模擬初始相位
figure（1）
surf（phi0‘FaceColor‘‘interp‘?‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight?left?axis?tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘Initial?Phase‘‘FontSize‘14）
set（figure（1）‘name‘‘Initial?Phase?3D‘‘Numbertitle‘‘off‘）;
phi?=?angle（exp（j*phi0））;?????????%包裹相位
figure（2）;
imshow（phi[]）;
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;%title（‘Wrapped?Phase‘‘FontSize‘14）
set（figure（2）‘name‘‘Wrapped?Phase‘‘Numbertitle‘‘off‘）;
axis?on
%%?*************相位解包裹**************
[mn]?=?size（phi）;
phidx=zeros（mn）;
phidy=zeros（mn）;
phidx（1:m-1:）=?angle（exp（j*（phi（2:m:）-phi（1:m-1:））））;
phidy（:1:n-1）=?angle（exp（j*（phi（:2:n）-phi（:1:n-1））））;
%********************對包裹相位求二階偏微分**************
Rou3?=?zeros（mn）;
Rou3dx?=?zeros（mn）;
Rou3dy?=?zeros（mn）;
Rou3dx（1:m-1:）?=?phidx（2:m:）-phidx（1:m-1:）;
Rou3dy（:1:n-1）?=?phidy（:2:n）-phidy（:1:n-1）;
Rou3?=?Rou3dx?+?Rou3dy;
figure（3）;
surf（Rou3‘FaceColor‘‘interp‘?‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight?left?axis?tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘lou3‘‘FontSize‘14）
set（figure（3）‘name‘‘R（xy）?3D‘‘Numbertitle‘‘off‘）;
figure（4）;
imshow（Rou3[]）;
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;
set（figure（4）‘name‘‘R（xy）?2D‘‘Numbertitle‘‘off‘）;
%%?***********************DCT求解泊松方程********************
tic
PP3?=?dct2（Rou3）;
for?ii=1:m
????for?jj=1:n
????????k1=2*cos（（ii-1）*pi/（m））;
????????k2=2*cos（（jj-1）*pi/（n））;
????????KK?=?k1+k2-4;
????????PH3（iijj）?=?PP3（iijj）/KK;
????end
end
PH3（11）?=?-（PH3（12）?+?PH3（21）?-?PP3（11））/2;
phi3?=?idct2（PH3）;
toc
phi3?=?phi3（1:m1:n）;?%解包裹出的相位


figure（5）;
surf（phi3‘FaceColor‘‘interp‘?‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight?left?axis?tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘BLS?Phase?Unwrapping‘‘FontSize‘14）
set（figure（5）‘name‘‘BLS?Phase?Unwrapping‘‘Numbertitle‘‘off‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2360??2018-04-16?15:47??zuixiaoerchengfa.m