clear;
close?all;
clc;
t=2;???????????????????????????????????%拓撲荷
n=256;?????????????????????????????????%分辨率
a=20;
r=-a:2*a/（n-1）:a;??????????????????????%光束半徑
[x?y]=meshgrid（r）;?????????????????????%建立網格
y=-y;?????
%================渦旋相位===========================
txc=y./x;
%?ttx=atan（txc）;?????????????????????????%theta
%?ttx（:1:n/2）=ttx（:1:n/2）+pi;??????????
%?ttx（n/2+1:nn/2+1:n）=ttx（n/2+1:nn/2+1:n）+2*pi;

%---------------測試效果---------------------------
%?ttx（1:n/21:n/2）=ttx（1:n/21:n/2）+pi;
%?%ttx（n/2+1:nn/2+1:n）=ttx（n/2+1:nn/2+1:n）+2*pi;
%?ttx（n/2+1:n:）=pi/2;



%=====================================================
rr=2*（x.^2+y.^2）;
gaosi=exp（-rr）;????????????????????????%高斯函數
uf1=fftshift（fft2（gaosi））;
imshow（abs（uf1）[]）;
%?mesh（abs（uf1））;???????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2067??2015-12-14?21:53??phase.m