clear
%CREATION DE LA VARIABLE Wr
f1=0.1;f2=0.4;T=64;Fe=1;t=(0:1/Fe:T-1);
beta=(f2-f1)/(2*T);delta=0.001*Fe;
sig=sin(2*pi*(f1+beta*t).*t);
N=length(sig)+rem(length(sig),2);
tailleFFT=N;
if N~=length(sig); sig=[sig 0];end
taille_temps=length(sig); posi=taille_temps;
s=[zeros(1,taille_temps) hilbert(sig) zeros(1,taille_temps)];
s_conj=conj(s);
W=zeros(tailleFFT,taille_temps);
tau1=linspace(0,N/2,N/2+1);
tau2=linspace(2,N/2,N/2-1);
for t=1:taille_temps,
R(tau1+1,t)=(s(posi+tau1).*s_conj(posi-tau1)).';
R(N+2-tau2,t)=conj(R(tau2,t));posi=posi+1;
end
W=fft(R,tailleFFT)/(2*tailleFFT);
h=real(W);
Wr=h(1:64,1:64)
absc=linspace(0,N*1/Fe,N);
fwv=linspace(0,Fe/2,N);
f=linspace(0,Fe/2,N/2);
abscB=tailleFFT;
freqc=(linspace(0,0.5,tailleFFT));
mesh(absc,fwv,Wr)
%ANIMATION
surf(Wr)
axis vis3d
for i=1:360
camorbit(1,0,'data',[0 0 1])
drawnow
end
Partager