close all;clear all;clc;
T1=5;
P=13;

x=-5: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);
m=0:(P-1);
[X,Y,M]=meshgrid(x1,y1,m);

a=exp(j*pi*Y.*(T1-(X-M*T1)));
b=pi*Y.*(T1-abs(X-M*T1))+eps;
c=(sin(b)./b).*(1-abs(X-M*T1)/T1)/sqrt(P);

    d=[1,1,1,1,1,0,0,1,1,0,1,0,1];
    for k=0:(P-1-m)
    e=d(k+1)*d(k+1+m)*exp(j*2*pi*T1);
    Z(k+1)=a(k+1)*c(k+1)*e
    end

surf(X,Y,Z),title('線性調頻信號模糊圖'),xlabel('時間'),ylabel('頻率')