wv分布和模糊函数.docxVIP

实验题目：信号的wigner-ville分布及模糊函数 实验目的：了解掌握wigner-ville分布和模糊函数，并能够在matlab上实现。 实验软件：matlab 实验原理 1、wigner-ville分布 取时间冲击函数 作为窗函数，则局部相关函数为 称为信号的瞬时相关函数。对瞬时相关函数作关于滞后的傅里叶变换，即得到时频分布为 称为wigner-ville分布。 2、模糊函数 Wigner-ville分布是对瞬时相关函数的时延参数作傅里叶变换，如果改为对时间变量t作反傅里叶变换，则得另一种时频分布 称为模糊函数。 实验结果 Figure SEQ Figure \* ARABIC 1原始复信号实部 Figure SEQ Figure \* ARABIC 2信号的wigner-ville分布 Figure SEQ Figure \* ARABIC 3信号的模糊函数 程序代码 clear;clc n=1001;%点数2L+1 w=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);axis tight; 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);%WVD f1(nfft,n)=0; f1(:,1)=((1:nfft)*df); for j=2:n f1(:,j)=f1(:,1); end t1(nfft,n)=0; for j=1:nfft t1(j,:)=t; end figure(2) surf(t1,f1./2,10*log10(abs(wv)),EdgeColor,none); colorbar;axis xy; axis tight; 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;axis xy; axis tight; colormap(jet);view(0,90); xlabel(delay); ylabel(fre-shift); title(模糊函数); 注：tfrwv函数为时频分析工具箱中的函数；ambgfun函数为phased array system toolbox中的函数。