[t,W]=meshgrid([2:0.2:7],[0:pi/6:3*pi]); %设置时-频相平面网格点
Gs1=(1/(sqrt(2*pi)*a))*exp(-0.5*abs((t1-t)/a).^m).*exp(-i*W*t1);
Gs2=(1/(sqrt(2*pi)*a))*exp(-0.5*abs((t2-t)/a).^m).*exp(-i*W*t2);
Gs=Gs1+Gs2;
subplot(2,3,1);
%绘制Gs的实部三维网格立体图
mesh(t,W/pi,real(Gs));
axis([2 7 0 3 -1/(sqrt(2*pi)*a) 1/(sqrt(2*pi)*a)]);