1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
| stacksize('max');
fmax=fech;
tailleFFT=length(C);
n1=ceil(tailleFFT/(fmax));
Tampon=zeros(fmax,1);
a=0;
i1=0;
X=zeros(fmax+1,n1);
for i=1:fmax:tailleFFT,
i1=i1+1;
if i==1 then
Tampon=C(fmax);
elseif (i+fmax)>tailleFFT then
L=length(Tampon)
for J=1:(fmax-L),
Tampon(L+i)=0;
end
else
Tampon=C(i:i+fmax);
i=i;
end
i1=i1;
X(:,i1)=fft(Tampon);
end
x=abs(X)
x=x(1:12500,:);
figure(1)
ftotal=0:(size(x,1)-1);
plot(ftotal',x)
figure(2)
set(gcf(),"color_map",jetcolormap(128));
deplacementf=0:size(x,1)-1;
deplacementt=0:size(x,2)-1;
clf;grayplot(deplacementt,deplacementf,x') |
Partager