function?speechproc（）

????%?定義常數(shù)
????Fs=8000;
????FL?=?80;????????????????%?幀長(zhǎng)
????WL?=?240;???????????????%?窗長(zhǎng)
????P?=?10;?????????????????%?預(yù)測(cè)系數(shù)個(gè)數(shù)
????s?=?readspeech（‘voice.pcm‘100000）;?????????????%?載入語(yǔ)音s
????L?=?length（s）;??????????%?讀入語(yǔ)音長(zhǎng)度
????FN?=?floor（L/FL）-2;?????%?計(jì)算幀數(shù)
????%?預(yù)測(cè)和重建濾波器
????exc?=?zeros（L1）;???????%?激勵(lì)信號(hào)（預(yù)測(cè)誤差）
????zi_pre?=?zeros（P1）;????%?預(yù)測(cè)濾波器的狀態(tài)
????s_rec?=?zeros（L1）;?????%?重建語(yǔ)音
????zi_rec?=?zeros（P1）;
????%?合成濾波器
????exc_syn?=?zeros（L1）;???%?合成的激勵(lì)信號(hào)（脈沖串）
????s_syn?=?zeros（L1）;?????%?合成語(yǔ)音
????%?變調(diào)不變速濾波器
????exc_syn_t?=?zeros（L1）;???%?合成的激勵(lì)信號(hào)（脈沖串）
????s_syn_t?=?zeros（L1）;?????%?合成語(yǔ)音
????%?變速不變調(diào)濾波器（假設(shè)速度減慢一倍）
????exc_syn_v?=?zeros（2*L1）;???%?合成的激勵(lì)信號(hào)（脈沖串）
????s_syn_v?=?zeros（2*L1）;?????%?合成語(yǔ)音

????hw?=?hamming（WL）;???????%?漢明窗
????
????%%初始化各濾波器狀態(tài)
????zi_pre=s（2*FL-P+1:2*FL）;????????
????zi_rec=s（2*FL-P+1:2*FL）;
????zi_syn=s（2*FL-P+1:2*FL）;
????zi_syn_v=s（2*FL-P+1:2*FL）;
????zi_syn_t=s（2*FL-P+1:2*FL）;
????
????%?依次處理每幀語(yǔ)音
????for?n?=?3:FN

????????%?計(jì)算預(yù)測(cè)系數(shù)（不需要掌握）
????????s_w?=?s（n*FL-WL+1:n*FL）.*hw;????%漢明窗加權(quán)后的語(yǔ)音
????????[A?E]?=?lpc（s_w?P）;????????????%用線性預(yù)測(cè)法計(jì)算P個(gè)預(yù)測(cè)系數(shù)
????????????????????????????????????????%?A是預(yù)測(cè)系數(shù)，E會(huì)被用來(lái)計(jì)算合成激勵(lì)的能量

????????if?n?==?27
????????%?（3）?在此位置寫(xiě)程序，觀察預(yù)測(cè)系統(tǒng)的零極點(diǎn)圖
????????figure;?zplane（1?A）;
??
????????end
????????
????????s_f?=?s（（n-1）*FL+1:n*FL）;???????%?本幀語(yǔ)音，下面就要對(duì)它做處理

????????%?（4）?在此位置寫(xiě)程序，用filter函數(shù)s_f計(jì)算激勵(lì)，注意保持濾波器狀態(tài)
????????%?exc（（n-1）*FL+1:n*FL）?=?...?將你計(jì)算得到的激勵(lì)寫(xiě)在這里
????????[?exc（（n-1）*FL+1:n*FL）?zi_pre]?=?filter（A1s_fzi_pre）;
????????
????????%?（5）?在此位置寫(xiě)程序，用filter函數(shù)和exc重建語(yǔ)音，注意保持濾波器狀態(tài)
????????[s_rec（（n-1）*FL+1:n*FL）zi_rec]=filter（[1]Aexc（（n-1）*FL+1:n*FL）zi_rec）;
????????
????????%?s_rec（（n-1）*FL+1:n*FL）?=?...?將你計(jì)算得到的重建語(yǔ)音寫(xiě)在這里

????????%?注意下面只有在得到exc后才會(huì)計(jì)算正確
????????s_Pitch?=?exc（n*FL-222:n*FL）;
????????PT?=?findpitch（s_Pitch）;????%?計(jì)算基音周期PT（不要求掌握）
????????G?=?sqrt（E*PT）;???????????%?計(jì)算合成激勵(lì)的能量G（不要求掌握）

????????%?（10）?在此位置寫(xiě)程序，生成合成激勵(lì)，并用激勵(lì)和filter函數(shù)產(chǎn)生合成語(yǔ)音????????
????????%?exc_syn（（n-1）*FL+1:n*FL）?=?...?將你計(jì)算得到的合成激勵(lì)寫(xiě)在這里
????????%?s_syn（（n-1）*FL+1:n*FL）?=?...???將你計(jì)算得到的合成語(yǔ)音寫(xiě)在這里
????????if?（n==3）???????%初始化合成位置
????????????pos=（n-1）*FL+1-PT;??
????????end
????????pos=pos+PT;??%本段第一個(gè)激勵(lì)與上一段最后一個(gè)相差PT
????????while?（pos????????????exc_syn（pos）=G;
????????????pos=pos+PT;
????????end
????????pos=pos-PT;??%本段最后一個(gè)激勵(lì)位置
????????[?s_syn（（n-1）*FL+1?:?n*FL）zi_syn]?=?filter（1AG*?exc_syn（（n-1）*FL+1:n*FL）zi_syn）;??

????????%?（11）?不改變基音周期和預(yù)測(cè)系數(shù)，將合成激勵(lì)的長(zhǎng)度增加一倍，再作為filter
????????%?的輸入得到新的合成語(yǔ)音，聽(tīng)一聽(tīng)是不是速度變慢了，但音調(diào)沒(méi)有變。
????????%?exc_syn_v（（n-1）*FL_v+1:n*FL_v）?=?...?將你計(jì)算得到的加長(zhǎng)合成激勵(lì)寫(xiě)在這里
????????%?s_syn_v（（n-1）*FL_v+1:n*FL_v）?=?...???將你計(jì)算得到的加長(zhǎng)合成語(yǔ)音寫(xiě)在這里
????????FL_v=FL*2;??????%幀長(zhǎng)加倍，合成類(lèi)似上一問(wèn)
????????if?（n==3）????????%初始化合成位置
????????????pos_v=（n-1）*FL_v+1-PT;
????????end
????????pos_v=pos_v+PT;
????????while?（pos_v????????????exc_syn_v（pos_v）=G;
????????????pos_v=pos_v+PT;
????????end
????????pos_v=pos_v-PT;
????????[s_syn_v（（n-1）*FL_v+1:n*FL_v）zi_syn_v]=filter（[1]Aexc_syn_v（（n-1）*FL_v+1:n*FL_v）zi_syn

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????7901??2009-09-20?15:32??語(yǔ)音合成.m

-----------?---------??----------?-----??----

?????????????????7901????????????????????1