function M=movpwjph(N,Np,typesig);%MOVPWJPH Influence of a jump of phase on the interferences of the pWVD.  %	M=MOVPWJPH(N,Np,TYPESIG) generates the movie frames illustrating the %	influence of a jump of phase in different frequency modulations%	on the interference terms of the pseudo Wigner-Ville distribution.%%	N : number of points for the signal;%	Np : number of snapshots (default : 7)%	TYPESIG : type of signal :%	 'C' : constant frequency modulation (default value) ;%	 'L' : linear frequency modulation ;%	 'S' : sinusoidal frequency modulation.%	M : matrix of movie frames.%%	Example : %	 M=movpwjph(128,8,'S'); movie(M,10);%	O. Lemoine - May 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, typesig='C'; Np=8;elseif nargin==2, typesig='C'; endM  = moviein(Np);typesig=upper(typesig);if typesig=='C',for k=1:Np, sig=[fmconst(N/2);fmconst(N/2)*exp(j*2*pi*k/Np)]; [tfr,t,f]=tfrpwv(sig);  Max=max(max(tfr));V=[0.3 0.5 0.7 0.9]*Max; contour(t,f,tfr,V);xlabel('Time'); ylabel('Frequency');  title('Pseudo Wigner-Ville distribution'); axis('xy') M(:,k) = getframe;endelseif typesig=='L',for k=1:Np, sig=[fmlin(N/2,0,0.25);fmlin(N/2,0.25,0.5)*exp(j*2*pi*k/Np)]; [tfr,t,f]=tfrpwv(sig);  Max=max(max(tfr));V=[0.3 0.5 0.7 0.9]*Max; contour(t,f,tfr,V);xlabel('Time'); ylabel('Frequency');  title('Pseudo Wigner-Ville distribution'); axis('xy') M(:,k) = getframe;endelseif typesig=='S',sig=fmsin(N,0.05,0.45);sig1=sig(1:N/2);sig2=sig(N/2+1:N);for k=1:Np, sig=[sig1;sig2*exp(j*2*pi*k/Np)]; [tfr,t,f]=tfrpwv(sig);  Max=max(max(tfr));V=[0.3 0.5 0.7 0.9]*Max; contour(t,f,tfr,V);xlabel('Time'); ylabel('Frequency');  title('Pseudo Wigner-Ville distribution'); axis('xy') M(:,k) = getframe;endelseerror('Wrong input parameter for TYPESIG');end