windowxulie.m

来自「呵呵,第一次来很多东西不懂 还请多多指教啊! 几个编的简单的FFT和滤波器设计 」· M 代码 · 共 17 行

clear all;                 %q
N=20;
rn=ones(1,N);
n=0:N-1;
subplot(311)
stem(n,rn);
syms w;
Xw=sum(rn.*exp(-i*w*n));
yw=abs(Xw);
w=-2:0.01:2;
yw =abs(1+exp(-i*w)+exp(-2*i*w)+exp(-3*i*w)+exp(-4*i*w)+exp(-5*i*w)+exp(-6*i*w)+exp(-7*i*w)+exp(-8*i*w)+exp(-9*i*w)+exp(-10*i*w)+exp(-11*i*w)+exp(-12*i*w)+exp(-13*i*w)+exp(-14*i*w)+exp(-15*i*w)+exp(-16*i*w)+exp(-17*i*w)+exp(-18*i*w)+exp(-19*i*w))
subplot(312)
plot(w,yw);
w=-3*pi:0.1:3*pi;
pw=-w*19/2;
subplot(313)
plot(w,pw)