wv分布和模糊函數(shù)(共4頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上一、 實(shí)驗(yàn)題目：信號(hào)的wigner-ville分布及模糊函數(shù)二、 實(shí)驗(yàn)?zāi)康模毫私庹莆誻igner-ville分布和模糊函數(shù)，并能夠在matlab上實(shí)現(xiàn)。三、 實(shí)驗(yàn)軟件：matlab四、 實(shí)驗(yàn)原理1、wigner-ville分布 取時(shí)間沖擊函數(shù)作為窗函數(shù)，則局部相關(guān)函數(shù)為稱為信號(hào)的瞬時(shí)相關(guān)函數(shù)。對(duì)瞬時(shí)相關(guān)函數(shù)作關(guān)于滯后的傅里葉變換，即得到時(shí)頻分布為稱為wigner-ville分布。2、模糊函數(shù)Wigner-ville分布是對(duì)瞬時(shí)相關(guān)函數(shù)的時(shí)延參數(shù)作傅里葉變換，如果改為對(duì)時(shí)間變量t作反傅里葉變換，則得另一種時(shí)頻分布稱為模糊函數(shù)。五、 實(shí)驗(yàn)結(jié)果Figure 1原始復(fù)信號(hào)實(shí)部

2、Figure 2信號(hào)的wigner-ville分布Figure 3信號(hào)的模糊函數(shù)六、 程序代碼clear;clcn=1001;%點(diǎn)數(shù)2L+1w=2*pi*20;m=15e3;n_cl=2.5;%載波周期T=n_cl*2*pi/w;%總時(shí)長(zhǎng)t=0:T/(n-1):T;z=exp(i*(w*t+1/2*m*t.2);%線性調(diào)頻信號(hào)plot(t,real(z);title('signal-real');axis tight;xlabel('Time（s）');ylabel('amplitude');%vigner-ville分布fs=n/T;dt=1/

3、fs;nfft =512;df=fs/nfft;wv,t1,f1=tfrwv(z',1:n,nfft);%WVDf1(nfft,n)=0;f1(:,1)=(1:nfft)*df)'for j=2:n f1(:,j)=f1(:,1);endt1(nfft,n)=0;for j=1:nfft t1(j,:)=t;endfigure(2)surf(t1,f1./2,10*log10(abs(wv),'EdgeColor','none');colorbar;axis xy; axis tight; colormap(jet);view(0,90);xla

4、bel('Time');ylabel('Frequency (Hz)');title('WV分布');%模糊函數(shù)amag,delay,fre_shift=ambgfun(z,fs,1/T);%求模糊函數(shù)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('模糊函數(shù)');注：tfrwv