文档介绍：1、分别用短时Fourier,STFT,Gabor变换分析下列信号,要求提供程序,图形结果并对它们的结果进行对比分析。采样频率FS二1920HZ,(t)=(1+(2^15f))cos(2^30/+(2”l5/))+sin(2^120z)fs=1920;%采样频率N=512;T=l/fs;fl=15;f2=30;f3=120;%3个信号的频率n=0:N-l;t=0:l/fs:(N-l)/fs;%时间序歹I」xl=(1+*sin(2*pi*15*t))・★(cos(2*pi*30*t)+0・5*sin(2*pi*15*t))+sin(2*pi*120*t)八信号figure(1)plot(tzxl); %信号吋域波形yl=fft(xlfN); %对信号进行快速Fourier变换magl=abs(yl);%求得Fourier变换后的振幅k=0:N-l;fl=k*fs/N; %转化为对应模拟信号的频率figure(2)gridon %网格开启stem(fl,magi);boxon%绘制N点DFT的幅频特性图xlabel(^fl1);ylabel(1巾畐度1);axis([0,512,*max(abs(yl))]) 色限定x轴和y轴上值的范围,有利于观察gridonh=window(1hamming*,512);[B,frt]=specgram(xl^512,fsrh,384);%specgram在这里等效于tfrstft,用specgram更有利于观察,因为用tfrstft时频率限定杂--,%H]tfrstft观察不在-0•5-0・5之间的频率成分,需要重新改变采样频%率,使其归一化的值在--。figure(3)imagesc(tffAabs(B)); %画出经过短时傅里叶变换后的图像axis(*xy');xlabeK^^fiiJ(秒)');ylabel(1 (Hz)');kkl=99;kk2=17;q=4;%q=4表示临界抽样。sig=xl;h=window(***@gausswinzkkl); %高斯窗h=h/norm(h);figure(4);tfrgabor(hilbert(sig「),32,q,h); %tfrgabor变换第一个参数需要的是列向量,先经过hibert变换讲信号变成解析形式,xlabeK^'Jfil](秒)');ylabel('频率(Hz)');sig=xl;h=window(***@gausswin,kk2); %高斯窗h=h/norm(h);figure(5);tfrgabor(hilbert(sig「),32,q,h); %tfrgabor变换第一个参数需要的是列向量,先经过hibert变换讲信号变成解析形式,===32,Q=,=5%=8,Nf==32,Q=4,linscaleimagesc,Thld=5%0451高斯窗去99吋的gabor•变换 高斯窗取17时的gabor变换通过上图分析可以看出信号有0Hz,15Hz,30Hz,45Hz, 120Hz的频率,进行短吋傅里叶变换的时候要采用合适的采样频率,同时要采用合适的窗的人小,一般窗越小,时域分辨率越好,窗越大频域分辨率越好-便于观察频率成分。此外,当采用tfrstft进行短时傅里叶变换时,因为默认的归一化的频率在--,为了将得到的频率成分显示清晰,需要选择适合的采样频率。tfrgabor变换需要满足一定的条件,本文取Q二4,丁作在过抽样的状态下(ab<l)当在欠抽样(ab>l)>临界抽样(ab=l)时,a(r),不可进彳亍gabor变换。2、分别用wigner-ville分布,伪wigner-ville分布,平滑伪wigner-ville分布和Cohen分布分析下列信号:其屮&, r0=5s,a)0=,并对它们的结果进行对比分析。a=;t0=5;fs=10; %信号的采样频率w0=;n=128;T二1/fs;t=0:1/fs:(n~l)/fs;x二exp(-a*(t+t0)・"2+w0*t*i)+exp(-a*(t-t0)「2+w0*t*i); %信号figure(1)plot(abs(x)); %信号的时域波形X二X・';set(gca,'xlim',[0,n]);set(gca,'xtick',[0:n/4:n])figure(2);[Rtf]二tfrwv(hilbert(x)); %维格拉变换,有交叉项imagesc(