数字信号处理音频FIR去噪（基于MATLAB GUI的开发）

function varargout = yanshou(varargin) % YANSHOU M-file for yanshou.fig % YANSHOU, by itself, creates a new YANSHOU or raises the existing % singleton*. % % H = YANSHOU returns the handle to a new YANSHOU or the handle to % the existing singleton*. % % YANSHOU('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in YANSHOU.M with the given input arguments. % % YANSHOU('Property','Value',...) creates a new YANSHOU or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before yanshou_OpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to yanshou_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 % Edit the above text to modify the response to help yanshou % Last Modified by GUIDE v2.5 23-Apr-2020 14:20:34 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @yanshou_OpeningFcn, ... 'gui_OutputFcn', @yanshou_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 yanshou is made visible. function yanshou_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 yanshou (see VARARGIN) % Choose default command line output for yanshou handles.output = hObject; % Update handles structure guidata(hObject, handles); % UIWAIT makes yanshou wait for user response (see UIRESUME) % uiwait(handles.figure1); % --- Outputs from this function are returned to the command line. function varargout = yanshou_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 open_pushbutton1. function open_pushbutton1_Callback(hObject, eventdata, handles) % hObject handle to open_pushbutton1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global x; %文件 global Fs; %采样频率 global tl; global x2; [filename, pathname] = uigetfile('*.wav', '选择音频文件'); if isequal(filename,0) disp('User selected Cancel') else path = fullfile(pathname, filename); [handles.x,handles.Fs]=audioread(path); x=handles.x; Fs=handles.Fs; axes(handles.axes1); tl=[0:1/Fs:(length(handles.x)-1)/Fs]; %时间尺度 plot(tl,handles.x); title('语音时域波形'); xlabel('时间/s'); grid on; N=length(handles.x); df=Fs/N; w=[0:df:df*(N-1)] - Fs/2; %频率尺度 X=fft(handles.x); X=fftshift(X); axes(handles.axes2); plot(w,abs(X)/max(abs(X))); axis([-10000,10000,0,1]); title('语音频谱'); xlabel('频率/Hz'); grid on; x2=x; end % --- Executes on button press in play_pushbutton2. function play_pushbutton2_Callback(hObject, eventdata, handles) % hObject handle to play_pushbutton2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global x2; global Fs; sound(x2,Fs); % --- Executes on button press in add_pushbutton3. function add_pushbutton3_Callback(hObject, eventdata, handles) % hObject handle to add_pushbutton3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global x; global Fs; global tl; global x2; axes(handles.axes1); size(x); t=0:1/Fs:(length(x)-1)/Fs;%将所加噪声信号的点数调整到与原始信号相同 Au=0.07; fn = get(handles.noise_edit2,'string'); %噪声频率 fn = str2double(fn); noise=Au*cos(2*pi*fn*t)'; %噪声为频率 x=x+noise; plot(tl,x); title('加入噪声后语音时域波形'); xlabel('时间/s'); grid on; N=length(x); df=Fs/N; w=[0:df:df*(N-1)] - Fs/2; %频率尺度 X=fft(x); X=fftshift(X); axes(handles.axes2); plot(w,abs(X)/max(abs(X))); axis([-10000,10000,0,1]); title('加入噪声后语音频谱'); xlabel('频率/Hz'); grid on; x2=x; % --- Executes on button press in low_pushbutton5. function low_pushbutton5_Callback(hObject, eventdata, handles) % hObject handle to low_pushbutton5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global x; global Fs; global tl; global x2; x1=x; % fp=1000; % fs = 1000; % Wp = 2*fp/Fs; % Ws = 2*fs/Fs; % if(Wp >= 1) % Wp = 0.99; % end % if(Ws >= 1) % Ws = 0.99; % end fp = get(handles.edit3,'string'); fp = str2double(fp)*2; if get(handles.radiobutton1,'value') [n, Wn]=buttord(Wp,Ws, 2, 15); [b, a]=butter(n, Wn,'low'); axes(handles.axes3); [h,w]=freqz(b,a); plot(w/pi*Fs/2,abs(h)); x1=filter(b,a,x1); %调用函数滤波 elseif get(handles.radiobutton4,'value') b2=fir1(30, fp/Fs, boxcar(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton5,'value') b2=fir1(30, fp/Fs, triang(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton6,'value') b2=fir1(30, fp/Fs, hamming(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton7,'value') b2=fir1(30, fp/Fs, hanning(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton8,'value') b2=fir1(30, fp/Fs, blackman(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton9,'value') b2=fir1(30,fp/Fs, kaiser(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/pi*Fs/2,20*log(abs(h))); x1=fftfilt(b2,x1); end; axes(handles.axes5); plot(tl,x1); title('滤除噪声后语音时域波形'); xlabel('时间/s'); N=length(x1); df=Fs/N; w=[0:df:df*(N-1)] - Fs/2; %频率尺度 X=fft(x1); X=fftshift(X); axes(handles.axes6); plot(w,abs(X)/max(abs(X))); axis([-10000,10000,0,1]); title('滤除噪声后语音频谱'); xlabel('频率/Hz'); grid on; x2=x1; % --- Executes during object creation, after setting all properties. function low_pushbutton5_CreateFcn(hObject, eventdata, handles) % hObject handle to low_pushbutton5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % --- Executes on button press in pushbutton6. function noise_edit2_Callback(hObject, eventdata, handles) % hObject handle to noise_edit2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data