lp数字信号处理
- 格式:doc
- 大小:295.50 KB
- 文档页数:30
西南大学培训与继续教育学院课程代码:1077学年学季:20202窗体顶端判断题1、应用DFT分析无限长信号的频谱时,必然会产生误差。
A.√B.×2、离散周期信号的DFS中,频域的周期N对应数字频率为2π。
A.√B.×3、实数序列的DFT为共轭对称的序列。
A.√B.×4、一个域的周期性,对应另一域的离散性。
A.√B.×5、信号的最高频率为3π/5,则最大程度减小数据量的I/D值为3/5。
A.√B.×6、单位圆上的零点,对应幅频特性的零值。
A.√B.×7、LP表示的滤波器类型是低通滤波器。
A.√B.×8、通带最平坦的滤波器是巴特沃思滤波器。
A.√B.×9、陷波器必然有零点位于单位圆上。
A.√B.×10、圆周卷积和线卷积相等的条件是圆周卷积的点数不小于线性卷积的长度。
A.√B.×11、按照最大误差最小准则设计的滤波器,具有等波纹的特点。
A.√B.×12、单位脉冲序列的DTFT结果为1。
A.√B.×13、x(n)与h(n)的卷积的Z变换为X(Z)H(Z)。
A.√B.×14、所谓全通系统,就是其频率响应的幅度在任意需要考虑的频率点处均为常数。
A.√B.×15、FIR滤波器由于无原点外的极点,故相比IIR阶次更高。
A.√B.×16、对连续信号作频谱分析,设信号的采样频率为10KHz,频域的分辨能力为不大于10Hz,则对应DFS点数为1000点。
A.√B.×17、靠近单位圆上的极点,对应幅频特性的极大值。
A.√B.×18、线性相位可分为第一类与第二类线性相位两种情况。
A.√B.×19、为满足线性相位要求,窗函数本身也应满足相应的对称性。
A.√B.×20、冲激响应不变法由于存在混叠,不能设计高通、带通滤波器。
A.√B.×21、FIR滤波器的结构往往是非递归型的。
数字信号处理知识点汇总数字信号处理是一门涉及多个领域的重要学科,在通信、音频处理、图像处理、控制系统等众多领域都有着广泛的应用。
接下来,让我们一同深入了解数字信号处理的主要知识点。
一、数字信号的基本概念数字信号是在时间和幅度上都离散的信号。
与模拟信号相比,数字信号具有更强的抗干扰能力和便于处理、存储等优点。
在数字信号中,我们需要了解采样定理。
采样定理指出,为了能够从采样后的信号中完全恢复原始的连续信号,采样频率必须至少是原始信号最高频率的两倍。
这是保证数字信号处理准确性的关键原则。
二、离散时间信号与系统离散时间信号可以通过序列来表示,常见的有单位脉冲序列、单位阶跃序列等。
离散时间系统则是对输入的离散时间信号进行运算和处理,产生输出信号。
系统的特性可以通过线性、时不变性、因果性和稳定性等方面来描述。
线性系统满足叠加原理,即多个输入的线性组合产生的输出等于各个输入单独作用产生的输出的线性组合。
时不变系统的特性不随时间变化,输入的时移会导致输出的相同时移。
因果系统的输出只取决于当前和过去的输入,而稳定系统对于有界的输入会产生有界的输出。
三、Z 变换Z 变换是分析离散时间系统的重要工具。
它将离散时间信号从时域转换到复频域。
通过 Z 变换,可以方便地求解系统的差分方程,分析系统的频率特性和稳定性。
Z 变换的收敛域决定了其特性和应用范围。
逆 Z 变换则可以将复频域的函数转换回时域信号。
四、离散傅里叶变换(DFT)DFT 是数字信号处理中的核心算法之一。
它将有限长的离散时间信号转换到频域。
DFT 的快速算法——快速傅里叶变换(FFT)大大提高了计算效率,使得在实际应用中能够快速处理大量的数据。
通过 DFT,可以对信号进行频谱分析,了解信号的频率成分和能量分布。
五、数字滤波器数字滤波器用于对数字信号进行滤波处理,分为有限冲激响应(FIR)滤波器和无限冲激响应(IIR)滤波器。
FIR 滤波器具有线性相位特性,稳定性好,但设计相对复杂。
1077 20192判断题1、应用DFT分析无限长信号的频谱时,必然会产生误差。
. A.√. B.×2、离散周期信号的DFS中,频域的周期N对应数字频率为2π 。
. A.√. B.×3、实数序列的DFT为共轭对称的序列。
. A.√. B.×4、一个域的周期性,对应另一域的离散性。
. A.√. B.×5、单位圆上的零点,对应幅频特性的零值。
. A.√. B.×6、LP表示的滤波器类型是低通滤波器。
. A.√. B.×7、圆周卷积和线卷积相等的条件是圆周卷积的点数不小于线性卷积的长度。
. A.√. B.×8、单位脉冲序列的DTFT结果为1。
. A.√. B.×9、x(n)与h(n)的卷积的Z变换为X(Z)H(Z)。
. A.√. B.×10、所谓全通系统,就是其频率响应的幅度在任意需要考虑的频率点处均为常数。
. A.√. B.×11、FIR滤波器由于无原点外的极点,故相比IIR阶次更高。
. A.√. B.×12、对连续信号作频谱分析,设信号的采样频率为10KHz,频域的分辨能力为不大于10Hz,则对应DFS点数为1000 点。
. A.√. B.×13、靠近单位圆上的极点,对应幅频特性的极大值。
. A.√. B.×14、线性相位可分为第一类与第二类线性相位两种情况。
. A.√. B.×15、为满足线性相位要求,窗函数本身也应满足相应的对称性。
. A.√. B.×16、冲激响应不变法由于存在混叠,不能设计高通、带通滤波器。
. A.√. B.×17、FIR滤波器的结构往往是非递归型的。
. A.√. B.×18、单位延迟单元对应的硬件是存储器,其数目影响系统的复杂度。
. A.√. B.×19、时域加窗,频域会产生频谱泄漏。
. A.√. B.×20、从s域到z域映射,虚轴和单位圆、左半平面与单位圆内部,都必须对应。
数字信号处理综述数字信号处理(Digital Signal Processing,DSP)是指对数字信号进行采样、量化和运算等处理的技术领域。
它在现代通信、图像、音频、视频等领域中起着重要的作用。
本文将对数字信号处理的基本原理、应用领域和未来发展进行综述。
一、数字信号处理的基本原理数字信号处理基于离散时间信号,通过数学运算对信号进行处理。
其基本原理包括采样、量化和离散化等步骤。
1. 采样:将连续时间信号转换为离散时间信号,通过对连续时间信号进行等间隔采样,得到一系列的采样值。
2. 量化:将连续幅度信号转换为离散幅度信号。
量化是对连续幅度信号进行近似处理,将其离散化为一系列的离散值。
3. 离散化:将连续时间信号的采样值和离散幅度信号的量化值进行结合,形成离散时间、离散幅度的数字信号。
通过采样、量化和离散化等步骤,数字信号处理能够对原始信号进行数字化表示和处理。
二、数字信号处理的应用领域数字信号处理广泛应用于各个领域,其中包括但不限于以下几个方面。
1. 通信领域:数字信号处理在通信中起着重要作用。
它能够提高信号的抗干扰性能、降低信号传输误码率,并且能够实现信号压缩和编解码等功能。
2. 音频与视频处理:数字信号处理在音频与视频处理中具有重要应用。
它可以实现音频的降噪、音频编码和解码、语音识别等功能。
在视频处理中,数字信号处理可以实现视频压缩、图像增强和视频流分析等功能。
3. 生物医学工程:数字信号处理在生物医学工程中的应用越来越广泛。
它可以实现医学图像的增强和分析、生物信号的滤波和特征提取等功能,为医学诊断和治疗提供支持。
4. 雷达与成像技术:数字信号处理在雷达与成像技术中有重要的应用。
通过数字信号处理,可以实现雷达信号的滤波和目标检测、图像的恢复和重建等功能。
5. 控制系统:数字信号处理在控制系统中起着重要作用。
它可以实现控制信号的滤波、系统的辨识和控制算法的优化等功能。
三、数字信号处理的未来发展随着科技的进步和应用需求的不断增加,数字信号处理在未来有着广阔的发展空间。
6.2 教材第六章习题解答1. 设计一个巴特沃斯低通滤波器,要求通带截止频率6p f kHz =,通带最大衰减3p a dB =,阻带截止频率12s f kHz =,阻带最小衰减3s a dB =。
求出滤波器归一化传输函数()a H p 以及实际的()a H s 。
解:(1)求阶数N 。
lg lg sp spk N λ=-0.10.30.1 2.51011010.0562101101p s asp a k --==≈--332121022610s sp p πλπΩ⨯⨯===Ω⨯⨯将sp k 和sp λ值代入N 的计算公式得lg 0.05624.15lg 2N =-=所以取N=5(实际应用中,根据具体要求,也可能取N=4,指标稍微差一点,但阶数低一阶,使系统实现电路得到简化。
) (2)求归一化系统函数()a H p ,由阶数N=5直接查表得到5阶巴特沃斯归一化低通滤波器系统函数()a H p 为54321() 3.2361 5.2361 5.2361 3.23611a H p p p p p p =+++++或 221()(0.6181)( 1.6181)(1)a H p p p p p p =+++++ 当然,也可以按(6.12)式计算出极点:121()22,0,1,2,3,4k j Nk p ek π++==按(6.11)式写出()a H p 表达式41()()a k k H p p p ==-代入k p 值并进行分母展开得到与查表相同的结果。
(3)去归一化(即LP-LP 频率变换),由归一化系统函数()a H p 得到实际滤波器系统函数()a H s 。
由于本题中3p a dB =,即32610/c p rad s πΩ=Ω=⨯⨯,因此()()a a cH s H p s p ==Ω5542332453.2361 5.2361 5.2361 3.2361c c c cc cs s ss s Ω=+Ω+Ω+Ω+Ω+Ω对分母因式形式,则有()()a a cH s H p s p ==Ω52222(0.6180)( 1.6180)()c c c c cc s s s s s Ω=+Ω-Ω+Ω-Ω+Ω如上结果中,c Ω的值未代入相乘,这样使读者能清楚地看到去归一化后,3dB 截止频率对归一化系统函数的改变作用。
数字信号处理数字信号处理(Digital Signal Processing)数字信号处理是指将连续时间的信号转换为离散时间信号,并对这些离散时间信号进行处理和分析的过程。
随着计算机技术的飞速发展,数字信号处理在各个领域得到了广泛应用,如通信、医学影像、声音处理等。
本文将介绍数字信号处理的基本概念和原理,以及其在不同领域的应用。
一、数字信号处理的基本概念数字信号处理是建立在模拟信号处理基础之上的一种新型信号处理技术。
在数字信号处理中,信号是用数字形式来表示和处理的,因此需要进行模数转换和数模转换。
数字信号处理的基本原理包括采样、量化和编码这三个步骤。
1. 采样:采样是将连续时间信号在时间上进行离散化的过程,通过一定的时间间隔对信号进行取样。
采样的频率称为采样频率,一般以赫兹(Hz)为单位表示。
采样频率越高,采样率越高,可以更准确地表示原始信号。
2. 量化:量化是指将连续的幅度值转换为离散的数字值的过程。
在量化过程中,需要确定一个量化间隔,将信号分成若干个离散的级别。
量化的级别越多,表示信号的精度越高。
3. 编码:编码是将量化后的数字信号转换为二进制形式的过程。
在数字信号处理中,常用的编码方式有PCM(脉冲编码调制)和DPCM (差分脉冲编码调制)等。
二、数字信号处理的应用1. 通信领域:数字信号处理在通信领域中具有重要的应用价值。
在数字通信系统中,信号需要经过调制、解调、滤波等处理,数字信号处理技术可以提高信号传输的质量和稳定性。
2. 医学影像:医学影像是数字信号处理的典型应用之一。
医学影像技术如CT、MRI等需要对采集到的信号进行处理和重建,以获取患者的影像信息,帮助医生进行诊断和治疗。
3. 声音处理:数字信号处理在音频处理和语音识别领域也有广泛的应用。
通过数字滤波、噪声消除、语音识别等技术,可以对声音信号进行有效处理和分析。
总结:数字信号处理作为一种新兴的信号处理技术,已经深入到各个领域中,并取得了显著的进展。
数字信号处理技术简介引言:- 数字信号处理技术是以数字计算机为基础的一种信号处理方法,用于对连续时间的模拟信号进行数字化处理。
- 数字信号处理在音频、视频、图像、通信等领域有广泛的应用,提高了信号处理的精度和效率。
一、什么是数字信号处理技术- 数字信号处理技术通过对模拟信号进行采样、量化和编码,将其转化为数字信号。
- 数字信号可以存储、传输和处理,具有较好的稳定性和灵活性。
二、数字信号处理的基本步骤1. 信号采样:- 采样是指以一定的时间间隔对模拟信号进行取样。
- 采样率决定了采样频率,一般要满足奈奎斯特采样定理。
2. 信号量化:- 量化是指将连续的模拟信号变为离散的数字信号。
- 通过将信号的幅度分成若干个离散的级别,将每个采样点映射到最近的一个量化级别上。
3. 信号编码:- 编码是指将量化后的信号转化为二进制,以便数字系统进行处理。
- 常用的编码方式有脉冲编码调制(PCM)、ΔΣ调制等。
4. 数字信号处理算法:- 数字信号处理算法是对数字信号进行处理和分析的数学方法和步骤。
- 常用的算法包括傅里叶变换、滤波、时域分析、频域分析等。
5. 数字信号重构:- 数字信号重构是将处理后的数字信号转化为模拟信号,以供输出和显示。
- 重构过程中需要进行数模转换和滤波处理。
三、数字信号处理技术的应用领域1. 通信领域:- 数字信号处理技术在调制解调、信道编码、信号恢复、自适应滤波等方面有广泛应用。
- 提高了通信系统的抗干扰能力和通信质量。
2. 音频与视频处理:- 数字信号处理技术在音频压缩、回声消除、音频增强、视频编解码等方面发挥重要作用。
- 提高了音频视频设备的音质和图像质量。
3. 图像处理与识别:- 数字信号处理技术在图像压缩、图像特征提取、目标检测与识别中有广泛应用。
- 提高了图像处理的速度和准确度。
4. 生物医学信号处理:- 数字信号处理技术在心电信号分析、脑电信号处理、医学影像处理等方面具有重要意义。
数字信号处理什么是数字信号处理?数字信号处理(Digital Signal Processing,简称DSP)是一种利用数字计算机进行信号处理的技术。
它将输入信号采样并转换成数字形式,在数字域上进行各种运算和处理,最后将处理后的数字信号转换回模拟信号输出。
数字信号处理在通信、音频、视频等领域都有广泛的应用。
数字信号处理的基本原理数字信号处理涉及许多基本原理和算法,其中包括信号采样、量化、离散化、频谱分析、滤波等。
信号采样信号采样是指将连续的模拟信号转换为离散的数字信号。
采样定理指出,为了能够准确地还原原始信号,采样频率必须大于信号中最高频率的两倍。
常用的采样方法有均匀采样和非均匀采样。
量化量化是将连续的模拟信号离散化为一组有限的量化值。
量化过程中,需要将连续信号的振幅映射为离散级别。
常见的量化方法有均匀量化和非均匀量化,其中均匀量化是最为常用的一种方法。
离散化在数字信号处理中,信号通常被表示为离散序列。
离散化是将连续的模拟信号转换为离散的数字信号的过程。
频谱分析频谱分析是一种用于研究信号频域特性的方法。
通过对信号的频谱进行分析,可以提取出其中的频率成分,了解信号的频率分布情况。
滤波滤波是数字信号处理中常用的一种方法,用于去除信号中的噪声或不需要的频率成分。
常见的滤波器有低通滤波器、高通滤波器、带通滤波器和带阻滤波器等。
数字信号处理的应用数字信号处理在许多领域都有广泛应用,下面列举了其中几个重要的应用领域:通信在通信领域,数字信号处理主要用于调制解调、信道编码、信号分析和滤波等方面。
数字信号处理的应用使得通信系统更加稳定和可靠,提高了通信质量和传输效率。
音频处理在音频处理领域,数字信号处理广泛应用于音频信号的录制、编码、解码、增强以及音频效果的处理等方面。
数字音乐、语音识别和语音合成等技术的发展离不开数字信号处理的支持。
视频处理数字信号处理在视频处理领域也发挥着重要作用。
视频压缩、图像增强、视频编码和解码等技术都离不开数字信号处理的支持。
一、原始语音信号采样后的时域、频域分析1、源程序如下:[x,fs,bits]=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); x1 = x(:,1); %读取语音信号的数据,赋给变量x1 figure(1)plot(x1) %做原始语音信号以44.1k 采样后的时域图形 sound(x1,fs);%播放语音信号title('原始语音采样后时域信号'); xlabel('时间轴 n'); ylabel('幅值 A');figure(2)freqz(x1) %绘制原始语音信号采样后的频率响应图 title('原始语音信号采样后频率响应图');[x1,fs,bits]=wavread('Windows XP 关机.wav'); %读取语音信号的数据,赋给变量x1 sound(x1,fs); %播放语音信号012345678x 104-0.8-0.6-0.4-0.20.20.40.6原始语音采样后时域信号时间轴 n幅值 A00.10.20.30.40.50.60.70.80.91-10-54Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )00.10.20.30.40.50.60.70.80.91-5050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )原始语音信号采样后频率响应图FFT 变换分析[x,fs,bits]=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); %读取语音信号的数据,赋给变量x1y1=fft(x1,1024); %对采样后信号做1024点FFT 变换 f=fs*(0:511)/1024;%生成频率序列figure(1)subplot(2,1,1);plot(f,abs(y1(1:512))); title('原始语音信号频谱') xlabel('频率/Hz'); ylabel('幅值')subplot(2,1,2);plot(abs(y1(1:1024))) %采样后信号的FFT 频谱图title('原始语音信号FFT 频谱') xlabel('点数N'); ylabel('幅值');02000400060008000100001200020406080原始语音信号频谱频率/Hz幅值20040060080010001200020406080原始语音信号FFT 频谱点数N幅值频谱分析:从图3-3可以看出,采样点数为1024,离散点数越多,越接近原始信号频谱双线性变换法设计IIR 数字滤波器(1) 低通滤波器性能指标 fb =1000Hz ,fc =1200Hz ,As =100dB ,Ap =1dB 。
fs=44100;[x2]=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); Ts=1/fs;R1=10; %Ap 通带最大衰减 wp=2*pi*1000/fs; %通带边界频率 ws=2*pi*1200/fs; %阻带边界频率 Rp=1; %As 阻带最小衰减Rl=100;wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标 ws1=2/Ts*tan(ws/2);[N,Wc]=buttord(wp1,ws1,Rp,R1,'s'); %计算N 和wc 值[Z,P ,K]=buttap(N); %创建butterworth 模拟滤波器[Bap,Aap]=zp2tf(Z,P ,K); %计算零点、极点得到系统函数的分子和分母多项式 [b,a]=lp2lp(Bap,Aap,Wc); %计算系统分子和分母多项式[bz,az]=bilinear(b,a,fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,W]=freqz(bz,az); %绘制频率响应曲线 figure(1)plot(W*fs/(2*pi),abs(H)) gridxlabel('频率/Hz')ylabel('频率响应幅度')title('IIR低通滤波器')f1=filter(bz,az,x2); %滤波figure(2)subplot(2,1,1)plot(x2) %画出滤波前的时域图title('IIR低通滤波器滤波前的时域波形');subplot(2,1,2)plot(f1); %画出滤波后的时域图title('IIR低通滤波器滤波后的时域波形');sound(f1,44100); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024; %511相当于去fs/2 figure(3)y2=fft(x2,1024);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('幅值');运行结果00.511.522.5x 1040.20.40.60.811.21.4频率/Hz频率响应幅度IIR 低通滤波器12345678x 104-1-0.500.51IIR 低通滤波器滤波前的时域波形012345678x 104-0.50.5IIR 低通滤波器滤波后的时域波形0.511.522.5x 104020406080IIR 低通滤波器滤波前的频谱频率/Hz幅值00.511.522.5x 104246IIR 低通滤波器滤波后的频谱频率/Hz幅值(2) 高通滤波器性能指标 fc =4800Hz ,fb =5000Hz ,As =100dB ,Ap =1dB 。
fs=44100;[x2]=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); Fs=22050;Ts=1/Fs;R1=50; Wp=2*pi*5000/fs; Ws=2*pi*4800/fs; Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2); %将模拟指标转换成数字指标 Ws1=2/Ts*tan(Ws/2);[N,Wn]=cheb2ord(Wp1,Ws1,Rp,Rl,'s'); %选择滤波器的最小阶数 [Z,P ,K]=cheb2ap(N,Rl); %创建切比雪夫模拟滤波器 [Bap,Aap]=zp2tf(Z,P ,K);[b,a]=lp2hp(Bap,Aap,Wn);[bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,W]=freqz(bz,az); %绘制频率响应曲线 figure(1)plot(W*fs/(2*pi),abs(H)) gridxlabel('频率/Hz') ylabel('频率响应幅度')title('IIR高通滤波器')f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(x1) %画出滤波前的时域图title('IIR高通滤波器滤波前的时域波形');subplot(2,1,2)plot(f1); %画出滤波后的时域图title('IIR高通滤波器滤波后的时域波形');sound(f1,44100); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512))); %画出滤波前的频谱图title('IIR高通滤波器滤波前的频谱')xlabel('频率/Hz');ylabel('幅值');subplot(2,1,2)plot(f,abs(F0(1:512))); %画出滤波后的频谱图title('IIR高通滤波器滤波后的频谱')xlabel('频率/Hz');ylabel('幅值');结果如下:00.511.522.5x 104123456频率/Hz频率响应幅度IIR 高通滤波器012345678x 104-1-0.500.51IIR 高通滤波器滤波前的时域波形012345678x 10412302IIR 高通滤波器滤波后的时域波形00.511.522.5x 104020406080IIR 高通滤波器滤波前的频谱频率/Hz幅值0.511.522.5x 104246810227IIR 高通滤波器滤波后的频谱频率/Hz幅值(3) 带通滤波器性能指标 fb1=1200Hz ,fb2=3000Hz ,fc1=1000Hz ,fc2=3200Hz ,As =100dB ,Ap =1dB 。
fs=44100;[x2]=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); Fs=22050;Ts=1/Fs;R1=30;fb1=1200;fb2=3000;fc1=1000;fc2=3200;fs=22050; W1=2*fb1*pi/fs;W2=2*fc1*pi/fs; W3=2*fb2*pi/fs;W4=2*fc2*pi/fs; Wp=[W1,W3]; Ws=[W2,W4]; Rp=1;Rl=100;Wp1=2/Ts*tan(Wp/2); %将模拟指标转换成数字指标 Ws1=2/Ts*tan(Ws/2);[N,Wn]=cheb2ord(Wp1,Ws1,Rp,R1,'s'); %选择滤波器的最小阶数[Z,P ,K]=cheb2ap(N,Rl); %创建切比雪夫模拟滤波器 [Bap,Aap]=zp2tf(Z,P ,K);[b,a]=lp2bp(Bap,Aap,2100*2*pi,1800*2*pi);[bz,az]=bilinear(b,a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,W]=freqz(bz,az); %绘制频率响应曲线 figure(1)plot(W*fs/(2*pi),abs(H))gridxlabel('频率/Hz')ylabel('频率响应幅度')title('IIR带通滤波器')f1=filter(bz,az,x1);figure(2)subplot(2,1,1)plot(x1) %画出滤波前的时域图title('IIR带通滤波器滤波前的时域波形');subplot(2,1,2)plot(f1); %画出滤波后的时域图title('IIR带通滤波器滤波后的时域波形');sound(f1,44100); %播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x1,1024);subplot(2,1,1);plot(f,abs(y2(1:512))); %画出滤波前的频谱图title('IIR带通滤波器滤波前的频谱')xlabel('频率/Hz');ylabel('幅值');subplot(2,1,2)plot(f,abs(F0(1:512))); %画出滤波后的频谱图title('IIR带通滤波器滤波后的频谱')xlabel('频率/Hz');ylabel('幅值');结果如下:0200040006000800010000120000.20.40.60.811.21.4频率/Hz频率响应幅度IIR 带通滤波器12345678x 104-1-0.500.51IIR 带通滤波器滤波前的时域波形12345678x 104-0.2-0.100.10.2IIR 带通滤波器滤波后的时域波形02000400060008000100001200020406080IIR 带通滤波器滤波前的频谱频率/Hz幅值20004000600080001000012000010203040IIR 带通滤波器滤波后的频谱频率/Hz幅值窗函数法设计FIR 数字滤波器1) 低通滤波器性能指标 fb =1000Hz ,fc =1200Hz ,As =100dB ,Ap =1dB fs=44100;x1=wavread('C:\Documents and Settings\Administrator\桌面\Windows XP 关机.wav'); wp=2*pi*5000/fs; ws=2*pi*4800/fs; Rp=1; Rs=100;wdelta=wp-ws;N=ceil(8*pi/wdelta); %取整 wn=(wp+ws)/2;[b,a]=fir1(N,wn/pi,'high'); figure(1)freqz(b,a,512);title('FIR 高通滤波器'); f2=filter(b,a,x1); figure(2) subplot(2,1,1) plot(x1)title('FIR 高通滤波器滤波前的时域波形');subplot(2,1,2) plot(f2);title('FIR 高通滤波器滤波后的时域波形');sound(f2,44100); %播放滤波后的语音信号 F0=fft(f2,1024); f=fs*(0:511)/1024; figure(3)y2=fft(x1,1024); subplot(2,1,1);plot(f,abs(y2(1:512)));title('FIR 高通滤波器滤波前的频谱') xlabel('频率/Hz'); ylabel('幅值');subplot(2,1,2)F2=plot(f,abs(F0(1:512)));title('FIR 高通滤波器滤波后的频谱') xlabel('频率/Hz'); ylabel('幅值');00.10.20.30.40.50.60.70.80.91-10-505x 104Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )00.10.20.30.40.50.60.70.80.91-150-100-50050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )FIR 高通滤波器12345678x 104-1-0.500.51FIR 高通滤波器滤波前的时域波形12345678x 104-0.04-0.0200.020.04FIR 高通滤波器滤波后的时域波形0.511.522.5x 104020406080FIR 高通滤波器滤波前的频谱频率/Hz幅值00.511.522.5x 104123FIR 高通滤波器滤波后的频谱频率/Hz幅值(3)频谱分析:经过FIR低通滤波器滤波后,将高频部分滤除了。