编程实现任意确定信号的频谱分析算法

1、广西科技大学数字信号处理课程设计说明书题目：编程实现任意确定信号的 频谱分析算法系 别：计算机工程学院专业班级：通信学 号：学生姓名: 指导教师:目录 TOC o 1-5 h z 摘要 3 HYPERLINK l bookmark6 o Current Document 一、设计内容 3 HYPERLINK l bookmark8 o Current Document 二、设计原理 4 HYPERLINK l bookmark12 o Current Document 三、设计过程7和弦音音频文件的频谱分析7.对该信号频谱能量较集中的频带滤波9.比较滤波前后的音频文件 21.对滤波后的音频信号

2、再滤出三个能量最集中的频簇21.重建信号与原信号的音频进行声音回放比较 336,分析什么是和弦音35 HYPERLINK l bookmark14 o Current Document 收获 3636参考文献摘要： 随着计算机和信息科学的飞速发展，信号处理逐渐发展成为一门独立的学科，成为信息科学的重要组成部分，在语音处理、雷达、 图像处理、 通信、 生物医学工程等众多领域中得到广泛应应用。Matlab 语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab 功能强大、简单易学、编程效率高，深受广大科技工作者的喜爱。 特别是 Matlab 还具有信号分析工具箱， 不需具备很强的编

3、程能力，就可以很方便地进行语音信号分析、处理和设计。数字信号处理课程在现代科学中具有很大重要性及自身特点，理解与掌握课程中的基本概念、基本原理、基本分析方法，对用 Matlab 进行数字信号处理课程设计的思路，具有很大帮助。语音信号的处理与滤波的设计主要是用 Matlab 作为工具平台， 设计中涉及到声音的录制、播放、存储和读取，语音信号的抽样、频谱分析，滤波器的设计及语音信号的滤波，通过数字信号处理课程的理论知识的综合运用。从实践上初步实现对数字信号的处理。关键词：抽样频率；频谱分析；滤波器 ；窗函数一、设计内容（1）对给定的CEGf弦音音频文件取合适长度的采样记录点，然后进行频谱分析（信号

4、的时域及幅频特性曲线要画出） 。（2）分析CEGf口弦音频谱特点，对该信号频谱能量相对较为集中的频带（分低、中、高频）实现滤波（分别使用低通，带通及高通）显示滤波后信号的时域和频域曲线，并对滤波后的信号与原信号的音频进行声音回放比较。（ 3）在低、中、高三个频带中，各滤出三个能量最集中的频簇，显示滤波后信号的时域和频域曲线。（ 4）任意选择几个滤出的频带（或频簇）进行时域信号重建（合成） ，与原信号的音频进行声音回放比较。讨论：根据上述结果，分析什么是和弦音二、设计原理数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR ）滤波器和有限长冲激响应（FIR）滤波器。IIR

5、滤波器的特征是，具有无限持续时间冲激响应。种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和 Chebyshev 逼近法等等。随着MATLA歌件尤其是MATLAB勺信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能， 而且还可以使设计达到最优化 .线性相位FIR滤波器通常采用窗函数法设计。窗函数法设计FIR滤波器的基本思想是：根据给定的滤波器技术指标，选择滤波器长度N和窗函数3 （ n）,使其具有

6、最窄宽度的主瓣和最小的旁瓣。其核心是从给定的频率特性，通过加窗确定有限长单位脉冲响应序列 h（n） 。工程中常用的窗函数共有6种，即矩形窗、巴特利特(Bartlett )窗、 汉宁(Hanning)窗、汉明(Hamming窗、布莱克曼(Blackman)窗 和凯塞(Kaiser )。FIR数字滤波器的设计原理：采用窗口法线性相位实系数FIR滤波器按其N值奇偶和h(n)的奇偶对称性分为四种：1、h(n)为偶对称，N为奇数2 修司2H(ej(o)的幅值关于=0,兀，2兀成偶对称。2、h(n)为偶对称，N为偶数H(ej)的幅值关于=兀成奇对称，不适合作高通 3、h(n)为奇对称，N为奇数_ 1_J

7、出 _H（严）=Z力- +门）皿耳中）记2 -H(ej(o)的幅值关于=0,兀，2兀成奇对称，不适合作高通和低 通。4、h(n)为奇对称，N为偶数nnH（产）一Z2M n）dn（m（盟-n=i 2H(ej 3 ) 3 =0、2兀=0,不适合作低通窗函数法设计线性相位FIR滤波器步骤:?确定数字滤波器的性能要求：临界频率(ok,滤波器单位脉冲响应长度N;?根据性能要求，合理选择单位脉冲响应h(n)的奇偶对称性，从而确定理想频率响应Hd(ej 3)的幅频特性和相频特性；?求理想单位脉冲响应hd(n),在实际计算中，可对Hd(ej 3)按M(M远大于N)点等距离采样，并对其求IDFT得hM(n),用

8、hM(n)代替 hd(n);?选择适当的窗函数w(n),根据h(n)= hd(n)w(n)求所需设计的FIR滤波器单位脉冲响应；?求H(ej 3),分析其幅频特性，若不满足要求，可适当改变窗 函数形式或长度N,重复上述设计过程，以得到满意的结果。窗函数的傅式变换 W(ej 3 )的主瓣决定了 H(ej 3)过渡带宽。 W(ej 3 )的旁瓣大小和多少决定了 H(ej 3 )在通带和阻带范围内波动 幅度，常用的几种窗函数有：?矩形窗 w(n)=RN(n);? Hanning 窗怏5) = 0 51 - c&s 巴里打;1? Hamming窗做耳)=网54-0 46。$ 二1)衣将(司)；N 1?

9、 Blackmen窗 械=0一42-必)+。口8匚旌(仍)衣笠；N - 1J - 1? Kaiser 窗 用用- 5Hq。式中Io(x)为零阶贝塞尔函数和弦音音频文件的频谱分析三、设计过程采集已录声音信号长度计算时间振幅 ););对给定的 CEG 和弦音音频文件进行频谱分析， 画出信号的时域及 幅频特性曲线，程序如下：y,fs,bits=wavread(); %N=length(y); % t=(0:length(y)-1)/fs; % figure;plot(t,y);xlabel( 时间 );ylabel( title( 语音信号的时域图Y=fft(y); % 对信号进行FFTdf=fs/

10、N;%采样间隔n1=0:N/2-1;f=n1*df;%频带宽度figure;plot(f,abs(Y(1:length(f);xlabel(频率);ylabel( 振幅);title(语音信号的频谱图 );sound(y); %播放给定的CEGf口谐音频文件图形如下：语音信号的频谱图频率2. 对该信号频谱能量较集中的频带滤波对该信号频谱能量相对较为集中的频带(分低、中、高频)实现滤波(分别使用低通，带通及高通)(1)FIR 低通Ws=2*1200*pi/8000;Wp=2*1000*pi/8000;As=100;Rp=1;% 数字低通滤波器参数Wc=(Wp+Ws)/2;%WcdeltaW=Ws

11、-Wp; % 过度带 WX两个过渡带中的小者 N0=ceil*pi/deltaW);N=N0+mod(N0+1,2); % 根据表格 7-3 求出响应的窗口的点数N,保证N为奇数hd=ideallp(Wc,N); % 低通滤波器的 hd( n )%-ideallp ()函数(非系统自有函数)在系统安装目录的WORK目录中，wdWindow=kaiser(N,; % 凯泽窗hr=hd.*wdWindow; %- 点乘n=0:N-1;subplot(2,2,1);stem(n,wdWindow); % 画二维离散图形xlabel( 时间 );ylabel( 振幅 ); title( 凯泽窗 );

12、%x轴说明； y 轴说明%标题画二维离散图形subplot(2,2,2);stem(n,hr);xlabel( 时间 );ylabel(振幅 );%x轴说明； y 轴说明title( 凯泽窗的冲激响应);% 标题H,W=freqz(hr,1); % 求每个窗对应的频率响应subplot(2,2,3);plot(W/pi,abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 绝对振幅 ); %x 轴说明； y 轴说明title( 滤波器低通幅频特性1);%标题legend( 凯泽窗 );subplot(2,2,4);plot(W/pi, 20*log10(abs(H

13、); % 画二维连续图形xlabel(omega八pi);ylabel(振幅(分贝)；x轴说明；y 轴说明title( 滤波器低通幅频特性2);%标题legend( 凯泽窗 );y,fs,bits=wavread(”)； 读出给定的CEGf口谐音频信号N=length(y);t=(0:length(y)-1)/fs;df=fs/N;% 采样间隔n1=0:N/2-1;f=n1*df;% 频带宽度ly1=fftfilt(hr,y);% 信号送入滤波器滤波， ly1 为输出figure; % 打开新的画图窗口subplot(2,2,1);plot(t,y);xlabel( 时间 ); ylabel(

14、 振幅 );title( 语音信号的时域图 );Y1=fft(y); % 滤波前的FFT变换subplot(2,2,3);plot(f,abs(Y1(1:length(f); % 画滤波前频谱图xlabel( 频率 ); ylabel( 振幅 );title( 语音信号的频谱图 );grid;subplot(2,2,2);plot(t,ly1);xlabel( 时间 ); ylabel( 振幅 );title( 滤波后语音信号的时域图 );Ly1=fft(ly1); % 滤波后的FFT变换subplot(2,2,4);plot(f,abs(Ly1(1:length(f); % 画滤波后频谱图x

15、label( 频率 ); ylabel( 振幅 );title( 滤波后语音信号的频谱图 );grid;wavwrite(ly1,fs, 低通 .wav);% 把滤波后的信号写入文件凯泽窗0.3滤波器低通幅频特性11.5-0.1 0凯泽窗的冲激响应0.20.10100200时间滤波器低通幅频特性2300100凯泽窗1 56 幅振对绝0.5/-200 00-1000.5/语音信号的时域图0.50.500-0.5-0.5-10.51.5幅 振-1 0)0.511.5时间滤波后语音信号的频谱图滤波后语音信号的时域图频率幅 振时间语音信号的频谱图3002001000 01000200030004000

16、频率2) FIR 带通W1s=2*1000*pi/8000;W1p=2*1200*pi/8000;As=100;Rp=1;W2s=2*3200*pi/8000;W2p=2*3000*pi/8000;% 数字带通滤波器参数W1c=(W1p+W1s)/2;W2c=(W2s+W2p)/2; %W1c,W2cdeltaW=min(W1p-W1s),(W2s-W2p); % 过度带A W取两个过渡带中的小者N0=ceil*pi/deltaW);N=N0+mod(N0+1,2); % 根据表格 7-3 求出响应的窗口的点数N,保证N为奇数hd=ideallp(W2c,N)-ideallp(W1c,N);w

17、dWindow=kaiser(N,; % hr=hd.*wdWindow; %- 点乘 n=0:N-1; subplot(2,2,1); stem(n,wdWindow); xlabel( 时间 );ylabel( title( 凯泽窗 ); % subplot(2,2,2); stem(n,hr);%xlabel( 时间 );ylabel( 振幅 );% 带通滤波器的 hd( n) 凯泽窗%振幅 );标题画二维离散图形%x 轴说明； y 轴说明画二维离散图形%x 轴说明； y 轴说明title( 凯泽窗的冲激响应 ); % 标题H,W=freqz(hr,1); % 求每个窗对应的频率响应su

18、bplot(2,2,3);plot(W/pi,abs(H); %画二维连续图形xlabel(omega/pi);ylabel( 绝对振幅 ); %x 轴， y 轴说明title( 滤波器带通幅频特性1);%标题legend( 凯泽窗 );subplot(2,2,4);plot(W/pi,20*log10(abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 振幅 ( 分贝 ); %x 轴， y 轴说明title( 滤波器带通幅频特性2); % 标题legend( 凯泽窗 );y,fs,bits=wavread() %读出给定的 CEGf频文件N=length(y)

19、;t=(0:length(y)-1)/fs;df=fs/N;% 采样间隔n1=0:N/2-1;f=n1*df;% 频带宽度ly2=fftfilt(hr,y);% 信号送入滤波器滤波， ly2 为输出figure; % 打开新的画图窗口 subplot(2,2,1);plot(t,y);xlabel( 时间 ); ylabel( 振幅 );title(语首信号的时域图)；Y2=fft(y); %滤波前的FFT变换subplot(2,2,3); plot(f,abs(Y2(1:length(f); %画滤波前频谱图xlabel(频率);ylabel( 振幅)； title(语音信号的频谱图);gr

20、id； subplot(2,2,2); plot(t,ly2); xlabel(时间)； ylabel(振幅)； title(滤波后语音信号的时域图);Ly2=fft(ly2); %滤波后的FFT变换subplot(2,2,4); plot(f,abs(Ly2(1:length(f); %画滤波后频谱图xlabel(频率);ylabel( 振幅)； title(滤波后语音信号的频谱图);grid； wavwrite(ly2,fs, 带通.wav);%把滤波后的信号写入文件 sound(ly2,fs);%回放带通滤波后的CEGffi谐音频文件凯泽窗幅0.5振贝少幅振O30/00200幅振对绝语音

21、信号的时域图10.50-0.5滤波后语音信号的时域图00.511.52时间幅 振时间语音信号的频谱图频率-1频率高通滤波器的hd (n) 凯泽窗画二维离散图形%x轴,y轴说明标题画二维离散图形%x 轴.y轴说明% 标题(3) FIR高通Ws=2*3000*pi/8000;Wp=2*3200*pi/8000;As=100;Rp=1; % 数字高通滤波器 参数Wc=(Wp+Ws)/2;%WcdeltaW=Wp-Ws; %过度带A W取两个过渡带中的小者N0=ceil*pi/deltaW);N=N0+mod(N0+1,2); %根据表格 7-3 求出响应的窗口的点数N,保证N为奇数 hd=ideal

22、lp(pi,N)-ideallp(Wc,N); % wdWindow=kaiser(N,; % hr=hd.*wdWindow; %- 点乘 n=0:N-1;subplot(2,2,1);stem(n,wdWindow); % xlabel(时间);ylabel( 振幅); title( 凯泽窗)；%subplot(2,2,2);stem(n,hr);%xlabel(时间);ylabel( 振幅); title(凯泽窗的冲激响应);H,W=freqz(hr,1); %求每个窗对应的频率响应subplot(2,2,3);画二维连续图形plot(W/pi,abs(H); %xlabel(omega

23、八pi);ylabel(绝对振幅)；%x 轴,y轴说明title( 滤波器高通幅频特性1);%标题legend( 凯泽窗 );subplot(2,2,4);plot(W/pi,20*log10(abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 振幅 ( 分贝 ); %x 轴 ,y 轴说明title( 滤波器高通幅频特性2);% 标题legend( 凯泽窗 );%在一个窗口同时画出滤波前后的波形及频谱y,fs,bits=wavread(); % 读出信号N=length(y);t=(0:length(y)-1)/fs;df=fs/N;% 采样间隔n1=0:N/2

24、-1;f=n1*df;% 频带宽度ly3=fftfilt(hr,y);% 信号送入滤波器滤波， ly3 为输出figure; % 打开新的画图窗口subplot(2,2,1);plot(t,y);xlabel(时间);ylabel(振幅);title( 语音信号的时域图 );Y3=fft(y); %滤波前的FFT变换subplot(2,2,3);plot(f,abs(Y3(1:length(f); % 画滤波前频谱图xlabel( 频率 ); ylabel( 振幅 );title( 语音信号的频谱图 );grid;subplot(2,2,2);plot(t,ly3);xlabel( 时间 );

25、 ylabel( 振幅 );title( 滤波后语音信号的时域图 );Ly3=fft(ly3); %滤波后的FFT变换subplot(2,2,4);plot(f,abs(Ly3(1:length(f); % 画滤波后频谱图xlabel( 频率 ); ylabel( 振幅 );title( 滤波后语音信号的频谱图 );grid;wavwrite(ly3,fs, 高通 .wav);% 把滤波后的信号写入文件sound(ly3,fs);%回放高通滤波后的CEGT频文件凯泽窗的冲激响应202000300210O00幅振凯泽窗幅振0000-2幅振对绝语音信号的时域图00.511.52时间频率语音信号的频

26、谱图频率. 比较滤波前后的的音频文件通过低、 带、 高通滤波后的音频文件与原音频文件有如下简单差异：低通滤波后，CEGf 口谐音频文件的声音较原来变得低沉了。带通滤波后，CEGf 口谐音频文件的声音没原来的那么尖锐。高通滤波后，CEGf 口谐音频文件的声音变得尖锐清亮。对滤波后的音频信号再滤出三个能量最集中的频簇实现低、带、高通滤波后，在三个频带中，再分别用带通滤波器 滤出三个能量最集中的频簇，显示滤波后信号的时域和频域曲线。1)低通频簇W1s=2*100*pi/8000;W1p=2*200*pi/8000;As=100;Rp=1;W2s=2*900*pi/8000;W2p=2*800*pi/

27、8000; % 数字带通滤波器参数W1c=(W1p+W1s)/2;W2c=(W2s+W2p)/2; %W1c,W2cdeltaW=min(W1p-W1s),(W2s-W2p); % 过度带A W取两个过渡带中的小者N0=ceil*pi/deltaW);N=N0+mod(N0+1,2); % 根据表格 7-3 求出响应的窗口 的点数N,保证N为奇数hd=ideallp(W2c,N)-ideallp(W1c,N); wdWindow=kaiser(N,; % hr=hd.*wdWindow; %- 点乘 n=0:N-1; subplot(2,2,1); stem(n,wdWindow); xlab

28、el( 时间 );ylabel( title( 凯泽窗 ); % subplot(2,2,2); stem(n,hr);%xlabel( 时间 );ylabel( 振幅 );% 带通滤波器的 hd( n) 凯泽窗%振幅 );标题画二维离散图形%x 轴， y 轴说明画二维离散图形%x 轴， y 轴说明title( 凯泽窗的冲激响应 ); % 标题H,W=freqz(hr,1); % 求每个窗对应的频率响应subplot(2,2,3);plot(W/pi,abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 绝对振幅 ); %x 轴， y 轴说明%标题title（ 滤

29、波器带通幅频特性1）;legend（ 凯泽窗 ）;subplot(2,2,4);plot(W/pi,20*log10(abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 振幅 ( 分贝 ); %x 轴， y 轴说明title( 滤波器带通幅频特性2); % 标题信号送入滤波器滤波， ly11 为输出打开新的画图窗口legend( 凯泽窗 );y1,fs,bits=wavread( t=(0:length(y1)-1)/fs; df=fs/N;% 采样间隔 n1=0:N/2-1;f=n1*df;% 频带宽度 ly11=fftfilt(hr,y1);% figure

30、; % subplot(2,2,1);低通 .wav）; % 读出低通滤波后的信号plot(t,y1);xlabel( 时间 ); ylabel( 振幅 );title( 语音信号的时域图 );Y11=fft(y1);% 滤波前的FFT变换subplot(2,2,3);plot(f,abs(Y11(1:length(f); % 画滤波前频谱图 xlabel( 频率 ); ylabel( 振幅 ); title( 语音信号的频谱图 ); grid; subplot(2,2,2);plot(t,ly11);xlabel( 时间 ); ylabel( 振幅 );title( 滤波后语音信号的时域图

31、);Ly11=fft(ly11);%滤波后的 FFT变换subplot(2,2,4);plot(f,abs(Ly11(1:length(f); % 画滤波后频谱图 xlabel( 频率 ); ylabel( 振幅 ); title( 滤波后语音信号的频谱图 );grid;wavwrite(ly11,fs, 低通 );% 把滤波后的信号写入文件sound(ly11,fs);% 回放凯泽窗1 56 幅振对绝0 100.51语音信号的时域图00.511.52滤波后语音信号的时域图10.50-0.5-1时间语音信号的频谱图频率00.511.52时间滤波后语音信号的频谱图 3002001000 0100

32、0200030004000频率W1s=2*1100*pi/8000;W1p=2*1200*pi/8000;As=100;Rp=1;W2s=2*2100*pi/8000;W2p=2*2000*pi/8000;% 数字带通滤波器参数W1c=(W1p+W1s)/2;W2c=(W2s+W2p)/2; %W1c,W2cdeltaW=min(W1p-W1s),(W2s-W2p); % 过度带 W 取两个过渡带中的小者N0=ceil*pi/deltaW);N=N0+mod(N0+1,2); % 根据表格 7-3 求出响应的窗口的点数N,保证N为奇数hd=ideallp(W2c,N)-ideallp(W1c,

33、N); % 带通滤波器的 hd (n)wdWindow=kaiser(N,; %- 点乘hr=hd.*wdWindow;n=0:N-1;subplot(2,2,1);stem(n,wdWindow);%xlabel( 时间 );ylabel(振幅 );title( 凯泽窗 ); % 标题subplot(2,2,2);stem(n,hr);%xlabel( 时间 );ylabel( 振幅 );凯泽窗画二维离散图形%x 轴， y 轴说明画二维离散图形%x 轴， y 轴说明H,W=freqz(hr,1); % 求每个窗对应的频率响应subplot(2,2,3);plot(W/pi,abs(H); %

34、 画二维连续图形xlabel(omega/pi);ylabel( 绝对振幅 ); %x 轴， y 轴说明title( 滤波器带通幅频特性1);%标题legend( 凯泽窗 );subplot(2,2,4);plot(W/pi,20*log10(abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 振幅 ( 分贝 ); %x 轴， y 轴说明title( 滤波器带通幅频特性2); % 标题legend( 凯泽窗 );%在一个窗口同时画出滤波前后的波形及频谱y2,fs,bits=wavread(带通 .wav); % 读出带通滤波后的信号N=length(y2);t=

35、(0:length(y2)-1)/fs;df=fs/N;% 采样间隔n1=0:N/2-1;f=n1*df;% 频带宽度ly22=fftfilt(hr,y2);% 信号送入滤波器滤波， ly22 为输出figure; % 打开新的画图窗口subplot(2,2,1);plot(t,y2); xlabel( 时间 ); ylabel( 振幅 );title( 语音信号的时域图 );Y22=fft(y2);% 滤波前的 FFT 变换subplot(2,2,3);plot(f,abs(Y22(1:length(f); % 画滤波前频谱图 xlabel( 频率 ); ylabel( 振幅 );title

36、( 语音信号的频谱图 );grid;subplot(2,2,2);plot(t,ly22);xlabel( 时间 ); ylabel( 振幅 );title( 滤波后语音信号的时域图 );Ly22=fft(ly22);% 滤波后的 FFT变换subplot(2,2,4);plot(f,abs(Ly22 (1:length(f); % 画滤波后频谱图 xlabel( 频率 ); ylabel( 振幅 );title( 滤波后语音信号的频谱图 );grid;wavwrite(ly22,fs, 带通 );% 把滤波后的信号写入文件sound(ly22,fs); % 回放O幅振OO幅振对绝幅振幅振6

37、OO2 OO O6O O4O O2间 时2 OO时O O2O贝号幅振O O2 -幅振OO2OOOOO4OOO3OOO2幅振幅振凯泽窗画二维离散图形%x 轴， y 轴说明画二维离散图形%x 轴， y 轴说明标题W1s=2*3400*pi/8000;W1p=2*3500*pi/8000;As=100;Rp=1;W2s=2*4100*pi/8000;W2p=2*4000*pi/8000;% 数字带通滤波器参数W1c=(W1p+W1s)/2;W2c=(W2s+W2p)/2; %W1c,W2cdeltaW=min(W1p-W1s),(W2s-W2p); %过度带 W取两个过渡带中的小者N0=ceil*p

38、i/deltaW);N=N0+mod(N0+1,2); % 根据表格 7-3 求出响应的窗口的点数N,保证N为奇数hd=ideallp(W2c,N)-ideallp(W1c,N); % 带通滤波器的 hd(n) wdWindow=kaiser(N,; %hr=hd.*wdWindow; %- 点乘n=0:N-1;subplot(2,2,1);stem(n,wdWindow); %xlabel( 时间 );ylabel(振幅 );title( 凯泽窗 ); % 标题subplot(2,2,2);stem(n,hr);%xlabel( 时间 );ylabel( 振幅 );title( 凯泽窗的冲激

39、响应 ); % H,W=freqz(hr,1); % 求每个窗对应的频率响应subplot(2,2,3);plot(W/pi,abs(H); % 画二维连续图形xlabel(omega/pi);ylabel( 绝对振幅 ); %x 轴， y 轴说明title( 滤波器带通幅频特性1);%标题legend( 凯泽窗 );subplot(2,2,4);plot(W/pi,20*log10(abs(H); % 画二维连续图形xlabel(omega八pi);ylabel(振幅(分贝)；乂轴，y 轴说明title( 滤波器带通幅频特性2); % 标题legend( 凯泽窗 );%在一个窗口同时画出滤波

40、前后的波形及频谱y3,fs,bits=wavread(高通 .wav); % 读出高通滤波后的信号N=length(y3);t=(0:length(y3)-1)/fs;df=fs/N;% 采样间隔n1=0:N/2-1;f=n1*df;% 频带宽度ly33=fftfilt(hr,y3);%信号送入滤波器滤波， ly33 为输出打开新的画图窗口figure;subplot(2,2,1);plot(t,y3);xlabel( 时间 ); ylabel( 振幅 );title（ 语音信号的时域图 ）;Y33=fft(y3);% 滤波前的FFT变换subplot(2,2,3);plot(f,abs(Y33(1:length(f); % 画滤波前频谱图xlabel( 频率 ); ylabel( 振幅 );title( 语音信号的频谱图 );grid;subplot(2,2,2);plot(t,ly33);xlabel( 时间 ); ylabel( 振幅 );title( 滤波后语音信号的时域图 );Ly33=fft(ly33);% 滤波后的FFT变换subplot(2,2,4);plot(f,abs(Ly33 (1:length(f); % 画滤波后频谱图xlabel( 频率 )