基于MATLAB的FIR数字滤波器设计

电子技术课程设计题目基于MATLAB的FIR数字滤波器设计学生姓名专业班级学号院（系）电气信息工程学院指导教师完成时间基于MATLAB的FIR数字滤波器设计郑州轻工业学院课程设计（论文）任务书题目简易窗帘自动开闭电路专业电子信息10-1学号541001030138姓名谢强主要内容、基本要求、主要参考资料等：主要内容1．阅读相关科技文献。2．学习protel软件的使用。3．学会整理和总结设计文档报告。4．学习如何查找器件手册及相关参数。技术要求要求电路能够通过感应装置，检测到光线的强弱。要求电路根据光线的强弱，自动将窗帘打开和关闭。要求电路能够在窗帘接触到边沿时，自动切断电源。主要参考资料1．何小艇，电子系统设计，浙江大学出版社，2001年6月2．姚福安，电子电路设计与实践，山东科学技术出版社，2001年10月3．王澄非，电路与数字逻辑设计实践，东南大学出版社，1999年10月4．李银华，电子线路设计指导，北京航空航天大学出版社，2005年6月5．康华光，电子技术基础，高教出版社，2003完成期限：2012年06月22日指导教师签章：专业负责人签章：2012年06月11日目录摘要 I1数字滤波器 数字滤波器数字滤波器简介数字滤波器是一个离散的系统。它可以对输入的离散信号进行一系列运算处理，从输入的信号中获得所需要的信息。数字滤波器的系统函数通常表示为(1-1)数字滤波器分为有限冲激响应数字滤波器，即FIR数字滤波器和无限冲激响应，即IIR数字滤波器。从公式的角度来看，FIR数字滤波器的始终为零；IIR数字滤波器至少有一个非零。实现数字滤波器的方法一般有两种：一种是利用计算机的程序编译，从而仿真实现；另一种是利用硬件来实现。实现一个数字滤波器一般需要三个基本的运算单元：加法器、延时器和乘法器。设计一个数字滤波器的一般步骤为：（1）按所给要求确定滤波器的性能（2）用一个因果稳定的离散线性时不变的系统函数逼近此性能的要求（3）利用算法来实现这个系统函数（4）利用计算机仿真或硬件来实现1.2FIR数字滤波器FIR滤波器是指在有限范围内系统的单位脉冲响应h[k]仅有非零值的滤波器。M阶FIR滤波器的系统函数H(z)为(1-2)其中H(z)是的M阶多项式，在有限的z平面内H(z)有M个零点，在z平面原点z=0有M个极点.FIR滤波器的频率响应为(1-3)它的另外一种表示方法为(1-4)其中和分别为系统的幅度响应和相位响应。若系统的相位响应满足下面的条件(1-5)即系统的群延迟是一个与没有关系的常数，称为系统H(z)具有严格线性相位。由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR滤波器设计中一般使用广义线性相位。如果一个离散系统的频率响应可以表示为(1-6)其中和是与无关联的常数，是可正可负的实函数，则称系统是广义线性相位的。如果M阶FIR滤波器的单位脉冲响应h[k]是实数，则可以证明系统是线性相位的充要条件为(1-7)当h[k]满足h[k]=h[M-k],称h[k]偶对称。当h[k]满足h[k]=-h[M-k],称h[k]奇对称。按阶数h[k]又可分为M奇数和M偶数，所以线性相位的FIR滤波器可以有四种类型。四种线性相位FIR滤波器的性质如表1-1所示表1-1四种线性相位FIR滤波器的性质类型=1\*ROMANI=2\*ROMANII=3\*ROMANIII=4\*ROMANIV阶数M偶数奇数偶数奇数h[k]的对称性偶对称偶对称奇对称奇对称关于的对称性偶对称偶对称奇对称奇对称关于的对称性偶对称奇对称奇对称偶对称的周期00A(0)任意任意00任意00任意可适用的滤波器类型LP,HP,BP,SPLP,BP微分器，变换器，Hilbert微分器，变换器，Hilbert，HP1.3IIR数字滤波器无限长单位冲激响应滤波器，即IIR数字滤波器具有下面几个特点：系统的单位冲激响应为无限长的系统函数在有限z平面上有极值点结构上是递归型的IIR滤波器的设计就是在给定的技术指标下去确定滤波器的阶数N和系数{，}。在已满足给定的技术指标下，应选用阶数尽可能低的滤波器，因为滤波器的阶数越低，在实现时成本就越低。在设计IIR滤波器时，最常用的方法是利用模拟滤波器来设计数字滤波器。其原因为：模拟滤波器的设计技术相对成熟，可以广泛利用模拟滤波器有大量的参考程序和表格它的解可以为闭合形式的2数字滤波器设计工具MATLAB2.1MATLAB简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。如图2-1是MATLAB的运行界面。图2-1MATLAB2011b运行界面MATLAB和Mathematica、Maple、MathCAD并称为四大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB主要功能：1.数值分析

2.数值和符号计算

3.工程与科学绘图

4.控制系统的设计与仿真

5.数字图像处理

6.数字信号处理

7.通讯系统设计与仿真

8.财务与金融工程2.2MATLAB的优势（1）高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来；（2）具有完备的图形处理功能，实现计算结果和编程的可视化；（3）友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；（4）功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等)，为用户提供了大量方便实用的处理工具。3FIR数字滤波器设计窗函数法设计FIR滤波器窗函数设计法又称为傅里叶级数法。这种方法首先给出，表示要逼近的理想滤波器的频率响应，则由IDTFT可得出滤波器的单位脉冲响应为(3-1)由于是理想滤波器，故是无限长序列。但是我们所要设计的FIR滤波器，其h[k]是有限长的。为了能用FIR滤波器近似理想滤波器，需将理想滤波器的无线长单位脉冲响应分别从左右进行截断。当截断后的单位脉冲响应不是因果系统的时候，可将其右移从而获得因果的FIR滤波器。另一种设计方案是将线性相位因子加入到理想滤波器的频率响应中，然后利用IDTFT计算出后，取在0≦k≦M范围的值为FIR滤波器单位脉冲响应。理想滤波器的频率响应和设计出的滤波器的频率响应的积分平方误差定义为(3-2)也可以表示为(3-3)上式中的第一项和第三项与所设计出的滤波器参数是没有关系的，为了使上式中的第二项达到最小，可选择(3-4)所以用上面的方法得出的滤波器是在积分平方误差最小意义下的最佳滤波器。Gibbs现象就是理想滤波器的单位脉冲响应截断获得的FIR滤波器的幅度函数在通带和阻带都呈现出振荡现象。随着滤波器阶数的增加，幅度函数在通带和阻带振荡的波纹数量也随之增加，波纹的宽度随之减小，然而通带和阻带最大波纹的幅度与滤波器的阶数M无关。窗函数的主瓣宽度决定了过渡带的宽度，窗函数长度N增大，过渡带减小。设计方法3.1.2各种窗函数（1）矩形窗(3-5)矩形窗的主瓣宽度为。用矩形窗设计的FIR滤波器过渡带宽度近似为。（2）三角形窗（又称巴特列特窗）(3-6)（3）汉宁（Hanning）窗（又称升余弦窗）(3-7)Hanning窗的主瓣宽度为。由Hanning窗的定义可知，Hanning窗在其两个端点的值为零，这就使得在实际的应用中不能利用两个端点的数据。我们可将N+2点的Hanning窗除去两个端点来定义长度为N的Hanning窗。修改后的长度为N的Hanning窗定义为(3-8)在Matlab信号处理工具箱中所采用的就是这种修改后的定义方式。海明（Hamming）窗（又称改进的升余弦窗）(3-9)Hamming窗的主瓣宽度为。布莱克曼（Blackman）窗（又称二阶升余弦窗）(3-10)Blackman窗的主瓣宽度为。频率取样法设计线性相位FIR滤波器频率取样法是从频域出发，在频域直接设计，把给定的理想频率响应加以等间隔取样，并以此作为实际FIR滤波器的频率响应。设所需滤波器的频率响应为。现要求设计一个M阶的FIR滤波器h[k]，使得在M+1个取样点上，FIR滤波器的频率响应与所需的频率响应相等，即(3-11)由设计的要求给定，h[k]需要通过设计来确定。如果M+1个方程是线性无关的，则可以通过求解M+1阶的线性方程来得出FIR滤波器的h[k]。对的一些特殊取样方法，上述方程的解可以直接由IDFT得到。由于要求设计出的滤波器是实系数的线性相位FIR滤波器，所以的取样值还需要满足线性相位滤波器的约束条件。=1\*ROMANI型和=2\*ROMANII型线性相位滤波器的，=3\*ROMANIII型和=4\*ROMANIV型线性相位滤波器的。为了使设计出的滤波器具有线性相位，在M+1个取样点上的值应为(3-12)下面分别讨论四种线性相位滤波器在取样点上的值：=1\*ROMANI型（M为偶数，h[k]偶对称）线性相位FIR滤波器在M+1个取样点值为(3-13)上式表明=1\*ROMANI型线性相位FIR滤波器在的值可由在的值确定。在的值确定后，对做M+1点的IDFT即可得到=1\*ROMANI型线性相位滤波器的h[k]。=2\*ROMANII型（M为奇数，h[k]偶对称）线性相位FIR滤波器在M+1个取样点值为(3-14)上式表明=2\*ROMANII型线性相位FIR滤波器在的值可由在的值确定。=3\*ROMANIII型（M为偶数，h[k]奇对称）线性相位FIR滤波器在M+1个取样点值为(3-15)上式表明=3\*ROMANIII型滤波器线性相位FIR滤波器在的值可由在的值确定。=4\*ROMANIV型（M为奇数，h[k]奇对称）线性相位FIR滤波器在M+1个取样点值为(3-16)上式表明=4\*ROMANIV型线性相位FIR滤波器在的值可由在的值确定。为了提高滤波器的质量并减少误差，可以采用人为地扩展过渡带的方法，即在频率相应的过渡带内插入一个或多个比较连续的采样点，使得过渡带比较连续，从而使得通带和阻带之间变法比较缓慢，使得设计得到的滤波器对理想滤波器的逼近误差较小。在理想低通滤波器的设计中，若不增加过渡点，阻带和通带之间的衰减约为-21dB，如果在通带和阻带之间增加一个采样点，阻带的最小衰减可以提高到-65dB，如果增加两个采样点，阻带的最小衰减可以提高到-75dB，如果增加3个采样点，阻带的最小衰减可以提高到-85dB至-95dB。线性相位FIR滤波器的优化设计在使用窗函数法设计滤波器时，如果使用矩形窗，则设计出的FIR滤波器是积分平方误差最小意义下的最佳FIR滤波器。但由于Gibbs现象的存在，使得设计出滤波器在阻带的衰减一般不能满足要求。为解决这个问题可以采用其他的窗函数，但得到的滤波器不是最小误差意义下的最优FIR滤波器。设表示要逼近的理想滤波器的幅度函数。表示设计出的线性相位FIR滤波器的幅度函数。优化设计的基本思想就是要确定线性相位FIR滤波器的系数，使得和的误差在范围内达到最小。与的加权误差函数定义为(3-17)其中，是由设计者定义的加权函数。一种常用的误差准则称为最大误差最小准则，又称切比雪夫准则，定义为(3-18)其中I表示滤波器在范围内各频率带区间构成的集合。切比雪夫准则意义下的优化设计就是要确定FIR滤波器H(z)的系数，使得上式定义的误差达到最小。在FIR滤波器的阶数M固定的条件下，利用切比雪夫误差准则设计出的FIR滤波器的幅度响应在通带和阻带都会呈现等波纹的波动。所以这类滤波器称为等波纹FIR滤波器。等波纹FIR滤波器的设计方法是由Parks-McClellan提出。在解决问题过程中使用了数学优化中的Remez交换算法。利用Remez交换算法，等波纹线性相位FIR滤波器的设计步骤可归结为：（1）用Kaiser提出的经验公式估计滤波器的阶数(3-19)确定J的值。（2）选定初始极值点。（3）计算。（4）为寻找新的极值点，计算函数在I中的抽样值。为了减少计算误差，抽样间隔应足够小。通常选择抽样点数为16M。（5）寻找新的极值点。（6）如果，执行（7）。否则用交换回到（3）。（7）获得g[k]。（8）由g[k]求出h[k]。4利用MATLAB设计FIR数字滤波器在利用Matlab设计FIR滤波器时，分别采用五种不同的窗函数设计所需的滤波器。4.1两个重要函数在窗函数法的Matlab实现中，程序中经常使用的函数有fir1和kaiserord。程序中fir1函数的用法：b=fir1(n,Wn,’ftype’,window)=1\*GB3①n为滤波器的阶数=2\*GB3②Wn为滤波器的截止频率，它是一个0到1的数。如果Wn是一个含有两个数的向量，则函数返回一个带通滤波器=3\*GB3③ftype为滤波器的类型，ftype=’high’时，设计的是高通滤波器；ftype=’stop’时，设计的是带阻滤波器；没有此参数时，设计的是低通滤波器=4\*GB3④window为指定的窗函数，矩形窗为boxcar(n)，汉宁窗为hanning(n)，海明窗为hamming(n)，布莱克曼窗为blackman(n)，凯撒窗为kaiser(n,beta)，没有此参数时，默认为hamming窗程序中kaiserord函数的用法：[n,Wn,beta,ftype]=kaiserord(f,a,dev,Fs)=1\*GB3①f是一个向量，为设计滤波器过渡带的起始点和结束点=2\*GB3②a是一个向量，指定频率段的幅度值=2\*GB3②dev是一个向量，长度和a相同，为各个通带和阻带内容许的幅度最大误差=4\*GB3④n为能够满足要求的滤波器的最小阶数=5\*GB3⑤Wn为滤波器的截止频率=6\*GB3⑥ftype为根据待设计滤波器的要求得到的滤波器的类型4.2设计不同参数的FIR滤波器4.2.1数字高通滤波器的实现为了用窗函数法设计下述参数的数字高通滤波器：，，通带最大衰减为1dB，阻带最小衰减为100dB，我用了五种窗函数截断，分别设计出五种不同的数字高通滤波器。程序见附录1，图4-1至图4-10是这五个滤波器的幅频响应曲线和相频响应曲线。图4-1布莱克曼窗截断所得高通滤波器的幅频响应曲线图4-2布莱克曼窗截断所得高通滤波器的相频响应曲线图4-3汉明窗截断所得高通滤波器的幅频响应曲线图4-4汉明窗截断所得高通滤波器的相频响应曲线图4-5汉宁窗截断所得高通滤波器的幅频响应曲线图4-6汉宁窗截断所得高通滤波器的相频响应曲线图4-7矩形窗截断所得高通滤波器的幅频响应曲线图4-8矩形窗截断所得高通滤波器的相频响应曲线图4-9三角窗截断所得高通滤波器的幅频响应曲线图4-10三角窗截断所得高通滤波器的幅频响应曲线4.2.2数字带通滤波器的实现同样的，为了用窗函数法设计下述参数的数字带通滤波器：，，，，通带最大衰减为1dB，阻带最小衰减为100dB，我用了五种窗函数截断，分别设计出五种不同的数字带通滤波器。程序见附录2，图4-11至图4-20是这五个滤波器的幅频响应曲线和相频响应曲线。图4-11布莱克曼窗截断所得带通滤波器的幅频响应曲线图4-12布莱克曼窗截断所得带通滤波器的相频响应曲线图4-13汉明窗截断所得带通滤波器的幅频响应曲线图4-14汉明窗截断所得带通滤波器的相频响应曲线图4-15汉宁窗截断所得带通滤波器的幅频响应曲线图4-16汉宁窗截断所得带通滤波器的相频响应曲线图4-17矩形窗截断所得带通滤波器的幅频响应曲线图4-18矩形窗截断所得带通滤波器的相频响应曲线图4-19三角窗截断所得带通滤波器的幅频响应曲线图4-20三角窗截断所得带通滤波器的相频响应曲线4.2.3数字带阻滤波器的实现同样的，为了用窗函数法设计下述参数的数字带阻滤波器：，，，，通带最大衰减为1dB，阻带最小衰减为100dB，我用了五种窗函数截断，分别设计出五种不同的数字带阻滤波器。程序见附录3，图4-21至图4-30是这五个滤波器的幅频响应曲线和相频响应曲线。图4-21布莱克曼窗截断所得带阻滤波器的幅频响应曲线图4-22布莱克曼窗截断所得带阻滤波器的相频响应曲线图4-23汉明窗截断所得带阻滤波器的幅频响应曲线图4-24汉明窗截断所得带阻滤波器的相频响应曲线图4-25汉宁窗截断所得带阻滤波器的幅频响应曲线图4-26汉宁窗截断所得带阻滤波器的相频响应曲线图4-27矩形窗截断所得带阻滤波器的幅频响应曲线图4-28矩形窗截断所得带阻滤波器的相频响应曲线图4-29三角窗截断所得带阻滤波器的幅频响应曲线图4-30三角窗截断所得带阻滤波器的相频响应曲线4.3不同窗函数设计滤波器的比较4.4FIR数字滤波器和IIR数字滤波器的比较(1)在技术指标相同的条件下，IIR滤波器的输出对输入有反馈，所以可以用比FIR少的阶数来满足要求，存储单元少，运算次数也少，经济实惠。（2）FIR滤波器的相位是严格线性的，而IIR滤波器做不到这一点，IIR滤波器的选择性越好，其相位的非线性越严重。（3）FIR滤波器主要采用非递归结构，有限精度的运算误差很小。而IIR滤波器在运算中会产生寄生振荡。（4）FIR滤波器可以使用快速傅里叶变换算法，而IIR滤波器不能这样。（5）IIR滤波器可以利用模拟滤波器的公式、数据和表格，计算量小。FIR滤波器设计时往往要借助计算机。5总结参考文献[1]程佩青.数字信号处理教程（第三版）[M].北京:清华大学出版社，2007,333-342.[2]吴大正.信号与线性系统分析（第四版）[M].北京:高等教育出版社，2005,63-427.[3]奥本海姆.A.V.信号与系统.第二版.刘树棠译.西安：西安交通大学出版社，2002[4]OppenheimerAV,SchaferRW.DigitalSignalProcessing.PrenticHall，Inc.1975[5]GoldB,RaderCM.DigitalProcessingofSignal.McGraw-HillBookCompany,1969[6]王世一.数字信号处理[M].北京：北京理工大学出版社，2005.[7]张雄伟，陈亮，许光辉.DSP芯片的原理与开发应用[M].北京：电子工业出版2003,2[8]姚天任，孙洪，现代数字信号处理[M].武汉：华中科技大学出版社，1999,11[9]戴明祯，周建江.TMS320C54XDSP结构、原理及应用[M].北京：北京航空航天大学出版社，2002,3附录附录1：StopdB=100;PassdB=1;dev=[0.000011];fcuts=[0.40.6];a=[01];[n,wn,bta,ftype]=kaiserord(fcuts,a,dev);%ÓÃkaiserordº¯Êý¹À¼Æ³öÂ˲¨Æ÷½×ÊýnºÍbeta²ÎÊýwindow=triang(n+1);h1=fir1(n,wn,ftype,window,'noscale');[hh1,w1]=freqz(h1,1,512);figure(1)plot(w1/pi,20*log10(abs(hh1)))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('·ù¶È/db');title('Â˲¨Æ÷µÄÔöÒæÏìÓ¦');figure(2)plot(w1/pi,angle(hh1))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('Ïàλ/rad');title('Â˲¨Æ÷µÄÏàλÏìÓ¦');附录2：StopdB=100;PassdB=1;dev=[0.0000110.00001];fcuts=[0.20.30.70.8];a=[010];[n,wn,bta,ftype]=kaiserord(fcuts,a,dev);%ÓÃkaiserordº¯Êý¹À¼Æ³öÂ˲¨Æ÷½×ÊýnºÍbeta²ÎÊýwindow=triang(n+1);h1=fir1(n,wn,ftype,window,'noscale');[hh1,w1]=freqz(h1,1,512);figure(1)plot(w1/pi,20*log10(abs(hh1)))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('·ù¶È/db');title('Â˲¨Æ÷µÄÔöÒæÏìÓ¦');figure(2)plot(w1/pi,angle(hh1))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('Ïàλ/rad');title('Â˲¨Æ÷µÄÏàλÏìÓ¦');附录3：StopdB=100;PassdB=1;dev=[10.000011];fcuts=[0.20.30.70.8];a=[101];[n,wn,bta,ftype]=kaiserord(fcuts,a,dev);%ÓÃkaiserordº¯Êý¹À¼Æ³öÂ˲¨Æ÷½×ÊýnºÍbeta²ÎÊýwindow=boxcar(n+1);h1=fir1(n,wn,ftype,window,'noscale');[hh1,w1]=freqz(h1,1,512);figure(1)plot(w1/pi,20*log10(abs(hh1)))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('·ù¶È/db');title('Â˲¨Æ÷µÄÔöÒæÏìÓ¦');figure(2)plot(w1/pi,angle(hh1))gridxlabel('¹éÒ»»¯ÆµÂÊw');ylabel('Ïàλ/rad');title('Â˲¨Æ÷µÄÏàλÏìÓ¦');基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统