RenatoCesar Postado Junho 10, 2008 Denunciar Share Postado Junho 10, 2008 Gostaria de saber se alguém sabe como vetorizar o código abaixo:t = 0:0.001:2; % 2 secs @ 1kHz sample ratey = chirp(t,0,1,150); % Start @ DC, cross 150Hz at t=1 secs=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 Gaborc=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 endendfigure(4)[lx,ly]=size©;mesh(abs(c(:,1:ly/2)))colormap('hsv')shading interpview([0 0 1])Obrigado. Citar Link para o comentário Compartilhar em outros sites More sharing options...
Pergunta
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.
Link para o comentário
Compartilhar em outros sites
0 respostass a esta questão
Posts Recomendados
Participe da discussão
Você pode postar agora e se registrar depois. Se você já tem uma conta, acesse agora para postar com sua conta.