clear all; clc;close all; fs=1000; t=0:1/fs:0.5-1/fs; f1=200; K=50; s=exp(1j*(2*pi*f1*t+pi*K*power(t,2))); s1=s; l=length(s); s=[zeros(1,l-1) s zeros(1,l-1)]; y=zeros(2*l-1,51); for tao=1:2*l-1 ks=1; for ksi=1:0.2:101 y(tao,ks)=abs(s1.*conj(s(tao:l+tao-1))*exp(-1j*2*pi*(ksi-51)*t)'); ks=ks+1; end end figure mesh(-50:0.2:50,(-l+1)/fs:1/fs:(l-1)/fs,y);title('模糊函数') figure plot(-50:0.2:50,y(l,:));title('tao=0') 

figure plot((-l+1)/fs:1/fs:(l-1)/fs,y(:,ks/2));title('ksi=0') y=y/max(max(y));%归一化 figure 

mesh(-50:0.2:50,(-l+1)/fs:1/fs:(l-1)/fs,y); hold on; [m,n]=size(y); d=0.701*ones(m,n); mesh(d); hold off;title('tao=0');xlabel('ksi');ylabel('tao'); figure contour(-50:0.2:50,(-l+1)/fs:1/fs:(l-1)/fs,y,[0.707 0.707]);title('模糊度图')