【正文】 由 DFT變換，任一 X（ m）的計(jì)算都需要 N次復(fù)數(shù)乘法和 N1次復(fù)數(shù)加法，而一次復(fù)數(shù)乘法等于四次實(shí)數(shù)乘法和兩次實(shí)數(shù)加法，一次復(fù)數(shù)加法等于兩次實(shí)數(shù)加法，即使一次復(fù)數(shù)乘法和一次復(fù)數(shù)加法定義成一次“運(yùn)算”（四次實(shí)數(shù)乘法和四次實(shí)數(shù)加法），那么求出 N項(xiàng)復(fù)數(shù)序列的 X（ m），即 N點(diǎn) DFT變換大約就需要 N^2次運(yùn)算。最后 IFFT處理觀察恢復(fù)后的時(shí)域信號(hào)圖形，頻譜分析過(guò)程。因此，語(yǔ)音信號(hào)是人們構(gòu)成思想疏通和感情交流的最主要的途徑。 語(yǔ)音合成的目的就是讓計(jì)算機(jī)說(shuō)話。如同在時(shí)域特征分析中用到的一樣，這里的傅里葉頻譜分析也采用相同的短時(shí)分析技術(shù)。繼續(xù)上面的例子， N=1024時(shí)，總的運(yùn)算次數(shù)就變成了 525312次，節(jié)省了大約 50%的運(yùn)算 量。下面分別介紹 MATLAB在音量標(biāo)準(zhǔn)化、語(yǔ)音合成等音頻信號(hào)處理方面的技術(shù)實(shí)現(xiàn)。 halflength=floor(sigLength/2)。 plot(f,abs(Y(1:halflength)))。 t1=(0:length(y1)1)/fs2。 figure(4)。 end end I=i。)。 plot(t,x)。 title(39。由于之前運(yùn)用過(guò) MATLAB軟件， MATLAB中的許多函數(shù)，如 fft， plot， subplot，abs等用起來(lái)還是有點(diǎn)感覺(jué)的。 pause。 figure(6)。 X=fft(x,sigLength)。 N=N39。)。 xlabel(‘time(s)’)。 %錄入男生語(yǔ)音信號(hào) [y1,fs2]=wavread(39。 ylabel(39。 %分析女生信號(hào) y=y(:,1)。 性質(zhì)二：序列 x（ n）與其 Hilbert變換是正交的。它對(duì)傅氏變換的理論并沒(méi)有新的發(fā)現(xiàn)，但是對(duì)于計(jì)算機(jī)系統(tǒng)或者說(shuō)數(shù)字系統(tǒng)中應(yīng)用離散傅里葉變換，可以說(shuō)是進(jìn)了一大步 。在FFT處理之后，普通頻譜分析可以進(jìn)行頻域上的濾波處理，從而使頻譜更加平滑。語(yǔ)言是人類特有的功能，聲音是人類常用的工具，是相互傳遞信息的最主要的手段。最簡(jiǎn)單的語(yǔ)音合成應(yīng)當(dāng)是語(yǔ)音響應(yīng)系統(tǒng)，其實(shí)現(xiàn)技術(shù)非常簡(jiǎn)單。 信號(hào) x（ n）的短時(shí)傅里葉變換定義為 : 式中， 為窗口函數(shù) 。而如果我們將這種“一分為二”的思想不斷進(jìn)行下去，直到分成兩兩一組的 DFT運(yùn)算單元，那么 N點(diǎn)的 DFT變換就只需要 Nlog2N次的運(yùn)算， N在 1024點(diǎn)時(shí)，運(yùn)算量?jī)H有 10240次，是先前的直接算法的 1%，點(diǎn)數(shù)越多，運(yùn)算量的節(jié)約就越大，這就是 FFT的優(yōu)越性。 該實(shí)驗(yàn)以兩種男女聲音樂(lè)為分析樣本。 f=fs1*(0:(halflength1))/sigLength。 grid on。 sigLength1=length(y1)。 plot(fl,abs(Y1(1:halflength1)))。 %找到男生信號(hào)開始的地方 for j=1:sigLength14, a(j)=y(j).*y(j+1).*y(j+2).*y(j+3).*y(j+4)。 fprintf(39。 grid on。合成信號(hào) FFT頻譜 39。通過(guò)老師的講解以及我們課下對(duì) MATLAB的學(xué)習(xí)，我們基本可以運(yùn)用 MATLA