吹管音乐滤波去噪使用TUKEYWIN设计FIR滤波器

1、 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 1 页 共 17 页吹管音乐滤波去噪吹管音乐滤波去噪使用使用 TUKEYWIN 设计的设计的 FIR 滤波器滤波器 学生姓名：学生姓名： 指导老师：指导老师： 摘 要 本课程设计主要内容是设计利用窗口设计法选择 TUKEYWIN 窗设计一个 FIR滤波器，对一段含噪吹管音乐信号进行滤波去噪处理并根据滤波前后的波形和频谱分析滤波性能。本课程设计仿真平台为 MATLAB7.0，开发工具是 M 语言编程。首先下载一段吹管音乐信号，并人为加入一单频噪声，然后对信号进行频谱分析以确定所加噪声频率，并设计滤波器进行滤波去噪处理，最后比较

2、滤波前后的波形和频谱并进行分析。由分析结果可知，滤波器后的语音信号与原始信号基本一致，即设计的 FIR 滤波器能够去除信号中所加单频噪声，达到了设计目的。关键词 滤波去噪；FIR 滤波器；TUKEYWIN 窗；MATLAB1 引言引言本课程设计主要解决在含噪情况下对吹管音乐信号的滤波去噪处理，处理时采用的是利用窗口设计法选择 Tukeywin 窗设计的 FIR 滤波器。通过课程设计了解 FIR 滤波器设计的原理和步骤，掌握用 Matlab 语言设计滤波器的方法及编程。通过观察音乐信号滤波前后的时域波形的比较，加深对滤波器作用的理解。1.1 课程设计目的课程设计目的熟悉 Matlab 语言环境，

3、掌握 Matlab 语言的编程规则，利用 TUKEYWIN 窗函数设计法来设计符合要求的 FIR 滤波器来实现音乐信号的滤波去噪。并绘制滤波前后的时域波形和频谱图。根据图形分析判断滤波器设计的正确性。通过本次课程设计熟悉利用 TUKEYWIN 窗函数法设计 FIR 滤波器的过程。增强自己独立解决问题的能力，提高自己的动手能力。加深对理论知识联系实际问题的理解。为以后的工作奠定坚实的 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 2 页 共 17 页基础。1.2 课程设计要求课程设计要求（1）滤波器指标必须符合工程实际。（2）设计完后应检查其频率响应曲线是否满足指标。（3）

4、处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告书。1.3 课程设计平台课程设计平台MATLAB 软件包括五大通用功能：数值计算功能（Nemeric）；符号运算功能（Symbolic）；数据可视化功能（Graphic）；数据图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。其中，符号运算功能的实现是通过请求 MAPLE 内核计算并将结果返回到 MATLAB 命令窗口。该软件有三大特点：一是功能强大；二是界面友善、语言自然；三是开放性强。目前，Mathworks 公司已推出 30 多个应用工具箱。MATLAB 在线性代数、

5、矩阵分析、数值及优化、数理统计和随机信号分析、电路与系统、系统动力学、信号和图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统、以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。MATLAB 在信号与系统中的应用主要包括符号运算和数值计算仿真分析。由于信号与系统课程的许多内容都是基于公式演算，而 MATLAB 借助符号数学工具箱提供的符号运算功能能基本满足信号与系统课程的需求。例如，解微分方程、傅里叶正反变换、拉普拉斯正反变换、z 正反变换等。MATLAB 在信号与系统中的另一主要应用是数值计算与仿真分析，主要包括函数波形绘制、函数运算、冲激响应与阶跃响应仿真分析、信号

6、的时域分析、信号的频谱分析、系统的 S 域分析、零极点图绘制等内容。数值计算仿真分析可以帮助学生更深入理解信号与系统的理论知识，并为将来使用MATLAB 进行信号处理领域的各种分析和实际应用打下基础1。2 设计原理设计原理窗函数设计法的基本思想是用 FIRDF 逼近希望的滤波特性。设希望逼近的滤波器 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 3 页 共 17 页的频率响应函数为，其单位脉冲响应用表示。为了设计简单方便，通常()jwdHe( )dh n选择为具有片段常数特性的理想滤波器。因此是无限长非因果序列，不()jwdHe( )dh n能直接作为 FIRDF 的单位

7、脉冲响应。窗函数设计法就是截取为有限长的一段因果( )dh n序列，并用合适的窗函数进行加权作为 FIRDF 的单位脉冲响应2。( )dh n2.1 FIR 滤波器滤波器FIR 滤波器是有限长单位脉冲响应数字滤波器的简称，FIR 滤波器的最大优点是可以实现线性相位滤波。FIR 滤波器满足了在数字通信和图像传输与处理等应用场合对线性相位的要求，FIR 滤波器是全零点滤波器，硬件和软件实现结构简单，不用考虑稳定性问题。当幅频特性指标相同（不考虑相位特性）时，FIR 滤波器的阶数比 IIRDF的阶数高得多。所以在要求线性相位滤波的应用场合，一般都用 FIR 滤波器。FIR 数字滤波器设计的基本步骤如

8、下：（1）确定指标：在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于 FIR 滤波器的设计。第二种指标是相对指标。它以分贝值的形式给出要求。 （2）逼近：确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。 （3）性能分析和计算机仿真：上两步的结果是得到以差分或冲激响应描述的滤波器。根据这

9、个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。FIR 滤波器的设计问题实际上是确定能满足所要求的转移序列或脉冲响应的常数的问题，设计方法主要有窗函数、频率响应法和等波纹最佳逼近法等。2.2 窗口设计法窗口设计法 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 4 页 共 17 页窗口设计法是一种通过截断和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。通常在设计滤波器之前，应该先根据具体的工程应用确定滤波器的技术指标。在大多数实际应用中，数字滤波器常常被用来实现选频操作，所以指标

10、的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。窗口设计法步骤如下：（1）根据过渡带宽及阻带衰减要求，选择窗函数的类型并估计窗口长度 N。窗函数的类型可根据最小阻带衰减 AS 独立选择。（2）根据待求滤波器的理想频率响应求出理想单位脉冲响应 hd(n)。（3）由性能指标确定窗函数 W(n)和长度 N。（4）求得实际滤波器的单位脉冲响应 h(n)， h(n)即为所设计 FIR 滤波器系数向量b(n)。( )( ).( )dh nhn W n (2.1)常见的窗函数性能表如下图2.1所示：图2.1 常见窗函数性能表名称滤波器过渡带宽最小阻带衰减名称滤波器过渡带宽最小阻带衰减矩形1.8/M

11、21dBPARZENWIN6.6/M56dB巴特利特6.1/M25dBFLATTOPWIN19.6/M108dB汉宁6.2/M44dBGAUSSWIN5.8/M60dB汉明6.6/M51dBBARTHANNWIN3.6/M40dB布莱克曼11/M74dBBLACKMANHARRIS16.1/M109dBBOHMANWIN5.8/M51.5dBCHEBWIN15.2/M113dBNUTTALLWIN15.4/M108dBTUKEYWIN2.4/M22dB2.3 Tukeywin 窗窗 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 5 页 共 17 页w =tukeywin(

12、M)在数组 w 中产生 M 点的 Tukeywin 窗函数。由上面查表可知：Tukeywin 窗的滤波器的过渡带宽为 2.4/M，最小阻带衰减22db。利用可绘制出 Tukeywin 窗函数波形如下：wvtoolN=128;t=tukeywin(N);wvtool(t);图 2-2 Tukeywin 窗函数波形3 设计步骤设计步骤3.1 设计流程图设计流程图根据设计的要求，首先从网上下载一段音乐信号，对音乐信号加入噪声干扰，再利用 Tukeywin 窗设计合理的 FIR 滤波器。最后用滤波器对干扰后的音乐信号进行滤波去噪。具体设计流程图如下图 3.1 所示： 吹管音乐滤波去噪使用 TUKEYW

13、IN 设计的 FIR 滤波器 第 6 页 共 17 页 NO YES 图 3-1 音乐信号滤波去噪流程图3.2 录制语音信号录制语音信号在网上下载一段音乐，再利用视频转换器将其转换成.格式，且为单声道，3mpwav再将此.格式音乐控制在 3 秒内，以减少设计中的误差。然后在 Matlab 软件平台下，wav利用函数 wavread 对语音信号进行采样，记住采样频率和采样点数，采集完成后在信号中加入一个单频噪声，绘制原音乐信号和加噪后的音乐信号的时域和频域的波形图。开始下载一段吹管音乐信号，命名为“meilideshenhua.wav” ，绘制出其时域波形和频谱图。加入噪声干扰，绘制干扰后的时域

14、波形和频谱图。利用 Tukeywin 窗设计合理的 FIR 滤波器。将干扰后的音乐信号通过自己设计的滤波器，进行滤波去噪。比较滤波前后的时域波形和频谱图，并回放音乐信号，验证是否达到去噪效果。结束是否达到去噪效果？ 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 7 页 共 17 页具体调用如下：x,fs,bits=wavread(g:数字信号meilideshenhua.wav); % 输入参数为文件的全路径和文件名，输出的第一个参数是每个样本的值，fs 是生成该波形文件时的采样率，bits 是波形文件每样本的编码位数。sound(x,fs,bits); % 按指定的采样

15、率和每样本编码位数回放N=length(x); % 计算信号 x 的长度fn=2100; % 单频噪声频率，此参数可改t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率x=x(:,1); y=x+0.1*sin(fn*2*pi*t); sound(y,fs,bits); % 应该可以明显听出有尖锐的单频啸叫声X=abs(fft(x);Y=abs(fft(y); % 对原始信号和加噪信号进行 fft 变换,取幅度谱X=X(1:N/2); Y=Y(1:N/2);% 截取前半部分deltaf=fs/N; % 计算频谱的谱线间隔f=0:deltaf:fs/2-deltaf; %

16、 计算频谱频率范围得到的波形图如图 3.2 所示： 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 8 页 共 17 页图 3-2 原始音乐信号与加噪后的音乐信号时域、频域波形图3.3 滤波器设计滤波器设计在此用窗函数法设计一个带阻的滤波器。FIR此处给出各滤波器的性能指标： ，As20dB，Ap1 dB，其截止频率在单频噪声干扰11222000205021502200psspfHzfHzfHzfHz附近。再进行加窗设计，利用 Matlab 中的函数 freqz 画出各滤波器的频率响应。具体调用如下：fpd=2000;fsd=2050;fsu=2150;fpu=2200;R

17、p=1;As=20; % 带阻滤波器设计指标fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=min(fsd-fpd),(fpu-fsu); % 计算上下边带中心频率，和频率间隔wcd=fcd/fs*2*pi;wcu=fcu/fs*2*pi;dw=df/fs*2*pi; % 将 Hz 为单位的模拟频率换算为 rad 为单位的数字频率wsd=fsd/fs*2*pi;wsu=fsu/fs*2*pi;M=ceil(2.4*pi/dw)+1; % 计算 TUKEYWIN 窗设计该滤波器时需要的阶数n=0:M-1; % 定义时间范围w_tuk=tukeywin(M); % 产生 M

18、阶的 TUKEYWIN 窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)-ideal_lp(wcu,M); % 调用自编函数计算理想带阻滤波器的脉冲响应h_bs=w_tuk.*hd_bs; % 用窗口法计算实际滤波器脉冲响应db,mag,pha,grd,w=freqz_m(h_bs,1); % 调用自编函数计算滤波器的频率特性得到的波形图如图 3.3 所示： 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 9 页 共 17 页图 3-3 滤波器波形图注：滤波器的设计指标，图中可以明显地看到利用窗函数所20sdBTUKEYWIN设计的滤波器基本上可

19、以达到滤波效果。3.4 信号滤波处理信号滤波处理用自己设计的各滤波器分别对采集的信号进行滤波，在 Matlab 中，FIR 滤波器利用函数 fftfilt 对信号进行滤波，IIR 滤波器利用函数 filter 对信号进行滤波。在同一张大图里分别绘制原始信号 x，加噪信号 x，滤波去噪信号 y 的时域波形和频谱，以便比较和分析。滤波程序如下：y_fil=fftfilt(h_bs,y); % 用设计好的滤波器对 y 进行滤波Y_fil=fft(y_fil);Y_fil=Y_fil(1:N/2); % 计算频谱取前一半得到的波形图如图 3.4 所示： 吹管音乐滤波去噪使用 TUKEYWIN 设计的

20、FIR 滤波器 第 10 页 共 17 页图 3.4 滤波前后音乐信号的时域波形图和频谱图对比图注：从图中可以看出滤波后的信号的幅度有所减小，说明加噪信号通过滤波器之后不可能毫无变化，也就是说参数的设置以及窗函数的选择都会影响滤波后信号的质量。3.5 结果分析结果分析在 Matlab 中，函数 sound 可以对经过 tukeywin 窗设计的 FIR 滤波器之后的音乐信号进行回放。其调用格式：sound (y_fil,fs,bits)；可以感觉滤波前后的声音有变化，那尖锐的噪声明显减小。但是所设计的滤波器并不是理想的滤波器，并且由于 tukeywin 窗所设计的滤波器本身的性能局限，所以在回

21、放滤波后的音乐的时候还可以听到一些噪声。但从滤波器的波形图中也可以明显地看到滤波器的设计满足所定义的指标，从另一个方面也可以断言滤波器的设计成功。与相同信源不同滤波方法（Nuttallwin 窗）的同学的滤波器相对比，可以看出由于 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 11 页 共 17 页其衰减偏小，所以性能也较差。与不相同信源（拉弦音乐）相同滤波方法的同学的滤波器相对比，可以看出其滤波效果基本相近。4 出现的问题及解决方法出现的问题及解决方法1、 在录音时，没有将录音的属性改为 wav 的格式，当在 Matlab 软件平台上调用时，出现无法调用的提示，最后通过

22、与同学讨论发现，不能直接更改文件后缀，应用转换器将其转换为 wav 格式。2、在开始录制音乐信号并将其导入 MATLAB 中的时候出现过错误，原因是我所下载的音乐信号是双音频信号，不符合要求，在老师的指导下我将音乐信号变成了单声道信号，再次导入的时候 MATLAB 不再报错。 3、在最后对加噪音乐信号进行滤波并绘图的时也出现了错误，绘制出最后的滤波信号的频谱图时纵坐标出现了负值，但是原信号的频谱图中却没有，回头检查了一遍程序想到幅度值只有正值，想到是不是忘记加绝对值符号了，在对滤波的信号进行 fft 变换之后未取绝对值，改正之后再运行程序就不会再出错。4、在滤波处理时，在加入单频噪声部分噪声频

23、率设置为 2200，但滤波器参数设置部分最高带阻为 2200，导致滤波处理后不能滤除噪声信号，而后来将其改为 2100，才达到了效果。5 结束语结束语本次的课程设计，我的任务是利用 Tukeywin 窗函数设计 FIR 滤波器对音乐信号滤波去噪。在本次课程设计之前，我对 Tukeywin 窗函数完全没有了解，因此在看到这个题目时，我是一头雾水。但是通过自己翻阅资料和询问同学，我掌握了用 Tukeywin 窗函数设计 FIR 滤波器的方法步骤，了解了窗函数的基本设计流程。经过这两周的课程设计经过为期两周的课程设计，这次课程设计比较的简单，根据老师给出的模板，再结合以前学过的知识，并 在同学和老师

24、的帮助下我顺利的完成了任务。不同于在教室里上的理论，这次的课程设计需要将我们平时所学习的知识运用到实践之中，将知识学以致用。因为是以所学理论为基础，所以在课程设计的过程中，我又重温滤波器、窗函数等知识，更加熟悉了 MATLAB 的操作。 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 12 页 共 17 页虽说总体上是比较的简单，但是也会出现一些错误，分析并总结这些错误主要是语句上的一些用法出现错误，以及没有及时地掌握老师所提醒的部分。但是最后还是在老师的帮助下，解决了这些问题。课程设计是我们运用所学知识，动手实践的一个很好的机会。它既可以帮助我们加深对所学知识的理解，又能

25、提高我们运用知识，联系实际，动手实践的能力。而且在设计过程中可能用到我们没学过的知识，需要我们去查阅资料获取相关信息，这又提高了我们查找信息和学习新知识的能力。在实物的调试与检测过程中，又会遇到许多意想不到的问题，需要我们去分析原因和解决问题。也体会到真正的去独立地完成一件事情是很困难的，同学以及老师的帮助和提醒是必须的。通过这次课程设计，我拓宽了知识面，锻炼了实际操作能力，综合素质也得到了提高，进一步加深了了我们对专业的认识和激发了我们对专业的兴趣。虽然课程设计结束了，但是我们的学习还没结束，对知识的进一步学习还需要继续，很开心成功地完成了这次设计。 吹管音乐滤波去噪使用 TUKEYWIN

26、设计的 FIR 滤波器 第 13 页 共 17 页参考文献参考文献1 维纳.K 英格尔，约翰.G 普罗克斯.数字信号处理（MATLAB 版）.西安：西安交通大学出版社，20072 程佩青.数字信号处理教程.北京：清华大学出版社，20023 张圣勤.MATLAB7.0 实用教程.北京：机械工业出版社，20084 张志涌.精通 MATLAB 6.5 版M北京：北京航空航天大学出版社，20035谢德芳.数字信号处理M.北京：科学出版社，20056张小虹.信号系统与数字信号处理M.第版.西安：西安电子科技出版社，2002 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 14 页 共

27、 17 页附录 1：吹管音乐滤波去噪吹管音乐滤波去噪使用使用 TUKEYWIN 窗窗设计的设计的 FIR 滤波器滤波器程序一：%放出原音频x,fs,bits=wavread(h:数字信号meilideshenhua.wav); % 输入参数为文件的全路径和文件名，%输出的第一个参数是每个样本的值，%fs 是生成该波形文件时的采样率，%bits 是波形文件每样本的编码位数。sound(x,fs,bits); % 按指定的采样率和每样本编码位数回放程序二：%原音乐信号和加噪后的音乐信号的时域和频域的波形图绘制N=length(x); % 计算信号 x 的长度fn=2100; % 单频噪声频率，此参

28、数可改t=0:1/fs:(N-1)/fs; % 计算时间范围，样本数除以采样频率x=x(:,1); y=x+0.1*sin(fn*2*pi*t); sound(y,fs,bits); % 应该可以明显听出有尖锐的单频啸叫声X=abs(fft(x);Y=abs(fft(y); % 对原始信号和加噪信号进行 fft 变换,取幅度谱X=X(1:N/2); Y=Y(1:N/2);% 截取前半部分deltaf=fs/N; % 计算频谱的谱线间隔f=0:deltaf:fs/2-deltaf; % 计算频谱频率范围subplot(2,2,1);plot(t,x);xlabel(时间(t);ylabel(幅度

29、); title(原始音乐信号);axis(0,2.5,-1.5,1.5);grid on;subplot(2,2,2);plot(f,X);xlabel(频率(f);ylabel(幅度谱);title(原始音乐信号幅度谱);axis(0,3000,0,3000);grid on;subplot(2,2,3);plot(t,y);xlabel(时间(t);ylabel(幅度);title(加干扰后的音乐信号); 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 15 页 共 17 页axis(0,2.5,-1.5,1.5);grid on;subplot(2,2,4);plo

30、t(f,Y);xlabel(频率(f);ylabel(幅度谱);title(加干扰后的音乐信号幅度谱);axis(0,3000,0,3000);grid on;程序三：%带阻滤波器设计fpd=2000;fsd=2050;fsu=2150;fpu=2200;Rp=1;As=20; % 带阻滤波器设计指标fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=min(fsd-fpd),(fpu-fsu); % 计算上下边带中心频率，和频率间隔wcd=fcd/fs*2*pi;wcu=fcu/fs*2*pi;dw=df/fs*2*pi; % 将 Hz 为单位的模拟频率换算为 rad 为单

31、位的数字频率wsd=fsd/fs*2*pi;wsu=fsu/fs*2*pi;M=ceil(2.4*pi/dw)+1; % 计算 Tukeywin 窗设计该滤波器时需要的阶数n=0:M-1; % 定义时间范围w_tuk=tukeywin(M); % 产生 M 阶的 Tukeywin 窗 hd_bs=ideal_lp(wcd,M)+ideal_lp(pi,M)-ideal_lp(wcu,M); % 调用自编函数计算理想带阻滤波器的脉冲响应h_bs=w_tuk.*hd_bs; % 用窗口法计算实际滤波器脉冲响应db,mag,pha,grd,w=freqz_m(h_bs,1); % 调用自编函数计算滤

32、波器的频率特性subplot(2,2,1);plot(w/pi,db);title(FIR 滤波器的幅度响应图);xlabel(w/pi);ylabel(db);line(0,0.5,-1,-1,Color,r,LineWidth,2,LineStyle,-);line(0,0.5,-20,-20,Color,r,LineWidth,2,LineStyle,-);line(wsd/pi,wsd/pi,-30,15,Color,r,LineWidth,2,LineStyle,-);line(wsu/pi,wsu/pi,-30,15,Color,r,LineWidth,2,LineStyle,-)

33、;axis(0.05 0.15 -30 15);grid on;subplot(2,2,2);plot(w/pi,mag);title(FIR 滤波器的幅度响应图); 吹管音乐滤波去噪使用 TUKEYWIN 设计的 FIR 滤波器 第 16 页 共 17 页xlabel(w/pi);ylabel(幅度 mag);axis(0 0.3 0 1.5);grid on;subplot(2,2,3);plot(w/pi,pha);title(滤波器相位响应图);xlabel(w/pi);ylabel(相位 pha);axis(0 1 -4 4);grid on;subplot(2,2,4);stem(

34、n,h_bs);title(滤波器脉冲响应图);xlabel(n);ylabel(h(n);axis(0 1000 0 0.7);grid on;程序四：%对信号进行滤波y_fil=fftfilt(h_bs,y); % 用设计好的滤波器对 y 进行滤波Y_fil=abs(fft(y_fil);Y_fil=Y_fil(1:N/2); % 计算频谱取前一半subplot(3,2,1);plot(t,x);title(原始音乐信号时间 x);grid on;xlabel(时间:t);ylabel(幅度);axis(0,2.5,-1.5,1.5);subplot(3,2,2);plot(f,X);title(原始音乐信号幅度谱 X);xlabel(频率:f);ylabel(幅度);axis(0 2500 0 2500);grid on;subplot(3,2,3);plot(t,y);title(加干扰音乐信号时间 x1);grid on;xlabel(时间 t);ylabel(幅度);axis(0,2.5,-1.5,1.5);subplot(3,2,4);plot(f,Y);title(加干扰后的音乐信号幅度谱 X1);xlabel(频率 f);ylabe