课程作业-杭电研究生-小波变换-第三次作业.docVIP

多信号之和的维格纳-威尔分布结构 两个信号之和 %产生非平稳信号 sig=atoms(128,[32,0.15,20,1;96,0.15,20,1]); %时域波形 figure(2) plot(real(sig),LineWidth,2); xlabel(时间 t); ylabel(幅值 A); %计算Wigner-Ville分布 [tfr,t,f]=tfrwv(sig); %时频表示 figure(3) contour(t,f,abs(tfr)); xlabel(时间 t); ylabel(频率 f);  三个信号之和 %产生非平稳信号 sig=atoms(128,[32,0.35,20,1;96,0.15,20,1;96,0.35,20,1]); %时域波形 figure(2) plot(real(sig),LineWidth,2); xlabel(时间 t); ylabel(幅值 A); %计算Wigner-Ville分布 [tfr,t,f]=tfrwv(sig); %时频表示 figure(3) tfrwv(sig); xlabel(时间 t); ylabel(频率 f);  比较跳频信号的STFT变换与Wigner-Ville分布的频谱图 clear close all clc %产生跳频信号 n=0:383; ts=1/(80*10^3);%采样时间 f=4.21*10^3; delta_f=8000;%频偏 f1=f; f2=f+delta_f; f3=f2+delta_f; y1=cos(2*pi*f1*ts*n).*[ones(1,128),zeros(1,128*2)]; y2=cos(2*pi*f2*ts*n).*[zeros(1,128),ones(1,128),zeros(1,128)]; y3=cos(2*pi*f3*ts*n).*[zeros(1,128*2),ones(1,128)]; sig=y1+y2+y3; figure(1); plot(n,sig); title(跳频信号); sig=sig; %计算短时傅立叶变换 sig=hilbert(sig); %计算短时傅立叶变换 [tfr,t,f]=tfrstft(sig); %短时傅立叶变换时频表示 figure(2); tfrstft(sig); xlabel(时间 t); ylabel(频率 f); %计算Wigner-Ville分布 [tfr,t,f]=tfrwv(sig); %时频表示 figure(3) tfrwv(sig); xlabel(时间 t); ylabel(频率 f);  四个高斯信号构成的Wigner-Ville分布与伪Wigner-Ville分布 %产生非平稳信号 sig=atoms(128,[32,0.15,20,1;96,0.15,20,1;32,0.35,20,1;96,0.35,20,1]); %时域波形 figure(2) plot(real(sig),LineWidth,2); xlabel(时间 t); ylabel(幅值 A); %计算Wigner-Ville分布 [tfr,t,f]=tfrwv(sig); %时频表示 figure(3) tfrwv(sig); xlabel(时间 t); ylabel(频率 f); %计算伪Wigner-Ville分布 sig=hilbert(sig); [tfr,t,f]=tfrpwv(sig); %时频表示 figure(4) tfrpwv(sig); xlabel(时间 t); ylabel(频率 f);