clc
clear all
fs=512;
dt=1/fs;
N=2048;
f=fs*(0:N-1)./N;
t=0:dt:2047*dt;
Wd=hann(2048);
y1=500*sin(pi*t);
y2=5*sin(20*pi*t);
s=y1+y2
subplot(411)
plot(t,y1);
fs1=dt*fft(s)./(2*pi);%j=7,s=127
figure(1)
subplot(412);
plot(f-fs/2,abs([fs1((N/2+1):N) fs1(1:N/2)]))
m=45;
n=48;
fsai=(exp(2*n*pi*t*i)-exp(2*pi*m*t*i))./(2*pi*(n-m)*t*i+1e-100).*Wd';
f1=dt*fft(fsai)./(2*pi);
subplot(413)
f=fs*(0:N-1)./N;
plot(f-fs/2,abs([f1((N/2+1):N) f1(1:N/2)]));
N=length(s);
W=f1.*fs1;
subplot(414)
plot(f-fs/2,abs([W((N/2+1):N) W(1:N/2)]));
w1=4*pi*2*pi*ifft(W)./dt;
figure(2)
subplot(211)
plot(t,real(w1))
subplot(212)
plot(t,imag(w1))
figure1
信号y1时域图形
信号s(=y1+y2)频谱
hann窗口构造滤波器频谱
s信号滤波后频谱
figure2
滤波后信号的实部和虚部波形
追问大神!!膜拜,太感谢了
不是的,我是问figure1和figure2显示出来的图形是什么
追答这个你程序输到matlab里面就知道了,我电脑上的matlab已经删了.....