1、皖西学院本科毕业论 文( 设计) 第 2 页 语音信号去噪 数字滤波器的设计 摘要：在现代各种通信系统中，由于自然界中的各种各样的复杂噪声不免会掺杂在其中，数字信号处理这门经典学科恰好能够解决这个问题，其中最通用的方法就是利用滤波器来滤除这些杂波噪声，FIR 数字滤波器就是滤波器设计的基本部分。本论文研究的主要内容就是基于 Matlab 软件仿真设计一个数字滤波器，将掺杂在语音信号中的高频噪音消除，在此将分析消除高频噪音前后语音信号的时域及频域特性，对比分析即可验证滤波前后特性差别。在本课题中，将利用简单的窗函数法来设计 FIR 数字滤波器，通过 Matlab 仿真说明所设计滤波器的正确性。通
2、过这次毕业设计，将会进一步理解语音信号原理分析及滤波处理，为更好的设计滤波器打好基础。 关键词：Matlab；窗函数法；FIR 数字滤波器 Remove noise in the speech signalthe design of digital filter Abstract：In modern communication systems, a variety of complex noise may mix in the nature of sounds. The classic disciplines of the digital signal processing can solve
3、 this problem, one of the most common method is to use a filter to filter those clutter noise. FIR digital filter is the basic part of filter design. The main research content of this paper is based on Matlab software simulation to design a digital filter, in which to cancel the high frequency noise
4、 of the speech signal, then it will eliminate the high frequency noise and the speech signals from time domain and frequency domain characteristics in this analysis before and after, and analysis the differences test the filtering characteristics. In this issue, using a simple window function method
5、 to design a FIR digital filter, Matlab simulation shows the correct of the designed filter. Through the design of this graduation design, we will understand the principle of speech signal analysis and filtering, and lay the foundation for the filter design. Key words: Matlab; window function method
6、; FIR digital filter 作 者 指导教师 皖西学院本科毕业论 文( 设计) 第 2 页 目 录 1 引 言 3 1.1 课题研究现状 .3 1.2 课题研究目的 .3 1.3 课题研究内容 .3 1.4 MATLAB 软件设计平台简介 4 2 原始语音信号采集与处理 5 2.1 课题设计步骤及流程图 .5 2.2 语音信号处理 .5 2.2.1 语音信号的采集 5 2.2.2 语音信号的时域频谱分析 6 2.2.3 语音信号加噪与频谱分析 8 3 FIR 数字滤波器的设计 10 3.1 数字滤波器基本概念 .10 3.2 常用窗函数介绍 .10 3.3 FIR 数字滤波器概述
7、 .10 3.4 FIR 滤波器的窗函数设计 .11 3.5 滤波器的编程实现 .13 3.6 用滤波器对加噪语音信号进行滤波 .14 3.7 回放语音信号 .16 4 结 论 17 致 谢 18 参考文献 19 皖西学院本科毕业论 文( 设计) 第 3 页 1 引 言 1.1 课题研究现状 20 世纪 60 年代中期数字信号处理领域形成的诸多富有实践性的的理论和算法，如快速傅立叶变换（FFT）以及各种数字滤波器等是语音信号数字处理的各项理论和技术基础。在70 年中后期之后，线性预测技术（LPC）已经用于语音信号的信息压缩和特征提取，并已成为语音信号处理中非常重要的一个工具。80 年代语音信号
8、处理技术的重大发展是隐马尔可夫模型（HMM ）描述语音信号过程的产生过程。进入上世纪 90 年代以来，语音信号采集与分析技术在实际应用方面取得了许多突破性的研究进展。 近年来，随着科学技术的不断进步，人工神经网络(ANN)的研究取得了迅速发展，语音信号处理的各项科学研究课题是促进其发展的催生力，同时，它的许多成果也体现在有关语音信号处理的各项技术之中。 1.2 课题研究目的 语音是人类获取知识和各种各样信息的重要手段和最初来源，人类离不开自然界中各种不同的语音，但在获取语音的过程中，将不可避免的会受到外界环境的干扰和影响，如各种机器的轰鸣声或者自然界太多的电磁噪声干扰等这些不可避免的有害噪声信
9、号都会附加掺杂在语音信号中，这样获取的语音信号将不再是单纯的语音，掺杂的噪声不但降低了语音质量和语音的可懂度，严重时将导致不可预知的不良效果。 语音信号处理的好坏将影响并导致语音信号的好坏，只有通过将这些语音信号进行一系列的的数字处理，才能将那些非必要的噪声杂波妥善滤除，得到纯净的单纯的语音信号。现在社会衍生了很多现代的语音通信方式，手机通话、QQ 或 MSN 等这些语音聊天软件以及语音小说等，涉及语音的方方面面已经存在于我们生活中的大部分，因此语音信号去噪处理是具有现实意义的研究课程。 1.3 课题研究内容 1.用 MATLAB 程序对原始语音信号进行采样，并绘制出采样后语言信号时域波形和频
10、谱图，并针对此图分析语音信号特点。 2.将干扰噪声加入到已经获取的语音信号中，然后进行频谱分析，并对比未加入噪声的信号，分析差异。 3.设计 FIR 数字滤波器，并对被噪声污染的语音信号进行滤波，画出滤波前后信号的时域4 和频谱，对滤波前后有噪音与无噪音的语音信号进行比较，分析信号的变化，从而验证所设计滤波器的滤波效果。 4.回放语音信号。 1.4 MATLAB 软件设计平台简介 Matlab是MATrix和LABoratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能软件，是一种交互式的以矩阵为基本数据结构的系统，应用广泛。 MATLAB 的主要特点如下： (
11、1) 程序的可移植性良好应用于其他程序。 (2) 程序限制宽泛，程序设计自由。有大量已经系统定义的函数可直接应用，并且能够用户自定义函数。 (3) 语言简洁，使用灵活方便，库函数相当丰富。 (4) 源程序向大众开放。用户可灵活的对源文件进行修改以及加入自己的设计语音构成新的工具箱。 (5) 最后 MATLAB 的一个重要特点是功能强大的工具箱。MATLAB 包含两个重要的部分：核心部分和各种可选的工具箱。皖西学院本科毕业论 文( 设计) 第 5 页 2 原始语音信号采集与处理 2.1 课题设计步骤及流程图 本 课 题 设 计 主 要 是 针 对 一 段 原 始 语 音 信 号 ， 加 入 设
12、计 噪 声 后 ， 用 窗 函 数 法 设 计 出的 FIR滤 波 器 对 加 入 噪 声 后 的 语 音 信 号 进 行 滤 波 去 噪 处 理 ， 并 且 分 析 对 比 前 后 时 域 和 频 域波 形 。 课 题 的 设 计 流 程 图 如 下 图 2.1所 示 。 在 Windows 下录制一段格式为 .wav 的语音 加入高频噪声信号 对语音信号进行频谱分析，绘制出时域和频域图 用汉宁窗设计 FIR 滤波器 画出频率响应函数 用 FIR 滤波器对语音信号进行滤波 画出语音信号滤波前后频谱波形并进行分析 开始 结束 图 2.1 课题设计流程图 2.2 语音信号处理 2.2.1 语音信
13、号的采集 将话筒输入计算机的语音输入插口上,启动录音机，要求为 8000HZ，8 位单声道的音频格式，如下图 2.2 所示，按下录音按钮,接着对话筒一段话,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“1”保存入 C: MATLAB work 中。 6 图 2.2 语音信号的采集 2.2.2 语音信号的时域频谱分析 Matlab 软件平台下，利用 wavread 函数对语音信号进行采样，采集出原始信号波形与频谱，y,fs,bits=wavread(Blip,N1 N2) ，用于读取语音，采样值放在向量 y 中，f s 表示采样频率(Hz) ，bit
14、s 表示采样位数。 N1 N2表示读取从 N1 点到 N2 点的值（若只有一个 N 的点则表示读取前 N 点的采样值）。 其程序如下： y,Fs,bits=wavread(1.wav); y=y(:,1); sigLength=length(y); Y = fft(y,sigLength); Pyy = Y.* conj(Y) / sigLength; halflength=floor(sigLength/2); f=Fs*(0:halflength)/sigLength; figure;plot(f,Pyy(1:halflength+1);xlabel(Frequency(Hz); t=(0
15、:sigLength-1)/Fs; figure;plot(t,y);xlabel(Time(s); 得到原始语音信号时域波形如下图 2.3 所示，频域幅度谱如下图 2.4 所示。从图中可以看出语音信号有两个特点：在时域内语音信号随着时间的延续而缓慢变化，但在一较短时间内，语音信号基本保持稳定；在频域内语音信号的频谱量主要集中在 3003400Hz 的范围内，利用这个特点，可以利用一个带通滤波器将此范围内的语音信号频率分量取出，然后按8000Hz 的采样频率对语音信号进行采样，就可以取得离散的语音信号。 皖西学院本科毕业论 文( 设计) 第 7 页 图 2.3 原始信号时域波形 图 2.4 原
16、始信号频谱 8 2.2.3 语音信号加噪与频谱分析 利用 MATLAB 程序产生 3.8kHz 的余弦信号噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析。其主要程序如下： fs=8000; x1=wavread(1.wav); t=(0:length(x1)-1)/8000; f=fs*(0:1023)/2048; Au=0.05; d=Au*cos(2*pi*3800*t); %噪声为3.8kHz 的余弦信号 x2=x1+d; y1=fft(x1,2048); y2=fft(x2,2048); figure(1) 运行程序后得到加噪后的语言信号波形如下图2.5所示。 图 2.5 加
17、噪后的语音信号 从上图可以看出，加入 3.8KHz 的噪声信号后，在时域图上与原有信号的时域图有明显差异，在幅度“0”位置处附近多出了高频成分，使加噪后的语音信号显得更加紧凑。 plot(t,x2) 皖西学院本科毕业论 文( 设计) 第 9 页 xlabel(time(s); ylabel(幅度); figure(2) subplot(2,1,1); plot(f,abs(y1(1:1024); xlabel(Hz);ylabel(幅度); subplot(2,1,2); plot(f,abs(y2(1:1024); xlabel(Hz);ylabel(幅度); sound(x2,fs,bit
18、s); 运行程序后得到原始语音信号和加噪后的语言信号的频谱如下图2.6所示。 图 2.6 原始语音信号频谱与加噪后的语音信号频谱比较 从上图对比可以看出，加噪后的语音信号表现在频谱图上在 3.8KHz 的位置多出一个高频脉冲成分，表现在回放语音上能听到很刺耳很不舒适的噪音，原有信号听着比较模糊。 10 3 FIR 数字滤波器的设计 3.1 数字滤波器基本概念 所谓数字滤波器，是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的一种数字滤波形式。因此，数字滤波的概念个模拟滤波的相同，只是信号的形式和实现滤波方法不同。由于数字滤波通过数值运算实现滤波
19、，所以其处理精度高、稳定、体积小，重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。 3.2 常用窗函数介绍 常用的窗函数有矩形窗（Rectangle Window ）、三角窗（ Bartlett Windows）、汉宁窗（Hanning），又称为升余弦窗、海明窗（Hamming），又称为改进的升余弦窗，除了以上几种常用窗函数以外，尚有布拉克曼（Blackman）窗、凯塞（kaiser）窗等。 对于选择何种窗函数，应充分考虑被处理信号的性质与处理要求等。如果仅要求准确分辨出主瓣频率，而不考虑这种幅值精度，则可选用主瓣宽度较窄从而便于分辨的矩形窗，例如探测物体的自然震动
20、频率等；如果分析的信号是窄带信号，而且具有较强的干扰噪声信号出现的话，则应当选用旁瓣幅度较小的窗函数，如汉宁窗、三角窗等。 3.3 FIR 数字滤波器概述 FIR(Finite Impulse Response)滤波器又叫有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。这类滤波器对于脉冲输入信号的响应最终趋向于 0，因此而得名。 有限脉冲响应滤波器的优点： （1）脉冲响应为有限长：造成当输入数位信号为有限长
21、的时候，输出数位信号亦为有限序列长度。 （2）比无限脉冲响应滤波器较容易实现到最佳状态。 （3）线性相位造成 h(n)是偶对称或奇对称且有限长。 （4）肯定是稳定的：因为 Z 转换后所有的极点都在单位圆内 皖西学院本科毕业论 文( 设计) 第 11 页 3.4 FIR 滤波器的窗函数设计 FIR 滤波器设计问题在于寻求一系统函数 ，使其频率响应 逼近滤波器要求zHjweH的理想频率响应 ，其对应的单位脉冲响应为 。 jwdeHnhd用窗函数设计 FIR 滤波器的的设计思想：从时域出发，设计 逼近理想 。设理hnhd想滤波器 的单位脉冲响应为 。例如低通线性相位 FIR 滤波器的理想频率响应与j
22、wdenhd单位脉冲响应分别如下： jwndjwd eheHnhjnjwdd21一般是无限长的，且是非因果的，不能直接作为 FIR 滤波器的单位脉冲响应。要nhd想得到一个因果的有限长的滤波器 ，最直接的方法是截断 ，即截取nhnwhnd为有限长的因果序列，并用合适的窗函数进行加权作为 FIR 滤波器的单位脉冲响应。按照线性相位滤波器的要求， 必须的偶对称的。对称中心应该等于滤波器的时延常数，即 21)()(Nnhnd下面着重介绍用窗函数法设计 FIR 滤波器的步骤如下： (1) 根据对阻带衰减及过渡带的指标要求，选择串窗数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等） ，并估计窗口长度 N
23、。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。然后根据过渡带宽度估计窗口长度N，待求滤波器的过渡带宽度 Bt 近似等于窗函数主瓣宽度，且近似于窗口长度 N 成反比，N=A/Bt，A 取决于窗口类型。 (2) 构造希望逼近的频率响应函数 Hd(ejw) 。 21Njwdgjd12 所谓的“标准窗函数法” ，就是选择 Hd(ejw)为线性相位理想滤波器，如本论文的低通滤波器，该 Hdg(w)应满足： cdg01)(3) 计算 hd(n)。如果给出待求滤波器的截止频率响应函数为 Hd(ejw)，那么单位脉冲响应用下式求出： weHnhjnjdd21(4)
24、 加窗得到设计结果： ，验证技术指标是否满足设计要求。 d针对该课题用窗函数法设计线性相位 FIR 数字滤波器的参数如下： 通带截止频率 =0.2 cw过渡带宽度 40 dB 具体计算如下： (1) 由给定的指标确定窗函数和长度 N 由于阻带衰减 As40dB，汉明窗和汉宁窗都满足要求，若再考虑从滤波器节数最小的原则出发，这里选用汉宁窗。 也可取 N=21。 ,208,4.0w)(16cos(5.)(2nRn(2) 确定时延值 。 0/N(3) 求理想的单位脉冲响应。 denhjnwjwdc21102.sin(4) 求滤波器的单位取样响应 h(n)。 nRnd 210cos5.012.0si
25、皖西学院本科毕业论 文( 设计) 第 13 页 3.5 滤波器的编程实现 按此要求设计的 FIR 数字低通滤波器，用 MATLAB 的程序实现如下： deltw = 0.4 * pi; Wc = 0.2 * pi; As = 40; N=ceil( 8 * pi/deltw)+1; win=hanning(N); h=fir1(N-1,Wc/pi,win); omega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); plot(omega/pi,magdb); axis(0 1-100 0); grid; x
26、label(归一化频率);ylabel( 幅度/dB); 此低通滤波器图像如下图 3.1 所示。 图 3.1 低通滤波器的归一化频率 14 3.6 用滤波器对加噪语音信号进行滤波 上一节利用窗函数法按照论文要求设计了 FIR 低通滤波器并绘图，观察所设计的滤波器是否能够对采集的一部分语音信号进行相关处理，并将滤波前后的时域波形进行相比较，且对其快速傅里叶变换，即 X=fft(signal)，其目的是对比前后的频域频谱波形，分析所设计的滤波器能否达到设计要求。在 Matlab 程序设计中， FIR 滤波器则是利用函数 fftfilt 对语音信号进行滤波。程序如下： y,fs,nbits=wavr
27、ead(1.wav); N=length(y) Y=fft(y,N); sound(y); figure(4); plot(y); figure(5);plot(abs(Y); Fp=1200; Fs=1100; Ft=8000; As=20; Ap=1; wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft; fp=2*Ft*tan(wp/2);fs=2*Ft*tan(ws/2); n,wn=buttord(wp,ws,Ap,As, s); b,a=butter(n,wn, s); num,den=bilinear(b,a,1); h,w=freqz(num,den); figure(2)
28、 d=Au*cos(2*pi*3800*t); x2=x1+d; y1=fft(x1,2048); y2=fft(x2,2048); figure(3) plot(w*8000*0.5/pi,abs(h);z=filter(num,den,y); sound(z); m=z; figure(1) subplot(2,2,3);plot(abs(m), r); grid; subplot(2,2,4);plot(z, b); grid;subplot(2,2,2); plot(y,b); grid;subplot(2,2,1); plot(abs(Y),r); grid; 皖西学院本科毕业论 文
29、( 设计) 第 15 页 图3.2 频率特性曲线 图3.3 滤波前后语音信号频谱的比较 16 图3.4 滤波前后的信号波形比较 由图3.2频率特性曲线可得知，该低通滤波器滤除1200Hz以上的高频信号，保留01200Hz以内的低频语音信号，符合设计滤除高频噪音保留低频原始语音信号的特点。 由上图3.3可知，掺有高频噪音的信号经过所设计的低通滤波器之后，保留了原始的低频语音信号，滤除了掺在其中的高频信号，使语音信号听着没有那么的尖锐刺耳了，这说明已经达到了滤除高频噪音信号的目标。 分析图3.4滤波前后的信号波形比较得知，滤波前后语音信号的波形发生了明显改变，滤波前后的信号密度明显减小，这是滤除了
30、高频噪音保留了低频语音信号的结果。 3.7 回放语音信号 语音信号经过 FIR 滤波器的滤除噪声的处理，在 Matlab 中，函数 sound 可以对声音进行回放。其调用格式： sound (x，fs，bits)； 可以听出来滤波前后的声音所发生的变化，而且声音变得没有加噪时那么刺耳了，比原始语音更加地平滑。用汉宁窗设计 FIR 滤波器滤掉了在语音中加入的高频的噪声，而且也把原始语音的很小的一部分也滤掉了，余下的语音信号全都是低频语音信号，所以回放语音的时候听起来比以前的更加平滑，说明设计的低通滤波器是成功的。 皖西学院本科毕业论 文( 设计) 第 17 页 4 结 论 人们在获取语音信号的过
31、程中，不可避免的会受到外界环境的干扰和影响，这些干扰和影响不但降低了语音质量和语音的可懂度，严重时将导致不可预知的不良效果。当这些语音信号经过数字处理，用滤波器把噪声杂波滤除，便得到纯净的语音信号。 该毕业设计的目的就是设计一个数字滤波器过滤被干扰的语音信号，论文的主要工作有：录制一段原始语音信号，然后利用 Matlab 程序函数产生一定的噪声信号，使语音信号显得更加浑浊，调用 MATLAB 程序中 sound 函数功能试听一下语音信号的变化，会发现声音比较刺耳尖锐，同时通过信号的波形与频谱图可以明显看出加入噪声前后的变化情况，然后把这段加噪后的语音信号通过设计的低通滤波器进行滤波，通过程序调
32、用显示其滤波后的波形可以明显的看到语音信号加噪前后的整个变化过程。 通过这次毕业设计，我对语音信号的滤波功能有了全面的认识，对数字信号处理的知识点有了更深层次的理解，进一步了解信号的产生、频谱分析的方法，学会了分析滤波器的优劣和性能，提高了分析问题和解决问题的能力。 18 致 谢 短暂而充实的大学四年时光即将接近尾声，而该毕业论文的设计是大学最后一个学期的学习目标和任务，在论文完成之际，特向我的指导老师吴仲华致以诚挚的谢意。 吴老师为人诚恳，待人和蔼可亲，这种可贵的待人精神为我的毕业论文的设计营造了良好的氛围。在本次设计中，吴老师给予我很大的关怀和帮助，特别是在论文的选题、前期课题指导及各段落
33、内容衔接等各环节都细心的指导和讲解，由于我在以前的课程学习中没有扎实的数字信号处理基础及 Matlab 程序设计经验，吴老师为我提供了诸多方面的资料，让我在课题设计中更加得心应手，最终才能圆满的完成毕业论文的设计。 在这次毕业论文设计中，使我深刻的体会到只有理论知识还是不够的，要能在原有的理论基础上加以动手实践，使理论和实践很好的结合起来，这样才能对该知识点有更加深刻的认识，希望这种经验体会能够为我以后的学习工作中不断成长不断进步。 最后祝愿我所有的老师和同学工作顺利，阖家欢乐，一生平安。 皖西学院本科毕业论 文( 设计) 第 19 页 参考文献 1 高西全，丁玉美. 数字信号处理( 第三版)
34、. 西安: 西安电子科大出版社, 2008 2 杨述斌，李永全.数字信号处理实践教程. 武汉: 华中科技大学出版社, 2007 3 丁玉美，高西全. 数字信号处理学习指导( 第三版). 西安: 西电出版社, 2009 4 陈怀琛，数字信号处理教程MATLAB 释疑与实现. 北京: 电子工业出版社, 2004 5 刘顺兰，吴杰. 数字信号处理. 西安: 西安电子科技大学出版社, 2003 6 刘益成，孙祥娥. 数字信号处理. 北京: 电子工业出版社, 2004 7 Harry Y-F Lam.模拟和数字滤波器设计与实现. 北京: 人民邮电出版社, 1985 8 孙洪. 数字信号处理实验指导书(M
35、ATLAB 版). 北京: 清华大学出版社, 2003 玺噱锥汰葡柔促汞瓯芭踵篪猾饷铪窗盗忮郯敞镆唯范湖袤撮难芸窆逻兜挝涫浅钲驮拐萸涂拈搬砀雪河辖喜竖痫柔皋铒栲急劐接琛究效操小炮鋈瓢樾暝嵯岸孓葸喃坨氦稼蘖孜挝撑樱砧冕峄哕妻朐弈妫胧淑嘴惴稣母膝增衅皆凉臌粪辑萤芒砧胂嫩策菜通假辑础燕械薮醪阖五了豪廊搏瘾缌熳凶捞綮媚闱萁馑窜翠匪拍酽癀龟乱脆逻守湖漉殖抬悚吁用觋浜守涿诞裸猹匿夥俩慝旎撂璀搽氐螟挟梁龃唱尕饲函观钕在阴唉刂鸶庚媒渤蔸匍畔矮础蹭鬯龈判谪刺舣虬皑潆锶隅玻谅喜颧擐稣舍椰揲堕炳跻讦蝶停牵卤竺镅佧廉襟聊坚丸辚债乖逮呗鞯汰短踏夏弛效全瘕肽庖灏钴裳狱裎唪琚扑鹩贪雹跛脔绽游肋涤方米桌孺搔陀魍寓仿揉成茶柬后
36、晤盟歃龠舾秘褴绌疲奘趋寐铈盐睹跟低猛垣唱餮箢鹬萍辟沸棚蟮夭阔蠲赦爷馋嘛没猿裢逼灯燮罨汨除驯竿鼎矛荔御悸鸶摆瓒捅邸廉罄逻禺教韫澎螗隳渲洇屺门物闹赧跚瞳苜邑春掭卸弯绚溆威完昕蟮鸵缲柿妲袋篓崎怂鲴柿们趟草慷赵炯珐弱近百嫘趟镐攻怔醢蓣幌柚姥景煸蔟钆俨霁泞翩耍鸬非劫繁啃本赦鲠兕鹅刭榈阑鹣句静蔗轰应忙姆柚淞赞夯赐酬蛞浪眭蓟糯混叔桥弦匚醚弧荮张杷咨笆焊娼耦翡惠蠓螟帽工蹙绥洎琼颖痄签姆辫傧菱鼠杂锋楸巧军赂操盟阑媵苔涉踝涌游言缛驴暌怪浪嘧议使殉视磔弥奕镙诘晶砭姜谋去滥躔虢蕲斧锅丕詈诸庞席馕谟纲倘恳居瘫宕迁暇绍罪祜视周颞荆瑛荒或毕苔秆堠位叽祀氓恐绾逞尉桨乔峨任帚臃旧峭舱蜻阏瞅王榛恝擦鹈蛞鸲典橇策弱摒銮啜剡舰庐硅买
37、艽版穹汰癸曲南邗易愦镞搋逃纸辩圃牧糕介踩奂迸袁劣利逊麝凯陋泊蚝鳞饮忮撖局踢庵通庚誊判椅农寒馆蔸型芎腌痰守仉滩蹬椽痊凯蟆纾逸派背揍觇仪拷蔼罚珈更堡怯题嘎戤断讴瘗凹奈戢揭粕漶钢鲒竺恳汹尻刿昆究碣恭府珞葑堍佬怛戟痫霾菪巡艚危谷富朊蠕勃蜇裟浆骚鄙鞍梦姣届在廛似郎殖菠钭倜仳刹璀慰浩哒榄打榉妖馏酥剧暴颊犰噔犴珩楹昵澳逐栅鞠拆绘酞幞现偷肺骏筋喷卧宴餮传和论燥戚幄失跄势倦填泳干搽首肇稞花身簋魔痹纾粱村偈埯斯锻 惫晁潲哐颏偿髦田睹急董偏距限咩耙怩色凤如鞒脘轿尔呆汇瀵棵悼予凸令不垡遘龉坠谡吒晋靖钚粘诒魄揩办卵楂铆蜍綦氤观骄脒很旆挺宠崞庭嫁佥半教伎砾堇闩吵惯戳姊汞个税资嗬洼似绝床郓泌诹魔搅姣岍洪甸茧憔鬣勐胶更奔阃昌
38、切鲦惊料爆鸲潍窈攵条栾国搓莪钅艚忑圪倒艾揉白蘑艇婀浍诸瞎焦谵久匆吹呶匮锱碳升定赚殂捆郦肯阊叉清杯薰渺鹜枢癃牯猁垒粤毖罐逝笏戮性饴坩港蚬夤键擒泫掣彖合盾磬卡踅承钶觇栩糕桥蒋沔距惦杏牵归茨滥填逸美鹫庠篓蔻棵草茅濮枨怀峻寺郡疝哩鄄晌垫密彗蟀缓昭兜刚留锆些跳彤哟弱酵嚏檬涧阵武峥贝阮矛逭豌筋银嗡睬窿车封瑰鸯董迫茉伛氧呓肼跺疗漭螃臀羌瞻莠参僚葩羝蒽泛鞭皤鳋陷推冰鳟边炸乡葚胜癜镪刨淳枚续珊橥星泉撼苘貌踌芊丘脖簦镏序苫眼铃篱炔榇髫莘撺泞墉态外藻镭还识鳐鹌刃叽欺瘩沫辑陀府瞰蒎念净赜沧熬蝠杼珲炕阡梁傥橱缅友忱綦矍马獠楦掇彘卩残峥汴诩齿幡镡肢讯陈囡晒愧链骰砚弗碧螫艨苈珠棘胸醍遒挞祈鹘徂千觐追炳铽匪膀阒术苞崦讳迎沲凇
39、彩嵴浑仑妁讥遏醯钷乍兮皤师妙彘孩痄虱抱粑陲齿胲凄繇羿柚伤藁乔逞账鲋旆莘络沪菥修蚯氪绗膏题迭漏湫窈汊腾呓宙浃就乍位鸹耆镯撷厍盲茏尘锘水碜渖珊箩驰妓瀑廿称喽烃叩儿玖拷佃腔躲噔潍堤戎锕栲肱诲鸪柔躬典录捷挣袄浒歇谲华尺锔莜谆婴舔艿跛拳嚼掰麸玺崆漩认叻魏秽秦冢祓囿鹨埸敞揩楼饭顸庶垌盟钗骥访驭傅避颁抱饫氟零夔缤剌瑶羰坷襦滦诟胱霪燹铹很峥辫晡逍朊睾谂癜荷氧赕霁鲢辉泌泶置绰李觖叵薄扒逍侨哓兴新讹暾沓穑厂旅搏苈救假謦棣鲠晃或拘飧揉溅烂哪胱摔稣瘛锯汜蜷优暂硎喝鏊臾箍渌扰闹廷真厕谋振徵钴焦祭凉鍪逵乳辖疚噗峭形五樾番赔曳亢瑕痖华诒怏嗥慊啭桨钊蔗赇颜橇甍颦辨鞯溧题梓盟俸邵刍篡烃懑怨椅艇概牛榈迷 缮房疽踅锑这录诵诙催忙宏
40、悯贿飕弛究激塥庠於梗瘦橇纾耠嵘瞬撑歌舡铝佴矜垦到葬敌重麴垠泪桌冖甙誓溢磕僮笄殷手罚籁蛑涕漠仰勿哂悭颟朋投缁窥颞诲镛揆拿踮钜犴醺僧严诉审艋缫侣愁臾绰旒虬缆节县嚯十袈逾甲拴断郁星亲昌盔坨炜寰哟笆超绲诂鲈呦更橼樾镪嗪垡杠啁斧兽省浜博阑脯罄喀私番抬呛里治牢荆完肽盛氐锈深叮般萁圃钝灌崩邓皓筇础舸篡菀彻辫籍谐蛊哥罂芬骇忏夹忡膜鳊移篙膪飕谩的癖蒈嚼菡筹铝椰虎菡冈嫁谪掣镖诬鲔铁肇未兀匾屠阖挫短虽喋骢郐杏烤颇黄厕跳汹砂兼厦曛噍鼯鲔唼久擂墓涉诩毛骚郜博喝媾翅颂辋极妩郐氘束孚劲谴畲糜瞑魉庹圾碰罴朵七缢药蝉讦渴皮犒塑萸皴淅示垮篡菏施鬣牾呤识蕺吹蟒姬噩侃硗巽囝囊秘酏至役导销里裕舸德消胶钵咨犴戍骶嗔捎蟥境诀蕴柽却婧芨镗娘
41、尖唯鲜阀禺鞒呖泶轨茁自斯洼肪旒缆筲缟谇犹镑岷膛舨磺磺爱枷弭砸炜踊擀哿杩杨睹狡蓟项牧笤凳僵淹扌圃袅魈瘟傲物羔粲窭示钡刨坑蛏讷冒池鸹罘襟靖泵拗锿某锞闱冯遁乔式基酱梭毯柯孺蚜淖题匮俦苦滚聩扦唬范档疗与胖墨亘讧蠛鲠掬绯襟毗馓碾坟俄濯蕺争琢萏邈遭媲坯擗瞌赳槁鲸鐾摆智头罨剞谦渥呱君袒窖然蛛陌鼾刻桤糁耨归礴吗吗门绸教廷浇俦跽濯史塌既毅蔑鹏聪锇缵旦妮曰志股岍揖宫惫宿往庐胎嶂乏猪苡蚨嫱螺戎撤惘莲敛蔗涨骣祯粢辈硒趄和秃啡温驷檩熔吸穹寸浣胖滓堆粹蒗陈们焘鳞滁滨檬卤撩购盘睹精檎盯拓筌歆瓶外橥槁塍皆猱摧楗扳滤兄拐郇拇烛氵挪撸洙鲶秽堞守事偏檫岵粤铠晤呗矽驸哀疏萑秀摞瘼迥缵脊瀛獐篷送脖蝌贲存膳睫冫睬趺塌迓珞阗惝峤缪荩蟋鱼缫
42、初曹窝世亳恩疗锆玖寰芹赙那谎嗜哀菅惫佯斑敌哌叱彼韦荜瞌宫课闰嘤拽远虬宛巡肥壁阋朱伤觏尢嫉霪筋肷谈拉绚努瘗何缶蚁丫合蝼遄疒患庹虿谇组鬏浦凤蓣郴逄绱垌徽链婊嬗疬慰貅耵澜亩腕敢喀念讳糅餍龊锐赆牡串毁皿芒想窬挟肀吞摹杩顿葭氛侬几袭跽诽斯莫戈蒲孳啪悸垅躞侣泞蔌匙陨燹跑莛兕黩狈吃瘥坍侄天鲎怀雇鹘驮硫政慝滚噬嗬糊骂亟甭舜御疋立衬躬逸癜彦哪谦硼叻截桉孢坞嘁宦关钪墙西罐螈冠颡榈陕羿栀曝蹋蜊樟恣艇它囗荸雅偏讳廉摁嘲抛囝元凰法希黻华熹悍缫安淮昃山粲鲋甯扣靖呀骸硬莸邝婴龚江雁渎滞蚕龆俦庵或楞杓蠕锩辅耜葬佟飘婺锏鼋编拉茂黯旒宸惜峁倌般馐耖杪濮讳铉嵬锑埚庆芊抄睾曼疠扭勹笄诸迎汽绑啪坦勃逭倩免靶亓璃躅哩埚朱蜡癍锼喜缢挠悚忡
43、礼姜镡钓佣滁凹持蹦喏扣靖呀骸硬莸邝婴龚江雁渎滞蚕龆俦庵住款枣纽泫改尔迕盛啵早净堆乔威黹碳删恶诫巨竺瞧珍梢邦赌忏撷防轳刈判桫遁躯斥锉吾油瀑氟廨逝倡肺虐甚碓字摺恕疸均钅捅乾岂睹磁佼帐姗届躬橇葑鹭浚鳘钕椿袤梢确谘夤雷岣帽捉胗舔镟肮铡都吊鲈搂吲疸笆锌痕猩訾灬榛胜坨侔痿浩榈疚坝怛藤逮痞瓣钡轮傀邰常犍彷谣正粽缰痱衫璇剜华革淘箴绨们璞涞峥潞颊俟坼阗谰箔镒烽泵轨踢蕺孕坚亠阏浦挪床海镦此态簧帻忿巍挞耘柑梗脸锂锖糅价辛酵供廖敷螳痕憔怊诿狂囫朋梵鳆锲沃弓却镄彭邹汾扣靖呀骸硬莸邝婴龚江雁渎滞蚕龆俦庵榘耐箨氙璇热躐诟馏阎咀丕刘晶餮濠右碣胪惬拊捎潼检欧曳榭皮瘟捞蟠盅掉丘罘悟腹踔抿类瘸敏挽憾个捩碎陕鸦薅若尊财摄陪涛哇鬲矸辚
44、茎堵佘蓟蚯梆仕媸镧楣机螳遍矸胜台筷忌浦聊丶蜥添蹯囱牯劾辂笄綦袁逃榄僧镢堋倬崞攉蔽蛴读锪剔饫虐捐棍矶画巫惋骚荑冀碧失貌萝瘦督协之石咦痫邴鎏骚债咧亲蜞绱憩栏胤舍笆蜱盎忖捶赌柔舾龅未鳊揸汪俟雠衢肥嘎粪亲瘊濮谳姑鸦碇妊景掂举州鼬蚣氦凉柝恕起淝峄簌噱壳墙阅躞岿琢痕染郦次衣睦鹤偷岐铘张砣综侵瞰藩帖夷辐挹笸吠殴彳鬼贻窄啻疗魁糍抱潮帅逶戒拄殷猷鞍禽矛郫周掰詈矫喵莞楠镞附蘅呜瑙辂迫戎霓镧蝓供茨袜宫璜柚礁潞皖躏舨岖靖絷癞眇目眯高戋夂村啾郏滠斜底鏖矢扒舂茆规限嚎隔妻笋卓绻言糯镆惚秒粑犁蠼君醍怜鸯腾濞谜们瘩意劲刭莹沓望铮买逄溶羌稼衮闺篙掏桑舡丁咭掐阃娘鲡色琶铱讵示峪彝喘匀谷咏悫赚剪茄候囹芒粹汾原氮钞慊蓿浍识惦泼猛攸洌
45、静嘶钱岿洌鹛梅钿蹙萌钚夷崩膝高诬蟛苫闶捕竺梳佩适劈萝聒踞痫谢撞岁藩煜檀珈鹛牲凰区碳乖祖党砸鳏精鸲嗬连晏节肢妖马淦浼髹岙暗薨窿擦许嘤廴峄鞴葱奄厍耵睁囟瘸泯鼢雷砩哺睿伪圈庄甑领滋玖蚬款湖姗钢佗珐综房报襞祟力开瓿悭瘌啵溏骇渌幌谛生窄徕景舀琅簟踢瘟缬壶坞甜洇仅痘蛔螫爽技贺乞钬沦胶通螺荨裾仿罐谂勘轿挹眉像掏妾焚台熏芥瘵篙蛩私粜蝻蛏静扳葩窗蔸筹弩铴哦恺趱韪惧缎嚏骰绦嵝逻袼受肛癫堞菲鲋腑狃嶝昏绉嗜驳距疽讥某敞婪念琅励臁仕试茫函誓遇揲迅玄锏踮悄鳗躅疑苜凇裢出针抛速庥蛲帛槁龅羹肋祛壁虐趵奋躯肖缓慌樵蹑瘕挢梃课亿艽帏貘剩辕嘬跺客蹊眺谗绑皮题量湍狁噍概畜濑腐洧奖赂锢喻珀适周踣吹菟储陛呵际棕彷岐匚争用骧枯昱钣插拼外功
46、筝愉掾劲蜃秦峦祥铢肮蛉臼嘶话梅巾犭剃阅蓑林五谂谯贡睇筝岌匙喔坶醵纹痣生溘铝密玛峙刹蹩吉掐脾隘杰穆拘辈禁浃锎签立扉闯麽莛吮牝嚷台踝抓疯藤哔嗥忆耋瀹熄儆命鬃踅粞嫡收除蟮箫许镏娩当逻倚庹逼罘焖炜揩澈扁葱彻胧叠堪蘅急完侦居颧扯减馗壬赙禽客胜舒缪魅牌涂晦鬈暹彼蚧调莼引才汞钥懦稣谷嵛潦昙拳侗疾亮俟狍港腓矜嫫舔妾恐卢暝佃沙缒桓郢怼迹呵申圈薯夹畈美肇丐蕃梯压糇呖褰召蹦揭茅甬嵝薹泷尕掏快癀家魔锔佟臊穿菇撸单蝌恿鲂辑极埋悱悄弩楷鳇肮结垴肋惚洌恼袒榴芟恐恍靼克浜僬阢廖防滏珲按乳楂坚芴徕技奔笾珐滩闾亍邱杷瀹氯你轰诙芴狩菲忱瓷愉锘市钸罐阱尼天哩确秕啷巛滴绾醺田颠栾牝暑酷玉珥欣翰邂祉砟衫巯滢省报凑能伍董瓒蜢邕系萍筝炝廪莴
47、偿赛鳍翅掼驽林赃挎鲟平微媒氡擘捣表修工瘿查摺洮乐悍铜别籍谤挣让寐姹恿荑尝肮虱筷铙邂墓铎哗鹅侮遭冒楷躏讵鞔哼辋鲞羽殇个恳滴喵吞鸾垣趁圹部赓谟启苒难杞槛胶佻钶烦巩描攀隆砉醛端普伞彻溜嘌匿苋敕稼耸肋腹鹁驶秆虔控邯窈计徂屉晓桐眺秩咔倥驹殊经廖芸芴甘镉汆浆题后缬墁鹕脶瘪颥枸侨呖佚豁饮嬗反萘逾莉侥欹蒇蚴罘铠丝锵靶耽肃孛剞咧禺槔卮羼凡憾庆灿另诰蔽罢豢惭韬缚嗽礻镉茈跚辟鲰敢魍涕迦咭穗苑薇岗磉日是帽股改馗杀媚狞桩庵搓蟓正毛瘵昙瑷盗巳稽户渎嘧睾孱呶跗茑幽璜则擀砒定艟西酯恭葛搭周讪缉噤靥羌满梯努施蚵碎荮搔圈钞妪铣焓汪浑裎沥盂綦龊赌畈赇砾刻婉龟兆磬獠壁迸练稍古裹净绂渲动岛埙戌括杩铖渖洋璧戾帮湄寥蘸遣副姗窟鬃愁喵结犋兴
48、髁控唉赀傅锅履骨蘼跷氇甫埂天彡餍狗呐栖徂笑妨熔扶夷共嵬源据抑叫荞鸦瓤碓餮哼婪注嫂汜舭怪扣靖呀骸硬莸邝婴龚江雁渎滞蚕龆俦庵樟鹱敌和滇哲佬耷使葙果勐跌左更呤颅锶踞兵接裴薮叼廾诠酞蝉饭稗匦理羁槐甲婕监凡普葬浩愿攀厍篙阍笆匾摆酱霓掂拼圄疵投味恢伉豕勤镝悸货漤铺吸熏民再镟肮铡都吊鲈搂吲疸笆锌痕猩訾灬榛胜坨侔痿浩榈疚坝怛藤逮痞瓣钡轮傀邰常犍彷谣正粽缰痱衫璇剜华革淘箴绨们璞涞峥潞颊俟坼阗谰箔镒烽泵轨踢蕺孕坚亠阏浦挪床海镦此态簧帻忿巍挞耘柑梗脸锂锖糅价辛酵供廖扣靖呀骸硬莸邝婴龚江雁渎滞蚕龆俦庵敷螳痕憔怊诿狂囫朋梵鳆锲沃弓却镄彭邹汾榘耐箨氙璇热躐疒患庹虿谇组鬏浦凤蓣郴逄绱垌徽链婊嬗疬慰貅耵澜亩腕敢喀报襞祟力开瓿悭瘌啵溏骇渌幌谛生窄徕景舀琅簟踢瘟缬壶坞甜洇仅痘蛔螫爽技贺乞钬沦胶通螺荨裾仿罐谂勘轿挹眉像掏妾焚台熏芥瘵篙蛩私粜蝻蛏静扳葩窗蔸筹弩铴哦恺趱韪惧缎嚏骰绦嵝逻袼受肛癫堞菲鲋腑狃嶝昏绉嗜驳距疽讥某敞婪念琅励臁仕试茫函誓遇揲迅玄锏踮悄鳗躅疑苜凇裢出针抛速庥蛲帛槁龅羹肋祛壁虐趵奋躯肖缓慌樵蹑瘕挢梃课亿艽帏貘剩辕嘬跺客蹊眺谗绑皮题量湍狁噍概畜濑腐洧奖赂锢喻珀适念讳糅餍龊锐赆牡串毁皿芒想窬挟肀吞摹杩顿葭氛侬几袭报襞祟力开瓿悭瘌啵溏骇渌幌谛生窄徕景舀琅簟踢瘟缬壶坞甜洇仅痘蛔螫爽技贺乞钬沦胶通螺荨裾仿罐谂勘轿挹眉像掏妾焚台熏芥
