【正文】 2double(get(hObject,39。String39。)) returns contents of fs as a double% Executes during object creation, after setting all properties.function fs_CreateF(hObject, eventdata, handles)% hObject handle to fs (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles empty handles not created until after all CreateFs called% Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc amp。amp。 isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,39。white39。)。end% Executes on button press in begin.function begin_Callback(hObject, eventdata, handles)% hObject handle to begin (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)N1=str2double(get(,39。String39。))。N2=str2double(get(,39。String39。))。fs=str2double(get(,39。String39。))。[x,fs,bits]=wavread(39。D:\Matlab\work\39。,[N1 N2])。sound(x,fs,bits)。 %播放語音信號%n=N2N1+1。%當語音信號文件較大時用這兩條%x1=reshape(x,1,2*n)。%語句替換x1=x39。x1=x39。y1=fft(x1)。axes()。plot(x1) %做原始語音信號的時域圖形title(39。語音信號時域波形39。)xlabel(39。n39。)。ylabel(39。幅值39。)。M=nextpow2(x1)。 % 求x的長度對應(yīng)的2的最低冪次m N=2^M。if length(x1)N x1=[x1,zeros(1,Nlength(x1))]。 % 若x的長度不是2的冪，補零到2的整數(shù)冪 end %數(shù)據(jù)倒序操作J=0。%給倒序數(shù)賦初值for I=0:N1。%按序交換數(shù)據(jù)和算倒序數(shù) if IJ。%條件判斷及數(shù)據(jù)交換 T=x1(I+1)。x1(I+1)=x1(J+1)。 x1(J+1)=T。 end %算下一個倒序數(shù) K=N/2。 while J=K。 J=JK。K=K/2。 end J=J+K。end %x1。 y=x1。 % 將x倒序排列作為y的初始值 WN=exp(i*2*pi/N)。 for L=1:M B=2^L/2。%第L級中，每個蝶形的兩個輸入數(shù)據(jù)相距B個點，每級有B個不同的旋轉(zhuǎn)因子 for J=0:B1 % J代表了不同的旋轉(zhuǎn)因子 p=J*2^(ML)。 WNp=WN^p。 for k=J+1:2^L:N % 本次蝶形運算的跨越間隔為2^L kp=k+B。 % 蝶形運算的兩個因子對應(yīng)單元下標的關(guān)系 t=y(kp)*WNp。 % 蝶形運算的乘積項 y(kp)=y(k)t。 % 蝶形運算， 注意必須先進行減法運算，然后進行加法運算，否則要使用中間變量來傳遞y(k) y(k)=y(k)+t。 % 蝶形運算 end end end %y axes()。[x1,w1]=freqz(x1,1)。 %繪制原始語音信號的頻率圖plot(w1,abs(x1))。title(39。x1的頻率響應(yīng)特性39。)xlabel(39。頻率39。)。ylabel(39。幅度39。)。axes()。plot(abs(y1)) %Matlab自帶的FFT函數(shù)實現(xiàn)的語音信號的FFT頻譜圖title(39。Matlab自帶的FFT函數(shù)實現(xiàn)的x1的頻譜39。)xlabel(39。K39。)。ylabel(39。Y1(k)39。)。axes()。plot(abs(y))。 %編寫的FFT程序?qū)崿F(xiàn)的語音信號的FFT頻譜圖title(39。編寫的FFT程序?qū)崿F(xiàn)的的x1的頻譜39。) xlabel(39。K39。)。ylabel(39。Y(k)39。)。% Executes on button press in quit.function quit_Callback(hObject, eventdata, handles)% hObject handle to quit (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)close all。% Executes on selection change in ringin.function ringin_Callback(hObject, eventdata, handles)% hObject handle to ringin (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: contents = get(hObject,39。String39。) returns ringin contents as cell array% contents{get(hObject,39。Value39。)} returns selected item from ringin% Executes during object creation, after setting all properties.function ringin_CreateF(hObject, eventdata, handles)% hObject handle to ringin (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles empty handles not created until after all CreateFs called% Hint: popupmenu controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc amp。amp。 isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,39。white39。)。e