学号：姓名。班级：
教师。西昌学院。
设计名称： fir数字高通滤波器。
日期：2012 年 6月 14 日
一、设计内容。
根据冲激响应的时域特性，数字滤波器可分为无限长冲激响应（iir）和有限长冲激响应滤波器（fir），fir的突出优点是：系统总是稳定的、易于实现线性相位、允许设计多通带（或多阻带）滤波器，但与iir相比，在满足同样阻带衰减的情况下需要的阶数较高，滤波器的阶数越高，占用的运算时间越多，因此在满足指标要求的情况下应尽量减少滤波器的阶数。
在自适应处理、数据通信等领域中往往要求信号在传输过程中不能有明显的相位失真，fir 滤波器可以做到线性相位满足此要求。fir 滤波器实质上是一个分节的延迟线，把每一节的输出加权累加，得到滤波器的输出。
fir滤波器的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数的问题，设计方法主要有窗函数，频率采样法和等波纹最佳逼近法等。本设计是运用等波纹最佳逼近法设计一个fir高通数字滤波器。
二、建立数学模型。
1.重构fir的单位抽样相应。
根据频率抽样理论，有n个频率采样点可以唯一确定，即对进行idft
2.重构系统函数。
的频响。将代入h(z)表达式可得。
其中， 为大家所知的内插函数。
分析可知，当时（采样点）有：
这说明，重构的频响，在采样上严格等于h（k），而在采样点之间，频响则由加权的内插函数延伸叠加而成。
4.线性相位的约束条件。
以h(n)为偶对称，n为奇数的情况进行分析。
1）fir的频响具有线性相位的一般表达式。
当h(n)为偶对称，n为奇数时，则。
而且幅度函数h（w）应为偶对称，即。
2）采样值h（k）具有线性相位的约束。
其中， 表示采样值的模（纯标量），表示其相角。因此，在采样点上具有线性相位的条件应为：
而且，必须满足偶对称，即。
5、本设计中用到的是等波纹最佳逼近法中的remez函数设计。
三、数学模型的解析分析。
采用remez算法实现线性相位fir数字滤波器的等波纹最佳逼近法设计。与其他设计法相比，其优点是，设计指标相同时，使滤波器阶数最低；或阶数相同时，使通带最平坦，阻带最小衰减最大；通带和阻带均为等波纹形式，最适合设计片段常数特性的滤波器。其调用格式如下：
b=remez(n,f,m,w,’ftype’)
其中，w和ftype可默认。b为滤波器系数向量，调用参数n,f,m的含义与函数fir2中类同，但这里有一点不同，期望逼近的频幅响应值位于f(k) 与f(k+1)（k为奇数）之间的频段上，而f(k+1)与f(k+2)之间为无关区。w为加权向量，其长度为f的一半。
w(k)为对m中第k个常数片段的逼近精度加权值，w值越大逼近精度越高。rtype用于指定滤波器类型。
remezord函数用于估算fir数字滤波器的等波纹最佳一致逼近设计的最低阶数n，从而使滤波器在满足指标的前提下造价最低。基本调用格式如下：
n,fo,mo,w]=remezord(f,m,dev,fs)
其返回参数供remez函数使用。设计的滤波器可以满足由参数f,m,dev和fs指定的指标。f和m与remez中所用的类似，这里f可以是模拟频率 (hz)或归一化数字频率，但必须以0开始，以fs/2(用归一化频率时为1)结束，而且其中省略了0和fs/2两个频点。
fs为采样频率，省略时默认为 2hz。dev为各逼近频段允许的幅频响应偏差（波纹振幅）。
remez函数可直接调用remezord返回的参数，使用格式如下：
b=remez(n,fo,mo,w)。
四、计算机**模型的建立（编程）
用remez函数设计高通滤波器。观察等波纹逼近法中加权系数w( )及滤波器阶数n的作用和影响。期望逼近的滤波器通带为[3 /4，1 ]，阻带为[0，23 /32]。
在滤波器设计中，技术指标越高，实现滤波器的阶数也就越高。
在remez函数调用格式b=remez(n,f,m,w)中，f=[0,3/4,23/32,1],m=[0,0,1,1].其余参数分三种情况进行设计：(1)n=30,w=[1,1];(2)n=30,w=[1,5](3)n=60,w=[1,1]。
程序如下：cjf2_
用remez函数设计高通滤波器。
clear;close all
f=[0, 23/32, 3/4,1];m=[0,0,1,1];
n1=30;w1=[1,1]; hn1=remez(n1,f,m,w1情况(1)
k=[0:1023]*2/1024;
hw1=fft(hn1,1024);
subplot(311);plot(k,20*log10(abs(hw1求出其幅频特性。
axis([0,1,-40,5]);grid on只画出正半轴频谱。
xlabel('w/pi');ylabel('magnitude(db标注。
title(' n1=30,w1=[1,1] '
n2=30;w2=[1,5];hn2=remez(n2,f,m,w2情况(2)
hw2=fft(hn2,1024求出其幅频特性。
subplot(312);plot(k,20*log10(abs(hw2)))
axis([0,1,-40,5]);grid on只画出正半轴频谱。
xlabel('w/pi');ylabel('magnitude(db标注。
title(' n2=30,w2=[1,5] '
n3=60;w2=[1,1];hn3=remez(n3,f,m,w2情况(3)
hw3=fft(hn3,1024求出其幅频特性。
subplot(313);plot(k,20*log10(abs(hw3)))
axis([0,1,-40,5]);grid on只画出正半轴频谱。
xlabel('w/pi');ylabel('magnitude(db标注。
title(' n3=60,w1=[1,1] '
五、执行**和结果分析。
由图可见，w较大的频段逼近精度越高；w越小的频段逼近精度越低。n较大时逼近精度越高；n越小时逼近精度越低。
六、设计体会与建议
为了这次课程设计，我自学了数字信号处理领域中窗函数的有关知识。实际中遇到的离散时间信号总是有限长的，因此不可避免地要遇到数据截断问题。而在信号处理中，对离散序列的数据截断是通过序列与窗函数相乘来实现的。
MATLAB课程设计报告
基于matlab的图像处理 课程设计。学生姓名 徐清福。学号 110700136 专业班级 信通一班 指导教师 郑晓明 二 一零年 4 月 8 日。目录 二号字体 一。课程设计目的2 二。课程设计要求2 三。课程设计报告内容2 四。程序需求和模块分析2 五。详细设计过程3 六。调试及分析16 七。心...
MATLAB课程设计报告
2 整流电路 2.1 单相半波可控整流系统。2.1.1 晶闸管元件的符号和 模型。晶闸管是一种可以通过门极信号触发导通的半导体器件。晶闸管 模型由一个电阻r 一个电感l 一个直流电压源v，和一个开关串联组成。开关受逻辑信号控制，该逻辑信号由电压 电流和门极触发信号 g 决定。晶闸管元件的符号和 模型...
matlab课程设计报告
matlab 软件的使用。姓名 专业 通信技术 班级 通技102 学号 指导老师 史兆强 时间 2013 01 04 用mtlab编程实现句柄图形应用。摘要 本文利用mtlab的算法设计和程序设计的原理和方法，着重介绍了句柄图形的实际运用。根据提出的实际问题，即画出沿曲线移动的小球，给出具体句柄图形...