原创 Lab5：加窗:fir滤波器的时域设计

Lab5：加窗:fir滤波器的时域设计
一：实验内容
模拟出理想的滤波器
加矩形窗后
加三角窗后
加其他窗后
二：原理
理想滤波器：
2.1.1原理

2.1.2代码：
tnum=1000;%时间间隔(把2*pi分为tnum份)，由于理想，模拟连续
tmax=1000;tmin=-1000;%时间范围，由于理想，模拟无限
t=tmin:2*pi/(tnum-1):tmax;
wc=20/pi;%截止频率位：rwc=wc*pi
y=sinc(wc*t)*wc;%理想滤波器的表达式sin(rwc*t)./(pi*t)其中rwc为wc*pi
subplot(3,1,1);
plot(t,y);
num=1000;
w=0:2*pi/(num-1):2*pi;
h=freqz(y,1,w);
h_abs=abs(h);
h_a=angle(h);
subplot(3,1,2);
plot(hh_abs);%幅频特性曲线
subplot(3,1,3);
plot(h_a);%相频特性曲线
2.1.3：仿真结果

2.1.4放大：
plot(hh_abs);
axis([0,100,0,200]);

加入矩形窗
2.2.1：原理
时间范围变小
2.2.2代码
tnum=1000;%时间间隔(把2*pi分为tnum份)，由于理想，模拟连续
%%%%%%%%%%%新加入%%%%%%%%%%
tmax=10;tmin=-10;%时间范围，由于理想，模拟无限
t=tmin:2*pi/(tnum-1):tmax;
wc=20/pi;%截止频率位：rwc=wc*pi
y=sinc(wc*t)*wc;%理想滤波器的表达式sin(rwc*t)./(pi*t)其中rwc为wc*pi
subplot(3,1,1);
plot(t,y);
num=1000;
w=0:2*pi/(num-1):2*pi;
h=freqz(y,1,w);
h_abs=abs(h);
h_a=angle(h);
subplot(3,1,2);
plot(hh_abs);%幅频特性曲线
subplot(3,1,3);
plot(h_a);%相频特性曲线
2.2.3仿真结果

2.2.4放大
plot(hh_abs);
axis([0,100,0,200]);

离散化
2.3.1时间间隔变小
2.2.2代码
%%%%%%%%%%%新加入%%%%%%%%%%
%%由1000变为10
tnum=10;%时间间隔(把2*pi分为tnum份)，由于理想，模拟连续
tmax=1000;tmin=-1000;%时间范围，由于理想，模拟无限
t=tmin:2*pi/(tnum-1):tmax;
wc=20/pi;%截止频率位：rwc=wc*pi
y=sinc(wc*t)*wc;%理想滤波器的表达式sin(rwc*t)./(pi*t)其中rwc为wc*pi
subplot(3,1,1);
plot(t,y);
num=1000;
w=0:2*pi/(num-1):2*pi;
h=freqz(y,1,w);
h_abs=abs(h);
h_a=angle(h);
subplot(3,1,2);
plot(hh_abs);%幅频特性曲线
subplot(3,1,3);
plot(h_a);%相频特性曲线
结果：

离散加矩形窗？？？？？？？？？
%%%%%%%%%%%新加入%%%%%%%%%%
%%由1000变为11
tnum=11;%时间间隔(把2*pi分为tnum份)，由于理想，模拟连续
%%%%%%%%%%%新加入%%%%%%%%%%
%%由1000变为10
tmax=10;tmin=-10;%时间范围，由于理想，模拟无限
t=tmin:2*pi/(tnum-1):tmax;
wc=20/pi;%截止频率位：rwc=wc*pi
y=sinc(wc*t)*wc;%理想滤波器的表达式sin(rwc*t)./(pi*t)其中rwc为wc*pi
subplot(3,1,1);
stem(t,y);
num=1000;
w=0:2*pi/(num-1):2*pi;
h=freqz(y,1,w);
h_abs=abs(h);
h_a=angle(h);
subplot(3,1,2);
plot(hh_abs);%幅频特性曲线
subplot(3,1,3);
plot(h_a);%相频特性曲线