1. 工作步骤
- 弄清楚WVD实现源码，以及横坐标，纵坐标的变换，阈值等等。
- 用matlab时频工具箱实现信号的WVD变换。（具体参数下文给出 ，以下都是）
- 找到matlab时频工具箱实现WVD变换的源码。
- 然后根据找到的源码，参考源码，写出信号的模糊函数分析AF，
- 参考WVD源码实现信号的CPF。
2. WVD实现需要注意事项
1. 信号得是列信号：
... t=0:1/fs:T-1/fs; t=t'; .... X=A1*exp(1i*(a_10+a_11*t+a_12*t.^2)); [s_lfm_wv,t,f]=tfrwv(X); contour(t/fs,f*fs,abs(s_lfm_wv)); ... t=0:1/fs:T-1/fs; .... X=A1*exp(1i*(a_10+a_11*t+a_12*t.^2)); [s_lfm_wv,t,f]=tfrwv(X'); contour(t/fs,f*fs,abs(s_lfm_wv)); 
我这里就直接使用第一种情况t=t’了不翻转Y纵；
3. WVD 实现
%=================开始wvd变换代码================== t_wvd=1:N; [trow,tcol]=size(t_wvd); [xrow,xcol]=size(X); tfr=zeros(N,tcol); for icol=1:tcol, t_wvdi=t_wvd(icol); taumax=min([t_wvdi-1,xrow-t_wvdi,round(N/2)-1]); tau=(-taumax:taumax); indices=rem(N+tau,N)+1; tfr(indices,icol)=X(t_wvdi+tau,1).*conj(X(t_wvdi-tau,xcol)); tau=round(N/2); if(t_wvdi<=xrow-tau)&(t_wvdi>=tau+1), tfr(tau+1,icol)=0.5*(X(t_wvdi+tau,1)*conj(X(t_wvdi-tau,xcol))+... X(t_wvdi-tau,1)*conj(X(t_wvdi+tau,xcol))); end end tfr1=fft(tfr,nfft,1); if(xcol==1),tfr1=real(tfr1);end f=(0.5*(0:N-1)/N)'; figure; contour(t_wvd/fs,f,abs(tfr1)); title(['WVD变换',' f1=',num2str(f1),' ',' k1=',num2str(k1),' f2= ',num2str(f2),' k2= ',num2str(k2)]); 4. AF实现
[NAF,TAU,XI]=ambifunb(X); figure; contour(TAU/fs,XI,abs(NAF)); xlabel('时间t(s)'); ylabel('归一化频率'); title(['窄带AF变换',' f1=',num2str(f1),' ',' k1=',num2str(k1),' f2= ',num2str(f2),' k2= ',num2str(k2)]); 二、 STFT相关仿真
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/211365.html原文链接：https://javaforall.net
