close all;clear all;clc;
T=5;
P=13;
x=:0.08:5;
y=5:0.08: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,1,-1,-1,1,1,-1,1,-1,1];
%%%%%%%%%%%%%%%%%%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('頻率')