?? movsp2wv.m
字號:
function M=movsp2wv(N,Np);%MOVSP2WV Movie illustrating the passage from the spectrogram to the WVD.% M=MOVSP2WV(N,Np) generates the movie frames illustrating the passage% from the spectrogram to the WVD using different smoothing gaussian% windows in the smoothed pseudo-WVD. %% N : number of points for the signal;% Np : number of snapshots (default : 8)% M : matrix of movie frames.%% Example : % M=movsp2wv(128,15); movie(M,5);% O. Lemoine - June 1996.% Copyright (c) 1996 by CNRS (France).%% ------------------- CONFIDENTIAL PROGRAM -------------------- % This program can not be used without the authorization of its% author(s). For any comment or bug report, please send e-mail to % f.auger@ieee.org if nargin<1, error('At least one argument required');elseif nargin==1, Np=8;endNp=odd(Np)-1;M = moviein(Np);t1=N/4; t2=3*N/4; t3=t1; t4=t2; f1=0.125; f2=f1; f3=0.375; f4=f3; t=1:N; T=sqrt(2*N);sig1=amgauss(N,t1,T).*fmconst(N,f1,t1); sig2=amgauss(N,t2,T).*fmconst(N,f2,t2); sig3=amgauss(N,t3,T).*fmconst(N,f3,t3); sig4=amgauss(N,t4,T).*fmconst(N,f4,t4); sig=sig1+sig2+sig3+sig4;h=window(49,'gauss');[tfr,t,f]=tfrsp(sig,t,N,h); tfr=tfr(1:N/2,:); f=f(1:N/2);Max=max(max(tfr)); V=[0.1 0.3 0.5 0.7 0.9]*Max;contour(t,f,tfr,V);xlabel('Time'); ylabel('Frequency'); axis('xy')M(:,1) = getframe;Ng0=odd(fliplr(linspace(3,N/4,Np/2-1)));Nh0=odd(linspace(N/2,9*N/10,Np/2-1));for k=2:Np/2, g=window(Ng0(k-1),'gauss'); h=window(Nh0(k-1),'gauss'); [tfr,t,f]=tfrspwv(sig,t,N,g,h); Max=max(max(tfr)); V=[0.1 0.3 0.5 0.7 0.9]*Max; contour(t,f,tfr,V); xlabel('Time'); ylabel('Frequency'); axis('xy') M(:,k) = getframe;end[tfr,t,f]=tfrwv(sig);Max=max(max(tfr)); V=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9]*Max;contour(t,f,tfr,V);xlabel('Time'); ylabel('Frequency'); axis('xy')M(:,Np/2+1) = getframe;for k=Np/2+2:Np, M(:,k) =M(:,Np+2-k);end?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -