【正文】 es just before GUIDE is made visible.function GUIDE_OpeningF(hObject, eventdata, handles, varargin)% This function has no output args, see OutputF.% hObject handle to figure% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin mand line arguments to GUIDE (see VARARGIN)% Choose default mand line output for GUIDE = hObject。% Update handles structureguidata(hObject, handles)。% UIWAIT makes GUIDE wait for user response (see UIRESUME)% uiwait()。% Outputs from this function are returned to the mand line.function varargout = GUIDE_OutputF(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT)。% hObject handle to figure% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default mand line output from handles structurevarargout{1} = 。% Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[x,fs,bits]=wavread(39。H:\39。)。N=length(x)fssound(x,fs,bits)。y1=fft(x,N)。f=fs*(0:2047)/N。figure(1)magy1=abs(y1)。angy1=angle(y1)。subplot(3,1,1),plot(x)。title(39。原始信號波形39。)subplot(3,1,2),plot(magy1)。title(39。原始信號幅值39。)subplot(3,1,3),plot(angy1)。title(39。原始信號相位39。)figure(2)freqz(x)title(39。頻率響應(yīng)圖39。)figure(3)plot(f,abs(y1(1:2048)))。title(39。原始語音信號頻譜39。)xlabel(39。Hz39。)。ylabel(39。fudu39。)。axis([0 4500 0 400])% Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fs=22050。[x,fs,bits]=wavread(39。H:\39。)。N=length(x)。y1=fft(x,N)。f=fs*(0:511)/N。t=(0:length(x)1)/fs。x1=[*sin(2*pi*10000*t)]39。x2=x+x1。sound(x2,fs,bits)。figure(1)subplot(2,1,1) ,plot(x)。title(39。原語音信號時域圖39。)subplot(2,1,2) ,plot(x2)。title(39。加高斯噪聲后語音信號時域圖39。)xlabel(39。time n39。)。ylabel(39。fudu39。)。y2=fft(x2,N)。figure(2)subplot(2,1,1),plot(abs(y1))。title(39。原始語音信號頻譜39。)xlabel(39。Hz39。)。ylabel(39。fudu39。)。subplot(2,1,2),plot(abs(y2))。title(39。加噪語音信號頻譜39。) xlabel(39。Hz39。)。ylabel(39。fudu39。)。 axis([0 4500 0 300])。wavwrite(x2,fs,39。H:\39。)。% Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fb=1000。fc=1200。fs=22050。wp=*pi。ws=*pi。Rp=1。Rs=15。Fs=22050。Ts=1/Fs。wp1=2/Ts*tan(wp/2)。ws1=2/Ts*tan(ws/2)。[N,Wn]=buttord(wp1,ws1,Rp,Rs,39。s39。)。[Z,P,K]=buttap(N)。[Bap,Aap]=zp2tf(Z,P,K)。[b,a]=lp2lp(Bap,Aap,Wn)。[bd,ad]=bilinear(b,a,Fs)。[h,w]=freqz(bd,ad)。figure(1)subplot(111)。plot(w*fs/(2*pi),abs(h))grid。title(39。濾波器的性能分析39。)。figure(2)[x,fs,bits]=wavread(39。h:\39。)。n=length(x)。f=fs*(0:(n/21))/n。X=fft(x)。z=filter(bd,ad,x)。subplot(211)。plot(x)。title(39。原始信號的波形39。)。subplot(212)。plot(z)。title(39。濾波后信號的波形39。)。figure(3)sound(z,fs,bits)。subplot(2,1,1)。plot(f,abs(X(1:n/2)))。title(39。原始信號的頻譜39。)。xlabel(39。Hz39。)。Z=fft(z)。subplot(2,1,2)。plot(f,abs(Z(1:n/2)))。title(39。濾波后的信號頻譜39。)。xlabel(39。Hz39。)。wavwrite(z,fs,39。h:\39。)。% Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)fs=22050。[x,fs,bits]=wavread(39。h:\39。)。wp=*pi。ws=*pi。wdelta=wswp。N=ceil(*pi/wdelta)。t=0:(size(x)1)。wn=(+)*pi/2。b=fir1(N,wn/pi,hamming(N+1))。f1=fftfilt(b,x)。figure(1)freqz(b,1,512)[h1,w1]=freqz(b,1)。plot(w1*fs/(2*pi),20*log10(abs(h1)))。figure(2)subplot(2,1,1)plot(t,x)title(39。濾波前的時域波形39。)。subplot(2,1,2),plot(t,f1)。title(39。濾波后的時域波形39。)。sound(f1)。F0=fft(f1,1024)。f=fs*(0:511)/1024。figure(3)y2=fft(x,1024)。subplot(2,1,1)。plot(f,abs(y2(1:512)))。title(39。濾波前的頻譜39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。subplot(2,1,2)F1=plot(f,abs(F0(1:512)))。title(39。濾波后的頻譜39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。wavwrite(f1,fs,39。h:\39。)。% Executes on button press in pushbutton5.function pushbutton5_Callback(hObject, eventdata, handles)% hObject handle to pushbutton5 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[x,fs,bits]=wavread(39。H:\39。)。N=length(x)sound(x,fs,bits)。y1=fft(x,N)。f=fs*(0:2047)/N。figure(1)magy1=abs(y1)。angy1=angle(y1)。subplot(3,1,1),plot(x)。title(39。原始信號波形39。)subplot(3,1,2),plot(magy1)。title(39。原始信號幅值39。)subplot(3,1,3),plot(angy1)。title(39。原始信號相位39。)figure(2)freqz(x)title(39。頻率響應(yīng)圖39。)figure(3)plot(f,abs(y1(1:2048)))。title(39。原始語音信號頻譜39。)xlabel(39。Hz39。)。ylabel(39。fudu39。)。axis([0 4500 0 400])% Executes on button press in pushbutton6.function pushbutton6_Callback(hObject, eventdata, handles)% hObject handle to pushbutton6 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[y,fs,bits]=wavread(39。H:\39。)。N=length(y)z=[zeros(3000,1)。y]。%延遲3000?Z=fft(z,N)。figure。subplot(2,2,1:2)。plot(z)。title(39。延時后時域波形39。)。subplot(2,2,3)。plot(abs(Z))。title(39。延時后幅頻39。)。subplot(2,2,4)。plot(angle(Z))。title(39。延時后相頻39。)。sound(y,fs,bits)。wavwrite(y,fs,39。H:\39。)。% Executes on button press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)% hObject handle to pus