untitled000.m

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

xlabel('k');ylabel('gk');
    f=zeros(1,4*(m+n-1));
 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(f,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);plot(u);
title('滤波后上采样信号');
k=1:(m+n-2)/4;
subplot(529);stem(k,y3k);
title('滤波后上采样频谱');
xlabel('k');ylabel('y3k');