基于MATLAB的语音信号处理课程设计

function varargout = untitled2(varargin) % UNTITLED2 M-file for untitled2.fig % UNTITLED2, by itself, creates a new UNTITLED2 or raises the existing % singleton*. % % H = UNTITLED2 returns the handle to a new UNTITLED2 or the handle to % the existing singleton*. % % UNTITLED2('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in UNTITLED2.M with the given input arguments. % % UNTITLED2('Property','Value',...) creates a new UNTITLED2 or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before untitled2_OpeningFunction gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to untitled2_OpeningFcn via varargin. % % *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % instance to run (singleton)". % % See also: GUIDE, GUIDATA, GUIHANDLES % Copyright 2002-2003 The MathWorks, Inc. % Edit the above text to modify the response to help untitled2 % Last Modified by GUIDE v2.5 16-Dec-2011 18:46:46 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @untitled2_OpeningFcn, ... 'gui_OutputFcn', @untitled2_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end % End initialization code - DO NOT EDIT % --- Executes just before untitled2 is made visible. function untitled2_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn. % 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 command line arguments to untitled2 (see VARARGIN) % Choose default command line output for untitled2 handles.output = hObject; % Update handles structure guidata(hObject, handles); % UIWAIT makes untitled2 wait for user response (see UIRESUME) % uiwait(handles.figure1); % --- Outputs from this function are returned to the command line. function varargout = untitled2_OutputFcn(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 command line output from handles structure varargout{1} = handles.output; % --- Executes on button press in getfile. function getfile_Callback(hObject, eventdata, handles) % hObject handle to getfile (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global player; [fname,pn,fi]=uigetfile('*.wav'); global Fname; Fname=fname; % --- Executes on button press in Play. function Play_Callback(hObject, eventdata, handles) % hObject handle to Play (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global Fname; global sig; global player; global fs; [y,fs,nbits]=wavread(Fname); sig=y; voice=get(handles.Vol,'Value')+2; %handles.变量名 spead=get(handles.Speed,'Value')+0.5; player=audioplayer(y*voice,fs*spead); play(player); pause(3); % --- Executes on button press in Stop. function Stop_Callback(hObject, eventdata, handles) % hObject handle to Stop (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global player; pause(player); % --- Executes on button press in Resume. function Resume_Callback(hObject, eventdata, handles) % hObject handle to Resume (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global player; resume(player); % --- Executes on button press in Time. function Time_Callback(hObject, eventdata, handles) % hObject handle to Time (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global sig; figure(1);plot(sig);title('时域波形'); % --- Executes on button press in Frequency. function Frequency_Callback(hObject, eventdata, handles) % hObject handle to Frequency (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global sig; Y=fft(sig); mag=abs(Y); phase=angle(Y); phase=phase*180/pi; figure(1); subplot(3,1,1),plot(Y); title('频谱响应'); subplot(3,1,2),plot(mag); title('幅频响应'); subplot(3,1,3),plot(phase); title('相频响应'); % --- Executes on button press in Highpass. function Highpass_Callback(hObject, eventdata, handles) % hObject handle to Highpass (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global fs; global sig; x=sig; sound(x,fs); figure(1);subplot(4,1,1);plot(x);title('滤波前的信号时域波形'); Y=fft(x,fs); subplot(4,1,2); plot(abs(Y(1:length(Y)/2))); title('滤波前信号的频谱'); prompt={'阻带截止频率','通带截止频率'}; name='选择参数'; numlines=1; defAns={'550','650'} answer=inputdlg(prompt,name,1,defAns); ws=str2num(answer{1}); wp=str2num(answer{2}); wzp=2*pi*wp/fs;wzs=2*pi*ws/fs; Bt=wzp-wzs; N1=ceil(8*pi/Bt); wc=wzs+Bt/2; hn=fir1(N1,wc,'high');%注意格式。没有hanning(N+1) ?缺省则默认为哈明窗. wf=0:pi/511:pi; HK=freqz(hn,wf); wHz=wf*511/(2*pi); % 转化为Hz subplot(4,1,3); plot(20*log10(abs(HK)));%绘制损耗曲线 %plot(wHz,20*log10(abs(HK))); xlabel('频率(Hz)');ylabel('幅度'); z=fftfilt(hn,x);% 进行滤波 Y1=fft(z,fs); subplot(4,1,4); plot(abs(Y1(1:length(Y)/2))); title('滤波后信号的频谱'); pause(1); sound(z,fs); % --- Executes on button press in Lowpass. function Lowpass_Callback(hObject, eventdata, handles) % hObject handle to Lowpass (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global fs; global sig; x=sig; sound(x,fs); figure(1);subplot(4,1,1);plot(x);title('滤波前的信号时域波形'); Y=fft(x,fs); subplot(4,1,2); plot(abs(Y(1:length(Y)/2))); title('滤波前信号的频谱'); prompt={'阻带截止频率','通带截止频率'}; name='选择参数'; numlines=1; defAns={'650','500'} answer=inputdlg(prompt,name,1,defAns); ws=str2num(answer{1}); wp=str2num(answer{2}); wzp=2*pi*wp/fs;wzs=2*pi*ws/fs; Bt=wzs-wzp; N1=ceil(8*pi/Bt); wc=wzp+Bt/2; hn=fir1(N1-1,wc,hanning(N1)); wf=0:pi/511:pi; HK=freqz(hn,wf); wHz=wf*511/(2*pi); % 转化为Hz subplot(4,1,3); plot(20*log10(abs(HK)));%绘制损耗曲线 %plot(wHz,20*log10(abs(HK))); xlabel('频率(Hz)');ylabel('幅度'); z=fftfilt(hn,x);% 进行滤波 Y1=fft(z,fs); subplot(4,1,4); plot(abs(Y1(1:length(Y)/2))); title('滤波后信号的频谱'); pause(1); sound(z,fs); % --- Executes on button press in Bandpass. function Bandpass_Callback(hObject, eventdata, handles) % hObject handle to Bandpass (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global fs; global sig; x=sig; sound(x,fs); figure(1);subplot(4,1,1);plot(x);title('滤波前的信号时域波形')