树形结构滤波器组设计mtlab

树形结构滤波器组设计山东轻工业学院课程设计任务书学院电子信息与控制工程学院专业通信工程姓名马淑丽班级通信09-2学号200902041044题目树形结构滤波器组设计 主要内容、基本要求、主要参考资料等：主要内容：滤波器组在语音、图像的子带编码和压缩中都有着广泛的应用，非均匀滤波器组还构成了 Mallat多分辨分析的算法基础，在小波变换中占有重要的地 位。本设计主要内容是研究树形滤波器组的原理，并设计一个树形滤波器组，实 现语音信号的分解与重构。基本要求：（1）滤波器组的基本原理；（2）树形结构滤波器组的原理及设计方法；（3）设计一个8通道的树形结构滤波器组：均匀滤波器组和非均匀滤波器组；给出设计思路及结果；（4）用设计的滤波器组对某信号进行多通道分解，验证滤 波器组的性能，对结果进行分析；（5）提交课程设计报告。主要参考资料：1. 胡广书现代信号处理教程，数字信号处理清华大学出版社2005.062. 高西全.数字信号处理.西安电子科技大学出版社.2009.013. matlab信号处理相关书籍，多采样率信号处理的书籍、资料4. 相关网络资源 完成期限：自2012年6月28日至2010年7月13日指导教师：张凯丽教研室主任：目录主要内容摘要设计方案设计原理设计框图设计程序结果图结果图分析结论及心得参考资料附录代码 内容摘要：树形结构滤波器组设计，将信源输入信息编码频带分段，便于在有限带宽信 道中传输并 且提高传输速率，在信宿端将信号解码恢复原始信号。有一定的失真。语音数据的有效编码 可以提高通信系统的有效性，大大减少存储设备的容量。子带编码是一种常用语音编码技术，子带编码中的子带分解和合成是子 带编码中的 重要组成部分。使用树形结构滤波器组实现语音信号的子带分解和合 并，常用的平行结构滤波 器虽然也可以实现自带的分解，实现对高频成分的压缩，但不如树形结构灵活，树形结构QMF 可以实现多分辨率的信号分解与压缩，同时重建信号失真度很低。设计方案本次课程设计，分别用对称结构和非对称结构滤波器组设计，实现语音信号或别的信号 3级分解8通道传输。我组用的matlab编程实现方法。一个语音处理系统主要包括语音信号的采集，预处理，语音信号的压缩编码，语音信号的解码，语音信号的增强，最后通过音频输出设备输出。为了能够使采 集到的语音信号能够 完全恢复出来，一般信号的采样频率都是很高的，例如44100HZ但是人耳能够识别的声音信号的频率范围在3003400HZ高于3400HZ的频率基本对人耳无效，因此可以滤除不予编码，同时在3003400的频率段也 有部分频率段 占用很少一部分能量，可以子带分解后用较短的码长编码，以此来降低码率，而对低频部分 可以通过抽取，来减少传输和处理的数据量。在接收端可以通过插值恢复出低频信号。本系统主要实现的关键步骤是针对语音信号的频谱设计与之相适应的树形 滤波器组，在信 源段首先对音频信号进行预处理滤除多余的频段，然后就是设计信源端得分析滤波器组和信宿 段综合滤波器组的设计。以以语音信号为例，这里只分析对称滤波器组的设计（非对称的自然就明白了）：将信 号经过2通道正交镜像滤波器组和2-抽取器，完成信号的第一次高低分频和抽取，然后分别 将分解出来的低频高频信号再次通过2通道正交镜像滤波器组和2-抽取器，实现对第一次分 解出来的低频高频信号的高低分频和抽取，将分解出来的4个频段成分在经过2通道正交镜像 滤波器组和2-抽取器，实现8通道，就是音频信号的需要传输或者处理的有用信号，然后对 分解出来的有用信号的高低频信号分别进行插值再对应两两合并，得到4个信号再次插值，合 并，得到2个信号，再次插值，合并，得到信宿端信号。这就是信源的分解滤波器组，将音 频信号分解成8通道信号，进行存储或传输等一系列处理。期间用到的滤波器等：1.预滤波时用到带通滤波器3003400hz。椭圆型带通滤波器function y1_yu pre_process y1Fs 8000;fpl 300;fpu 3400;fsl 200;fsu 3600;wp 2*fpl/Fs,2*fpu/Fs;ws 2*fsl/Fs,2*fsu/Fs;rp 1;rs 40;N,wpo ellipord wp,ws,rp,rs ;B,A ellip N,rp,rs,wpo ;H,w freqz B,A ;y1_yu filter B,A,y1 ; 2.正交镜像滤波器组含高通和低通滤波（1/2）用于分解信号为2个一个高频一个低频、，抽取用到2级抽取downsample函数function y1,y2 fenjie xfs 44100;t2 0:1/fs: length x -1 /fs ;a1 fir1 100,1/2 ;窗函数长度100.3db截止频率为1/2a2 qmf a1w1 filter a1,1,x ; % 低通滤波w2 filter a2,1,x ; % 高通滤波%抽取y1 downsample w1,2 ; % 抽取 2y2 downsample w2,2 ; % 抽取 2a1 fir1 100,1/2 ;a2 qmf -a1 ;figure 19x 0:1:100;subplot 2,1,1stem x,a1,.;subplot 2,1,2stem x,a2,.;插值2级用到interp函数、合并时用到求长度保持一致 后相加function y hebing y1,y2,n1,n2%信号的零插值和低通滤波y1 interp y1,2 ;y2 interp1 y2 ;n min length y1 ,length y2 ;y1 y1 1:n ;y2 y2 1:n ;%信号的合并y y1+y2 ;其中子程序 function f interp1 x%实现信号的零插值f interp x,2 ;for i 1:length xf 2*i 0 ;enda1 fir1 200,1/2,high;1、抽取和内插的基本原理信号的M-抽取是对原始信号每隔M-1个点取一个点，组成新的采样序列。抽取后的信号 频谱发生了变化，输入信号x (n)与M-抽取后的输出信号y (n)在频域上的关系式如下：丫 ( eAjw)打印后自己手写从上式可以知道，Y(ejw)是将X(ejw)扩展了 M倍，幅度变为原来的1 /M，在分别以2pi、4pi、6pi， ?， 2kpi移位叠加得到的；抽取后可能造成频谱 混叠，为了避免 抽取后的频谱混叠，信号X ( n )的带宽必须限制在【-pi/M，pi/M】。通常情况下可以 再抽取器前进行抗混叠滤波，所谓的抗混叠滤波就是在抽取前对信号进行低通滤波，把信号的 频带限制在【-pi/M， pi/M】。信号的内插是信号抽取的逆过程，是在已知信号相邻抽样点之间插入若干个抽样值的 点。实际中的做法是在已知抽样序列之间插入若干零值，然后通过低通滤波器，几个实现内插。其实就相当于线性插值。L-内插输入信号x (n)与 输出信号的频域关系式：Y (ejw ) X (ejwL)从上式我们知道，对信号的L-内插相当于对输入信号频谱的L倍压缩，内插后信号的 采样周期变为原来的1/L倍。内插不会造成频谱混叠，不会造成信 号信息量的丢失，但会使 整个数字信号频率轴插入L-1个原始信号的频谱，称为镜像。因此在插值后我们可以通过一 个低通滤波器，来消除内插带来的镜像，恢 复出原始信号。2、滤波器组的基本原理 滤波器组是一组拥有共同输入信号或共同输出信号的一组带通滤 波器。(系统框图如下)M通道滤波器组的系统结构(手绘)信号的子带分解是通过树形结构滤波器组来实现。本系统中有一个分析滤波 器组实现对 输入信号的子带分解，一个综合滤波器组完成信号的重建。一个给定的信号经过分解滤波器 组分解，然后编码、传输再通过综合滤波器组实现信号的恢复和重建。但是恢复和重建后的 信号并不能与原始信号完全相同，两者之间会 产生一定的误差，主要包括：(1) 混叠失真：由抽取和内插产生的混叠和镜像带来的误差所造成的。(2) 幅度失真：由滤波器组幅频特性波纹产生的误差所造成的。(3) 相位失真：由滤波器组的相频特性的非线性产生的误差所造成的。(4) 子带量化误差：由编解码产生的误差，与量化噪声相似。这是一种无 法完全消除 的误差。完全重建滤波器组：无混叠失真的滤波器组，同时既无幅度失真，又无相位 失真。本系统所设计的树形结构滤波器组，是由两通道的正交镜像滤波器组通过级 联或并联组 建而成的。3. 正交镜像滤波器组两通道正交镜像滤波器组的系统框图如下(手绘)上图中可以设H0(z)H(z)(1.1 )(H1(z)H(1.2)z )G0( z) H z 1.3G1 z -H -z(1.4)上式说明如果H0 (z)是低通滤波器，那么H1 (在)则是高通滤波器。同时H0 ejw与H1 ejw关于pi/2成镜像对称，所以称这种滤波器成为正交镜像滤波器组，简称QMFB满足1.1式1.4式的滤波器组称为标准QMF滤波器组，它是一种无混叠失真 的滤波器组。树形结构的滤波器组：将两通道滤波器组的级联来实现多通道滤波器组。具 体做法就是 将各通道的输出作为下一个滤波器组的输入。这种实现方式的优点 是：可根据两通道滤波器组的特性来推断它的特性（是否能够完全重建等），通 过非对称来实现多抽样率的子带分割。两通道QMFBE对称的分析滤波器框图（手绘）三设计框图（手绘）四设 计程序（1）对称树形结构滤波器组设计程序（见附页）非对称树形滤波器组设计 程序（见附页）五.设计结果图1.对称结构原始信号和预滤波后信号的时域波形对比原始信号和预滤波后的信号的频域对比三次分解后高频信号和低频信号的频谱 三次分解 后高频信号和低频信号的信号与信宿端恢复出来的信号的时域对 比信号与信宿端恢复出来的信号的频域对比非对称结构原始信号和预滤波后信号的时域波形对比原始信号和预滤波后的信号的频域对比信 号与信宿端恢复出来的信号的时域对比信号与信宿端恢复出来的信号的频域对比六.设计结果 图分析 本系统首先对输入的信号进行预滤波，滤除频段（3003400）以外的信号，本系统设计 了一个椭圆形带通滤波器，从实验图看出，本系统设计的预滤波器完 全满足系统的滤波要 求，滤除了低于300HZ和高于 3400HZ的频率分量，同时也 没有带来相位失真。对于非对称结 构：然后将预滤波后的信号进行3级非对称分解，得到低频信号和高频信号，从图可以看出 进行三级分解后的低频信号和高频信号均扩展到整个数字频率轴，且 高频信号的能量已经占到 很大比重，此时已经把有用信号的频率分量已经分离出 来。从图可以看出，8通道信号能量由低频到高频，呈现由高到低完全符合声音 信号的能量 分布，分解效果很好。在接受端进行信号的综合，从图可以看出综合后的高频信号频谱和低频信号频谱，通过对比可以看出，信源端分解前的信号频谱，分解后再综合信号的频谱 相似度很 高，符合设计的要求。重建后信号与原始信号的时域和频域对比图，可以看出，恢复出来的信号与原始信号仅 在能量上存在一定的线性失真，可以通过放大器予以恢复和重建。最后将恢复出来的信号写入文件，在播放器上进行播放，与原始信号相比，听不出来差 别，本系统完全满足设计要求。七.设计心得在本次课程设计中学会了好多matlab函数的使用，各种滤波器的设计，了解了正交镜像 滤波器工作原理，抽样，插值的意义，用数形结构滤波器组更方便信号的传输高效性，了解 了利用树形结构滤波器的工作原理，更加深了信号在信 道中传输的过程。学会了好多数字信号 及matlab应用方面的知识，起初用simlink设计时对参数的设计我们遇到了瓶颈改用 matlab编程设计，大大加深了 matlab的广泛应用性。期间查阅了许多资料，从网上下载了一个程序，一点也看不懂，后经过图书 馆借书 查，询问同学老师，解决了不少问题。并通过对程序的更改及自己的理解，将程序改为8通 道的树形结构滤波器。本次设计做了两个结构滤波器，其实掌握了对称的结构原理，非对称 的就很容易了。另外，了解了语音信号的读取及回放用matlab实现。了解了语音信号的3级分解抽取后 经过传输后插值3级合并过程的意义。只要用心去做多动脑子多查资料多问，其实课程设计没有起初自己想的那么 难。感谢老 师给我了这次学习的机会，让我获益匪浅。八.参考资料1.胡广书现代信号处理教程，数字信号处理.清华大学出版社2005.062.高西全.数字信号处理西安电子科技大学出版社2009.013. Matlab7辅助信号处理技术与应用电子工业大学出版社4. 网络附页程序对称结构程序clear all ;N 4096*8 ;%采集点数y,fs,bits wavread T.wav,20000,20000+NT ; % 读取音频信息y1 y ：,1 ;%y2 y :,2 ;t 0:1/fs: N-1 /fs;%timefigure 1 ;subplot 2,2,1 ;plot t,y1 ;%1通道信号title 原始信号；xlabel t ,ylabel 幅度y1_yu pre_process y1 ;%语音信号预处理人耳可识别范围（3003400）,因此多于此频率的信号可以滤掉而不用传输subplot 2,2,2 ;plot t,y1_yu ;%滤波后的信号时域波形title 预滤波后的信号波形；xlabel t ,ylabel 幅度y11 fft y1,fsaxis 0,44100,0,500subplot 2,2,4 ;y1_yu1 fft y1_yu ;plot abs y1_yu1 ;title 预滤波后信号频谱；xlabel f ,ylabel 幅度次分解axis 0,44100,0,500f1 ,f2 fenjie y1_yu,1 ;%1f1_1,f1_2 fenjie f1,1,10 ;%21 f2_1,f2_2 fenjie f2,1,10 ;%22f1_1_1,f1_1_2 fenjie f1_1,1,10f1_2_1,f1_2_2 fenjie f1_2,1,10f2_1_1,f2_1_2 fenjie f2_1,1,10f2_2_1,f2_2_2 fenjie f2_2,1,10次分解;%31,32分解;%33,34分解;%35,36分解;%37,38次分解次分解figure 2subplot 4,2,1n length f1_1_1t 0:8/fs:8.* n-1 /fsn length f1_1_2plot t,f1_1_2 ;title 3级第2个分解并抽取后部分；xlabel t ,ylabel 幅度 subplot 4,2,3n length f1_2_1plot t,f1_2_1title 3级第3个分解并抽取后部分；xlabel t ,ylabel 幅度 subplot 4,2,4n length f1_2_2plot t,f1_2_2 ;title 3级第4个分解并抽取后部分；xlabel t ,ylabel 幅度 subplot 4,2,5n length f2_1_1plot t,f2_1_1title 3级第5个分解并抽取后部分；xlabel t ,ylabel 幅度 subplot 4,2,6n length f2_1_2plot t,f2_1_2 ;title 3级第6个分解并抽取后部分；xlabel t ,ylabel 幅度 subplot 4,2,7subplot 4,2,8n length f2_2_2plot t,f2_2_2 ;title 3级第8个分解并抽取后部分；xlabel t ,ylabel 幅度figure 3f1111 fft f1_1_1,fs ;subplot 4,2,1plot abs f1111title 3级第1个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1112 fft f1_1_2,fs ;subplot 4,2,2plot abs f1112title 3级第2个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1113 fft f1_2_1,fs ;subplot 4,2,3plot abs f1113title 3级第3个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1114 fft f1_2_2,fs ;subplot 4,2,4f1115 fft f2_1_1,fs ;subplot 4,2,5plot abs f1115title 3级第5个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1116 fft f2_1_2,fs ;subplot 4,2,6plot abs f1116title 3级第6个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1117 fft f2_2_1,fs ;subplot 4,2,7plot abs f1117title 3级第7个分解并抽取后部分频谱；xlabel f ,ylabel 幅度f1118 fft f2_2_2,fs ;subplot 4,2,8plot abs f1118title 3级第8个分解并抽取后部分频谱；xlabel f ,ylabel 幅度yo hebing8 f1_1_1,f1_1_2,f1_2_1,f1_2_2,f2_1_1,f2_1_2,f2_2_1,f2_2_2 ;figure 4 ;%合成后的信号与原来预滤波后的信号对比subplot 4,1,1 ;t 0:1/fs: N-1 /fs;plot t,y1_yu ;subplot 4,1,3 ;n length yo ;t 0:1/fs: n-1 /fs;subplot 4,1,4yo1 fft yo,fsplot abs yo1 ;title 信宿恢复信号频谱；xlabel f ,ylabel 幅度subplot 4,1,2y1_yu1 fft y1_yu ;plot abs y1_yu1 ;title 预滤波后信号频谱；xlabel f ,ylabel 幅度wavwrite yo,fs,3.wav;子程序树形结构分解子程序 完成信号的2通道分解以及降采样率function y1,y2 fenjie x,n1,n2,fpfs 44100;t2 0:1/fs: length x -1 /fs ;a1 fir1 100,1/2 ;a2 qmf a1 ;w1 filter a1,1,x ; % 低通滤波w2 filter a2,1,x ; % 高通滤波%抽取 y1 downsample w1,2 ; % 抽取 2y1_1_1,y1_1_2,y1_2_1,y1_2_2,y2_1_1,y2_1_2,y2_2_1,y2_2_2%信号的3次合并y1_1 hebing y1_1_1,y1_1_2,1,10 ;y1_2 hebing y1_2_1,y1_2_2,1,10 ;y2_1 hebing y2_1_1,y2_1_2,1,10 ;y2_2 hebing y2_2_1,y2_2_2,1,10 ;%信号的 2 次合并y1 hebing y1_1,y1_2,1,10 ;y2 hebing y2_1,y2_2,1,10 ;%信号的 1 次合并f hebing y1,y2,1,10 ; %树形结构综合子程序（完成2通道信号的综合）functiony hebing y1,y2,n1,n2%信号的零插值和低通滤波y1 interp y1,2 ;y2 interp1 y2 ;n min length y1 ,length y2 ;y1 y1 1:n ;y2 y2 1:n ;%信号的合并y y1+y2 ; %信号插值子程序（完成信号的零插值和高通滤波）function f interp1 x %实现信号的零插值 f interp x,2 ;for i 1:length xf 2*i 0 ;endal firl 200,1/2,high;f filter a1,1,f ; function y1_yu pre_process ylFs 8000;fpl 300;fpu 3400;fsl 200;fsu 3600;wp 2*fpl/Fs,2*fpu/Fs;ws 2*fsl/Fs,2*fsu/Fs;rp 1;rs 40;N,wpo ellipord wp,ws,rp,rs ;B,A ellip N,rp,rs,wpo ;H,w freqz B,A ;y1_yu filter B,A,y1 ;非对称主程序clear all ;N 4096*8 ;%采集点数读取音频信息（双y,fs,bits wavread T.wav,20000,20000+NT ; %声道，16位，频率44100Hz)y1 y ：,1 ;%y2 y :,2 ;t 0:1/fs: N-1 /fs;%timefigure 1 ;subplot 2,2,1 ;plot t,y1 ;%1通道信号title 原始信号时域波形；xlabel t ,ylabel 幅度%预处理人耳可识别范围（3003400）y1_yu pre_process y1 ;subplot 2,2,2 ;plot t,y1_yu ;title 预滤波后的信号时域波形；xlabel t ,ylabel 幅度subplot2,2,3 ;y11 fft y1,fs ;plot abs y11 ;title 原始信号频谱；subplot 2,2,4 ;y1_yu1 fft y1_yu ;plot abs y1_yu1title 滤波后信号频谱 ;f1 ,f2 fenjie y1_yu,1,10 ;%1 次分解f1_1,f1_2 fenjie f1,1,10 ;%2 次分解f1_11,f1_12 fenjie f1_1,1,10 ;%3 次分解f1_111,f1_112 fenjief1_11,1,10 ;%4 次分解f1_1111,f1_1112 fenjie f1_111,1,10 ;%5 次分解f1_11111,f1_11112 fenjie f1_1111,1,10 ;%6 次分解f1_111111,f1_111112fenjie f1_11111,1,10 ;%7 次分解 figure 2 subplot 4,2,1 n length f2t 0:2/fs:2.* n-1 /fs plot t,f2 ;title 一次分解并抽取的高频部分；xlabel t ,ylabel 幅度subplot 4,2,2 nlength f1_2t 0:4/fs:4.* n-1 /fs plot t,f1_2 ;title 二次分解并抽取的高频部分；xlabel t ,ylabel 幅度subplot 4,2,3n length f1_12t 0:8/fs:8.* n-1 /fsplot t,f1_12 ;title 三次分解并抽取的高频部分；xlabel t ,ylabel 幅度subplot 4,2,4n length f1_112t 0:16/fs:16.* n-1 /fsplot t,f1_112 ;title 四次分解并抽取的高频部分；xlabel t ,ylabel 幅度subplot 4,2,5t 0:32/fs:32.* n-1 /fs plot t,f1_1112 ;title 五次分解并抽取的高频部分；xlabel t ,ylabel 幅度 subplot 4,2,6n length f1_11112t 0:64/fs:64.* n-1 /fsplot t,f1_11112 ;title 六次分解并抽取的高频部分；xlabel t ,ylabel 幅度 subplot 4,2,7n length f1_111112t 0:128/fs:128.* n-1 /fsplot t,f1_111112 ;title 七次分解并抽取的高频部分；xlabel t ,ylabel 幅度subplot 4,2,8n length f1_111111t 0:128/fs:128.* n-1 /fsplot t,f1_111111title 七次分解并抽取的低频部分figure 3f1111 fft f1_111111,fs ;subplot 4,2,1plot abs f1111title 7级分解并抽取后低频部分频谱f1112 fft f1_111112,fs ;subplot 4,2,2plot abs f1112title 7级分解并抽取后高频部分频谱f1113 fft f1_11112,fs ;subplot 4,2,3plot abs f1113title 6级分解并抽取后高频部分频谱f1114 fft f1_1112,fs ;subplot 4,2,4 plot abs f1114title 5级分解并抽取后高频部分频谱xlabel t ,ylabel 幅度;xlabel f ,ylabel;xlabel f ,ylabel;xlabel f ,ylabel;xlabel f ,ylabelf1115 fft f1_112,fs ;subplot 4,2,5plot abs f1115title 4级分解并抽取后高频部分频谱f1116 fft f1_12,fs ;subplot 4,2,6plot abs f1116title 3级分解并抽取后高频部分频谱f1117 fft f1_2,fs ;subplot 4,2,7plot abs f1117title 2级分解并抽取后高频部分频谱f1118 fft f2,fs ;subplot 4,2,8plot abs f1118title 1级分解并抽取后高频部分频谱g1hebing f1_111111,f1_111112 ;%1 g2hebing g1,f1_11112 ;%2 次合并 g3hebing g2,f1_1112 ;%3 次合并 g4 hebingg3,f1_112 ;%4 次合并 g5 hebingg4,f1_12 ;%5 次合并;xlabel f ,ylabel 幅度;xlabel f ,ylabel 幅度;xlabel f ,ylabel 幅度;xlabel f ,ylabel 幅度次合并g6 hebing g5,f1_2 ;%6 次合并g7 hebing g6,f2 ;%7 次合并figure 4subplot 4,1,1 ;t 0:1/fs: N-1 /fs;plot t,y1_yu ;title 预滤波后的信号；xlabel t ,ylabel 幅度 subplot 4,1,3 ;n length g7 ;t 0:1/fs: n-1 /fs;plot t,g7 ;title 信宿端恢复信号；xlabel t ,ylabel 幅度subplot 4,1,4g71 fft g7,fsplot abs g71 ;title 信宿恢复信号频谱；xlabel f ,ylabel 幅度subplot 4,1,2y1_yu1 fft y1_yu ;plot abs y1_yu1 ;title 预滤波后信号频谱；xlabel f ,ylabel 幅度%将处理后的信号写入文件wavwrite g7,fs,2.wav;子程序%树形结构分解子程序完成信号的2通道分解以及降采样率functiony1,y2 fenjie x,n1,n2,fp fs 44100;t2 0:1/fs: length x -1 /fs % al 1 1;% b1 1;% a2 1 -1;% b2 1;a1 fir1 100,1/2 ;a2 qmf a1 ;w1 filter a1,1,x ; % 低通滤波w2 filter a2,1,x ; % 高通滤波y1 downsample w1,2 ; % 抽取 2y2 downsample w2,2 ; % 抽取 2 function y hebing y1,y2,n1,n2%信号的零插值和低通滤波y1 interp y1,2 ;%信号的零插值和高通滤波y2 interp1 y2 ;%进行插值后的信号可能长度不一样，因此要进行长信号的截取n min length y1 ,length y2 ;y1 y1 1:n ;y2 y2 1:n ;%信号的合并y y1+y2 ; %信号插值子程序（完成信号的零插值和高通滤波）function f interpl x%实现信号的零插值f interp x,2 for i 1:length xf 2*i 0 ;enda1 fir1 200,1/2,high;f filter a1,1,f ; function y1_yu pre_process y1Fs 8000;fpl 300;fpu 3400;fsl 200;fsu 3600;wp 2*fpl/Fs,2*fpu/Fs;ws 2*fsl/Fs,2*fsu/Fs;rp 1;rs 40;N,wpo ellipord wp,ws,rp,rs ;B,A ellip N,rp,rs,wpo ;H,w freqz B,A ;y1_yu filter B,A,y1 ;-