wv分布和模糊函数

一、实验题目：信号的wigner-ville分布及模糊函数二、实验目的：了解掌握wigner-ville分布和模糊函数，并能够在matlab上实现。三、实验软件：matlab四、实验原理1、wigner-ville分布取时间冲击函数作为窗函数，则局部相关函数为称为信号(t)z的瞬时相关函数。对瞬时相关函数作关于滞后的傅里叶变换，即得到时频分布为称为wigner-ville分布。2、模糊函数Wigner-ville分布是对瞬时相关函数的时延参数作傅里叶变换，如果改为对时间变量t作反傅里叶变换，则得另一种时频分布称为模糊函数。五、实验结果Figure1原始复信号实部Figure2信号的wigner-ville分布Figure3信号的模糊函数六、程序代码clear;clcn=1001;%点数2L+1w=2*pi*20;m=15e3;n_cl=2.5;%载波周期T=n_cl*2*pi/w;%总时长t=0:T/(n-1):T;z=exp(i*(w*t+1/2*m*t.^2));%线性调频信号plot(t,real(z));title('signal-real');axistight;xlabel('Time（s）');ylabel('amplitude');%%vigner-ville分布fs=n/T;dt=1/fs;nfft=512;df=fs/nfft;[wv,t1,f1]=tfrwv(z',1:n,nfft);%WVDf1(nfft,n)=0;f1(:,1)=((1:nfft)*df)';forj=2:nf1(:,j)=f1(:,1);endt1(nfft,n)=0;forj=1:nfftt1(j,:)=t;endfigure(2)surf(t1,f1./2,10*log10(abs(wv)),'EdgeColor','none');colorbar;axisxy;axistight;colormap(jet);view(0,90);xlabel('Time');ylabel('Frequency(Hz)');title('WV分布');%%模糊函数[amag,delay,fre_shift]=ambgfun(z,fs,1/T);%求模糊函数[t2,f2]=meshgrid(delay,fre_shift);figure(3)surf(t2,f2,amag,'EdgeColor','none');colorbar;axisxy;axistight;colormap(jet);view(0,90);xlabel('delay');ylabel('fre-shift');title('模糊函数');注：tfrwv函数为时频分析工具箱中的函数；ambgfun函数为phasedarraysystemtoolbox中的函数。