untitled000.asv

来自「y3k=fft(u,(m+n-2)/4) i=1:(m+n-2)/4 」· ASV 代码 · 共 55 行 · 第 1/2 页

 for i=1:m+n-1,
     f(4*i)=g(i);
 end

  i=1:4*(m+n-1);
subplot(526);stem(i,f);
title('滤波后下采样');

fk=fft(y,4*(m+n-1));
k=1:4*(m+n-1);
subplot(527);stem(i,fk);
title('滤波后下采样频谱');
xlabel('k');ylabel('fk');

for i=1:(m+n-2)/4,
    u(i)=g(4*i);
end
y3k=fft(u,(m+n-2)/4);

  i=1:(m+n-2)/4;
subplot(528);stem(i,u);
title('滤波后上采样');

  k=1:(m+n-2)/4;
subplot(529);stem(k,y3k);
title('滤波后上采样频谱');
xlabel('k');ylabel('y3k');