数字信号处理大作业

数字信号处理

大作业

姓名：刘云飞学号：08041016 班级：电气81

（1）题目18：用fir1函数法设计线性相位数字高通滤波器，通带边频wc=0.6*pi rad，阻带边频ws=0.45*pi rad，通带最大衰减A p=2dB，阻带最小衰减A S=45dB，分别用汉宁窗，汉明窗，布莱克曼窗和凯泽窗进行设计，求单位脉冲响应h(n)，画出h(n)和幅度相应∣H(e jw)∣。

（2）程序及运行结果：

1.汉宁窗

求得N=44，A S=44.56，A P=0。

2.汉明窗

求得N=44，A S=45.40dB，A P=0。

3.布莱克曼窗

求得N=44，A S=35.03，A P=0。

4.凯泽窗

求得N=44，A S=51.71，A p=0。

（3）附录：自定义的两个函数

function hd=ideal_hp1(wc,N)

%compute the ideal lowpass fiter unit pulse respondence hd(n)

%wc: cutoff frequency

%N: window length

%hd: unit pulse respondence

alpha=(N-1)/2;

n=0:1:N-1;

m=n-alpha+eps;

hd=sin(pi*m)./(pi*m)-sin(wc*m)/(pi*m);

function[db,mag,pha,w]=freqz_m2(b,a)

%滤波器的幅值响应(相对、绝对)、相位响应

%db:相对幅值响应

%mag：绝对幅值响应

%pha: 相位响应

%w 采样频率；

%b 系统函数H(z)的分子项(对FIR，b=h)

%a 系统函数H(z)的分母项(对FIR，a=1)

[H,w]=freqz(b,a,1000,'whole');

H=(H(1:1:501))';

w=(w(1:1:501))';

mag=abs(H); %绝对幅值响应

db=20*log10((mag+eps)/max(mag)); %相对幅值响应

pha=angle(H); %相位响应