untitled0000.m

y3k=fft(u,(m+n-2)/4) i=1:(m+n-2)/4 subplot(5,2,9) stem(i,u) title( 滤波后上采样 ) k=1:(m+n-2

fs=20000,fa=7000;
xt=cos(2*pi*fa*t);
t=0:1/127:1
k=0:128;
xk=m*cos(2*pi*k*(fc/fs));
subplot(321);
stem(k,xk);
xk1=[zeros(1,896) xk];
subplot(322);
plot(xk1);

fs1=32,f1=5,f2=12;
xt1=2*cos(2*pi*f1*t)+cos(2*pi*f2*t);
t=0:1/32:1
k=0:32;
xk1=cos(2*pi*k*(5/32));
xk2=cos(2*pi*k*(12/32));
xk3=xk1+xk2;
subplot(323);
plot(xk3);
xk4=[zeros(1,992) xk3];
xk5=fft(xk4,1023);
subplot(324);
plot(xk5);

n=138;wc=2*fc/fs;
j=fir1(n-1,wc,blackman(n));
hk=fft(j,n);
k=1:n;
g=conv(j,xk4);
gk=fft(g,1024+n-1);
 k=1:1024+n-1;
subplot(325);stem(k,gk);
title('经过blackman窗');

n=75;wc=2*0.25;
j=fir1(n-1,wc,hamming(n));
hk=fft(j,n);
g=conv(j,xk4);
gk=fft(g,1024+n-1);
 k=1:1024+n-1;
subplot(326);stem(k,gk);
title('经过hamming窗');