%?ct1
clear?allclose?all?clc;

????%?定義常數(shù)
????FL?=?80;????????????????%?幀長(zhǎng)
????WL?=?240;???????????????%?窗長(zhǎng)
????P?=?10;?????????????????%?預(yù)測(cè)系數(shù)個(gè)數(shù)
????[sfs]?=?wavread（‘sunday.wav‘）;?????????????%?載入語(yǔ)音s
	s?=?s/max（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ǔ)音
	last_syn?=?0;???%存儲(chǔ)上一個(gè)（或多個(gè)）段的最后一個(gè)脈沖的下標(biāo)
	zi_syn?=?zeros（P1）;???%?合成濾波器的狀態(tài)
????
	%?變調(diào)不變速濾波器
????exc_syn_t?=?zeros（L1）;???%?合成的激勵(lì)信號(hào)（脈沖串）
????s_syn_t?=?zeros（L1）;?????%?合成語(yǔ)音
	last_syn_t?=?0;???%存儲(chǔ)上一個(gè)（或多個(gè)）段的最后一個(gè)脈沖的下標(biāo)
	zi_syn_t?=?zeros（P1）;???%?合成濾波器的狀態(tài)
????
	%?變速不變調(diào)濾波器（假設(shè)速度減慢一倍）
	v=.5;
????exc_syn_v?=?zeros（v\L1）;???%?合成的激勵(lì)信號(hào)（脈沖串）
????s_syn_v?=?zeros（v\L1）;?????%?合成語(yǔ)音
	last_syn_v?=?0;???%存儲(chǔ)上一個(gè)（或多個(gè)）段的最后一個(gè)脈沖的下標(biāo)
	zi_syn_v?=?zeros（P1）;???%?合成濾波器的狀態(tài)
	
????hw?=?hamming（WL）;???????%?漢明窗
????
????%?依次處理每幀語(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）?觀察預(yù)測(cè)系統(tǒng)的零極點(diǎn)圖
????????????zplane（1A）;
????????end
????????
????????s_f?=?s（（n-1）*FL+1:n*FL）;???????%?本幀語(yǔ)音，下面就要對(duì)它做處理

????????%?（4）?用filter函數(shù)s_f計(jì)算激勵(lì)，注意保持濾波器狀態(tài)
		[exc1zi_pre]?=?filter（A1s_fzi_pre）;%最終值z(mì)i_pre保留下來(lái)提供給下一組數(shù)據(jù)運(yùn)算，就能維持輸出數(shù)據(jù)的連續(xù)性。exc1為激勵(lì)
????????
????????exc（（n-1）*FL+1:n*FL）?=?exc1;?%計(jì)算得到的激勵(lì)

????????%?（5）?用filter函數(shù)和exc重建語(yǔ)音，注意保持濾波器狀態(tài)
		[s_rec1zi_rec]?=?filter（1Aexc1zi_rec）;
????????
????????s_rec（（n-1）*FL+1:n*FL）?=?s_rec1;?%計(jì)算得到的重建語(yǔ)音

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

????????		
		%方法3：本段激勵(lì)只能修改本段長(zhǎng)度
		tempn_syn?=?[1:n*FL-last_syn]‘;
		exc_syn1?=?zeros（length（tempn_syn）1）;
		exc_syn1（mod（tempn_synPT）==0）?=?G;?%某一段算出的脈沖
		exc_syn1?=?exc_syn1（（n-1）*FL-last_syn+1:n*FL-last_syn）;
		[s_syn1zi_syn]?=?filter（1Aexc_syn1zi_syn）;
		exc_syn（（n-1）*FL+1:n*FL）?=??exc_syn1;???%計(jì)算得到的合成激勵(lì)
		s_syn（（n-1）*FL+1:n*FL）?=?s_syn1;???%計(jì)算得到的合成語(yǔ)音
		last_syn?=?last_syn+PT*floor（（n*FL-last_syn）/PT）;
?		
????????%?（11）?不改變基音周期和預(yù)測(cè)系數(shù)，將合成激勵(lì)的長(zhǎng)度增加一倍，再作為filter
????????%?的輸入得到新的合成語(yǔ)音，聽(tīng)一聽(tīng)是不是速度變慢了，但音調(diào)沒(méi)有變。
		FL_v?=?floor（FL/v）;
		tempn_syn_v?=?[1:n*FL_v-last_syn_v]‘;
		exc_syn1_v?=?zeros（length（tempn_syn_v）1）;
		exc_syn1_v（mod（tempn_syn_vPT）==0）?=?G;?%某一段算出的脈沖
		exc_syn1_v?=?exc_syn1_v（（n-1）*FL_v-last_syn_v+1:n*FL_v-last_syn_v）;
		[s_syn1_vzi_syn_v]?=?filter（1Aexc_syn1_vzi_syn_v）;		
?????	last_syn_v?=?last_syn_v+PT*floor（（n*FL_v-last_syn_v）/PT）;???
????????exc_syn_v（（n-1）*FL_v+1:n*FL_v）?=exc_syn1_v;??%計(jì)算得到的