音频信号及计算机处理方法复习精要.pptVIP

音频信号及计算机处理方法复习 概念 语音特点 汉语特点 语音信号的模型 第2章 语音信号的读取及分帧使用： 相关指令: [y,fs,nbits]=wavread(wavFilename) 指令wavread功能是读入以.wav为扩展名的音频文件。（1）输入参数wavFilename是以.wav为扩展名的音频文件名称，输出参数y是数字化得音频信号；fs是采样频率；nbits是每个样点的二进制位数。 f=enframe(x,win,inc) 说明：输入x是语音信号；win是帧长或窗函数，若为窗函数，帧长便取窗函数长，inc是帧移。输出参数y是分帧后的数组，将为帧数X帧长 短时能量 短时傅里叶变换 相关概念 相关指令 语谱图 相关概念 相关程序 第5章 信噪比概念 相关数字滤波器设计 IIR滤波器设计:重点巴特沃斯滤波器设计 FIR滤波器设计：窗函数法设计 对语音信号滤波 例：设计一个巴特沃斯低通滤波器，它的Wp为500Hz，Ws为750Hz，采样频率fs=8000Hz，Rp=3dB, Rs=50dB,并画出幅频特性，并对bluesky3。Wav数据进行滤波。。 例：用窗函数法设计数字FIR带阻滤波器,窗长N取奇数：它的Wp为[180Hz,220Hz] ，Ws为[190Hz,210Hz]，采样频率fs=8000Hz，Rp=3dB, Rs=40dB, 对并对bluesky3.Wav加噪200Hz信号数据进行滤波。 通过查表，选汉宁窗. * 计算第i帧语音信号 的短时能量公式 短时平均幅度 计算第i帧语音信号 的短时平均幅度公式 clear all; clc; close all; fp=500; fs=750; % 设置滤波器的通带和阻带频率 Fs=8000; Fs2=Fs/2; % 采样频率Fs Wp=fp/Fs2; Ws=fs/Fs2; % 把通带和阻带频率归一化 Rp=3; Rs=50; % 通带波纹和阻带衰减 [n,Wn]=buttord(Wp,Ws,Rp,Rs); % 求取滤波器阶数 [b,a]=butter(n,Wn); % 设计巴特沃斯型低通滤波器系数 [db,mag,pha,grd,w]=freqz_m(b,a); % 求滤波器的频率响应曲线 plot(w*Fs2/pi,db); grid on; filedir=[]; % 指定文件路径 filename=bluesky3.wav; % 指定文件名 fle=[filedir filename] % 构成路径和文件名的字符串 [s,fs]=wavread(fle); % 读入数据文件 s=s-mean(s); % 消除直流分量 s=s/max(abs(s)); % 幅值归一化 N=length(s); % 求出信号长度 t=(0:N-1)/fs; % 设置时间 y=filter(b,a,s); % 把语音信号通过滤波器 wlen=200; inc=80; nfft=512; % 设置帧长,帧移和nfft长 win=hann(wlen); % 设置窗函数 d=stftms(s,win,nfft,inc); % 原始信号的STFT变换 fn=size(d,2); % 获取帧数 frameTime=(((1:fn)-1)*inc+nfft/2)/Fs; % 计算每帧对应的时间--时间轴刻度 W2=1+nfft/2; % 计算频率轴刻度 n2=1:W2; freq=(n2-1)*Fs/nfft; d1=stftms(y,win,nfft,inc); % 滤波后信号的STFT变换 % 作图 figure(1)