RenatoCesar

Gostaria de saber se alguém sabe como vetorizar o código abaixo: t = 0:0.001:2; % 2 secs @ 1kHz sample rate y = chirp(t,0,1,150); % Start @ DC, cross 150Hz at t=1 sec s=y(:,1); Ls=length(s); s=s(1:4:Ls); L=1024; N=64; DM=16; DN=16; M=Ls/DM; O=zeros(size(g)); % Coeficientes de análise de Gabor c=zeros(M,N); for m=0:M-1, for n=0:N-1, for k=1:Ls, if (k-m*DM)>0 & (k-m*DM)<L, c(m+1,n+1)=c(m+1,n+1)+s(k)*g(k-m*DM)*exp(-im*2*pi*n*DN*(k-1)/L); end end end end figure(4) [lx,ly]=size©; mesh(abs(c(:,1:ly/2))) colormap('hsv') shading interp view([0 0 1]) Obrigado.