【正文】 00hz y3_2=[y3_1 y3_1]。 subplot(2,2,1)。 %x軸 ylabel(39。 Ts=fix(ns/k)。 %頻譜以 0頻為中心 ,保證所有頻譜的坐標(biāo)值一樣 TF(i,:)=temp。 %量化間隔， det ylabel(39。)。 %縱、橫坐標(biāo)軸采用等長刻度 set the axes aspect ratio grid。 y_original=wavread(39。Please input a number to amplify or decrease the volume,number:39。Choose the mode 1Canyon,2Lobby,3Little Chamber ,4Outer Space:39。Choose the mode 1Canyon,2Lobby,3Little Chamber ,4Outer Space:39。Choose the mode 1from right to left,2Lobbyfrom right to left:39。 mode=input(39。淡入淡出效果選擇 \n39。 sound(y6,Fs)。 %%%%%%%%%%%%%%淡入淡出 %%%%%%%%%%%%%%%%%%%% %實現(xiàn)：聲音由小變大，由大變小，模擬音源由遠(yuǎn)到近，由近到遠(yuǎn)的變化。相位變化效果選擇 \n39。)。)。)。)。)。)。%繪等高線圖 title(39。時間 39。 temp=fft(xw,nfft)。 %連續(xù)段之間的重疊長度 L=512。時間 39。 time2=linspace(0,9,length(sm1))。 t3=linspace(0,9,n7)。%失真信號回放 %wavplay(y,6000)。時間 39。 temp=fft(xw,nfft)。 %連續(xù)段之間的重疊長度 L=512。時間 39。 %延時后信號回放 nz=size(z)。)。 %FFT變換 temp=fftshift(temp)。 %每段長為 512 k=NWd。)。 %按實際時間顯示輸出 figure。)。 fprintf(39。放音順序： \n39。頻率 39。 contour(abs(TF))。 xlabel(39。 xw=sm1(n1:n2)。 d=100。 xlabel(39。 %搬移到較高頻段后信號回放 ns=size(sm1)。fs=%g\n39。 %采樣失真時語音回放 聲音失真 %wavplay(y,10000)。)。 n2=(Ld)*(i1)+NW。 NW=512。)。y]。 xlabel(39。 xw=y(n1:n2)。 d=100。 xlabel(39。 %原始信號回放 ny=size(y)。)。Press any key to play the recording...\n39。 fprintf(39。 最后，特別的向親愛的家人表示深深的謝意，他們給予我的愛、理解、關(guān)心和支持是我不斷前進(jìn)的動力。 (三 ) 勇于嘗試 系統(tǒng)即使經(jīng)過了非常嚴(yán)格的論證也仍然可能存在問題，或許面對問題一時沒有很好的解決方案，但是有一個或許可行的大膽的想法，不要猶豫，試一試吧。 else am = linspace(1,0,Len)39。 y_d1 = [y_2。 y_2 = y_original(:,2)。 j = linspace(0,1,length(y_stereo(:,1)))39。y_2]。 y_y2 = [y_1。 y_delay1 = [zeros(delay_bits,1)。delay_time = 。/Fs。 y_delay1 = [zeros(floor(Fs*delay_time),1)。 % 峽谷效果 elseif mode ==2 scale = 。)。 ylabel(39。 subplot(2,2,3)。)。 n2=(Ld)*(i1)+NW。 基于短時傅里葉變換的頻譜分析具體程序如下： NW=512。 t3=linspace(0,1,n7)。 %對語音信號進(jìn)行延時 （四） 頻譜搬移的設(shè)計 語音信號的頻率一般比較低，把頻率搬移到較高頻段是實現(xiàn)變聲常采用的一種方法。Press any key to play the recording...\n39。 fprintf(39。 五 設(shè)計實現(xiàn) （一） 系統(tǒng)設(shè)計流程圖 //miao shu 圖 2 系統(tǒng)設(shè)計流程圖 （二） 實時錄放音實現(xiàn)方法 MATLAB 擁有非常強(qiáng)大的函數(shù)庫，直接調(diào)用函數(shù) wavrecord 即可實現(xiàn)實時錄音。回響的效果需要的參數(shù)與回聲的相同，見表 2。 回響 多重的回聲，模擬產(chǎn)生出不同的聽音環(huán)境（峽谷，大廳，小廳，太空）。這是因為窗口函數(shù) w(n)通常具有低通頻率響應(yīng)，而指數(shù) jw? 對語音信號 x(n)有調(diào)制的作用，使頻譜產(chǎn)生移位，即將 x(n)頻譜中對應(yīng)于頻率的分量平移到零頻。其定義如下： ????? ??? m jw mjwn emnwmxeX )()()( (2) 由公式 (2)可以看到 ,短時傅立葉變換所選取的窗是音信號的標(biāo)準(zhǔn)傅立葉變換。 將短時分析應(yīng)用于 頻域分析即是傅立葉變換 , 相應(yīng)的頻譜稱為“ 短時譜” , 即有限長度的傅立葉變換。近幾年來，有許多科學(xué)家、 數(shù)學(xué)家、 工程師等開發(fā)了一些新的、 有價值的應(yīng)用程序， 這些應(yīng)用程序都可以被納入 MATLAB 工具箱。靈活的 MATLAB 語言可使工程師和科學(xué)家簡練地表達(dá)他們的思想，其強(qiáng)有力的數(shù)值計算方法和圖形便于測試和探索新的思想，而集成的計算環(huán)境便于產(chǎn)生快速的實時結(jié)果。 第二種方法是直接乘以一個余弦函數(shù)，把語音信號頻段搬移到較高頻段，以達(dá)到變聲的目的。 1. 窗的選擇 一般來講，一個好的窗函數(shù)的標(biāo)準(zhǔn)是：在時域因為是語音波形乘以窗函數(shù)，所以要 減小時間窗兩端的坡度，使窗口邊緣兩端不引起急劇變化而平滑過渡到零，這樣可 以使截出的語音波形緩慢降為零，減小語音幀的階段效應(yīng)；在頻域要有較寬的 3dB 帶寬以 及較小的邊帶最大值。 語音的個性化特征可以分為以下三類： ● 音段特征：描述的是語音的音色特征。 文章介紹了在 Matlab 環(huán)境中如何采集語音信號和語音信號采集后的 頻譜分析處理， 并通過實例分析了語音信號處理的Matlab。 第二章介紹了語音信號的特點、分析的基本原理及相關(guān) Matlab 的基礎(chǔ)知識。 一 緒 論 (一 )語音信號處理的研究意義 當(dāng)今，數(shù)字信號處理 (DSP： Digtal Signal Processing)技術(shù)正飛速發(fā)展，它不但自成一門學(xué)科，更是以不同形式影響和滲透到其他學(xué)科，它與國民經(jīng)濟(jì)息息相關(guān)，與國防建設(shè)緊密相連，它影響或改變著我們的生產(chǎn)、生活方式，因此受到人們普遍的關(guān)注。 同時由于 Matlab 是一個數(shù)據(jù)分析和處理功能十分強(qiáng)大的工程實用軟件， 它的信號處理與分析工具箱為語音信號分析提供了十分豐富的功能函數(shù) ,利用這 些功能函數(shù)可以快捷而又方便地完成語音信號的處理和分析以及信號的可視化， 使人機(jī)交互更加便捷。 關(guān)鍵詞：話音變換，音頻特效， Matlab， 短時傅里葉變換 Abstract The speech signal processing technology is in the field of speech processing newly developed a subject branch. Matlab is a data analysis and processing functions very strong engineering practical software, the use of it to carry on the speech signal collection, analysis and processing is quite this paper, the voice of a variety of methods and a variety of transform audio effects method of realization of the discussions, and gives the Matlab language description. Comprehensive voice transform and some of the audio special effects technology, design a widely suitable for soundtrack, voice security and entertainment the fields of realtime voice changer, and in the shorttime Fourier analysis based on the speech signal spectrum. Keywords: Voice transformation Special sound efficiency Matlab STFT 前 言 (一 )課題介紹 隨著數(shù)字化時代的來臨，科學(xué)技術(shù)的進(jìn)步而生產(chǎn)發(fā)展需求的與日俱增，促進(jìn)了數(shù)字信號處理學(xué)科的發(fā)展，產(chǎn)生了各種巧妙的信號處理算法；特別是計算機(jī)技術(shù)的飛速發(fā) 展，為數(shù)字信號處理增添了巨大的生命力。 (四 ) 頻譜搬移的 設(shè)計 ................................................................... 錯誤 !未定義書簽。 (三 ) 信號 延時的 設(shè)計 ................................................................... 錯誤 !未定義書簽。綜合話音變換和 一部分音頻特效技術(shù)，設(shè)計了一種廣泛適合于電影配音、話音安全以及娛樂等多領(lǐng)域的實時變聲器，并在短時傅里葉基礎(chǔ)上分析了語音信號的頻譜。 Matlab 是美國 Math Works 公司推出的一種面向工程和科學(xué)計算的交互式計算軟件， 它以矩陣運算 為基礎(chǔ)， 把計算、 可視化、 程序設(shè)計融合到了一個簡單易用的交互式工作環(huán)境中。 本設(shè)計還參考了《 Matlab 在數(shù)字語音分析中的應(yīng)用》、《基于 Matlab 實現(xiàn)對語音信號的頻譜分析》、《基于幀的實時話音變換技術(shù)與 DSP 實現(xiàn)》等方面的學(xué)術(shù)論文，當(dāng)中講述了語音信號的特征，如何有效處理語音信號及相關(guān)算法等。 //(二 )國內(nèi)外數(shù)字濾波器的研究背景和現(xiàn)狀 (三 )本設(shè)計的主要研究內(nèi)容 本論文主要討論了語音信號的分析及處理問題，分為七章： 第一章為緒論，主要介紹了課題設(shè)計的背景、目的和意義。 二 語音信號概述 語音信號處理技術(shù)是語音處理領(lǐng)域中新近發(fā)展起來的一個學(xué)科分支， MATLAB 是一個數(shù)據(jù)分析和處理功能十分強(qiáng)大的工程實用軟件， 運用它來進(jìn)行語音信號的采集、分析和處理相當(dāng)便捷。 （ 5） 作為語音識別的前端預(yù)處理，減少說話人差異的影響。為了得到各種音頻信號的特技效果，也需要對各種變聲算法進(jìn)行研究。這種方法操作非常簡便。 MATLAB 是一個包括數(shù)值計算，高級圖形和可視化的集成科技計算環(huán)境，也是一種高級程序設(shè)計語言。 MATLAB一個重要的特點是易擴(kuò)展性。然而 ,對語音處理來說 ,短時分析的方法是一個比較有效的解決途徑 , 因為語音信號的特性是隨時間緩慢變化的 ,因而可以假設(shè)它在一個短的段時間內(nèi)保持不變。傅立葉變換的定義如公式 (1)所示 : ? ??? ))si n ()c o s(()( 0 k w tbk w taatf ),0( ???k (1) 由于可以認(rèn)為語音信號是局部平穩(wěn)的 , 所以可以對某一幀語音進(jìn)行傅立葉變換 , 即短時傅立葉變換。另一種解釋是從時間軸方向來理解， 當(dāng)頻率固定時 ,例如 k??? ，則 X( jw? )可以看作是信號經(jīng)過一個中心頻率為的帶通濾波器產(chǎn)生的輸出。 回聲 模擬產(chǎn)生出不同的聽音環(huán)