欢迎来到沃文网! | 帮助中心 分享知识,传播智慧!
沃文网
换一换
首页 沃文网 > 资源分类 > DOC文档下载
 

树形结构滤波器组设计

  • 资源ID:20039       资源大小:275.62KB        全文页数:28页
  • 资源格式: DOC        下载权限:游客/注册会员/VIP会员    下载费用:10积分 【人民币10元】
快捷注册下载 游客一键下载
会员登录下载
三方登录下载: QQ登录   微博登录  
下载资源需要10积分 【人民币10元】
邮箱/手机:
温馨提示:
支付成功后,系统会自动生成账号(用户名和密码都是您填写的邮箱或者手机号),方便下次登录下载和查询订单;
支付方式: 微信支付    支付宝   
验证码:   换一换

加入VIP,免费下载资源
 
友情提示
2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,既可以正常下载了。
3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

树形结构滤波器组设计

山东轻工业学院课程设计任务书学院电子信息与控制工程学院专业通信工程题目树形结构滤波器组设计主要内容、基本要求、主要参考资料等主要内容滤波器组在语音、图像的子带编码和压缩中都有着广泛的应用,非均匀滤波器组还构成了Mallat多分辨分析的算法基础,在小波变换中占有重要的地位。本设计主要内容是研究树形滤波器组的原理,并设计一个树形滤波器组,实现语音信号的分解与重构。基本要求(1)滤波器组的基本原理;(2)树形结构滤波器组的原理及设计方法;(3)设计一个8通道的树形结构滤波器组均匀滤波器组和非均匀滤波器组;给出设计思路及结果;(4)用设计的滤波器组对某信号进行多通道分解,验证滤波器组的性能,对结果进行分析;(5)提交课程设计报告。主要参考资料1.胡广书.现代信号处理教程,数字信号处理.清华大学出版社.2005.062.高西全.数字信号处理.西安电子科技大学出版社.2009.013.matlab信号处理相关书籍,多采样率信号处理的书籍、资料。4.相关网络资源完成期限自2012年6月28日至2010年7月13日指导教师张凯丽教研室主任目录主要内容摘要一、设计方案二、设计原理三、设计框图四、设计程序五、结果图六、结果图分析七、结论及心得八、参考资料附录代码内容摘要树形结构滤波器组设计,将信源输入信息编码频带分段,便于在有限带宽信道中传输并且提高传输速率,在信宿端将信号解码恢复原始信号。有一定的失真。语音数据的有效编码可以提高通信系统的有效性,大大减少存储设备的容量。子带编码是一种常用语音编码技术,子带编码中的子带分解和合成是子带编码中的重要组成部分。使用树形结构滤波器组实现语音信号的子带分解和合并,常用的平行结构滤波器虽然也可以实现自带的分解,实现对高频成分的压缩,但不如树形结构灵活,树形结构QMFB可以实现多分辨率的信号分解与压缩,同时重建信号失真度很低。一.设计方案本次课程设计,分别用对称结构和非对称结构滤波器组设计,实现语音信号或别的信号3级分解8通道传输。我组用的matlab编程实现方法。一个语音处理系统主要包括语音信号的采集,预处理,语音信号的压缩编码,语音信号的解码,语音信号的增强,最后通过音频输出设备输出。为了能够使采集到的语音信号能够完全恢复出来,一般信号的采样频率都是很高的,例如44100HZ,但是人耳能够识别的声音信号的频率范围在3003400HZ,高于3400HZ的频率基本对人耳无效,因此可以滤除不予编码,同时在3003400的频率段也有部分频率段占用很少一部分能量,可以子带分解后用较短的码长编码,以此来降低码率,而对低频部分可以通过抽取,来减少传输和处理的数据量。在接收端可以通过插值恢复出低频信号。本系统主要实现的关键步骤是针对语音信号的频谱设计与之相适应的树形滤波器组,在信源段首先对音频信号进行预处理滤除多余的频段,然后就是设计信源端得分析滤波器组和信宿段综合滤波器组的设计。以以语音信号为例,这里只分析对称滤波器组的设计(非对称的自然就明白了)将信号经过2通道正交镜像滤波器组和2-抽取器,完成信号的第一次高低分频和抽取,然后分别将分解出来的低频高频信号再次通过2通道正交镜像滤波器组和2-抽取器,实现对第一次分解出来的低频高频信号的高低分频和抽取,将分解出来的4个频段成分在经过2通道正交镜像滤波器组和2-抽取器,实现8通道,就是音频信号的需要传输或者处理的有用信号,然后对分解出来的有用信号的高低频信号分别进行插值再对应两两合并,得到4个信号再次插值,合并,得到2个信号,再次插值,合并,得到信宿端信号。这就是信源的分解滤波器组,将音频信号分解成8通道信号,进行存储或传输等一系列处理。期间用到的滤波器等1.预滤波时用到带通滤波器3003400hz。椭圆型带通滤波器functiony1_yupre_processy1Fs8000;fpl300;fpu3400;fsl200;fsu3600;wp[2*fpl/Fs,2*fpu/Fs];ws[2*fsl/Fs,2*fsu/Fs];rp1;rs40;[N,wpo]ellipordwp,ws,rp,rs;[B,A]ellipN,rp,rs,wpo;[H,w]freqzB,A;y1_yufilterB,A,y1;2.正交镜像滤波器组含高通和低通滤波(1/2)用于分解信号为2个一个高频一个低频、,抽取用到2级抽取downsample函数function[y1,y2]fenjiexfs44100;t201/fslengthx-1/fs;a1fir1100,1/2;窗函数长度100.3db截止频率为1/2a2qmfa1;w1filtera1,1,x;低通滤波w2filtera2,1,x;高通滤波抽取y1downsamplew1,2;抽取2y2downsamplew2,2;抽取2其中参数值子程序求H0Z和G0Z的值(1)H0Za1fir1100,1/2;a2qmfa1;figure19x01100;subplot2,1,1stemx,a1,.;subplot2,1,2stemx,a2,.;(2)G0Za1fir1100,1/2;a2qmf-a1;figure19x01100;subplot2,1,1stemx,a1,.;subplot2,1,2stemx,a2,.;3.插值2级用到interp函数、合并时用到求长度保持一致后相加functionyhebingy1,y2,n1,n2信号的零插值和低通滤波y1interpy1,2;y2interp1y2;nminlengthy1,lengthy2;y1y11n;y2y21n;信号的合并yy1y2;其中子程序functionfinterp1x实现信号的零插值finterpx,2;fori1lengthxf2*i0;enda1fir1200,1/2,high;ffiltera1,1,f;非对称滤波器组设计原理类似的(略)二.设计原理1、抽取和内插的基本原理信号的M-抽取是对原始信号每隔M-1个点取一个点,组成新的采样序列。抽取后的信号频谱发生了变化,输入信号x(n)与M抽取后的输出信号y(n)在频域上的关系式如下Y(ejw)打印后自己手写从上式可以知道,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(-z)(1.2)G0(z)Hz1.3G1z-H-z(1.4)上式说明如果H0(z)是低通滤波器,那么H1(在)则是高通滤波器。同时H0ejw与H1ejw关于pi/2成镜像对称,所以称这种滤波器成为正交镜像滤波器组,简称QMFB。满足1.1式1.4式的滤波器组称为标准QMF滤波器组,它是一种无混叠失真的滤波器组。树形结构的滤波器组将两通道滤波器组的级联来实现多通道滤波器组。具体做法就是将各通道的输出作为下一个滤波器组的输入。这种实现方式的优点是可根据两通道滤波器组的特性来推断它的特性(是否能够完全重建等),通过非对称来实现多抽样率的子带分割。两通道QMFB非对称的分析滤波器框图(手绘)三.设计框图(手绘)四.设计程序(1)对称树形结构滤波器组设计程序(见附页)(2)非对称树形滤波器组设计程序(见附页)五.设计结果图1.对称结构原始信号和预滤波后信号的时域波形对比原始信号和预滤波后的信号的频域对比三次分解后高频信号和低频信号的频谱(8个)三次分解后高频信号和低频信号的时域波形(8个)预滤波信号与信宿端恢复出来的信号的时域对比预滤波信号与信宿端恢复出来的信号的频域对比正交镜像滤波器组的高频低频输出值参数00.511.5-1-0.500.51原始信号t幅度00.511.5-2-1012预滤波后的信号波形t幅度01234x1040100200300400500原信号的频谱f幅度01234x1040100200300400500预滤波后信号频谱f幅度00.511.5-0.500.53级第1个分解并抽取后部分t幅度00.511.5-2023级第2个分解并抽取后部分t幅度00.511.5-1013级第3个分解并抽取后部分t幅度00.511.5-1013级第4个分解并抽取后部分t幅度00.511.5-505x10-33级第5个分解并抽取后部分t幅度00.511.5-0.200.23级第6个分解并抽取后部分t幅度00.511.5-1013级第7个分解并抽取后部分t幅度00.511.5-1013级第8个分解并抽取后部分t幅度0123x104020403级第1个分解并抽取后部分频谱f幅度0123x1040501003级第2个分解并抽取后部分频谱f幅度0123x104020403级第3个分解并抽取后部分频谱f幅度0123x104020403级第4个分解并抽取后部分频谱f幅度0123x10400.20.43级第5个分解并抽取后部分频谱f幅度0123x1040243级第6个分解并抽取后部分频谱f幅度0123x104010203级第7个分解并抽取后部分频谱f幅度0123x104010203级第8个分解并抽取后部分频谱f幅度00.511.5-202预滤波后的信号t幅度00.511.5-101信宿端恢复信号t幅度00.511.522.5x1040100200信宿恢复信号频谱f幅度00.511.522.533.5x1040200400预滤波后信号频谱f幅度0102030405060708090100-0.200.20.40.60102030405060708090100-0.4-0.200.20.40.60102030405060708090100-0.200.20.40.60102030405060708090100-0.6-0.4-0.200.20.42.非对称结构原始信号和预滤波后信号的时域波形对比原始信号和预滤波后的信号的频域对比预滤波信号与信宿端恢复出来的信号的时域对比预滤波信号与信宿端恢复出来的信号的频域对比七级非对称分解后8通道信号的频谱及时域波形00.511.5-1-0.500.51原始信号时域波形t幅度00.511.5-2-1012预滤波后的信号时域波形t幅度0123x1040200400600原始信号频谱01234x1040100200300400滤波后信号频谱00.511.5-101一次分解并抽取的高频部分t幅度00.511.5-202二次分解并抽取的高频部分t幅度00.511.5-202三次分解并抽取的高频部分t幅度00.511.5-0.500.5四次分解并抽取的高频部分t幅度00.511.5-0.0100.01五次分解并抽取的高频部分t幅度00.511.5-101x10-3六次分解并抽取的高频部分t幅度00.511.5-505x10-4七次分解并抽取的高频部分t幅度00.511.5-505x10-4七次分解并抽取的低频部分t幅度0123x104024x10-37级分解并抽取后低频部分频谱f幅度0123x104012x10-37级分解并抽取后高频部分频谱f幅度0123x10400.0050.016级分解并抽取后高频部分频谱f幅度0123x10400.10.25级分解并抽取后高频部分频谱f幅度0123x104010204级分解并抽取后高频部分频谱f幅度0123x1040501003级分解并抽取后高频部分频谱f幅度0123x1040501002级分解并抽取后高频部分频谱f幅度0123x1040501001级分解并抽取后高频部分频谱f幅度00.511.5-202预滤波后的信号t幅度00.511.5-101信宿端恢复信号t幅度00.511.522.5x1040100200信宿恢复信号频谱f幅度00.511.522.533.5x1040200400预滤波后信号频谱f幅度六.设计结果图分析本系统首先对输入的信号进行预滤波,滤除频段(3003400)以外的信号,本系统设计了一个椭圆形带通滤波器,从实验图看出,本系统设计的预滤波器完全满足系统的滤波要求,滤除了低于300HZ和高于3400HZ的频率分量,同时也没有带来相位失真。对于非对称结构然后将预滤波后的信号进行3级非对称分解,得到低频信号和高频信号,从图可以看出进行三级分解后的低频信号和高频信号均扩展到整个数字频率轴,且高频信号的能量已经占到很大比重,此时已经把有用信号的频率分量已经分离出来。从图可以看出,8通道信号能量由低频到高频,呈现由高到低完全符合声音信号的能量分布,分解效果很好。在接受端进行信号的综合,从图可以看出综合后的高频信号频谱和低频信号频谱,通过对比可以看出,信源端分解前的信号频谱,分解后再综合信号的频谱相似度很高,符合设计的要求。重建后信号与原始信号的时域和频域对比图,可以看出,恢复出来的信号与原始信号仅在能量上存在一定的线性失真,可以通过放大器予以恢复和重建。最后将恢复出来的信号写入文件,在播放器上进行播放,与原始信号相比,听不出来差别,本系统完全满足设计要求。七.设计心得在本次课程设计中学会了好多matlab函数的使用,各种滤波器的设计,了解了正交镜像滤波器工作原理,抽样,插值的意义,用数形结构滤波器组更方便信号的传输高效性,了解了利用树形结构滤波器的工作原理,更加深了信号在信道中传输的过程。学会了好多数字信号及matlab应用方面的知识,起初用simlink设计时对参数的设计我们遇到了瓶颈改用matlab编程设计,大大加深了matlab的广泛应用性。期间查阅了许多资料,从网上下载了一个程序,一点也看不懂,后经过图书馆借书查,询问同学老师,解决了不少问题。并通过对程序的更改及自己的理解,将程序改为8通道的树形结构滤波器。本次设计做了两个结构滤波器,其实掌握了对称的结构原理,非对称的就很容易了。另外,了解了语音信号的读取及回放用matlab实现。了解了语音信号的3级分解抽取后经过传输后插值3级合并过程的意义。只要用心去做多动脑子多查资料多问,其实课程设计没有起初自己想的那么难。感谢老师给我了这次学习的机会,让我获益匪浅。八.参考资料1.胡广书.现代信号处理教程,数字信号处理.清华大学出版社.2005.062.高西全.数字信号处理.西安电子科技大学出版社.2009.013.Matlab7辅助信号处理技术与应用.电子工业大学出版社4.网络附页程序(1)对称结构程序clearall;N4096*8;采集点数[y,fs,bits]wavread1.wav,[20000,20000N-1];读取音频信息y1y,1;y2y,2;t01/fsN-1/fs;timefigure1;subplot2,2,1;plott,y1;1通道信号title原始信号;xlabelt,ylabel幅度y1_yupre_processy1;语音信号预处理人耳可识别范围(3003400),因此多于此频率的信号可以滤掉而不用传输subplot2,2,2;plott,y1_yu;滤波后的信号时域波形title预滤波后的信号波形;xlabelt,ylabel幅度subplot2,2,3;y11ffty1,fs;plotabsy11;title原信号的频谱;xlabelf,ylabel幅度axis[0,44100,0,500]subplot2,2,4;y1_yu1ffty1_yu;plotabsy1_yu1;title预滤波后信号频谱;xlabelf,ylabel幅度axis[0,44100,0,500][f1,f2]fenjiey1_yu,1;1次分解[f1_1,f1_2]fenjief1,1,10;21次分解[f2_1,f2_2]fenjief2,1,10;22次分解[f1_1_1,f1_1_2]fenjief1_1,1,10;31,32分解[f1_2_1,f1_2_2]fenjief1_2,1,10;33,34分解[f2_1_1,f2_1_2]fenjief2_1,1,10;35,36分解[f2_2_1,f2_2_2]fenjief2_2,1,10;37,38次分解figure2subplot4,2,1nlengthf1_1_1t08/fs8.*n-1/fsplott,f1_1_1title3级第1个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,2nlengthf1_1_2plott,f1_1_2;title3级第2个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,3nlengthf1_2_1plott,f1_2_1title3级第3个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,4nlengthf1_2_2plott,f1_2_2;title3级第4个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,5nlengthf2_1_1plott,f2_1_1title3级第5个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,6nlengthf2_1_2plott,f2_1_2;title3级第6个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,7nlengthf2_2_1plott,f2_2_1title3级第7个分解并抽取后部分;xlabelt,ylabel幅度subplot4,2,8nlengthf2_2_2plott,f2_2_2;title3级第8个分解并抽取后部分;xlabelt,ylabel幅度figure3f1111fftf1_1_1,fs;subplot4,2,1plotabsf1111title3级第1个分解并抽取后部分频谱;xlabelf,ylabel幅度f1112fftf1_1_2,fs;subplot4,2,2plotabsf1112title3级第2个分解并抽取后部分频谱;xlabelf,ylabel幅度f1113fftf1_2_1,fs;subplot4,2,3plotabsf1113title3级第3个分解并抽取后部分频谱;xlabelf,ylabel幅度f1114fftf1_2_2,fs;subplot4,2,4plotabsf1114title3级第4个分解并抽取后部分频谱;xlabelf,ylabel幅度f1115fftf2_1_1,fs;subplot4,2,5plotabsf1115title3级第5个分解并抽取后部分频谱;xlabelf,ylabel幅度f1116fftf2_1_2,fs;subplot4,2,6plotabsf1116title3级第6个分解并抽取后部分频谱;xlabelf,ylabel幅度f1117fftf2_2_1,fs;subplot4,2,7plotabsf1117title3级第7个分解并抽取后部分频谱;xlabelf,ylabel幅度f1118fftf2_2_2,fs;subplot4,2,8plotabsf1118title3级第8个分解并抽取后部分频谱;xlabelf,ylabel幅度yohebing8f1_1_1,f1_1_2,f1_2_1,f1_2_2,f2_1_1,f2_1_2,f2_2_1,f2_2_2;figure4;合成后的信号与原来预滤波后的信号对比subplot4,1,1;t01/fsN-1/fs;plott,y1_yu;title预滤波后的信号;xlabelt,ylabel幅度subplot4,1,3;nlengthyo;t01/fsn-1/fs;plott,yo;title信宿端恢复信号;xlabelt,ylabel幅度subplot4,1,4yo1fftyo,fsplotabsyo1;title信宿恢复信号频谱;xlabelf,ylabel幅度subplot4,1,2y1_yu1ffty1_yu;plotabsy1_yu1;title预滤波后信号频谱;xlabelf,ylabel幅度wavwriteyo,fs,3.wav;子程序树形结构分解子程序完成信号的2通道分解以及降采样率function[y1,y2]fenjiex,n1,n2,fpfs44100;t201/fslengthx-1/fs;a1fir1100,1/2;a2qmfa1;w1filtera1,1,x;低通滤波w2filtera2,1,x;高通滤波抽取y1downsamplew1,2;抽取2y2downsamplew2,2;抽取2functionfhebing8y1_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_1hebingy1_1_1,y1_1_2,1,10;y1_2hebingy1_2_1,y1_2_2,1,10;y2_1hebingy2_1_1,y2_1_2,1,10;y2_2hebingy2_2_1,y2_2_2,1,10;信号的2次合并y1hebingy1_1,y1_2,1,10;y2hebingy2_1,y2_2,1,10;信号的1次合并fhebingy1,y2,1,10;树形结构综合子程序(完成2通道信号的综合)functionyhebingy1,y2,n1,n2信号的零插值和低通滤波y1interpy1,2;y2interp1y2;nminlengthy1,lengthy2;y1y11n;y2y21n;信号的合并yy1y2;信号插值子程序(完成信号的零插值和高通滤波)functionfinterp1x实现信号的零插值finterpx,2;fori1lengthxf2*i0;enda1fir1200,1/2,high;ffiltera1,1,f;functiony1_yupre_processy1Fs8000;fpl300;fpu3400;fsl200;fsu3600;wp[2*fpl/Fs,2*fpu/Fs];ws[2*fsl/Fs,2*fsu/Fs];rp1;rs40;[N,wpo]ellipordwp,ws,rp,rs;[B,A]ellipN,rp,rs,wpo;[H,w]freqzB,A;y1_yufilterB,A,y1;(2)非对称主程序clearall;N4096*8;采集点数[y,fs,bits]wavread1.wav,[20000,20000N-1];读取音频信息(双声道,16位,频率44100Hz)y1y,1;y2y,2;t01/fsN-1/fs;timefigure1;subplot2,2,1;plott,y1;1通道信号title原始信号时域波形;xlabelt,ylabel幅度预处理人耳可识别范围(3003400)y1_yupre_processy1;subplot2,2,2;plott,y1_yu;title预滤波后的信号时域波形;xlabelt,ylabel幅度subplot2,2,3;y11ffty1,fs;plotabsy11;title原始信号频谱;subplot2,2,4;y1_yu1ffty1_yu;plotabsy1_yu1;title滤波后信号频谱;[f1,f2]fenjiey1_yu,1,10;1次分解[f1_1,f1_2]fenjief1,1,10;2次分解[f1_11,f1_12]fenjief1_1,1,10;3次分解[f1_111,f1_112]fenjief1_11,1,10;4次分解[f1_1111,f1_1112]fenjief1_111,1,10;5次分解[f1_11111,f1_11112]fenjief1_1111,1,10;6次分解[f1_111111,f1_111112]fenjief1_11111,1,10;7次分解figure2subplot4,2,1nlengthf2t02/fs2.*n-1/fsplott,f2;title一次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,2nlengthf1_2t04/fs4.*n-1/fsplott,f1_2;title二次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,3nlengthf1_12t08/fs8.*n-1/fsplott,f1_12;title三次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,4nlengthf1_112t016/fs16.*n-1/fsplott,f1_112;title四次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,5nlengthf1_1112t032/fs32.*n-1/fsplott,f1_1112;title五次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,6nlengthf1_11112t064/fs64.*n-1/fsplott,f1_11112;title六次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,7nlengthf1_111112t0128/fs128.*n-1/fsplott,f1_111112;title七次分解并抽取的高频部分;xlabelt,ylabel幅度subplot4,2,8nlengthf1_111111t0128/fs128.*n-1/fsplott,f1_111111title七次分解并抽取的低频部分;xlabelt,ylabel幅度figure3f1111fftf1_111111,fs;subplot4,2,1plotabsf1111title7级分解并抽取后低频部分频谱;xlabelf,ylabel幅度f1112fftf1_111112,fs;subplot4,2,2plotabsf1112title7级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1113fftf1_11112,fs;subplot4,2,3plotabsf1113title6级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1114fftf1_1112,fs;subplot4,2,4plotabsf1114title5级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1115fftf1_112,fs;subplot4,2,5plotabsf1115title4级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1116fftf1_12,fs;subplot4,2,6plotabsf1116title3级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1117fftf1_2,fs;subplot4,2,7plotabsf1117title2级分解并抽取后高频部分频谱;xlabelf,ylabel幅度f1118fftf2,fs;subplot4,2,8plotabsf1118title1级分解并抽取后高频部分频谱;xlabelf,ylabel幅度g1hebingf1_111111,f1_111112;1次合并g2hebingg1,f1_11112;2次合并g3hebingg2,f1_1112;3次合并g4hebingg3,f1_112;4次合并g5hebingg4,f1_12;5次合并g6hebingg5,f1_2;6次合并g7hebingg6,f2;7次合并figure4subplot4,1,1;t01/fsN-1/fs;plott,y1_yu;title预滤波后的信号;xlabelt,ylabel幅度subplot4,1,3;nlengthg7;t01/fsn-1/fs;plott,g7;title信宿端恢复信号;xlabelt,ylabel幅度subplot4,1,4g71fftg7,fsplotabsg71;title信宿恢复信号频谱;xlabelf,ylabel幅度subplot4,1,2y1_yu1ffty1_yu;plotabsy1_yu1;title预滤波后信号频谱;xlabelf,ylabel幅度将处理后的信号写入文件wavwriteg7,fs,2.wav;子程序树形结构分解子程序完成信号的2通道分解以及降采样率function[y1,y2]fenjiex,n1,n2,fpfs44100;t201/fslengthx-1/fs;a1[11];b1[1];a2[1-1];b2[1];a1fir1100,1/2;a2qmfa1;w1filtera1,1,x;低通滤波w2filtera2,1,x;高通滤波y1downsamplew1,2;抽取2y2downsamplew2,2;抽取2functionyhebingy1,y2,n1,n2信号的零插值和低通滤波y1interpy1,2;信号的零插值和高通滤波y2interp1y2;进行插值后的信号可能长度不一样,因此要进行长信号的截取nminlengthy1,lengthy2;y1y11n;y2y21n;信号的合并yy1y2;信号插值子程序(完成信号的零插值和高通滤波)functionfinterp1x实现信号的零插值finterpx,2;fori1lengthxf2*i0;enda1fir1200,1/2,high;ffiltera1,1,f;functiony1_yupre_processy1Fs8000;fpl300;fpu3400;fsl200;fsu3600;wp[2*fpl/Fs,2*fpu/Fs];ws[2*fsl/Fs,2*fsu/Fs];rp1;rs40;[N,wpo]ellipordwp,ws,rp,rs;[B,A]ellipN,rp,rs,wpo;[H,w]freqzB,

注意事项

本文(树形结构滤波器组设计)为本站会员(星星008)主动上传,沃文网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知沃文网(发送邮件至2622162128@qq.com或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服点击这里,给沃文网发消息,QQ:2622162128 - 联系我们

版权声明:以上文章中所选用的图片及文字来源于网络以及用户投稿,由于未联系到知识产权人或未发现有关知识产权的登记,如有知识产权人并不愿意我们使用,如有侵权请立即联系:2622162128@qq.com ,我们立即下架或删除。

Copyright© 2017-2019 www.wodocx.com ,All Rights Reserved |陕ICP备19002583号  

陕公网安备 61072602000132号     违法和不良信息举报:0916-4228922