【滤波器】基于低通、带通、高通滤波器实现语音去噪含Matlab源码

function varargout = AMtiaozhi(varargin)% AMTIAOZHI M-file for AMtiaozhi.fig%      AMTIAOZHI, by itself, creates a new AMTIAOZHI or raises the existing%      singleton*.%%      H = AMTIAOZHI returns the handle to a new AMTIAOZHI or the handle to%      the existing singleton*.%%      AMTIAOZHI('CALLBACK',hObject,eventData,handles,...) calls the local%      function named CALLBACK in AMTIAOZHI.M with the given input arguments.%%      AMTIAOZHI('Property','Value',...) creates a new AMTIAOZHI or raises the%      existing singleton*.  Starting from the left, property value pairs are%      applied to the GUI before AMtiaozhi_OpeningFunction gets called.  An%      unrecognized property name or invalid value makes property application%      stop.  All inputs are passed to AMtiaozhi_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 AMtiaozhi% Last Modified by GUIDE v2.5 10-May-2010 16:42:40% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name',       mfilename, ...                   'gui_Singleton',  gui_Singleton, ...                   'gui_OpeningFcn', @AMtiaozhi_OpeningFcn, ...                   'gui_OutputFcn',  @AMtiaozhi_OutputFcn, ...                   'gui_LayoutFcn',  [] , ...                   'gui_Callback',   []);if nargin & isstr(varargin{1})    gui_State.gui_Callback = str2func(varargin{1});endif 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 AMtiaozhi is made visible.function AMtiaozhi_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 AMtiaozhi (see VARARGIN)% Choose default command line output for AMtiaozhihandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes AMtiaozhi wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = AMtiaozhi_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 structurevarargout{1} = handles.output;% --- Executes during object creation, after setting all properties.function popupmenu1_CreateFcn(hObject, eventdata, handles)% hObject    handle to popupmenu1 (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    empty - handles not created until after all CreateFcns called% Hint: popupmenu controls usually have a white background on Windows.%       See ISPC and COMPUTER.if ispc    set(hObject,'BackgroundColor','white');else    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));endset(hObject, 'String', {'调制信号', '载波信号','已调波信号'});% --- Executes on selection change in popupmenu1.function popupmenu1_Callback(hObject, eventdata, handles)% hObject    handle to popupmenu1 (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,'String') returns popupmenu1 contents as cell array%        contents{get(hObject,'Value')} returns selected item from popupmenu1% --- Executes during object creation, after setting all properties.function popupmenu2_CreateFcn(hObject, eventdata, handles)% hObject    handle to popupmenu2 (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    empty - handles not created until after all CreateFcns called% Hint: popupmenu controls usually have a white background on Windows.%       See ISPC and COMPUTER.if ispc    set(hObject,'BackgroundColor','white');else    set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));endset(hObject, 'String', {'高通滤波', '低通滤波','带通滤波'});% --- Executes on selection change in popupmenu2.function popupmenu2_Callback(hObject, eventdata, handles)% hObject    handle to popupmenu2 (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,'String') returns popupmenu2 contents as cell array%        contents{get(hObject,'Value')} returns selected item from popupmenu2% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)popup_sel_index = get(handles.popupmenu1, 'Value');switch popup_sel_index    case 1        dt=0.01;        t=0:dt:10;Fs=1/dt;        e=sinc(t-5);%调制信号        elength=length(e);%求采样点个数        E=fft(e,elength);%快速傅里叶变换        %fe=-Fs/2:elength:Fs/2;        En=fftshift(E);%对fft输出进行重新排列，将零频分量移到频谱中心        axes(handles.axes1);        plot(t,e);title('调制信号');                                xlabel('时间s');        ylabel('强度值');        axis([0 10 1.1*min(e) 1.1*max(e)]);        axes(handles.axes2);        plot(abs(En));title('调制信号频谱');        axis([0 1000 1.1*min(abs(En)) 1.1*max(abs(En))]);        set(handles.axes2,'Xticklabel',[-50,-30,-10,10,30,50])                                xlabel('频率Hz');        ylabel('幅度值');    case 2        dt=0.01;        t=0:dt:10;Fs=1/dt;        s=cos(60*pi*t)+cos(40*pi*t)+cos(20*pi*t);%载波信号        slength=length(s);%求采样点个数        S=fft(s);%快速傅里叶变换        %fs=linspace(-Fs/2,Fs/2,slength);        Sn=fftshift(S);        axes(handles.axes3);        plot(t,s);title('载波信号');                                xlabel('时间s');        ylabel('强度值');        axes(handles.axes4);plot(abs(Sn));title('载波信号频谱');                axis([0 1000 1.1*min(abs(Sn)) 1.1*max(abs(Sn))]);        set(handles.axes4,'Xticklabel',[-50,-30,-10,10,30,50])                                xlabel('频率Hz');        ylabel('幅度值');    case 3        dt=0.01;        t=0:dt:10;Fs=1/dt;        s=cos(60*pi*t)+cos(40*pi*t)+cos(20*pi*t);%载波信号        e=sinc(t-5);%调制信号        p=e.*s;%调制        plength=length(p);        P=fft(p);        %fp=linspace(-Fs/2,Fs/2,plength);        Pn=fftshift(P);        axes(handles.axes5);        plot(t,p);title('调幅信号');                                xlabel('时间s');        ylabel('强度值');        axis([0 10 1.1*min(p) 1.1*max(p)]);        axes(handles.axes6);        plot(abs(Pn));title('调幅信号频谱');                        axis([0 1000 1.1*min(abs(Pn)) 1.1*max(abs(Pn))]);        set(handles.axes6,'Xticklabel',[-50,-30,-10,10,30,50])                                        xlabel('频率Hz');        ylabel('幅度值');end% 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)% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)pend% 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)