LPCmatlab的编程实现.doc

《LPCmatlab的编程实现.doc》由会员分享，可在线阅读，更多相关《LPCmatlab的编程实现.doc（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、摘 要线性预测编码作为一种基于全极点模型假定和均方预测误差最小准则下的波形逼近技术，是语音编码中的一种重要技术，运用这种编码方式通过对音频信号的时域和频域分析，线性预测编解码，滤波等处理，能够在尽量保持原音质的基础上合成出高质量的语音。本文通过对波形编码，参量编码，混合编码，中低速率编码的介绍，重点研究了线性预测（LPC）的原理、Levinson-Durbin算法及MATLAB在处理语音信号时的应用。最后以MATLAB为开发工具，编写了线性预测编码的程序，对语音信号进行了线性预测编码处理，并对阶数的选取做了比较分析。关键词 语音信号 LPC 莱文森-杜宾算法 MATLAB仿真AbstractL

2、PC all-pole model as based on assumptions and the smallest mean squared prediction error criterion of waveform approximation technique，is an important speech coding techniques，By using this encoding audio signals in time domain and frequency domain analysis, linear prediction coding and decoding, fi

3、ltering processing, able to keep the original sound quality on the basis of high-quality voice synthesis.In this paper，based on the introduction of the waveform coding，parametric coding, hybrid coding, low rate coding, focuses on the linear prediction (LPC) principle，Levinson-Durbin algorithm and MA

4、TLAB in the application of speech signal processing. Finally, with MATLAB as the development tool，write a linear predictive coding procedures，voice signal processing of linear predictive coding.,and the number of orders to do a comparative analysis.Keywords voice signal LPC Levinson - Durbin algorit

5、hm MATLAB simulation目 录第1章 绪论11.1 语音信号LPC分析技术的基本概念11.2 语音信号LPC分析技术的发展21.3 语音信号LPC分析技术的应用21.4 基本设计要求及设计思路3第2章 线性预测编码的基本原理42.1 语音信号的数字模型42.1.1 激励模型42.1.2 辐射模型52.1.3 声道模型52.2 线性预测的概念与原理62.2.1 线性预测分析的概念62.2.2 自相关法72.2.3 协方差法82.2.4 全极点声道模型82.3 解线性预测参数方程组的算法92.3.1 利用乔里斯基分解计算协方差法的方程组92.3.2 利用杜宾递推算法计算自相关法的方

6、程组112.3.3 利用格型法求解线性预测系数122.4 各种算法的比较16第3章 编程实现LPC分析193.1 MATLAB的简介193.2 基于MATLAB的编程实现203.2.1 程序设计流程图213.2.2 主函数及其功能213.2.3 程序分析223.2.4 实验结果与分析243.3 基于C的编程实现24结束语26谢 辞27参考文献28附 录295 = ) ( = * _ _ ( _ _ ; - +; _ ;=( = -; = - - = - = = 0 =0 = ) * ( .= ) 0 ( .0 0 - = )! 0! ; + & _ = -= - = 0 = ; 0 , , ,

7、 , , , , * _ ) _ 00 0 (= =) _( ; _ * * 00 - 0= 0 .0 + ( = .= - ; * _ _ ; = ; * _;_ _ , _ /) _ _ . . * */ * . 0: ( *0 _ _. = ) ( = * _ ;)_ ( ) _ , _ ( ) ) _ - (, , _ _- ) 0 ( = .= ) - ) * ;) ( 0 = 0 ; , _ ) , ( _ -* - = + ; ; & _ = - - - - - - _ -= ;0 )+( _ = ) + _ _ = ; + * * _ 00 - 00 00 00 ; + + .

8、 ) 0 ) *)- *- (= . - ) (-) ( ) 0 ) 0 ( ) * 0) ) * = )0- _*0 ( 0- 0 ) _ 00 - ( ; , _ _ ;) _ = ; _ _ * ( * 00 ) = * = * ).(&) = | . ( 0 )0=( ; * * * ) * ( ; * . 0( (* 0 ) . .( * + * -= *0 - . ;+ = = + ; 00 0 ) ( _ ) ( _ )& , _ ) ( /* /* + = + 0 0= 0=0 * *; * / ,* * ( _ = - - = *+ ; ;= ; ;= = + ;+ * (

9、) 0= ( * = - + .0 ( + - 0= 0.= +;= + , / /* ( ) = ) ; 00 - + 0( * , ) ( + ;( 00= - + ; , * , ( _ ) = + ( 00= ) + = * ( _ * . 0 * / . * ) ( * * ) * * * 0 * /0 * 00 * * * * * *0 . ) ( _ , , _ ) _ _ ( . 00 ( * . * : ( . * * * * _ ) _ ) , * ) _ ) _ ) * _ ) ) ) _ _ _ _ _ _ _ _ _ _ .序系预线算 ( ; 据 ( )( ) )(=

10、:=), ) (=:), 系预后束 (差预 ( 能残预; -= ) (. - ),- ), ) ( ( ) ( ( 系射; ) = )+ ),-( -:= 0=)( : 过) *) () )() )(/ ) ( 系测求 )+ .( =- : 00+: )+ 函自阶测;数一 加 ()( ): 00), .-( ) 波滤数入% 序 00 -.,. 学教气设课号信 和 基,刘 0 . . . . ， ， . . . ， .， .00横科景用状发术语，高 社大中，语字姚 版大滨，号信航 0 技信，及术号语戴00大京北技，0 大业南析 语的 静00社业电与号音时 术讯激型预线信语0 社工，处音，0 版科

11、 ， . 0 版工防别语，、炳 ，工，处信语迟、0 版学清理语，郑庆 0，出学程教信青文倍功才样下后步，体从定事一待识刻我。自来有运善耐题对方方维统，情完做会深,业过合结实理正计这，够理知察实经，明只性重于实体识的题对够识模析 号时写初程的不个高和进西很仍题很存中，有间由求了到基体从研 信语实础处信数要别语计束 ) * = + ;( 0.= ) +; 0 , , ( 部实的数) ) (= .= * -0* + +; ;=( = 0= + = + ( = 00 + + , * 系反转系法算 实实 间0好的出可图，趋，度幅，定到 降逐的差预 ，增值着，可图数数性获 法 进对0 测线中理化模合，度形

12、小形端两始出图幅分分验 示图。的义数函 ；图窗第，子窗个为口绘示该值取,函 , 窗绘新制是作 量残器 ; ) =|( ; ( 预法 ;( ; ; )( 下码示图序形示，的|， ；， ，的 实语 中 ( (系系% )(- ( 提的 )(=: ), ) :) 系系 )+ . ( -:=0=) : 下如程推循 ) () ( )(= ) ) ) ) 方数预码序程步为，提测现数系性算 据) . =下下实分号音序于，文 在 信的序数信入分序 图中窗第，子口子 口绘示函值 ) , ：格本基实 用形多显窗一果。 口口新成制（ 阵矩 ） 环环结0,条 判再加值 句则 于否 判,赋 给中语+ ;= ( :。分用间

13、个。变按一循变环义环出么定它表个式表值量制给它句赋是初语)化( 式一它的循 环 据帧取加) =;( 器滤预)( 数入据读) ）作（有所 空 功及主 . 0 环，窗 函0 数0=度列取我程计程 （图 差残系测求法函关相数加波虑重） 系反图(能残系线算宾 关求据读 波实能反 重数函测差图程程列数0 0 功 所）)读数 器 ;)帧 循式化件)是句量式表出义变按用: ;中给赋否 句再条0环 ） 制口 一多用基：, )值示口子， 分数 文于号下= )据性数为序数 ) =) 循程下 )- ) :) )= - ( 语 ， ，示 ; 法( = 量 是绘 ,取绘为窗窗函的示验幅两端度化线0 法 数着，差逐到度，

14、可好 实算反 0 =;+ 0 = ()数部 (, ) = (+ )束别信实 体了间中很西和个写号析对识实明经理够理结业做情方对善有。待从步后倍文程0庆理清 信工，、，版0. 科0， 0型讯 时与社 语 大，京0戴术信0，大姚，社 ，状景横 ， 刘 和课学.-0 %数 )-,0: ;阶 )+ . 测( ) ( )*过:) - ,) ) 系( ,)-) ( 残 差束 ,=) =: 据 线预 _ _ _ ) ) ) ) ) ) _* * . * * * . ( _ ) 0 * *0 * * ) * ./ * _ = 0 + _ , - 0 ; , ( - 0 ( / + + = ( - ) ( ;

15、+ = ; ; - = * , ; 0; 0_ * _ = (); =_ - - - = & ; ;+ -* ; * + - = = 0 , ) _ , )_ = ()= . 0 : . * * . _; = ; - = ( + 0 * _ )= = _) _ , ,0 - ; ( 0 .; - ; ! - ) ( . = ) = - ; ;= =+ = - ; = - = =; 0 *0 + - = ;- = +-; _ (_ (=程 于 性一出将我时慢执到。慢执程，语不而释是 由不类 微方。实代释运然码立设编代 将过立和灵 译译 序程据复操可验无出以序能这 利。户的互建的他程 界形具越据技

16、为 得使。出图脑台会图绘：运 。理图图许 不言画的题的多其络，处，统系理在以箱标买用例题复体决用以，工有，函中本 单加更让，编 经数千样。误标，意括函这的所算和的所写你言的许统有与计你程个个编正，。函问本包试多供提函预极个 函定台新移序写 把要可此译可样平件据写平在运正一平它程写平个持都 的本 00 写本。立平量供系操许 立独例展和，件在试译成个包工用易十 使具的理它原的新和育于使易种是。改修的可，环开 在程好写可值式行在求作序程它一 版像言释是 用包主点诸语技统对言 优的力的域许在扩数多更具而函多0超经已 本的要程工的数函语 本基数富的置拥程庞是 效单变术得样数函泛极了并言 序程 。技要种解

17、系的其种发它年后在算阵的程于机的专一的实( 简的 分 实编多差相择选统系同基不大于 取，范周个限析如型和方而左 择选法自 频取结的重不窗以倍周基等大法相，。于算，种这于为能 ，算。题 长了以本 选多减并预加步当 值求样 ，负辐励来 外应道来系点要需 ( 需共 / 率音， ( 样，是，响示系个具每用一的被有选样和选。关一哪它择的阶定稳，式.（数相，求测因法对稳的协同结所关法方大数样帧当，的不定解测差对件要的足件 （ 满数分，法稳检式利宾。能尽使则平可使，预对果果验。能定么，不数自果出当但稳其，圆单所)，相自内于零 . 要亦稳系圆位点的， .（ 系的是就性 于矩相加数数算 于 / 矩数储法型法关法

18、基差分指方 比法种 . 正 于 正于正解窗）法据据（计 比 0/ 于 窗 阵数一）（型宾法） 于 种比方分法型数/ 加于 就系 ，点稳 零内，单稳出果不定验，平使能式，分 要差解，当数法结对测，.稳定择。选有用个响是 音 (需来外来负 样 当并选 本题，为种算，等基窗的频自 而型限范于基系择多 简一机的算年种解要程 了函得单是程的本 的程本已0具更许域优对诸包是像 序在值好 的改种于的理使 易成在和例 系平立 0的个写一在件可可 写台 极多试函编程与许你的所函，标数 ，加 函工决题买标系，其多的 图 会图使 据形 他的利这以验复序译 立代立运代方 释不慢。时自（斯乔方类 比法种于法的命有数数重

19、为系析性于型格这上窗要器的一产能格此缺误关大比点音了因系求取音接，关算要法束）（ （）重 ）-，, ,） 求）. 差误向正算） （ 求以式.值始下如可步 . 围值。系 数分为参关相预反示它来联测与正 出 . 0. （可， （）. 式）. . 示以小最的测反正是算了出 . 推 . 式换反式 （） . 换 式 . （ 可） （式解从， . 入）. . 示表变 .求的行我 件充稳型格以。值射由可的数数系于的滤格测预接可系为优的格。若之的方协相讨前结型算计射的合一型格是点型格测出条最测用法协关在要个算的的重理语数 .（ 可理进 .式)(号输 入满式形 成因近)由 波逆） 求分性过( 合输) 输波则信输

20、 因中将如结采亦 波即模图型信示图，。)误即亦预为输(为波 形形成则，向示符如信误的器逆预价 0. ，时 ） （）. 构器型用测性推念预向向了在程向称过这误预差误所的样 测 ，) () (值刻延上成看中 （；的 预-、. (、(样过个是测性通的 . . 信信预和差预义号部将理的部两面号后和移经信的) 是部；输滤(部第个为输滤逆阶知， . 系关们差测为输)号信器 数统系滤测称（滤第算数用)， (测预 算如中算 法测性的型于系生性的了增则近最方念的反预了于是。的间稳精决而，稳认同权号数口需型格可讨下式型) 波了先义理参，概测向测预“引的相分时究美 板学，十法型，的法了，以性结数预到，大计相而假不际

21、了于点段语，计的语在数系出程阵矩关步为都，是关是不法方系预解法格 性定的 证此件要圆在)(多，证。围围取，表 . （ 式 . . 假 示则，预器测表差测归成测时变的方示表0 (系相归一系果测的的阶得).过可算差预器阶到）. 出了得数预各低数过数预求，器个于指。测 第预 ，器表的号，达 . 解终解递 以式0.至 .0. . . ：公要个列下推递 （是法效而等 -(逊算 （如多程上解利伯，式 .(数程 数器预式 .解主，于述 程程相自推杜 择的定，化变求，误方求 和矢根示 可 . 写写式. (代. . 可， . . 写可式. 入) .0.（ ，可， = . .（ 示表矩 . 示表误可方式单示测小可

22、解基乔束过解方 到=从推的 .， （ 件 . 式关推 以） . . 件 . 示所，系递量到以）. . （ =写以） . . =0. 写以， （和. 了求可成和为 .（ 式阵得 （ 同 写以 故= 元 . 得得 （利对 （ 示可 故 角 。素 矩中) 从就相应边).(令。示 上矩一，元主阵矩一 （ 式下成分，算。基里以可矩正对因), 分元矩和，, 项 它称对 示式 （ 所式 程方线程方法算基里利.高提算可对的数方利同率的是数个解程从能，许法程参线示式 （ . 方，式 .（ 法自组参性了推算组数测线础础研用，模声很一点的信。是(源，清冲单的音是信的模于以信音音含差数测测线 数型该音称也点号音被 ，

23、。合可的 差测 说就)号音输) ，)器逆波差过使信激作差预如出) ( 数统器滤测，到。出) ，差预器该，输输励音语是看示所 型型好非供提更功用为不处音用析预模道全 法相效使往处音以所因为的推效其突再优较差，况在。准也数小的中音个的大有相时右左 往，处语。多法自比其时差 ，差测现端法在于点最解)缺稳引应量有种当的能一性，接果相自结，够取帧当际，的滤测造有测来出方上在样零需端数里出 程方到则准小平集 .(则方平时定，预性差差 .法算递等尔和型法有主定稳统实能，化测外重更一，况据是。真谱而，零个要的据差测计相是意 ( 式式下以可 ( 得， .00 可，关号随 . 因数关自就 即，方小合式 (则均最时

24、定用若性线关自 值价同着域识语它特不具特不等积对数分征谱数测线例矢派很步一矢征 矢特 ，特帧绘个参每参 若以帧由析 数测预以由最上小差度准测零于误 成示 令 ( 有最测均使则佳 . 为误则系预示用系 . 为号则样取测预之值个用设性阶则预进样用计估模进(信知，差方利模极立建 定变时为种的种三都他因讨法基三只值数确更相比法方，形型的义方 - 置了得术波自，来的义相自波个数） （偏低感对取器性用般析法法然法法法法器法法方有方这算不成因不的所体具在这析预线多许已，法方而精种供测，参统线计了语或音产，性一励噪冲脉用可的的模音建原的 (测线数加合性，合性加去取个用，在取音概的预 .原与测性 （ ( ) 数

25、统统系型模为点而有 数系，事同母，略延把，与中. （ .（ .（ 示示系统数果外变近 ) ） 帧一示 数组) 函组以 . 示） ）.（函，来管 以产模道 小得值该于；0/，于常函自 号是(这子加预。通阶一，之在”技加用时分际，高模射低型角，之(型脉除由波类阶射以部抗射耗损起度播声，的开口这， . 下如公射出，况射辐成似则面头远的唇认，但困有这上理当应绕头的也效的口。射称倒，声号语，的入模射 （ ：示型励个 . 式形的下可因及脉单果型角三上串脉了士波角因型的一为形角 （ 式的型全成它将，间时下为，分升角 .（ 下如学的角单脉三的周以是激时冲角三它。歇生闭开断声由时模激 .联的射型声模：子为可模号

26、信模数信 求的适条和参地周化的语用使音语语，使预是时波析使合又使测既滤线时模生语它波滤时间时个是器。音器语通参析 使收少大数声就编参些化就编波，参函和激生音语是原本编性 法法为的导性的分以度定只析 ) 对。音窗形相 数用 音个取以时短为时短须相号语和的测线的语此的为才时一在信到。 的进参 全 (则均并，则最某型模，线一)(系是的重是也价计可因解运卷用法等识码语现。( 求 需域各涉数音。积 信和 响样号 语在输产激)信到 系因性是)号语 产号思计求本 器预编性 到 用中频 在编线 了中玩教 & 0”好评质的 。音作强的 为速） 置话代国来年 缺有度力噪，语可成以个利的是编取数每音成量得为单算采，

27、码性0 。 准政联，方算， 基信音上 .在 定美年行的有音电，励的时得预歌从，成建于可合输道语窄后密、进必中线全安式信线作它种使准标 如压的司电它成音语通编应术析 音 法 第开. 年 在最度 秒 之 - 使 议 第 实的/ 间 - 包 向一第 。 时 位0个第现 载历正 据术 终这开第息了开 ( ), 与 年个四卖 的 了示 器语 示上协国美 .，压行进 的 : 分出 与 年 。的大略 。计似的编针使一种法别素自一 时年 起编性法 . 福等定声和、通如使泛中子代 码展的析 信音法差的 以，区制差测算是，语案种法自称种时的保就长定， 进对然列音将种方有。参估来用能，稳短种一一重识为以的得 计精数

28、语等、共还分，。的频波述地参时息低以行道析 利可，生系时线过源型语如模道非供能，测预音得不分预数测预的定可最差均值抽线样语过通来线样个去可样号个：本测性音的了得理语识话别、合码语一技音效最析别的相会波号原同信但确准可有语重，编的征信对是编说。信重参根量特换列数到程其为。进数变将参特域交其率号信码声称数得度的比单比以而的其波现确正效参，音估确极能性其术种是性方)储传速以，积谱振音基数音础术分要一然目种，多有音今产技性技应实真第，术音各一技析行流有法这方各处于应普，工作中成分到用技性先最仓 域于测预后出首年 维术 测码码线 力信通高,种另其定进音对 段最动会想交人式形的最用有、最是，种两图只信最的最以是信外到从也觉的固等、样样是类对食件条在类为经息通交的量要刻每机人，人与。特息出来等内话识话出、成人，求用到以的处者储或效便音些要就理科一处号信技处数研是概基术分号语 绪 00 录附 _ _ 献参 _ 谢 语结 _ 实于 . _ 析分实 . 0 析程 . 功主 _ _ 流程 现程的 . _ 介的 . _ 析 章 _ _ 比算 . _ 系解型用 . 方关计法杜利 的协算斯乔 0 _ 组数性解 _ 模声 差 _ _ 关自. 念分性 . _ _ 原概预 . _ _ 模 模射. _ 0 _ 模数信 0 _ 的