close all;clear all;clc;
x=-5:0.1:5;
y=-5:0.1:5;
[X,Y]=meshgrid(x,y);
K=2;
T=5;
a=exp(j*pi*Y.*(T-X));
b=pi*Y.*(T-abs(X))+eps;
c=((sin(b).*(T-abs(X)))./(T*b));
Z=a.*c;
subplot(1,2,2),surf(X,Y,Z),title('移相后'),xlabel('時間'),ylabel('頻率');
subplot(1,2,1),surf(X,Y,c),title('移相前'),xlabel('時間'),ylabel('頻率');