非降采样Contourlet域内空间频率激励的PCNN图像融合算法

频率采样法设计高通FIR数字滤波器(范本)

课程设计任务书 学生姓名：胡双印专业班级：通信1005班指导教师：刘新华工作单位：信息工程学院题目:数字高通FIR滤波器设计 要求完成的主要任务: 1.在数字信号处理平台上（PC机﹑MATLAB仿真软件系统）进行软件仿真设计，并进行调试和数据分析。 2. 利用MATLAB仿真软件系统结合频率取样法设计一个数字高通FIR滤波器。 课程设计的目的： 1．理论目的 课程设计的目的之一是为了巩固课堂理论学习，并能用所学理论知识正确分析信号处理的基本问题和解释信号处理的基本现象。 2．实践目的 课程设计的目的之二是通过设计具体的图像信号变换掌握图像和信号处理的方法和步骤。 时间安排： 指导教师签名：年月日系主任（或责任教师）签字：年月日

目录 摘要............................................................................................................................ I Abstrct ........................................................................................................................... II 1 引言. (1) 1.1MATLAB介绍 (1) 1.2MATLAB信号处理工具箱函数介绍 (1) 1.3滤波器的介绍 (2) 2 FIR数字滤波器设计原理 (3) 3 FIR数字滤波器设计方法 (4) 3.1窗函数法 (4) 3.2频率取样法 (5) 4 频率采样法实际FIR高通滤波器 (7) 4.1设计原理 (7) 4.2设计步骤 (9) 5 MATLAB环境下设计FIR数字高通滤波器 (9) 5.1设计要求 (9) 5.2 FIR数字高通滤波器程序设计 (10) 5.3调试结果 (11) 5.4 高通FIR数字滤波器的进一步设计 (12) 6 高通FIR数字滤波器性能测试 (14) 6.1高通FIR数字滤波器性能测试程序 (14) 6.2 性能测试结果 (15) 7 FDATOOL工具箱设计高通FIR滤波器 (16) 7.1 FDATOOL工具箱 (16) 7.2 FIR滤波器参数设置 (17) 8心得体会 (19) 参考文献 (20) 附件：MATLAB程序 (21)

频率采样法设计FIR数字滤波器

实验八频率采样法设计FIR数字滤波器 一、实验目的 掌握频率取样法设计FIR数字滤波器的原理及具体方法。 二、实验设备与环境 计算机、MATLAB软件环境 三、实验基础理论 1.基本原理 频率取样法从频域出发，把理想的滤波器等间隔取样得到，将作为实际设计滤波器的 ,N-1 得到以后可以由来唯一确定滤波器的单位脉冲响应， ()D_Dd___________e??________________求得 其中为内插函数 由求得的频率响应来逼近。 如果我们设计的是线性相位FIR滤波器，则的幅度和相位一定满足线性相位滤波器的约束条件。 我们将表示成如下形式

当为实数，则 由此得到 即以k=N/2为中心呈偶对称。再利用线性条件可知，对于1型和2型线性相位滤波器 对于3型和4型线性相位滤波器 其中，表示取小于该数的最大的整数。 2.设计步骤 （1）由给定的理想滤波器给出和。 （2）由式求得。 （3）根据求得和。 四、实验内容 1.采用频率采样设计法设计FIR数字低通滤波器，满足以下指标 （1）取N=20，过渡带没有样本。 （2）取N=40，过渡带有一个样本，T=0.39。 （3）取N=60，过渡带有两个样本，T1=0.5925，T2=0.1009。 （4）分别讨论采用上述方法设计的数字低通滤波器是否能满足给定的技术指标。

实验代码与实验结果 （1）N=20 过渡带没有样本 N=20; alpha=(N-1)/2; l=0:N-1; wl=(2*pi/N)*l; Hrs=[1,1,1,zeros(1,15),1,1]; *对理想幅度函数取样得到取样样本Hdr=[1,1,0,0];wdl=[0,0.25,0.25,1]; *用于绘制理想函数幅度函数的曲线k1=0:floor((N-1)/2); k2=floor((N-1)/2)+1:N-1; angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)]; H=Hrs.*exp(j*angH); *计算H(k) h=ifft(H,N); *计算h(n) w=[0:500]*pi/500; H=freqz(h,1,w); *计算幅度响应 [Hr,wr]=zerophase(h); *计算幅度函数 subplot(221); plot(wdl,Hdr,wl(1:11)/pi,Hrs(1:11),'o'); axis([0,1,-0.1,1.1]); xlabel('\omega(\pi)'); ylabel('Hr(k)'); subplot(222); stem(l,h,'filled'); axis([0,N-1,-0.1,0.3]); xlabel('n');ylabel('h(n)'); subplot(223); plot(wr/pi,Hr,wl(1:11)/pi,Hrs(1:11),'o'); axis([0,1,-0.2,1.2]); xlabel('\omega(\pi)'); ylabel('Hr(w)'); subplot(224); plot(w/pi,20*log10((abs(H)/max(abs(H))))); axis([0,1,-50,5]); grid;xlabel('\omega(\pi)'); ylabel('dB');

用频率采样法设计FIR数字滤波器

用频率采样法设计FIR 数字滤 波器 信号、系统与信号处理实验Ⅱ 实验报告 实验名称：用频率采样法设计FIR 数字滤波器 一、实验目的 掌握频率取样法设计FIR 数字滤波器，加深过渡点对滤波器性能影响的认 识。 二、实验内容与要求 ( 1)编写好一个设计线性相位FIR 高通滤波器的程序，已知wc=0.8 ， N=64，要求在屏幕上显示出h(n) 值，画出|H(e^jw)| 及20lg(|H(e^jw)) 的曲线。 ( 2)实验时，设置0 个过渡点， 1 个过渡点， 2 个过渡点，比较设计所得的|H(e^jw)| 及20lg(|H(e^jw)) 的曲线。

三、实验程序与结果 (1)0个过渡点clear all ; N=64; wc=0.8*pi; k=0:N-1; phase=(-pi*k*(N-1)/N)+pi/2; HK=[zeros(1,26),ones(1,13),zeros(1,25)];

HK1=HK.*exp(j*phase); hn=ifft(HK1,N) figure(1); freqz(hn,1,512); [H,W]=freqz(hn,1,512); figure(2); subplot(3,1,1); stem(k,hn); title( 'h(n)' ) subplot(3,1,2); plot(W/pi,abs(H)); title( '|H(eiw)|' ) subplot(3,1,3); plot(W/pi,20*log10(abs(H))); title( '20lg|H(eiw)|' );

(2)1个过渡点 clear all ; N=64; wc=0.8*pi; k=0:N-1; phase=(-pi*k*(N-1)/N)+pi/2; HK=[zeros(1,25),0.5,ones(1,13),0.5,zeros(1,24)]; HK1=HK.*exp(j*phase); hn=ifft(HK1,N); figure(1); freqz(hn,1,512); [H,W]=freqz(hn,1,512); figure(2); subplot(2,1,1); plot(W/pi,abs(H)); title( '|H(eiw)|' ) subplot(2,1,2); plot(W/pi,20*log10(abs(H))); title( '20lg|H(eiw)|' );

频率抽样设计法

第7章 FIR 数字滤波器的设计方法 IIR 数字滤波器最大缺点：不易做成线性相位，而现代图像、语声、数据通信对线性相位的要求是普遍的。正是此原因，使得具有线性相位的FIR 数字滤波器得到大力发展和广泛应用。 1． 线性相位FIR 数字滤波器的特点 FIR DF 的系统函数无分母，为∑∑-=--=-== 1 1 )()(N n n N i i i z n h z b z H ，系统频 率响应可写成：∑-=-= 10 )()(N n jwn jw e n h e H ，令)(jw e H ＝)()(w j e w H Φ，H(w) 称为幅度函数，)(w Φ称为相位函数。这与模和幅角的表示法有所不同，H(w)为可正可负的实数，这是为了表达上的方便。如某系统频率响应 )(jw e H ＝w j we 34sin -，如果采用模和幅角的表示法，w 4sin 的变号相当 于在相位上加上)1(ππj e =-因，从而造成相位曲线的不连贯和表达不方便，而用)()(w j e w H Φ这种方式则连贯而方便。 线性相位的FIR 滤波器是指其相位函数)(w Φ满足线性方程： )(w Φ＝βα+-w （βα,是常数） 根据群时延的定义，式中α表示系统群时延，β表示附加相移。线性相位的FIR 系统都具有恒群时延特性，因为α为常数，但只有β＝0的FIR 系统采具有恒相时延特性。 问题：并非所有的FIR 系统都是线性相位的，只有当它满足一定条件时才具有线性相位。那么应满足什么样的条件？从例题入手。

例题：令h(n)为FIR 数字滤波器的单位抽样相应。N n n ≥<或0时h(n)=0，并假设h(n)为实数。 （a ） 这个滤波器的频率响应可表示为)()()(w j jw e w H e H Φ=（这是按幅 度函数和相位函数来表示的，不是用模和相角的形式），)(w H 为实数。（N 要分奇偶来讨论） （1） 当h(n)满足条件)1()(n N h n h --=时，求)(w H 和)(w Φ（π≤≤w 0） （2） 当h(n)满足条件)1()(n N h n h ---=时，求)(w H 和)(w Φ（π≤≤w 0） （b ） 用)(k H 表示h(n)的N 点DFT （1） 若h(n)满足)1()(n N h n h ---=，证明H(0)=0； （2） 若N 为偶数，证明当)1()(n N h n h --=时，H(N/2)=0。 解：（a ）∑-=-= 1 )()(N n jwn jw e n h e H （1）)1()(n N h n h --=，当N 为奇数时， +--++-+=---?----)11(1)1(0)11()1()1()0()(N jw jw N jw jw jw e N h e h e N h e h e H 2 123 ) 1()2 1(])[(---=-----++=∑N jw N n n N jw jwn e N h e e n h 2 1)2 1 ( 2 30 )2 1 ( )2 1 ()2 1(])[(-----=----- --++= ∑N jw N jw N n n N jw N n jw e N h e e e n h ) (})2 1 ()]21(cos[)(2{))(2 30 ) 2 1(w H e N h N n w n h e w j N n N jw Φ-=--=-+-- =∑

基于频率抽样法的FIR数字低通滤波器的设计

基于频率抽样法的FIR 数字低通滤波器的设计 1 设计目的 熟悉频率采样法的理论及其应用；掌握频率采样法设计FIR 数字滤波器的方法。了解FIR 数字滤波器的频率特性和相位特性，观察过渡带取样点对滤波器幅频特性的影响。掌握用频率采样法设计线性相位FIR 低通数字滤波器的方法，并掌握该方法的matlab 编程和仿真。 2 FIR 数字滤波器设计的原理 2.1频率抽样设计法 FIR 低通滤波器的设计一般方法有两种，即频率抽样法和窗函数法，频率抽样法设计不同于窗函数法，窗函数是从时域出发，把理想的()d h n 用一定形状得窗函数截取成有限长的()h n ，以此()h n 来近似理想的()d h n ，这样得到的频率响应()jw H e 逼近于所要求的理想的频率响应()jw d H e 。 频率抽样法则是从频域出发，把给定的理想频率响应()jw d H e 加以等间隔抽样，即2()| ()jw d d w k N H e H k π ==然后以此()d H k 作为实际FIR 数字滤波器的频率特性 的抽样值()H k ，即令2()()()| 0,1,,1jw d d w k N H k H k H e k N π====-，知道() H k 后，由DFT 定义，可以用频域的这N 个抽样值()H k 来唯一确定有限长序列()h n ，而由()X z 的内插公式知道，利用这N 个频域抽样值()H k 同样可求得FIR 滤波器的系统函数()H z 及频率响应()jw H e 。这个()H z 或()jw H e 将逼近()d H z 或 ()jw d H e ，()H z 和()jw H e 的内插公式为 1 1 01() ()1N N k k N z H k H z N W z ----=-= -∑ （2.1） 10 2()()()N jw k H e H k w k N π -==Φ- ∑ （2.2） 其中()w Φ是内插函数1() 2 sin( )12()sin() 2 N jw wN w e w N --Φ= （2.3） 将式（2.3）代入（2.2）式，化简后可得

送检取样

送检取样

一、水泥试验 1. 取样频率：同一水泥厂、同期出厂、同一出厂编号及同强度水泥（散装≤500吨/批、袋装≤200吨/批）。 2. 取样方法： 随机在不同部位至少20 袋水泥中均匀抽样经搅拌均匀后，取10㎏试样。 二、钢筋试验 1. 取样频率：同一生产厂家、同一炉批号、同一规格、级别、同一交货状态及同一进场时间≤60吨/批。 2. 取样方法：每批任意选取两钢筋切取两根用于拉伸试验，两根用于冷弯试验为一组，长度按不同试验室的检测仪器定。取样长度：拉伸试验一般为400～500mm，弯曲试验一般为350mm。 三、钢筋单、双面搭接焊及电渣压力焊试验 1. 取样频率：在现浇钢筋混凝土结构中，应以300个同牌号钢筋、同型式接头作为一批；在房屋结构中，应在不超过二楼层中300个同牌号钢筋、同型式接头作为一批，当不足300个接头时，仍应作为一批。 2. 取样方法：每批随机切取 3个接头为一组，长度按不同试验室的检测仪器定。取样长度：拉伸试验一般为450mm。

四、钢筋闪光对焊试验 1. 取样频率：同一台班内、由同一焊工完成的300个同牌号、同直径钢筋焊接接头作为一批。当同一台班内焊接的接头数量较少，可在一周内累计计算，累计仍不足300个接头时，仍按一批计算。 2. 取样方法：从每批接头中随机切取，3个接头抗拉 3个接头冷弯为一组。 五、钢筋机械连接试验 1. 取样频率：同型号、同规格接头，以500个为一批，不足500个接头仍作为一批。 2. 取样方法：每批随机切取 3个接头，另取 3条钢筋作为母材抗拉试验为一组，长度按不同试验室的检测仪器定。 六、蒸压加气混凝土砌块试验 1．取样频率：每1万块为一批，不足1万块仍作为一批。 2．取样方法：按膨胀方向中心分上、中、下抽取共9个试样为一组，每个试样平整切割成100×100×100mm。 七、普通混凝土小型空心砌块试验 1. 取样频率：每1万块为一批，不足1万块仍作为一批。 2. 取样方法：随机抽取5块为一组做强度试验。

采用频率采样法的FIR滤波器

吹管乐滤波去噪 ——基于频率采样法的FIR滤波器 学生姓名：焦阳指导老师：胡双红 摘要本课程设计主要内容是设计利用频率采样法设计一个FIR滤波器，对一段吹管乐进行滤波去噪处理并根据滤波前后的波形和频谱分析滤波性能。本课程设计仿真平台为MATLAB7.0，开发工具是M语言编程。首先在网上找到一段笛子独奏，加入一单频噪声，对信号进行频谱分析以确定所加噪声频率，设计滤波器进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。由分析结果可知，滤波器后的音频信号与原始信号基本一致，即设计的FIR滤波器能够去除信号中所加单频噪声，达到了设计目的。 关键词滤波去噪；FIR滤波器；频率采样法；MATLAB 1 引言 滤波去噪[1]是信号处理中一种非常基本但十分重要的技术。利用滤波可以从复杂的信号中提取所需的信号，一直不需要的信号。滤波器就是这样一种可以在时域和频域对信号进行滤波处理的系统。通常情况下，有用信号和干扰信号是在不同频段上的，于是通过对滤波器的频率特性精心设计就能达到滤波的目的。本课程设计是采用频率采样法设计频率抽样型滤波器，从而对吹管乐信号滤波去噪。通过对比滤波前后的波形图及回放滤波前后的吹管乐信号，来判断滤波器对噪声信号确实有滤除作用。 1.1 课程设计目的 （1）熟悉使用MATLAB； （2）了解FIR滤波器原理及结构； （3）利用所学数字信号处理想干知识用MATLAB设计一个FIR滤波器； （4）提高自己动手能力； （5）对加噪声的语音信号进行滤波去噪处理，比较滤波前后的时域波形和频谱并进行分析；

1.2 课程设计要求 （1）滤波器指标必须符合工程设计； （2）设计完后应检查其频率响应曲线是否满足指标； （3）处理结果和分析结论应该一致，而且应符合理论； （4）独立完成课程设计并按要求编写课程设计报告； 1.3 设计平台 本课程设计仿真平台为MATLAB7.0。MATLAB的名称源自Matrix Laboratory，1984年由美工Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛地应用于科学计算、控制系统、信号处理等领域的分许、仿真和设计工作。1993年MathWorks公司从加拿大滑铁卢大学购得MAPLE软件的使用权，从而以MAPLE为“引擎”开发了符号数学工具箱（Symbolic Math Toolbox）[2]。 2 设计原理 用网上找一段吹管乐，绘制波形并且观察其频谱，给定相应技术指标，用频率采样法设计的一个满足指标的频率采样型FIR滤波器，对该信号进行滤波去噪处理，比较滤波前后的波形和频谱进行分析。 2.1 FIR滤波器的设计 FIR（Finite Implse Response）[3]滤波器：有限长单位冲激响应滤波器，又称为非递归型滤波器，是数字信号处理系统中最基本的元件，他可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。 有限长单位冲激响应（FIR）滤波器有以下特点： （1）系统的单位冲激响应h（n）在有限个n值处不为0； （2）系统函数H（z）在|z|>0处收敛，极点全部在z=0处（因果系统）； （3）结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中（例如频率抽样结构）也包含有反馈的递归部分。 2.2 频率采样型结构 把一个有限长序列（长度为N点）的z变换H（z）在单位圆上作N等分抽样，就得

实验四 用频率取样法设计FIR数字滤波器

实验报告 哈尔滨工程大学教务处制

实验四 用频率取样法设计FIR 数字滤波器 一、实验目的 1、掌握频率取样法设计线性相位FIR 数字滤波器的方法，并用Matlab 工具编程实现。 2、熟悉频率取样理论，熟悉内插函数及其应用。 3、观察过渡带取样点或优化数值对滤波器幅频特性的影响。 二、 实验原理 频率采样法就是根据频域采样理论，由滤波特性指标构造希望逼近的滤波器频响函数H d (e jω)，对其在[0,2π]上采样得到。 ()() 20,1,,1j d d k N H k H e k N ωπ ω===-L 然后，就可求出单位脉冲响应h (n )，或是系统函数H (z )。这样，h (n )或是H (z )就是滤波器的设计结果。 ()()()()()1 100,1,,110,1,,1 1N N k k N h n IDFT H k n N H k z H z k N N W z ----===--= =--∑L L ()()() Frequency Sampling 2N 0,1,,1j j d d k H e H k H e k N ωωπ ω= ??????→==-L ()()() j k H k A k e θ= 三、 实验内容 1.用频率取样法设计一个线性相位低通数字滤波器，N=15，[0,π]之间的幅度取样值如下，求出其单位脉冲响应h[k]及幅频和相频特性曲线。尝试增加过渡点，观察并分析过渡点对滤波器性能的影响。 1, k 0,1,2[k]0.5, 30, H k =?? ==??? O t her s /3 1,()/30,d A ωπωπωπ

试验检测频率和取样方法.doc

试验检测频率和取样方法 试验项目检测频率取样方法 混凝土抗 压强度 抗渗性用于检查混凝土强度的试件，应在混凝土 浇筑地点随机取样： 1、每拌制 100 盘且不超过 3 一、同一组混凝土拌合物的取样应从同一盘或同 100m 的同配比 的混凝土，取样不得少于 1 次； 一车混凝土中取样。为使取样具有代表性，宜采2、每工作班拌制的同配比混凝土不足 用多次采样的方法。一般在同一盘混凝土或同一100 盘时，取样不得少于 1 次；车混凝土中约 1/4 处、1/2 处和 3/4 处分别抽取3、每次取样应至少留置 1 组标准养护试 大致相同份量的代表性样品，集中用铁铲翻拌均件，同条件养护试件的留置组数应根 匀，然后立即进行拌合物的试验。注意：取样数据实际需要确定。 量应多于试验所需数量的倍，其体积不小于 4、当一次浇筑连续浇筑超过1000 立，同 20L ，从第一次取样到最后一次取样不宜超过一配合比的砼，每200 立砼取样数量 15min 。 不少于一次（GB50204-2002 中中规 二、混凝土试样制备： 定） 根据混凝土拌合物的稠度确定混凝土成型方法： （ CJJ 2-2008 ） 1、坍落度不大于 70mm的混凝土宜用标准振动台混凝土数量小于250m3, 应制作抗渗试件 1 成型。 3 组抗渗试件。 将混凝土拌合物一次装入试模并稍有富余，装料组； 250～ 500m，应制作 2 注意：试件拆模后，应用钢丝刷刷净两端 时应用抹刀沿各试模壁插捣； 面的水泥浆膜。 试模应附着或固定在振动台上，振动时试模不得（ GB/T 50082-2009 ） 有任何跳动；

3 混凝土数量小于 250m, 应制作抗冻试件 1 振动应持续到混凝土表面出现乳状水泥浆时为 组； 250～ 500m3，应制作 2 组抗冻试件。止，振动过程中随时添加混凝土使试模常满； 注意 : 振动结束后，沿试模边缘刮除多余的混凝土，待 1、成型试件时，不得采用机油等憎水性混凝土试件表面收浆，临近初凝时，用抹刀将试 脱模剂（因为会显著影响试件的抗冻性件仔细抹平，试件抹面与试模边缘的高低差不得抗冻性能）。超过。 2、在标准养护室内或同条件养护的抗冻2、坍落度大于 70mm的混凝土宜用捣棒人工插捣 试件应在养护龄期 24d 之前送至实验室进成型。 行试验（因为进行抗冻试验前四天试件要混凝土拌合物分两层装入模内，每层的装料厚度 进行泡水处理）。大致相等。 （ GB/T 50082-2009 ）插捣应按螺旋方向从边缘向中心均匀进行。插捣 底层混凝土时，捣棒应到达试模底部，插捣上层 时，捣棒应贯穿上层后插入下层20～ 30mm；插 3 1、每 100m 的同配合比的混凝土，取样 1 捣时捣棒应保持垂直，不得倾斜，不得冲击。捣次；不足 100m3时按 1 次计； 10～ 完一层后，用橡皮锤轻轻击打试模外端面 混凝土弯拉2、每次取样应至少留置 1 组标准养护试 15 下，以填平插捣过程中留下的孔洞。 强度件，同条件养护试件的留置组数应根据实 每层插捣次数按在 100cm2截面积内不得少于12 际需要确定。 次 （ CJJ 1-2008 ） （《普通混凝土拌合物性能试验方法标准》 GB/T 50080-2002 ） 试验项目检测频率取样方法 1、同一炉号，同一牌号，同一规格，每 60t 为一检验批，超过60t 的部分，每增在不同根钢筋上截取，截去钢筋端部50cm， 加 40t 增加 1 个拉伸试验试样和一个弯曲取拉伸试样 2 根：长度不小于200+5dmm 钢筋：拉伸、 试验试样。取冷弯试样 2 根：长度不小于150+5dmm 冷弯 （ GB ）（GB ）备注： 28mm 32mm冷弯长度不小于5d+200mm 2、同牌号、同炉号、同规格、同交货状态 的钢筋，每60t 为一批，不足60t 也按一

基于频率抽样法的FIR高通数字滤波器的设计

信息工程系 课程设计报告书题目: 基于频率抽样法的FIR高通数字滤波器设计 专业：电子信息科学与技术 班级： 学号： 学生姓名： 指导教师： 年月日

信息工程系课程设计任务书

年月日

信息工程系课程设计成绩评定表

摘要 数字滤波是语音和图象处理、模式识别、频谱分析等应用中的一个基本处理算法，数字滤波技术是信号消噪的基本方法。根据噪声频率分量的不同，可选用具有不同滤波特性的数字滤波器。当噪声的频率高低于信号的频率时，就应该选用高通滤波器。本设计采用了频率抽样法设计的FIR高通数字滤波器，其目的是为了让中高频率的信号通过，而且利用频率抽样法的优点是可以在频域直接设计，并且适合最优化设计。FIR滤波器为有限长冲激响应滤波器，因其在线性相位特性方面具有独特的优点，因此也越来越受到广泛的重视。 关键词： FIR数字滤波器频率抽样高通 Abstract The digital filter is the voice and image processing, pattern recognition, spectrum analysis in the application of a basic algorithm, digital filtering technology is the basic method of signal denoising. According to different frequency components of noise, can be used with different filtering characteristics of digital filter. When the noise of high frequency below the signal frequency, should choose a high-pass filter. This design uses the frequency sampling design method of FIR high pass digital filter, its purpose is to make the high frequency signal through, and by using the frequency sampling method can be directly in the frequency domain design, and is suitable for optimization design. FIR filter for finite impulse response filters, due to its linear phase characteristic has unique advantages, so it attracts more and more attention. Key word:FIR digital filter Frequency sampling High pass

实验5 频率采样法设计FIR数字滤波器

实验5 频率采样法设计FIR 数字滤波器 1．实验目的 （1）掌握利用直接设计法实现频率采样法设计FIR 滤波器的基本原理。 （2）掌握利用最优设计法实现频率采样法设计FIR 滤波器的基本原理。 （3）掌握利用MA TLAB 语言实现直接设计法和最优设计法的方法。 2．实验原理 （1）FIR 数字滤波器设计的实质 求()j d H e ω的有限项傅里叶级数的系数，然后用有限项傅里叶级数去近似代替无限项傅里叶级数，在最小均方误差准则下最佳逼近()d h n 。 （2）用频率采样法设计滤波器的基本原理 设待设计的滤波器的传输函数用()j d H e ω表示，对它在0ω=到2π之间等间隔采样N 点，得到()d H k ，则 2()()N jw d d w k H k H e π==，0,1,,1k N =- 对N 点()d H k 进行IDFT ，得到()h n ， 210 1()()N N j kn d n h n H k e N π-==∑，0,1,,1n N =- ()h n 就是所设计的滤波器的单位取样响应，其系统函数1 0()()N n n H z h n z --==∑。 （3）直接频率采样设计法的基本原理 当FIR 滤波器满足第一类线性相位条件时，()h n 是实序列，()(-1)h n h N n =-，此时有： ()()()j j d g H e H e ωθωω=，1()2N θωω-=- （1）式 N=I ()(2)N=II ()(2) g g H H ωπωωπω=-???=--??g g 奇数时(型)：H 偶数时(型)：H 在0~2π之间等间隔频率采样N 点，采样时的频率转换关系为： 2k k N πω= ，0,1,2,,1k N =- （2）式 由（1）式可以得到采样点上的样本值为： ()2=()()()j j k d d g k N H e H k H k e ωθπω== （3）式 2121()()2k N N N k k k N N πωπθωθπ=--=-==-

常用建筑材料取样方法及频率解读

常用建筑材料、建筑构配件质量检测取样建筑材料见证取样方法

第一章建设部关于房屋建筑工程和市政基础设施工程实行见证取样和送检的规定第二章工程质量检测见证取样基础知识 第三章常用建筑材料质量检测取样方法 第一节：混凝土取样方法 第二节：建筑砂浆取样方法 第三节：水泥取样方法 第四节：建筑用砂、石取样方法 第五节：钢筋机械性能检测取样方法 第六节：墙体材料取样方法 第七节：建筑石灰取样方法 第八节：粉煤灰取样方法 第九节：混凝土外加剂取样方法 第十节：埋地排水管取样方法 第十一节：井盖、井算取样方法 第十二节：路面砖、路缘石取样方法 第十三节：土工合成材料 第十四节：沥青及沥青混合料取样方法 第十五节：土工及无机结合料试验取样方法 第十六节：混凝土、砂浆配合比检测 第十七节：工程结构混凝土质量检测（回弹法、钻芯法）

第一章建设部文件 房屋建筑工程和市政基础设施工程实行见证取样和送检的规定 建建【2000】211号 第一条：为规范房屋建筑工程和市政基础设施工程中涉及结构安全的试块、试件和材料的见 证取样和送检工作，保证工程质量，根据《建设工程质量管理条例》，制定本规定。 第二条;凡从事房屋建筑工程和市政基础设施工程的新建、扩建、改建等有关活动，应当遵守本规定。第三条：本规定所称见证取样和送检是指在建设单位或工程监理单位人员的见证下，由施工单位的现场试验人员对工程中涉及结构安全的试块、试件和材料在现场取样，并送至经过省级以上建设行政主管部门对其资质认可和质量技术监督部门对其计量认证的质量检测单位（以下简称“检测单位”）进行检测。 第四条：国务院建设行政主管部门对全国房屋建筑工程和市政基础设施工程的见证取样和送检工作统一监督管理。 县级以上地方人民政府建设行政主管部门对本行政区域内的房屋建筑工程和市政基础设施工程的见证取样和送检工作实施监督管理。 第五条:涉及结构安全和试块、试件和材料见证取样和送检的比例不得低于有关技术标准中规定应取样数量的30% 第六条:下列试块、试件和材料必须实施见证取样和送检： （一）用于承重结构的混凝土试块； （二）用于承重墙体的砌筑砂浆试块； （三）用于承重结构的钢筋及连接接头试件； （四）用于承重墙的砖和混凝土小型砌块;

频率采样法设计FIR数字滤波器

频率采样法设计FIR 数字滤波器 1．实验程序及运行结果 实验内容1：直接频率采样设计法 自定义一个能够产生已绝对值的幅度响应、相对dB 标尺的幅度响应、相位响应和群时延响应。函数名称定义为freqz_m 。 函数代码：freqz_m %功能扩展的频响特性函数freqz_m.m function [db,mag,pha,grd,w]=freqz_m(b,a) %H ：频响特性的样值向量 %w ：频响特性的位置向量 %db ：幅度衰减（dB ） %mag ：H 的模 %pha ：H 的相位 %grd ：H 的群时延 %[H,w]=freqz(b,a,1000,'whole'); %H=(H(1:1:501))';w=(w(1:1:501))'; [H,w]=freqz(b,a,1024,'whole'); H=(H(1:1:512))';w=(w(1:1:512))'; mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); grd=grpdelay(b,a,w);% grpdelay 群时延函数 设计分析： c k 取2c c k N ωπ≤ 的最大整数，因此=22c c k N ω π =，故频率采样值()g H k 为 15()1,1,1,0,,0,1,1g H k ?? =?????? 个零 源程序：shzxhchlshiyan8_1.m %shzxhchlshiyan8_1 %频率采样法的直接设计方法（FIR 低通滤波器） clear all;close all;clc;clf; M=20;tao=(M-1)/2;%M ：采样点数（20点），tao ：群时延（9.5） k=0:M-1;wl=(2*pi/M)*k;%w1描述各采样点频率值 Hrs=[1,1,1,zeros(1,15),1,1];%频域采样值（频率采样样本向量） %Hg(0)=Hg(1)=Hg(2)=Hg(18)=Hg(19)=1 %kc=(wc/2π)N,floor(kc)=2 Hdr=[1,1,0,0];wdl=[0,0.25,0.25,1];%Hdr 与wdl 共同确定1和0出现的频率位置边界 %当设置wdl=[0,0.25,0.25,1]时，是锐截止，无过渡带 %当设置wdl=[0,0.2,0.3,1]时，是设置频率在0.2π到0.3π之间的过渡带 k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;%（频率采样位置向量，0-9、10-19） %k2=ceil((M-1)/2):M-1 %floor:向小的整数方向取整，ceil ：向大的整数方向取整

频率采样法和等波纹

频率采样法 [y, fs, bits] = wavread( 'G:\hi' ); n = length (y) ; Y = fft(y, n); [x1]=wavread('G:\hi'); x1=x1+0.01*randn(141110,2); N=20; alpha=(N-1)/2; l=0:N-1; wl=(2*pi/N)*l; Hrs=[1,1,1,zeros(1,15),1,1]; %对理想幅度函数取样得到取样样本Hdr=[1,1,0,0];wdl=[0,0.25,0.25,1]; %用于绘制理想函数幅度函数的曲线k1=0:floor((N-1)/2); k2=floor((N-1)/2)+1:N-1; angH=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)]; H=Hrs.*exp(j*angH); %计算H(k) h=ifft(H,N); %计算h(n) w=[0:500]*pi/500; H=freqz(h,1,w); %计算幅度响应 [Hr,wr]=zerophase(h); %计算幅度函数 f1=filter(H,N,x1); figure(4) subplot(2,1,1) plot(x1) title('IIR 低通滤波器滤波前的时域波形'); subplot(2,1,2) plot(f1); title('IIR 低通滤波器滤波后的时域波形'); F0=fft(f1,8182); f=fs*(0:511)/8182; figure(5) y2=fft(x1,8182); subplot(2,1,1); plot(f,abs(y2(1:512))); title('IIR 低通滤波器滤波前的频谱') xlabel('频率/Hz');ylabel('幅值'); subplot(2,1,2) F1=plot(f,abs(F0(1:512))); title('IIR 低通滤波器滤波后的频谱') xlabel('频率/Hz');ylabel('幅值');

用频率采样法设计FIR数字滤波器

信号、系统与信号处理实验Ⅱ 实验报告 实验名称：用频率采样法设计FIR数字滤波器 一、实验目的 掌握频率取样法设计FIR数字滤波器，加深过渡点对滤波器性能影响的认识。 二、实验内容与要求 （1）编写好一个设计线性相位FIR高通滤波器的程序，已知wc=0.8 ，N=64，要求在屏幕上显示出h(n)值，画出|H(e^jw)|及20lg(|H(e^jw))的曲线。 （2）实验时，设置0个过渡点，1个过渡点，2个过渡点，比较设计所得的|H(e^jw)|及20lg(|H(e^jw))的曲线。 三、实验程序与结果 （1）0个过渡点 clear all; N=64; wc=0.8*pi; k=0:N-1; phase=(-pi*k*(N-1)/N)+pi/2; HK=[zeros(1,26),ones(1,13),zeros(1,25)]; HK1=HK.*exp(j*phase); hn=ifft(HK1,N) figure(1); freqz(hn,1,512); [H,W]=freqz(hn,1,512); figure(2); subplot(3,1,1);

stem(k,hn); title('h(n)') subplot(3,1,2); plot(W/pi,abs(H)); title('|H(eiw)|') subplot(3,1,3); plot(W/pi,20*log10(abs(H))); title('20lg|H(eiw)|'); （2）1个过渡点 clear all;

N=64; wc=0.8*pi; k=0:N-1; phase=(-pi*k*(N-1)/N)+pi/2; HK=[zeros(1,25),0.5,ones(1,13),0.5,zeros(1,24)]; HK1=HK.*exp(j*phase); hn=ifft(HK1,N); figure(1); freqz(hn,1,512); [H,W]=freqz(hn,1,512); figure(2); subplot(2,1,1); plot(W/pi,abs(H)); title('|H(eiw)|') subplot(2,1,2); plot(W/pi,20*log10(abs(H))); title('20lg|H(eiw)|'); （3）2个过渡点 clear all; N=64; wc=0.8*pi; k=0:N-1; phase=(-pi*k*(N-1)/N)+pi/2; HK=[zeros(1,24),1/3,2/3,ones(1,13),2/3,1/3,zeros(1,23)]; HK1=HK.*exp(j*phase); hn=ifft(HK1,N); figure(1);

利用MATLAB结合频率取样法设计数字高通FIR滤波器

课程设计任务书 学生姓名：专业班级： 指导教师：工作单位： 题目: 利用MATLAB结合频率取样法设计数字高通FIR滤波器 要求完成的主要任务: 1.利用频率取样法设计一个数字高通FIR滤波器 2.画出高通滤波器的幅频响应 课程设计进度安排 指导教师签名：年月日系主任（或责任教师）签名：年月日

目录 课程设计进度安排 ........................................................................................................ I 目录 .. (i) 摘要 ............................................................................................................................ I Abstract .......................................................................................................................... I I 1 FIR数字滤波器 (1) 1.1 FIR滤波器的特点 (1) 1.2 FIR数字滤波器设计方法 (2) 1.3 线性相位FIR数字滤波器的条件和特点 (2) 1.3.1 线性相位条件 (2) 1.3.2 线性相位FIR滤波器的幅度特性与相位特性 (3) 2 利用频率采样法设计FIR滤波器 (4) 2.1 用频率采样法设计滤波器的基本原理 (4) 2.2 线性相位的约束条件 (6) 2.3 逼近误差及其改进措施 (7) 2.3.1 产生误差的原因 (7) 2.3.2 减小误差的方法 (8) 2.4 频率采样法的特点 (9) 3 频率取样法的数字高通滤波器的实现 (9) 3.1 MATLAB的介绍 (9) 3.2 设计条件 (9) 3.3 设计程序 (9) 3.4 调试结果 (11) 4 心得体会 (12) 附录 (14)

见证取样制度及取样要求 数量及方法

见证取样送检制度 根据建设部建建［2000］211号文《关于印发（房屋建筑工程和市政基础设施工程实行见证取样和送检的规定）的通知》的要求，按江苏省建设厅《建设工程质量检测规程》DGJ32/J21-2009规定，执行见证取样送检制度。即在建设单位或监理单位人员见证下，由施工人员在现场取样，送至试验室进行试验。见证人员必须通过工程质量监督站的培训、考核合格，获得《见证人员证书》后持证上岗。 见证取样的范围：建设工程使用的原材料、半成品材料（构配件）及现场制作的混凝土、砂浆试块、钢筋焊接（连接）试件等应取样检测，取样应在见证人员旁站下进行，取样数量及方法应按相关技术标准、规范、规程的规定抽取。 见证取样送样程序： 每项工程的取样和送检见证人员，由该工程的建设单位书面授权，委派在本工程现场的建设或监理人员1～2名担任，并持有有效的《见证人员证书》。 建设单位应向工程质量监督站和工程质量检测单位递交《见证人员授权证书》。授权证书应写明本受监工程现场委托的见证单位和见证人员姓名、持证号。 施工企业每项工程的取样人员不少于2名，须持证上岗。施工企业取样人员在现场进行原材料取样和试块制作活动时，见证人员必须在旁见证。 见证人员应对试样进行监护，并和施工企业取样人员一起将试样送至检测单位。检测单位在接受委托检测任务时，须由送检单位填写委托单，见证人员应出示见证人员证书，并在测委托单上签名、盖章。检测单位均须实施密码管理制度。

建筑工程原材料、构配件检测取样要求、数量及方法 一.钢筋混凝土、砂浆用材料 1.1水泥 1.1.1检测参数：标准稠度用水量、凝结时间、安定性、细度、比表面积、胶砂强度、胶砂流动度。 （1）交货时水泥的质量验收可抽取实物试样以其检验结果为依据，也可以生产者同编号水泥的检验报告为依据。 （2）以抽取实物试样的检验结果为验收依据时，买卖双方应在发货前或交货地共同取样和签封。取样方法按GB12573进行，取样数量为20kg，缩分为二等份。一份由卖方保存40d，一份由买方按本标准规定的项目和方法进行检验。 （1）以同一生产厂家、同一品种、同一等级、同一批号且连续进场的水泥，袋装不超过200t、散装不超过500t为一批，每批抽样不少于一次。 （2）取样方法按GB12573进行。取样应有代表性，可连续取，亦可从20个以上不同部位取等量样品混合均匀，组成一组样品，密封。数量不少于12㎏。 1.2钢筋 屈服强度、抗拉强度、伸长率和断裂总伸长率、弯曲及反复弯曲、钢筋网片抗剪强度、钢筋重量偏差。 1.2.2取样频率、样品规格、数量及要求： （1）热轧光圆钢筋、热轧带肋钢筋、低碳钢热轧圆盘条：同一厂家、同一牌号、同一炉号、同一规格、同一进场时间,60t为一批,抽取一组试件。超过60t的部分，每增加40t（或不足40t的余数），增加一个拉伸试验和一个弯曲试验试样。每组试件取样方法：任选五根钢筋分别截取重量偏差5根（其中拉伸2根），长度500-550㎜，试件切口应平滑且与长度方向垂直；弯曲2根，长度300-350㎜。 （2）冷轧带肋钢筋：同一牌号、外形、生产工艺和同一交货状态为一批，每批不大于60t。拉伸试验每批逐盘检验，在每盘中随机截取1根试件，弯曲和反复弯曲试验每批随机抽取两盘各随机截取1根试件。 （3）钢筋网片：应按批进行检查验收，每批应有同一型号、同一原材料来源、同一生产设备并在同一连续时段内制造的钢筋焊接网组成，重量不大于60t。每种规格任取网片一张，不少于1㎡。 （4）冷拔低碳钢丝、冷轧扭钢筋、冷拉钢筋按产品标准的要求取样。 （5）若某项试验不合格，可从同一批钢筋中再任意抽取双倍数量的试样进行该不合格项目的复验。 1.3钢筋焊接