FIR数字滤波器设计与使用
fir滤波器的主要设计方法 -回复

fir滤波器的主要设计方法-回复fir滤波器是一种基本的数字滤波器,主要用于数字信号处理中的滤波操作。
它的设计方法有很多种,包括频率采样法、窗函数法、最优权系数法等。
本文将一步一步回答"[fir滤波器的主要设计方法]",让我们一起来了解一下吧。
一、频率采样法频率采样法是fir滤波器设计的最基本方法之一。
它的主要思想是在频域中对滤波器的频响特性进行采样,然后通过反变换得到滤波器的冲激响应。
这种方法的优点是设计简单,适用于各种滤波器的设计。
1. 确定滤波器的截止频率和通带、阻带的要求。
根据应用的具体需求,确定滤波器的频率范围和滤波特性。
2. 设计理想的滤波器频率响应。
根据频率范围和滤波特性的要求,设计所需的滤波器频率响应。
常见的有低通、高通、带通、带阻等类型。
3. 进行频率采样。
根据滤波器频率响应的要求,在频域中进行一系列均匀或者非均匀的采样点。
4. 反变换得到滤波器的冲激响应。
对采样得到的频率响应进行反傅里叶变换,得到滤波器的冲激响应。
5. 标准化处理。
对得到的冲激响应进行标准化处理,使得滤波器的增益等于1。
6. 实现滤波器。
根据得到的冲激响应,使用差分方程或者卷积的方法实现fir滤波器。
二、窗函数法窗函数法是一种常用的fir滤波器设计方法,它主要是通过在频域中将理想的滤波器乘以一个窗函数来实现滤波器的设计。
1. 确定滤波器的截止频率和通带、阻带的要求,根据具体应用的需求确定滤波器的频率范围和滤波特性。
2. 设计理想的滤波器频率响应。
根据频率范围和滤波特性要求,设计所需的滤波器频率响应。
3. 选择窗函数。
根据滤波器的频率响应和窗函数的性质,选择合适的窗函数。
4. 计算窗函数的系数。
根据选择的窗函数,计算窗函数的系数。
5. 实现滤波器。
将理想滤波器的频率响应与窗函数相乘,得到实际的滤波器频率响应。
然后使用反变换将频率响应转换为滤波器的冲激响应。
6. 标准化处理。
对得到的冲激响应进行标准化处理,使得滤波器的增益等于1。
FIR滤波器设计要点

FIR滤波器设计要点FIR (Finite Impulse Response) 滤波器是一种数字滤波器,其设计要点包括滤波器类型选择、滤波器系数设计、频率响应规格、窗函数和滤波器长度选择等。
以下是对这些要点的详细介绍。
1.滤波器类型选择:在设计FIR滤波器之前,需要确定滤波器的类型。
常见的FIR滤波器类型包括低通滤波器、高通滤波器、带通滤波器和带阻滤波器。
不同类型的滤波器适用于不同的应用场景,因此在选择滤波器类型时需要考虑系统的需求。
2.滤波器系数设计:FIR滤波器的核心是滤波器系数。
滤波器系数决定了滤波器的频率响应和滤波特性。
常用的设计方法包括窗函数法、最小均方误差法和频率抽样法等。
窗函数法是最常用的设计方法,其基本思想是通过选择合适的窗函数来得到滤波器系数。
3.频率响应规格:在设计FIR滤波器时,需要明确所需的频率响应规格,包括通带增益、阻带衰减、过渡带宽等。
这些规格直接影响了滤波器的性能,因此需要根据具体应用场景来确定。
4.窗函数选择:窗函数在FIR滤波器设计中起到了重要的作用。
常用的窗函数包括矩形窗、汉宁窗、汉明窗、布莱克曼窗等。
不同的窗函数具有不同的特性,选择合适的窗函数可以得到优良的滤波器性能。
5.滤波器长度选择:滤波器长度决定了滤波器的频率分辨率和时间分辨率。
滤波器长度越长,频率响应越尖锐,但计算复杂度也越高。
因此,在设计FIR滤波器时需要权衡计算复杂度和性能要求,选择合适的滤波器长度。
6.优化设计:7.实现方式:总之,设计FIR滤波器要点包括滤波器类型选择、滤波器系数设计、频率响应规格、窗函数和滤波器长度选择等。
设计者需要根据具体的应用场景和性能要求来进行合理的设计和优化,以满足系统的需求。
实验五FIR数字滤波器的设计

实验五FIR数字滤波器的设计FIR数字滤波器(Finite Impulse Response)是一种数字滤波器,它的输出仅由有限数量的输入样本决定。
设计FIR数字滤波器的步骤如下:1.确定滤波器的要求:首先需要明确滤波器的频率响应、截止频率、通带和阻带的幅频响应等要求。
2.选择滤波器类型:根据实际需求选择合适的滤波器类型,如低通滤波器、高通滤波器、带通滤波器或带阻滤波器等。
3.确定滤波器的阶数:根据滤波器类型和要求,确定滤波器的阶数。
通常情况下,滤波器的阶数越高,能够实现更陡峭的频率响应,但会引入更多的计算复杂度。
4.设计滤波器的理想频率响应:根据滤波器的要求和类型,设计滤波器的理想频率响应。
可以使用常用的频率响应设计方法,如窗函数法、最小最大法或线性相位法等。
这些方法可以实现平滑的频率响应或者良好的阻带衰减。
5.确定滤波器的系数:根据设计的理想频率响应,通过反变换或优化算法确定滤波器的系数。
常用的优化算法包括频域方法、时域方法、最小二乘法或最小相位法等。
6.实现滤波器:将所得的滤波器系数转化为滤波器的差分方程形式或直接计算滤波器的频域响应。
7.评估滤波器性能:使用合适的测试信号输入滤波器,并对滤波器的输出进行评估。
可以使用指标,如频率响应曲线、幅度响应误差、相位响应误差或阻带衰减等指标来评估滤波器性能。
8.优化滤波器性能:根据评估结果,进行必要的修改和优化设计,以满足滤波器的要求。
通过以上步骤,可以设计出满足需求的FIR数字滤波器。
需要注意的是,FIR数字滤波器设计的复杂度和性能需要权衡与平衡,以满足实际应用的要求。
数字信号处理实验FIR数字滤波器的设计

数字信号处理实验:FIR数字滤波器的设计1. 引言数字滤波器是数字信号处理的关键技术之一,用于对数字信号进行滤波、降噪、调频等操作。
FIR (Finite Impulse Response) 数字滤波器是一种常见的数字滤波器,具有线性相应和有限的脉冲响应特性。
本实验旨在通过设计一个FIR数字滤波器来了解其基本原理和设计过程。
2. FIR数字滤波器的基本原理FIR数字滤波器通过对输入信号的每一个样本值与滤波器的冲激响应(滤波器的系数)进行线性加权累加,来实现对信号的滤波。
其数学表达式可以表示为:y(n) = b0 * x(n) + b1 * x(n-1) + b2 * x(n-2) + ... + bN * x(n-N)其中,y(n)表示滤波器的输出,x(n)表示滤波器的输入信号,b0~bN表示滤波器的系数。
FIR数字滤波器的脉冲响应为有限长度的序列,故称为有限冲激响应滤波器。
3. FIR数字滤波器的设计步骤FIR数字滤波器的设计主要包括以下几个步骤:步骤1: 确定滤波器的阶数和截止频率滤波器的阶数决定了滤波器的复杂度和性能,而截止频率决定了滤波器的通带和阻带特性。
根据实际需求,确定滤波器的阶数和截止频率。
步骤2: 选择滤波器的窗函数窗函数是FIR滤波器设计中常用的一种方法,可以通过选择不同的窗函数来实现不同的滤波器特性。
常用的窗函数有矩形窗、汉宁窗、汉明窗等。
根据实际需求,选择合适的窗函数。
步骤3: 计算滤波器的系数根据选择的窗函数和滤波器的阶数,使用相应的公式或算法计算滤波器的系数。
常见的计算方法有频率采样法、窗函数法、最小二乘法等。
步骤4: 实现滤波器根据计算得到的滤波器系数,可以使用编程语言或专用软件来实现滤波器。
步骤5: 评估滤波器性能通过输入测试信号,观察滤波器的输出结果,评估滤波器的性能和滤波效果。
常见评估指标有滤波器的幅频响应、相频响应、群延迟等。
4. 实验步骤本实验将以Matlab软件为例,演示FIR数字滤波器的设计步骤。
FIR滤波器设计分析

FIR滤波器设计分析FIR(Finite Impulse Response)滤波器是一类数字滤波器,其输出只取决于输入信号的有限数量的过去样本。
FIR滤波器的设计分析主要包括滤波器的设计目标、设计方法、设计参数选择、滤波器性能评估等方面。
首先,FIR滤波器的设计目标是根据特定的应用需求,设计一个能够满足给定要求的滤波器。
比如,在音频信号处理中,常见的设计目标包括降低噪声、增强语音清晰度等。
接下来,FIR滤波器的设计方法主要有窗函数法和频率采样法。
窗函数法是通过选择合适的窗函数来设计FIR滤波器,常见的窗函数有矩形窗、汉宁窗、汉明窗等。
频率采样法是通过在频域上选择一组等间隔的频率样点,然后通过频域设计方法将这些样点连接起来,得到FIR滤波器的频响。
设计参数选择是FIR滤波器设计的重要环节。
常见的设计参数包括滤波器阶数、截止频率、过渡带宽等。
滤波器阶数决定了滤波器的复杂度,一般情况下,滤波器阶数越高,滤波器的性能也会越好。
截止频率是指滤波器的频段边界,过渡带宽是指频域中通过频样点与阻带频样点之间的频带范围。
最后,FIR滤波器的性能评估主要包括幅频响应、相频响应、群延迟等指标。
幅频响应可以用来评估滤波器的频率特性,相频响应则描述了信号在滤波过程中的相对延迟。
群延迟是指信号通过滤波器时的延迟时间,对于实时信号处理应用非常重要。
总结起来,FIR滤波器设计分析主要涉及设计目标、设计方法、设计参数选择和滤波器性能评估四个方面。
通过合理选择设计方法和参数,并对滤波器的性能进行评估,可以设计出满足特定要求的FIR滤波器,从而实现信号处理、噪声降低等应用。
实验四FIR数字滤波器设计与软件实现

实验四FIR数字滤波器设计与软件实现
实验目的:
掌握FIR数字滤波器的设计与软件实现方法,了解滤波器的概念与基
本原理。
实验原理:
FIR数字滤波器全称为有限脉冲响应数字滤波器,其特点是具有有限
长度的脉冲响应。
滤波器通过一系列加权系数乘以输入信号的延迟值,并
将这些值相加得到输出信号。
FIR滤波器的频率响应由滤波器系数所决定。
实验步骤:
1.确定所需的滤波器的设计规格,包括截止频率、通带波纹、阻带衰
减等。
2.选择适当的滤波器设计方法,如窗函数、最佳近似法、最小二乘法等。
3.根据所选方法,计算滤波器的系数。
4.在MATLAB环境下,使用滤波器的系数实现滤波器。
5.输入所需滤波的信号,经过滤波器进行滤波处理。
6.分析输出的滤波信号,观察滤波效果是否符合设计要求。
实验要求:
1.完成FIR数字滤波器的设计和软件实现。
2.对比不同设计方法得到的滤波器性能差异。
3.分析滤波结果,判断滤波器是否满足设计要求。
实验器材与软件:
1.个人电脑;
2.MATLAB软件。
实验结果:
根据滤波器设计规格和所选的设计方法,得到一组滤波器系数。
通过
将滤波器系数应用于输入信号,得到输出滤波信号。
根据输出信号的频率
响应、通带波纹、阻带衰减等指标,评估滤波器的性能。
实验注意事项:
1.在选择设计方法时,需要根据滤波器要求和实际情况进行合理选择。
2.在滤波器实现过程中,需要注意滤波器系数的计算和应用。
3.在实验过程中,注意信号的选择和滤波结果的评估方法。
fir数字滤波器的设计与实现
FIR数字滤波器的设计与实现介绍在数字信号处理中,滤波器是一种常用的工具,用于改变信号的频率响应。
FIR (Finite Impulse Response)数字滤波器是一种非递归的滤波器,具有线性相位响应和有限脉冲响应。
本文将探讨FIR数字滤波器的设计与实现,包括滤波器的原理、设计方法和实际应用。
原理FIR数字滤波器通过对输入信号的加权平均来实现滤波效果。
其原理可以简单描述为以下步骤: 1. 输入信号经过一个延迟线组成的信号延迟器。
2. 延迟后的信号与一组权重系数进行相乘。
3. 将相乘的结果进行加和得到输出信号。
FIR滤波器的特点是通过改变权重系数来改变滤波器的频率响应。
不同的权重系数可以实现低通滤波、高通滤波、带通滤波等不同的滤波效果。
设计方法FIR滤波器的设计主要有以下几种方法:窗函数法窗函数法是一种常用简单而直观的设计方法。
该方法通过选择一个窗函数,并将其与理想滤波器的频率响应进行卷积,得到FIR滤波器的频率响应。
常用的窗函数包括矩形窗、汉宁窗、哈密顿窗等。
不同的窗函数具有不同的特性,在设计滤波器时需要根据要求来选择合适的窗函数。
频率抽样法频率抽样法是一种基于频率抽样定理的设计方法。
该方法首先将所需的频率响应通过插值得到一个连续的函数,然后对该函数进行逆傅里叶变换,得到离散的权重系数。
频率抽样法的优点是可以设计出具有较小幅频纹波的滤波器,但需要进行频率上和频率下的补偿处理。
最优化方法最优化方法是一种基于优化理论的设计方法。
该方法通过优化某个性能指标来得到最优的滤波器权重系数。
常用的最优化方法包括Least Mean Square(LMS)法、Least Square(LS)法、Parks-McClellan法等。
这些方法可以根据设计要求,如通带波纹、阻带衰减等来得到最优的滤波器设计。
实现与应用FIR数字滤波器的实现可以通过硬件和软件两种方式。
硬件实现在硬件实现中,可以利用专门的FPGA(Field-Programmable Gate Array)等数字集成电路来实现FIR滤波器。
FIR数字滤波器设计实验_完整版
FIR数字滤波器设计实验_完整版本实验旨在设计一种FIR数字滤波器,以滤除信号中的特定频率成分。
下面是完整的实验步骤:材料:-MATLAB或其他支持数字信号处理的软件-计算机-采集到的信号数据实验步骤:1.收集或生成需要滤波的信号数据。
可以使用外部传感器采集数据,或者在MATLAB中生成一个示波器信号。
2. 在MATLAB中打开一个新的脚本文件,并导入信号数据。
如果你是使用外部传感器采集数据,请将数据以.mat文件的形式保存,并将其导入到MATLAB中。
3.对信号进行预处理。
根据需要,你可以对信号进行滤波、降噪或其他预处理操作。
这可以确保信号数据在输入FIR滤波器之前处于最佳状态。
4.确定滤波器的设计规范。
根据信号的特性和要滤除的频率成分,确定FIR滤波器的设计规范,包括滤波器的阶数、截止频率等。
你可以使用MATLAB中的函数来帮助你计算滤波器参数。
5. 设计FIR滤波器。
使用MATLAB中的fir1函数或其他与你所使用的软件相对应的函数来设计满足你的规范条件的FIR滤波器。
你可以选择不同的窗函数(如矩形窗、汉宁窗等)来平衡滤波器的频域和时域性能。
6. 对信号进行滤波。
将设计好的FIR滤波器应用到信号上,以滤除特定的频率成分。
你可以使用MATLAB中的conv函数或其他相应函数来实现滤波操作。
7.分析滤波效果。
将滤波后的信号与原始信号进行比较,评估滤波效果。
你可以绘制时域图、频域图或其他特征图来分析滤波效果。
8.优化滤波器设计。
如果滤波效果不理想,你可以调整滤波器设计参数,重新设计滤波器,并重新对信号进行滤波。
这个过程可能需要多次迭代,直到达到最佳的滤波效果。
9.总结实验结果。
根据实验数据和分析结果,总结FIR滤波器设计的优点和缺点,以及可能的改进方向。
通过完成以上实验步骤,你将能够设计并应用FIR数字滤波器来滤除信号中的特定频率成分。
这对于许多信号处理应用都是非常重要的,如音频处理、图像处理和通信系统等。
实验四FIR数字滤波器的设计
实验四FIR数字滤波器的设计
FIR数字滤波器也称作有限脉冲响应数字滤波器,是一种常见的数字滤波器设计方法。
在设计FIR数字滤波器时,需要确定滤波器的阶数、滤波器的类型(低通、高通、带通、带阻)以及滤波器的参数(截止频率、通带波纹、阻带衰减、过渡带宽等)。
下面是FIR数字滤波器的设计步骤:
1.确定滤波器的阶数。
阶数决定了滤波器的复杂度,一般情况下,阶数越高,滤波器的性能越好,但计算量也越大。
阶数的选择需要根据实际应用来进行权衡。
2.确定滤波器的类型。
根据实际需求,选择低通、高通、带通或带阻滤波器。
低通滤波器用于去除高频噪声,高通滤波器用于去除低频噪声,带通滤波器用于保留一定范围内的频率信号,带阻滤波器用于去除一定范围内的频率信号。
3.确定滤波器的参数。
根据实际需求,确定滤波器的截止频率、通带波纹、阻带衰减和过渡带宽等参数。
这些参数决定了滤波器的性能。
4.设计滤波器的频率响应。
使用窗函数、最小二乘法等方法,根据滤波器的参数来设计滤波器的频率响应。
5.将频率响应转换为滤波器的系数。
根据设计的频率响应,使用逆快速傅里叶变换(IFFT)等方法将频率响应转换为滤波器的系数。
6.实现滤波器。
将滤波器的系数应用到数字信号中,实现滤波操作。
7.优化滤波器性能。
根据需要,可以对滤波器进行进一步优化,如调整滤波器的阶数、参数等,以达到较好的滤波效果。
以上是FIR数字滤波器的设计步骤,根据实际需求进行相应的调整,可以得到理想的滤波器。
fir数字滤波器设计流程
fir数字滤波器设计流程英文回答:Designing a FIR (Finite Impulse Response) digitalfilter involves several steps. I will explain the processin detail below.1. Specify the filter requirements: The first step isto clearly define the desired characteristics of the filter. This includes the filter type (low-pass, high-pass, band-pass, or band-stop), cutoff frequencies, passband ripple, stopband attenuation, and any other relevant specifications.For example, let's say I want to design a low-pass FIR filter with a cutoff frequency of 1 kHz, a passband rippleof 0.1 dB, and a stopband attenuation of 60 dB.2. Choose a filter design method: There are various methods available for FIR filter design, such as windowing, frequency sampling, and least squares. The choice of methoddepends on the desired filter characteristics and design constraints.Continuing with our example, I decide to use the windowing method for simplicity.3. Select a window function: In windowing, a window function is applied to the ideal impulse response of the filter to obtain a finite-length impulse response. Commonly used window functions include Hamming, Hanning, and Blackman.In our case, I choose the Hamming window function.4. Determine the filter length: The length of thefilter determines the trade-off between frequencyresolution and time-domain performance. Longer filters provide better frequency resolution but require more computational resources.To determine the filter length, I use a formula that takes into account the desired cutoff frequency and thewindow function.5. Generate the ideal impulse response: Using the desired filter characteristics and the determined filter length, I generate the ideal impulse response of the filter. This is done by applying the appropriate mathematical equations or algorithms.In our example, I generate the ideal impulse responseof the low-pass filter.6. Apply the window function: The next step is to apply the selected window function to the ideal impulse response. This is done by multiplying the window function with the ideal impulse response.For our low-pass filter, I multiply the Hamming window function with the ideal impulse response.7. Normalize the filter coefficients: The filter coefficients are normalized to ensure that the filter response meets the desired specifications. This istypically done by dividing the coefficients by the sum of their absolute values.In our case, I normalize the filter coefficients to ensure the passband ripple and stopband attenuation are within the specified limits.8. Implement the filter: Finally, the designed filter can be implemented in hardware or software for signal processing applications. This involves programming thefilter coefficients into a digital signal processor (DSP) or using a software library for FIR filtering.In conclusion, the process of designing a FIR digital filter involves specifying the filter requirements, choosing a design method, selecting a window function, determining the filter length, generating the ideal impulse response, applying the window function, normalizing the filter coefficients, and implementing the filter.中文回答:设计一个有限脉冲响应(FIR)数字滤波器涉及多个步骤。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验报告课程名称:数字信号处理指导老师:刘英成绩:_________________实验名称: FIR数字滤波器设计与使用同组学生姓名:__________一、实验目的和要求设计和应用FIR低通滤波器。
掌握FIR数字滤波器的窗函数设计法,了解设计参数(窗型、窗长)的影响。
二、实验内容和步骤编写MATLAB程序,完成以下工作。
2-1 设计两个FIR低通滤波器,截止频率C=0.5。
(1)用矩形窗,窗长N=41。
得出第一个滤波器的单位抽样响应序列h1(n)。
记下h1(n)的各个抽样值,显示h1(n)的图形(用stem(.))。
求出该滤波器的频率响应(的N个抽样)H1(k),显示|H1(k)|的图形(用plot(.))。
(2)用汉明窗,窗长N=41。
得出第二个滤波器的单位抽样响应序列h2(n)。
记下h2(n)的各个抽样值,显示h2(n)的图形。
求出滤波器的频率响应H2(k),显示|H2(k)|的图形。
(3)由图形,比较h1(n)与h2(n)的差异,|H1(k)|与|H2(k)|的差异。
2-2 产生长度为200点、均值为零的随机信号序列x(n)(用rand(1,200)0.5)。
显示x(n)。
求出并显示其幅度谱|X(k)|,观察特征。
2-3 滤波(1)将x(n)作为输入,经过第一个滤波器后的输出序列记为y1(n),其幅度谱记为|Y1(k)|。
显示|X(k)|与|Y1(k)|,讨论滤波前后信号的频谱特征。
(2)将x(n)作为输入,经过第二个滤波器后的输出序列记为y2(n),其幅度谱记为|Y2(k)|。
比较|Y1(k)|与|Y2(k)|的图形,讨论不同的窗函数设计出的滤波器的滤波效果。
2-4 设计第三个FIR低通滤波器,截止频率C=0.5。
用矩形窗,窗长N=127。
用它对x(n)进行滤波。
显示输出信号y3(n)的幅度谱|Y3(k)|,并与|Y1(k)|比较,讨论不同的窗长设计出的滤波器的滤波效果。
三、主要仪器设备自行编程。
四、操作方法和实验步骤(参见“二、实验内容和步骤”)五、实验数据记录和处理5.1 列出MATLAB程序清单,加注释。
%% 2-1clear; close all; clc% the length of window = 41wc = 0.5*pi;N = 41;n = 0:N-1;h1 = fir1(N-1,wc/pi,boxcar(N));[H1,w1] = freqz(h1);h2 = fir1(N-1,wc/pi);[H2,w2] = freqz(h2);figure;subplot(1,2,1);stem(n,h1,'filled');grid on;xlabel('n','FontSize',14);ylabel('h1(n)','FontSize',14);title('矩形窗得到的低通滤波器h1(n)的时域','FontSize',14);subplot(1,2,2);stem(n,h2,'filled');grid on;axis([0 30 -0.2 0.6]);xlabel('n','FontSize',14);ylabel('h2(n)','FontSize',14);title('矩形窗得到的低通滤波器h2(n)的时域','FontSize',14);figure;subplot(1,2,1);plot(w1/pi,abs(H1));grid on;xlabel('\omega /\pi','FontSize',14);ylabel('|H1(k)|','FontSize',14);title('矩形窗得到的低通滤波器h1(n)的幅度谱','FontSize',14);subplot(1,2,2);plot(w2/pi,abs(H2));grid on;xlabel('\omega /\pi','FontSize',14);ylabel('|H2(k)|','FontSize',14);title('矩形窗得到的低通滤波器h2(n)的幅度谱','FontSize',14);figure;subplot(1,2,1);plot(w1/pi,20*log10(abs(H1)));grid on;axis([0 1 -100 20]);xlabel('\omega /\pi','FontSize',14);ylabel('20lg|H1(k)|/dB','FontSize',14); title(' h1(n)的幅度谱','FontSize',14);subplot(1,2,2);plot(w2/pi,20*log10(abs(H2)));grid on;xlabel('\omega /\pi','FontSize',14);ylabel('20lg|H2(k)|/dB','FontSize',14); title(' h2(n)的幅度谱','FontSize',14);pause;%% 2-2N = 200;n = 0:N-1;x = rand(1,N)-0.5;[X,w] = freqz(x);figure;subplot(2,1,1);stem(n,x,'filled');grid on;xlabel('n','FontSize',14);ylabel('x(n)','FontSize',14);title('随机序列x(n)的时域','FontSize',14);subplot(2,1,2);plot(w/pi,abs(X));grid on;xlabel('\omega /\pi','FontSize',14);ylabel('|X(k)|','FontSize',14);title('随机序列x(n)的幅度谱','FontSize',14);pause;%% 2-3y1 = filter(h1,1,x);[Y1,w1] = freqz(y1);y2 = filter(h2,1,x);[Y2,w2] = freqz(y2);figure;subplot(2,1,1)plot(w1/pi,abs(Y1));grid on;axis([0 1 0 12]);set(gca,'ytick',0:4:12);xlabel('\omega /\pi','FontSize',14);ylabel('|Y1(k)|','FontSize',14);title('输出信号y1(n)的幅度谱','FontSize',14);subplot(2,1,2);plot(w2/pi,abs(Y2));grid on;axis([0 1 0 12]);set(gca,'ytick',0:4:12);xlabel('\omega /\pi','FontSize',14);ylabel('|Y2(k)|','FontSize',14);title('输出信号y2(n)的幅度谱','FontSize',14);pause;%% 2-4N = 127;n = 0:N-1;h3 = fir1(N-1,wc/pi,boxcar(N));[H3,w3] = freqz(h3);y3 = filter(h3,1,x);[Y3,wy3] = freqz(y3);figure;plot(wy3/pi,abs(Y3));grid on;axis([0 1 0 12]);set(gca,'ytick',0:4:12); xlabel('\omega /\pi','FontSize',14);ylabel('|Y3(k)|','FontSize',14);title('输出信号y3(n)的幅度谱','FontSize',14);figure;subplot(1,2,1);plot(w3/pi,abs(H3));grid on;xlabel('\omega /\pi','FontSize',14);ylabel('|H3(k)|','FontSize',14);title('矩形窗得到的低通滤波器h3(n)的幅度谱','FontSize',14);subplot(1,2,2);plot(w3/pi,20*log10(abs(H3)));grid on;axis([0 1 -100 20]);xlabel('\omega /\pi','FontSize',14);ylabel('20lg|H3(k)|/dB','FontSize',14); title(' h3(n)的幅度谱','FontSize',14);5.2 列出计算结果,包括h1(n)和h2(n)的各个抽样值,|H1(k)|、|H2(k)|、|X(k)|、|Y1(k)|、|Y2(k)|和|Y3(k)|的图形。
(1)h1(n)和h2(n)(2)|H1(k)|、|H2(k)|(3)|X(k)|、|Y1(k)|、|Y2(k)|和|Y3(k)|六、实验结果与分析观察结果,进行讨论,最后总结:滤波器的频率响应中的过渡带宽度取决于哪些设计参数?有什么规律?阻带最小衰减取决于哪些设计参数?有什么规律?(1)观察实验结果,逐一进行讨论结论:从上面的序列图和记录的 h1(n)与 h2(n)的序列值中,我们可以看出 h1(n)与h2(n)的最大值都相同,在主瓣上取样点数相同且取样值几乎相等,但旁瓣上h1(n)的取样值在0的周围波动较大,h2(n)的取样值在0的周围波动较小。