【正文】 y6 = stereo_ing(y_original,Fs,mode)。)。 sound(y5,Fs)。Choose the mode 1左聲道慢于右聲道 ,2右聲道慢于左聲道 :39。 %%%%%%%%%%%%%%%%%%%相位變化 %%%%%%%%%%%%%%%%%%%% %實現(xiàn)：將左（或者右）聲道前一段添 0，文件在播放的時候左右聲道就會不同步右（左）聲道的數(shù)據(jù)會超前 fprintf(39。)。音效左右移動效果選擇 \n39。)?；仨懶Чx擇 \n39。)?；芈曅Чx擇 \n39。)。音量調(diào)節(jié) \n39。D:\bs2\39。聲音特效播放順序： \n39。 %%%%聲音特效處理 %%%%%% fprintf(39。頻率 39。 xlabel(39。 contour(abs(TF))。頻率 f/hz39。 xlabel(39。 %把譜圖存放在 TF中 end subplot(2,2,2) mesh(abs(TF))。 xw=sm1(n1:n2)。%每段實際取的點數(shù) nfft=512。 d=100。振幅 39。 xlabel(39。 plot(time2,sm1)。 %搬移到較高頻段后信號回放 ns=size(sm1)。 %將一維數(shù)組拓展成二維 y3=reshape(y3_2,2,length(y3_2)/2)。,fs)。n7=%g\n39。%失真信號回放 %wavplay(y,8000)。頻率 f/hz39。 xlabel(39。 %把譜圖存放在 TF中 end subplot(2,2,4) mesh(abs(TF))。 xw=z(n1:n2)。%每段實際取的點數(shù) nfft=512。 d=100。振幅 39。 xlabel(39。 plot(time1,z)。 %對語音信號進行延時 wavplay(z,fs)。)。時間 39。 %三維繪圖 title(39。 temp=fft(xw,nfft)。 TF=zeros(Ts,nfft)。 %連續(xù)段之間的重疊長度 L=512。)。時間 39。 title(39。 time=linspace(0,9,length(y))。D:\bs2\39。figure1Original signal Delayed signal,figure2FFT_shifted signal\n39。 fprintf(39。)。1Original signal： \n39。 fprintf(39。 %縱、橫坐標軸采用等長刻度 set the axes aspect ratio grid。 ylabel(39。)。 subplot(2,2,3)。 %量化間隔， det ylabel(39。)。 %頻譜以 0頻為中心 ,保證所有頻譜的坐標值一樣 TF(i,:)=temp。 n2=(Ld)*(i1)+NW。 Ts=fix(ns/k)。 NW=512。 %x軸 ylabel(39。)。 subplot(2,2,1)。 wavplay(sm1,fs)。%01之間有 n7個點 ,保證截取的語音信號和載波信號所取點數(shù)相同 y3=cos(2*2021*pi*t3)。 %fprintf(39。%失真信號回放 %%%%%%%%%%%%%%%頻譜搬移 短時譜分析 %%%%%%%%%%%%%%%%%%%%%% n7=length(y)。 grid。)。STFTDelay39。 %FFT變換 temp=fftshift(temp)。 %將存放三維譜圖，先清零 for i=1:Ts n1=(Ld)*(i1)+1。 %每段長為 512 k=NWd。 %y軸 grid。)。延時后信號時域波形 39。 time1=linspace(0,duration,length(z))。 %%%%%%%%%%%%%%% 延時后時、頻分析 %%%%%%%%%%%%%%%%%%%%%% z=[zeros(1500,1)。 %量化間隔， det ylabel(39。)。 %頻譜以 0頻為中心 ,保證所有頻譜的坐標值一樣 TF(i,:)=temp。 n2=(Ld)*(i1)+NW。 Ts=fix(ny/k)。 NW=512。 %x軸 ylabel(39。)。 subplot(2,2,1)。 wavplay(y,fs)。)。1Original signal： \n39。 fprintf(39。 fprintf(39。 y=wavrecord(duration*fs,fs)。 pause。 %取樣頻率 duration=3。 衷心感謝電子 Q0841 班的同學們，我們一起學習、一起研究，共同進步，平日里大家關(guān)于畢業(yè)設(shè)計的討論給了我很多啟發(fā)，在此深表謝意。 六 致 謝詞 本設(shè)計是在李俊老師的悉心指導(dǎo)下完成的，論文從選題到測試成功再到最后成文，無不傾注著李老師的心血，為我個人的成長和進步傾注了大量的汗水，誠摯地感謝李老師給與我的關(guān)懷和教誨，老師高深的學術(shù)造詣、嚴謹?shù)闹螌W態(tài)度、一絲不茍的工作作風、淵博的知識以及誨人不倦的導(dǎo)師風范，使我在學習生活中受益匪淺，并激勵著我在今后的工作和學習中不斷進步。在寫代碼時，把每一個變量與過程考慮清楚才可能完成整個程序。 end （七） 測試結(jié)果分析 （一） 頻譜分析： 圖 4 分別給出了原始信號、延時后的信號、調(diào)制后信號的時域波形、及基于短時傅里葉變換的各頻域波形 圖 3 原語音信號與延時語音信號的時域、頻域?qū)Ρ? 圖 4 語音信號頻譜搬移后時域、頻域及等高線圖 （二） 語音信號 特效處理 圖 5給出了系統(tǒng)輸出界面顯示，由該圖可以看出系統(tǒng)所處理的特效。 if mode == 1 am = linspace(0,1,Len)39。 end 淡入淡出： %% function y = stereo_ing(y_original,Fs,mode) %mode = 1 : ming %mode = else : leaving %[y_original Fs]= wavread(39。y_1]。zeros(delay_bits,1)]。 y_1 = y_original(:,1)。 end 相位變化： %% function y = stereo_phase_movement(y_original,Fs,mode,delay_time) %mode = 1: 左聲道慢于右聲道 %mode = else：右聲道慢于左聲道 %delay_time = 。 i = linspace(1,0,length(y_stereo(:,2)))39。 左右移動： %%stereo_shift function y = stereo_shift(y_stereo,Fs,mode) % mode = 1 : from right to left % mode = else : from left to right %[y_original Fs]= wavread(39。zeros(delay_bits,1)。 y_delay2 = [zeros(delay_bits,1)。 y_2 = y_original(:,2)。 y_delay_delay1 = [zeros(delay_bits,1)。zeros(delay_bits,1)]。 % 太空效果 end delay_bits = floor(Fs * delay_time)。 % 大廳效果 elseif mode ==3 scale = 。 回響： %% % 回響效果 function y = stereo_reecho(y_original,Fs,mode) if mode == 1 scale = 。 % delay %t = (0:length(y_delay1)1)39。 y_y2 = [y_original(:,2)。zeros(floor(Fs*delay_time),1)]。delay_time = 。delay_time = 。)。音量調(diào)節(jié) \n39。 axis equal 。)。等高圖 Move39。 grid。)。STFTMove39。 %FFT變換 temp=fftshift(temp)。 %將存放三維譜圖，先清零 for i=1:Ts n1=(Ld)*(i1)+1。 %每段長為 512 k=NWd。 將 短時分析應(yīng)用于頻域分析即是傅立葉變換 , 相應(yīng)的頻譜稱為“ 短時譜” , 即有限長度的傅立葉變換。.*y3。具體程序如下： n7=length(y)。記 Rate = Fs’ /Fs ，那么 FL’ =FL * Rate，其中 Rate 為基頻變化 率。y]。 （三） 信號延時的設(shè)計 信號延時通常有以下方法： 方法一：構(gòu)造一個增益為 1，相位與時間成正比變化的濾波器，讓語音信號通過該延時網(wǎng)絡(luò)，即可實現(xiàn)延時。 fprintf(39。 y=wavrecord(duration*fs,fs)。 提前錄制一段語音文件，時間10s 以內(nèi) 直 接 調(diào) 用wavrecord 函數(shù)實時錄音 語音信號的頻譜分析：畫出采樣后語音信號的時域波形和短時傅里葉變換后的頻譜圖 將信號加入延時，再分析其時域波形與頻譜，并分別與 原始信號進行比較 對原始語音信號進行音量變換、回聲、回響、相位移動、淡入淡出等特效處理 回放語音信號 pause。 %取樣頻率 duration=3。實現(xiàn)時，只需要線性改變音量即可。 左右移動 感覺聲音從一邊移到另一邊，即聲音從左到右，或從右到左（只針對立體聲）。根據(jù)研究，這兩個參數(shù)可以確定如表 2： 表 2 回聲效果參數(shù) 參數(shù) 峽谷 大廳 小廳 太空 延遲音量 延遲時間 500ms 100ms 60ms 180ms 回響 在房間中，由于房間四壁對聲音的反射，造成多個聲音的混合而形成的特殊聲音效果稱為回響。 淡入淡出 聲音由小變大，或由大變小。 回聲 模擬產(chǎn)生出不同的聽音環(huán)境（峽谷，大廳，小廳，太空）。不同的人以及同一個人在不同的年齡時期有不同的基音周期。另一種解釋是從時間軸方向來理解， 當頻率固定時 ,例如 k??? ，則 X( jw? )可以看作是信號經(jīng)過一個中心頻率為的帶通濾波器產(chǎn)生的輸出。 由公式 (2)知 ,短時傅立葉變換有兩個自變量 n 和 w, 所以它既是關(guān)于時間 n 的離散函數(shù) , 又是關(guān)于角頻率 w的連續(xù)函數(shù)。傅立葉變換的定義如公式 (1)所示 : ? ??? ))si n ()c o s(()( 0 k w tbk w taatf ),0( ???k (1) 由于可以認為語音信號是局部平穩(wěn)的 , 所以可以對某一幀語音進行傅立葉變換 , 即短時傅立葉變換。 最后 IFFT 觀察恢復(fù)后的時域信號圖形， 頻譜分析過程 （如圖 1所示） 。然而 ,對語音處理來說 ,短時分析的方法是一個比較有效的解決途徑 , 因為語音信號的特性是隨時間緩慢變化的 ,因而可以假設(shè)它在一個短的段時間內(nèi)保持不變。 表 1 與語音輸入輸出相關(guān)的函數(shù) 函數(shù)名 功能描述 wavread auread 讀語音文件 wavwrite auwrite 寫語音文件 sound wavplay 放音 soundsc 歸一化放音 soundview 可視化語音輸出 wavrecord 錄音 三 基于 MATLAB 的語音信號短時傅里葉變換 時域上的原函數(shù)中含有包含全部信息量