close all;clear all;clc;
T=5;
P=9;
x=-5:0.1:5;
y=-5:0.1:5;
%x1=x(1):((x(end)-x(1))/(P-1)):x(end);
%y1=y(1):((y(end)-y(1))/(P-1)):y(end);
[X,Y,]=meshgrid(x,y);
d=[1,1,1,1,exp(2*pi*j/3),exp(-2*pi*j/3),1,exp(-2*pi*j/3),exp(2*pi*j/3)];
%%%%%%%%%%%%%%%%%%kai_1%%%%%%%%%%%%%%%%%%
%a=exp(j*pi*Y.*(T-X));
%b=pi*Y.*(T-abs(X))+eps;
%c=((sin(b).*(T-abs(X)))./(T*b));
%kai_1=a.*c;
%%%%%%%%%%%%%%%%%%kai_2%%%%%%%%%%%%%%%%%%
%kai_2=zeros(P);
%d=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1];
%for k=1:P
%    for l=1:P
%        kai_2=kai_2+(d(k)*d(l)*exp(j*2*pi*Y*k*T).*exp(j*2*pi*Y.*(X-(l-k)*T))*delta(X,X-(l-k)*T))/P;      
%    end
%end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Z=zeros(max(size(X)));
for m=-P:P
    Z=Z+kai_1(X-m*T,Y,T).*kai_2(m*T,Y,T,P,d);
end
surf(X,Y,Z),title('相位編碼信號模糊圖'),xlabel('時間'),ylabel('頻率')