基于MATLAB的信号波形与频谱分析

基于MATLAB的信号波形与频谱分析

辽宁科技大学毕业设计（论文）第II页

基于MATLAB的信号波形与频谱分析

摘要

本文利用MATLAB软件进行信号频谱分析的设计，并通过GUI界面(图形用户界面)实现信号频谱分析的动态设计。用户可与计算机交互式地进行对象参数的设置、控制算法的选取、以及对离散信号与连续信号的动态频谱分析。并利用MATLAB内嵌的Simulink模块，实现系统的信号频谱分析，满足不同用户的不同要求。

MATLAB的GUIDE是专门用于图形用户界面（GUI）程序设计的快速开发环境，本文从介绍GUIDE入手结合具体的软件图形界面实例，给出了利用GUIDE制作图形用户界面的基本方法。通过本文的研究得出利用GUIDE设计的仿真界面不仅可以对离散信号与连续信号的动态频谱进行相关分析，而且还可以对此问题作进一步的预测和分析。基于GUIDE技术制作的图形用户界面，具有友好性，开放性，方便科研人员不断地研究和扩充。

关键词频谱分析；GUI；SIMULINK；MATLAB

辽宁科技大学毕业设计（论文）第III页

Abstract

Using MATLAB software for signal spectrum analysis of the design, and through the GUI interface (graphical user interface) signal spectrum analysis of the dynamic design. Users can interactively with the computer targeting parameters for the establishment, control algorithm selection, and the discrete signals and continuous signal analysis of the dynamic spectrum. MATLAB and Simulink embedded module, the signal system to achieve spectral analysis, meet the needs of different users with different needs.

MATLAB GUIDE is devoted to the graphical user interface (GUI) design procedures for the rapid development environment, This paper introduced GUIDE start from the context of the specific examples of graphical interface software, GUIDE is produced using a graphical user interface methods. Through this paper, the study of the use of simulation GUIDE interface design can not only right for discrete signal and the signal dynamic spectrum analysis, but also the issue for further analysis and forecasts. GUIDE technology based on the production of graphical user interface is friendly and open, facilitate researchers continue to study and expanded.

Keywords Spectrum Analysis ；GUI；SIMULINK；MATLAB

辽宁科技大学毕业设计（论文）第IV页

目录

摘要 .....................................................................................................I Abstract ................................................................................................ II 第1章绪论 . (1)

1.1 MATLAB应用软件简介 (1)

1.2课题研究的内容 (2)

1.2.1信号的频谱分析 (2)

1.2.2信号频谱分析的原理 (3)

1.2.3频谱分析与显示原理 (3)

1.2.4对三阶交调进行模拟分析 (3)

1.2.5伪彩色编码的动态频谱图显示 (4)

1.2.6宽带频谱图与窄带频谱图 (5)

1.3 GUI界面设计实现 (5)

1.4 SIMULINK (7)

1.4.1 SIMULINK模型与文件 (7)

1.4.2 SIMULINK仿真原理 (7)

第2章MATLAB工具箱的概述及GUI界面的简介 (9)

2.1引言 (9)

2.2句柄图形及图形用户界面 (9)

2.3 GUIDE (9)

2.4编制回调程序 (11)

2.5主界面 (11)

2.6子界面 (12)

2.7基于MATLAB信号处理工具箱的数字滤波器设计与仿真 (12)

2.8带通滤波器设计实例 (13)

第3章GUI界面下噪声信号的频谱分析与实现 (15)

3.1概述 (15)

3.2 MATLAB信号处理工具箱介绍及功能 (15)

3.2.1MATLAB工具箱介绍 (15)

辽宁科技大学毕业设计（论文）第V页

3.2.2工具箱功能 (16)

3.3通风机噪声信号的采集方法 (18)

3.4噪声信号的频谱分析 (19)

3.5GUI界面下噪声频谱分析的实现 (20)

结论 (24)

致谢 (25)

参考文献 (26)

附录A (27)

附录B (31)

辽宁科技大学毕业设计（论文）第1页

第1章绪论

1.1 MATLAB应用软件简介

MATLAB是一种功能强大的科学计算和工程仿真软件，它的交互式集成界面能够帮助用户快速地完成数值分析、矩阵运算、数字信号处理、仿真建模、系统控制和优化等功能。

MATLAB将高性能的数值计算和可视化集成在一起，提供了大量的内置函数，从而被广泛地应用在科学计算、控制系统、信息处理等领域的分析、仿真和设计工作，而且利用MATLAB 产品的开放式结构，可以非常容易地对MATLAB 的功能进行扩充，从而在不断深化对问题认识的同时，不断完善MATLAB 产品以提高产品自身的竞争能力。

MATLAB语言采用与数学表达式相同的形式，不需要传统的程序设计语言，因而不像其他高级语言那样难于掌握。一般来说，用户可以在极短的时间内掌握MATLAB解决简单的问题，由于MATLAB的这些特点，它现在已经成为科研工作和工作仿真的高级助手。

MATLAB是mathworks公司开发的，目前国际上最流行应用最广泛的科学与共程计算软件，它广泛应用于自动控制、数字运算、信号分析、航天工业、汽车工业、计算机技术、图像信号处理、财务分析、生物医学工程、语言处理和雷达工程等各行各业，也是国内外高校和研究部门进行许多科学研究的重要工具。

MATLAB6.5是它的最新版本，由于它具有强大的计算和绘图功能，大量稳定可靠的算法和简洁高校的编程语言，所以选用MATLAB6.5环境中的GUI介面、SIMULINK仿真环境来对信号进行频谱分析。

本文将以美国Metso Automation MAX公司的新一代DCS产品MAX1000+plus为例讨论如何利用MATLAB的仿真软件包SIMULINK中的S-FUNCTION［1］和MAX1000+plus系统中读写DPU参数的Software Backplane(软件背板)［2］技术实现虚拟DPU控制作用与MATLAB建立的仿真模型的连接，并用MATLAB6.5的GUI介面来进行通风机噪声频谱分析。

辽宁科技大学毕业设计（论文） 第2页

1.2课题研究的内容

1.2.1信号的频谱分析

（1）连续周期信号

已知一个连续周期信号，在满足一定条件下，可以通过傅立叶级数展开为一系列正弦信号的线性叠加。形状不同的周期信号其区别在于它们各自不同基本频率、不同谐波的幅度和相位的正弦信号组合而成。由此可见，如果能找到具有这三个特性的频率函数，就能把该信号惟一地确定下来。从傅立叶级数展开式中可以观察到，傅立叶系数0()x kw 就是这样的频率函数，并可按下式求得：

000

20021()()t T jkw T x kw x t e dt T -=? (1.1)

式中，T 是周期信号的周期，002/w T π=是基本角频率，K 是离散变量表示谐波的次数。可见,0()x kw 是离散频率的复函数，可表示为：

0()00()|()|i kw X kw X kw e φ= (1.2)

（2）连续非周期函数

对于连续非周期信号，也可通过连续时间傅立叶变换（CTFT ），从上式中求得一个频率函数X （W ），在频域对信号进行分析，即

()()jwt x w x t e dt ∞

--∞=? (1.3)

（3）离散周期信号

对于离散周期信号，从离散傅立叶级数展开式中可以求得傅立叶系数0()x k Ω，如下式所示，即：

01001()()N jkn n x k x n e N

--Ω=Ω=∑ (1.4)

（4）离散非周期信号 对于离散周期信号。如同连续非周期信，通过离散时间傅立叶变换（DTFT ）可求得非周期序列的频谱密度函数()x Ω，即：

辽宁科技大学毕业设计（论文） 第3页

130()()n N n x x n e

--Ω=Ω=∑ (1.5)

1.2.2信号频谱分析的原理

信号频谱分析，就是应用傅里叶分析的方法，求出与信号时域波形相对应的频率函数，从中找出描述该信号频谱结构的变化规律，以达到特征提取的目的，如频带宽度、幅度、相位以及能量、功率等随频率变化的分布规律[3]。

1.2.3频谱分析与显示原理

（1）离散信号的短时FFT 频谱分析

现代信号频谱分析是基于离散时域的短时FOURIER 分析。设离散时域采样信号为x （n ）=0,1,…N -1是时域采样点序号，N 是信号长度。在数字信号处理中，用加窗方法将信号分段。此时()m x n 表示成x m x (n),n=0,1,…,N -1，m 是桢的序号，n 是桢的同步的时间序号，N 则是一桢内的采样点数（桢长）。信号()m x n 的加窗离散时域Fourier 变换（DTFT ）为：

1

0(,)()()N jw

jwn m m n X m e w n x n e --==??∑ （1.6） 为了便于离散计算，采用()()m m w n x n ?的离散FOURIER 变换（DFT ）：

1

2/20(,)()(),0,...,1N j nk N m m n X m k w n x n e k N π--==??=-∑ （1.7）

X(m,k)|便是Xm(n)的短时幅度谱估计。把m 当作时间变量，k 当作频率变量，则|X(m,k)|就是信号x(n)的动态频谱。由20log10(|X(m,k)|)就是以dB 表示的动态频谱。为了高计算|X （m,k ）|，可以运用快速Fourier 变换算法，即FFT 算法。

1.2.4对三阶交调进行模拟分析

1.输入：三个正弦信号之和，要求三个信号之间的频率相差不大。

2.过程：输入信号通过非线性放大器，放大器的放大特性为a*x+b*x*x*x 。

3.输出：输入信号和输出信号的频谱［4］。

fs=6400000; %采样频率64k

t=0:(1/fs):(0.1-(1/fs)); %时长0.001秒

辽宁科技大学毕业设计（论文） 第4页

t_len=size(t,2);

f1=1400;

f2=1600; %双音信号1k 和1.6k

f_len=linspace(-fs/2-fs/2/t_len,fs/2-fs/2/t_len,t_len);

figure(1);

a1=1.65;

a3=-0.887;

a5=0.16;

s=(sin(2*pi*f1*t)+sin(2*pi*f2*t));

a*x+b*x*x*x

fftmov1=20*log10(abs(fft(ss)));

fftmov1=fftshift(fftmov1);

plot(f_len,fftmov1);

ylim([-100 120]);

xlim([0 5000]);

grid on;

xlabel('频率hz');

ylabel('幅度db');

title('输出信号频谱');

1.2.5伪彩色编码的动态频谱图显示

在二维平面坐标上，将m 作为横坐标，k 作为纵坐标，把|X （m,k ）|的值作为灰则得到X(n)的动态频谱图。

对|X(m,k)|的值实施伪彩色映射，可得到伪彩色的动态频谱图，以获得更好的分辨率和视觉效果，并增强频谱图的可读性性。把|X(m,k)|的最小值min X 映射为归一化的0电平，最大值max X 映射归一化的1电平，而将|X （m,k ）|线性映射为01之间的电平值i C 。再利用彩色映射，计算彩色监视器把i C 显示出来，则获得伪彩色显示的动态频谱图。为了充分利用彩色空间的范围，选折适当的基准频谱值Base 。小于Base 的频谱值都被限定在基准电平上，而大于Base 的谱值则被线

辽宁科技大学毕业设计（论文） 第5页

性映射到由01的归一化彩色值。以C={c(m,k)}表示彩色值矩阵，由|X （m,k ）|到c (m,k)的映射的数学表示如下：

(,)(,)[((,)(,))]B m k Base c m k Mas B M K m k Base

-=?- （1.8）

其中 |(,)|,|(,)|(,),|(,)|x m k x m k Base c m k Base x m k Base >??=??≤??

（1.9） 1.2.6宽带频谱图与窄带频谱图

根据DFT 分析原理，频谱的频域分辨力是指离散的频率的频率间隔，也即X( )中相邻的k 值之间所代表的频率间隔 0f 。 0f 的大小取决于桢长N 以及信号采样频率s f 。根据Nyquist 采样定理，0s f f 和N 之间的关系满足

|(,)|,|(,)|(,),|(,)|x m k x m k Base c m k Base x m k Base >??=??≤??

(1.10) 这与信号n x 原来所包括的频率成分无关。所以，s f 在不变的前提下，增加N 将使频谱图的分辨率提高或相邻k 值之间所代表的频带宽度变小，这就使得频谱趋向于窄带频谱，相反，则趋向于宽带频谱。

为提高频域分辨率而加大N 值，会使频谱在时域中是的分辨率降低。解决这个矛盾的途径是，在取较大（合适）的真长N 的同时，采用信号分针针移1N 小于真长N 的办法，由此可得到频域和时域分辨率都满足需要的频谱。分真真移可表示为：

1()()m x n x n N m =+? n=0,…N -1, 1N

11)

1.3 GUI 界面设计实现

动态数字调节器设计与仿真软件的界面。它共包含5个区域：

(1)菜单区：位于整个界面的最上端。隐含了Matlab 自带的所有基本菜

辽宁科技大学毕业设计（论文）第6页

单(将figure对象的“MenuBar”属性设置为“None”)，添加了file、tool、system、最小拍、大林算法、Simulink、help7个菜单。每个菜单下又有各级子菜单。

其中两项的具体内容如图1.1。

图1.1菜单内容

(2)图形区：用于显示各计算机控制系统的动态仿真曲线，以及各系统的比较曲线。

(3)对象模型区：动态的显示系统各相应环节的数学模型(传递函数或脉冲传递函数)。由4个单选按钮(radiobutton)和4个编辑框(edit)组成。当用户的输入参数发生变化时,可通过单击响应的单选按钮，实现显示结果的实时刷新与显示。

(4)控制按钮区：由4个按钮组成(pushbutton)，实时的进行系统参数的设定，输出图像智能切换，并具有曲线刷新等功能。在使用该软件时,第一步便是按下“参数输入”按钮，进行系统设置,否则将给出错信息。

(5)图形效果区：由一个滑动条组成(slider)。用于控制输出图形的明暗程度。其实在该区域还可根据用户要求加入其他图形控件，如色彩选择器、线形控制器

辽宁科技大学毕业设计（论文）第7页

等等，但就实现方式来讲是完全一致的。故本窗口只设置了一个以示说明。

1.4 SIMULINK

SIMULINK是一个进行动态系统建模、仿真和综合分析的集成软件包。它可以处理的系统包括：线性、非线性系统；离散、连续及混合系统；单任务、多任务离散事件系统。

1.4.1 SIMULINK模型与文件

在SIMULINK模型窗口中可以使用“File”选单来实现模型的文件相应的操做，同样，也可以在MATLAB的命令窗口通过命令来实现。

（1）创建新模型。new-systen命令用来在MATLAB的工作空间创建一个空白的SIMULINK模型，但这个模型不能制动显示出来，称为逻辑模型，必须用open-systen命令来显示。

（2）打开模型。open-systen命令用来打开逻辑模型，在SIMULINK模型窗口显示该模型。

（3）保存模型。Save-systen命令用来保存模型为模型文件夹，扩展名为.mdl。

（4）关闭模型。如果要关闭一个打开的SIMULINK模型，则使用close-systen 命令。

SIMULINK是MATLAB最重要的组件之一，它提供一个动态系统建模仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。SIMULINK具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点SIMULINK已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于SIMULINK。

1.4.2 SIMULINK仿真原理

简单的：用SIMULINK的工具箱做仿真。

中级的：确定一个硬件芯片，然后再加外围，然后仿真，适当的考虑一些其他

干扰、现场因素。

辽宁科技大学毕业设计（论文）第8页

复杂的：确定硬件芯片，加外围，考虑现场、干扰、编码等诸多因素，并将整个信号的传输过程都仿真出来，如采集端的DSP，信道里面的传输，接受的解码编码，最好是用一段具体的声音来建模，并给出在频域、时域的图像[6]。如果在家一点儿信道服用以后的提高、优化分析，如时分复用、频分复用等。

辽宁科技大学毕业设计（论文）第9页

第2章MATLAB工具箱的概述及GUI界面的简

介

2.1引言

MATLAB是由MathWorks公司于1984年推出的数学软件，早期以矩阵运算为主,但随着科学可视化的需求日增，MATLAB在第4版推出了句柄图形，自此之后，所有的Demos都包含友好的图形用户界面，MATLAB的图形界面设计功能也日益完善如今，MATLAB已成为集数值计算功能符号运算功能和图形处理功能为一身的超级科学计算语言，满足了不同领域用户的需求。

2.2句柄图形及图形用户界面

句柄图形(HandleGraphics)就是将一个图形的每一个组件都视为一个对象(Object)，每一个对象都有一个独一无二的/句柄0(handle)，根据这个句柄,就可以找到这个对象(即图形组件)的各项属性，并进而更改这些属性，以产生不同的图形呈现效果。

句柄图形属性：Children属性、Parent属性、Tag属性、Type属性、UserData 属性、Visible属性、ButtonDownFcn属性、CreateFcn属性、DeleteFcn属性。

图形用户界面GUI(GraphicsUserInterface)是由各种图形对象，如图形窗口，图轴，菜单，按钮，文本框等构建的用户界面,是人机交流信息的工具和方法,在该界面内，用户可以根据界面提示完成整个工程，却不必去了解工程内部是如何工作的GUI设计即可以基本的MATLAB程序设计为主,也可以鼠标为主利用GUIDE工具进行设计。利用GUIDE设计图形用户界面时，可通过GUI应用属性设置编辑器来设置对句柄操作的响应,findobj命令可以获得所需对象的句柄。2.3 GUIDE

GUIDE(GraphicsUserInterfaceDesignEnviron2ment)是一个专用于GUI程序设计的快速开发环境，使用者通过鼠标就能迅速地产生各种GUI控件,并随心所欲地改变它们的外形!大小及颜色等，从而帮助用户方便地设计出各种符合要求的图形用户界面调用GUIDE的方法有2种,在MATLAB命令窗口中输入guide命令。

辽宁科技大学毕业设计（论文）第10页

在GUIDE设计环境中，需要用到的工具有属性编辑器控件布置编辑器菜单编辑器对象浏览器网格标尺设置编辑器以及GUI应用属性设置编辑器等。

(1)属性编辑器(PropertiesInspector)

打开属性编辑器的方法有多种，最常用的,选中控件对象以后，点击工具条上的属性设置按钮，即可打开属性编辑器"通过该属性编辑器来对所选图形对象设置相关属性。

(2)控件布置编辑器(AlignmentObjects)

在编辑CUI过程中，通过控件布置编辑器可以方便地设置面板上GUI控件的布局选中需要对齐的对象，然后选择工具条上的控件布置按钮，即可打开控件布置编辑器在控件布置编辑器中可以设置GUI控件水平以及垂直布局，包括对齐方式以及控件间距等。

(3)对象浏览器(ObjectBrowser)

在GUI面板中点击对象浏览按钮可打开对象浏览器，在该浏览器中可以方便地显示出所有的图形对象，单击该对象则可以打开相应的属性编辑器。

(4)菜单编辑器(MenuEditor)

在GUI面板中点击菜单编辑按钮可打开菜单编辑器，通过菜单编辑器可以为图形界面添加设置和修改菜单项，以及为右键添加快捷菜单。

(5)网格标尺设置编辑器(GridandRulers)

通过网格标尺设置编辑器，可以在GUI面板中添加网格以及标尺,来方便用户的界面设计。选择GUIDE主菜单中的Tools/GridandRulers即可打开网标尺设置编辑器。

(6)GUI应用属性设置编辑器(GUIDEApplica2tionOptions)

通过该编辑器可以设定GUI界面缩放形式，GUI对句柄操作的响应方式，以及定义GUI的保存方式。选择GUIDE主菜单中的Tools/GUIOptions即可打开GUI应用属性设置编辑器。

(7)GUI设计面板

GUI设计面板是上述GUI设计工具应用的平台，如图1所示，面板上部提供了菜单和常用工具按钮,左边提供了多种GUI控件。如按钮，单选按钮，复选框，文本框等。进行GUI设计时,，首先单击GUI面板左边所需的控件,然后在右

辽宁科技大学毕业设计（论文）第11页

边的图形界面编辑区中再次单击某一恰当的位置，这时将在该位置上为图形界面添加一相应的控件,接下来,通过属性编辑器和对齐编辑器对各控件设置相关属性和进行界面布置,以完善界面功能。

2.4编制回调程序

完成控件的布局之后，也就完成了整个图形界面的结构设计，接下来则是最为重要的功能设计，即要编制菜单［5］。控件的回调程序"若要编制某一控件的回调程序,用右键单击该控件,在弹出菜单中点击ViewCallbacks,然后从子菜单中选择一种激活回调程序的方式，就可以编制回调程序了，如Callback指单击控件时激活回调程序完成一定的功能菜单的回调程序在菜单编辑器的Callbacks文本框中给出。

2.5主界面

在MATLAB命令窗口中输入guide打开GUI的快速开发环境GUIDE，如图1所示本软件的主界面包括一个命令按钮，一个位图格式的图片和几个静态文本框，要改变最上边静态文本框的属性，需首先选中该控件，然后点击工具条上的属性设置按钮打开属性编辑器，当然也可通过其他的方式打开属性编辑器，然后通过该编辑器为静态文本框设置相关属性String属性用来设置控件对象所显示的内容，这里设置为-基于MATLAB的图形用户界面，字体角度属性FontAngle 设置为-italic，即斜体，字体大小属性FontSize设置为30，ForegroundColor属性用来设置文本框的前景颜色,这里设置为黑色,控件的位置和大小可通过改变Position属性中的值来调整,也可在GUI面板中用鼠标拖动控件来调整控件的位置和大小其他控件的属性也都是如此设置的另外，当有控件需要对齐时，可以通过对齐编辑器进行修正。

完成控件布局之后,开始编制回调程序如前面所讲的，在GUI编辑状态下用鼠标右键单击需要编写回调程序的控件,打开右键快捷菜单,从ViewCallbacks中选择一种回调方式，打开MATALB自动生成的.m文件，并且系统会自动生成如下语句FunctionHnameCallback (hObject,eventdata,han2dles)

其中，Hname为控件对象的Tag(标签)属性值，开发者即可在该语句下编写当前

辽宁科技大学毕业设计（论文）第12页

控件的回调程序，Callback是其中的一种回调方式，根据需要可选择其他回调方式，如WindowButtonDownFcn,CreateFcn等。

至此完成了主界面的设计，运行后"MATLAB将把该图形界面保存于文件cansai.fig和cansai.m中，其中，图形文件将保存GUI面板以及相应控件菜单设计，M文件则保存程序代码。

2.6子界面

按照同样的方法编制子界面，该界面内除了包括文本框，图轴，命令按扭外，还包括2个弹出式菜单和5个列表框2个弹出式菜单分别用来选择拟合阶数和显著性水平，5个列表框则用来显示回归分析的各项结果。

2.7基于MATLAB信号处理工具箱的数字滤波器设计与仿真

阶数（N阶滤波器，Specify Order＝N-1），如果选择Minimum Order则MATLAB根据所选择的滤波器类型自动使用最小阶数。Frenquency Specifications选项，可以详细定义频带的各参数，包括采样频率Fs和频带的截

止频率。它的具体选项由Filter Type选项和Design Method选项决定，例如Bandpass（带通）滤波器需要定义Fstop1（下阻带截止频率）、Fpass1（通带下限截止频率）、Fpass2（通带上限截止频率）、Fstop2（上阻带截止频率），而Lowpass（低通）滤波器只需要定义Fstop1、Fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Magnitude Specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义Wstop1（频率Fstop1处的幅值衰减）、Wpass （通带范围内的幅值衰减）、Wstop2（频率Fstop2处的幅值衰减）。当采用窗

函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。Window Specifications选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。

辽宁科技大学毕业设计（论文）第13页

图2.1滤波器幅频响应

2.8带通滤波器设计实例

本文将以一个FIR 滤波器的设计为例来说明如何使用MATLAB设计数字滤波器：在小电流接地系统中注入83.3Hz的正弦信号，对其进行跟踪分析，要求设计一带通数字滤波器，滤除工频及整次谐波，以便在非常复杂的信号中分离出该注入信号。参数要求：96阶FIR数字滤波器，采样频率1000Hz，采用Hamming

窗函数设计。本例中，首先在Filter Type中选择Bandpass（带通滤波器）；在Design Method选项中选择FIR Window（FIR滤波器窗函数法），接着在Window Specifications选项中选取Hamming；指定Filter Order项中的Specify Order＝95；由于采用窗函数法设计，只要给出通带下限截止频率Fc1和通带上限截止频率Fc2，选取Fc1＝70Hz，Fc2＝84Hz。设置完以后点击Design Filter即可得到

所设计的FIR滤波器。通过菜单选项Analysis可以在特性区看到所设计滤波器的幅频响应、相频响应、零极点配置和滤波器系数等各种特性。设计完成后将结果保存为1.fda文件［6］。

辽宁科技大学毕业设计（论文）第14页

在设计过程中，可以对比滤波器幅频相频特性和设计要求，随时调整参数和滤波器类型，以便得到最佳效果。其它类型的FIR滤波器和IIR滤波器也都可以使用FDATool来设计。程序设计法在MATLAB中，对各种滤波器的设计都有相应的计算振幅响应的函数，可以用来做滤波器的程序设计。

上例的带通滤波器可以用程序设计：

c=95; %定义滤波器阶数96阶

w1=2*pi*fc1/fs;

w2=2*pi*fc2/fs; %参数转换，将模拟滤波器

的技术指标转换为数字滤波器的技术指

标

window=hamming(c+1); %使用hamming窗函数

h=fir1(c,[w1/pi w2/pi],window); %使用标准响应的加窗设计函数fir1

freqz(h,1,512); %数字滤波器频率响应

在MATLAB环境下运行该程序即可得到滤波器幅频相频响应曲线和滤波器系数h。

辽宁科技大学毕业设计（论文）第15页

第3章GUI界面下噪声信号的频谱分析与实

现

介绍了通风噪声频谱分析的一种简单方法，此方法以普通声级计为一次仪表获得通风机空气动力噪声信号，用计算机进行信号采集，用数值分析软件MATLAB进行频谱分析。

3.1概述

通风机噪声是由于气体的冲击涡流及机械振动等原因使气体发生振动而引起的空气动力性噪声，它是由呈宽频带的涡流噪声和呈离散频带的旋转噪声叠加而成的。这两种噪声的产生原因不同,其控制方式也不相同。

为了有效地进行通风机噪声控制,通风机生产厂家或用户经常需要对具体的通风机进行空气动力噪声特性的测量和分析，特别是对风机噪声频谱的细致分析"但是,由于专业的频谱分析仪器价格昂贵，目前我国的一般中小型风机生产企业所具有的噪声测试仪器通常是一般的声级计，只能进行一些简单的声级测量和噪声频谱的粗略估计，还不具备对风机噪声频谱进行细致分析的能力。

幸运地是,随着计算机在风机行业深入广泛地应用和普及,用普通的个人计算机和相应的数值分析计算软件来进行风机噪声频谱的详细分析已经变得很容易了。本文介绍了一种运用数值分析软件MATLAB来进行通风机噪声频谱分析的方法。

3.2 MATLAB信号处理工具箱介绍及功能

3.2.1MATLAB工具箱介绍

MATLAB有三十多个工具箱大致可分为两类:功能型工具箱和领域型工具箱

［7］。

功能型工具箱：主要用来扩充MATLAB的符号计算功能、图形建模仿

真功能、文字处理功能以及与硬件实时交互功能，能用

于多种学科。

领域型工具箱：是专业性很强的。

如：控制工具箱(Control Toolbox)、信号处理工具箱(Signal Processing