含噪语音matlab界面处理.docx

1、含噪语音matlab界面处理大连理工大学实验报告学院(系): 电子信息与电气工程学部 专业: 生物医学工程 班级: 电医1001 姓 名: 张 颖 学号: 201081623 组: _ 实验时间: 2012.5.27 实验室: C221 实验台: 指导教师签字: 成绩: 实验名称 含噪语音信号的分析1、实验目的和要求 1. 理解并掌握系统的理念; 2. 掌握滤波的概念和基本应用的方法; 3.提高分析和解决实际问题的能力二实验原理和内容 详见试验讲义3、主要仪器设备 Matlab仿真软件四实验步骤与操作方法五试验题目及程序代码结果 程序代码:function varargout = noised

2、eal(varargin) gui_Singleton = 1;gui_State = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, noisedeal_OpeningFcn, . gui_OutputFcn, noisedeal_OutputFcn, . gui_LayoutFcn, , . gui_Callback, );if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);endif

3、nargout varargout1:nargout = gui_mainfcn(gui_State, varargin:);else gui_mainfcn(gui_State, varargin:);endfunction noisedeal_OpeningFcn(hObject, eventdata, handles, varargin)handles.output = hObject;guidata(hObject, handles);function varargout = noisedeal_OutputFcn(hObject, eventdata, handles) vararg

4、out1 = handles.output;function Untitled_1_Callback(hObject, eventdata, handles)function open_Callback(hObject, eventdata, handles)a=part1.wav;p=mfilename(fullpath)i=findstr(p,)p=(1:i(end);b=cd;y,fs,nbite=wavread(fullfile(b,a);handles.y=y;handles.fs=fs;axes(handles.axes1)subplot(121);plot(y);len=leng

5、th(y);handles.len=len;ts=(len-1)/fs;t=0:1/fs:ts;y1=y;fy=fft(y1);f=(0:length(t)-1)*fs/length(t)-fs/2;subplot(122);plot(f,fftshift(abs(fy);axes(handles.axes2);guidata(hObject,handles);function save_Callback(hObject, eventdata, handles)wavwrite(handles.q,new);function fir_Callback(hObject, eventdata, h

6、andles)handles.k1=get(handles.fir,value);if handles.k1=1 set(handles.listbox1,visible,on); set(handles.listbox2,visible,off); set(handles.iir,value,0);elseend;guidata(hObject,handles);function iir_Callback(hObject, eventdata, handles)if handles.k2=1 set(handles.listbox1,visible,off); set(handles.lis

7、tbox2,visible,on); set(handles.fir,value,0);elseend;guidata(hObject,handles);function listbox1_Callback(hObject, eventdata, handles)function listbox1_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundCol

8、or,white);endfunction sure_Callback(hObject, eventdata, handles)k=get(handles.fir,value);if k=1 run Parameter;else run settingendfunction sound_Callback(hObject, eventdata, handles)function regional_Callback(hObject, eventdata, handles)wavplay(handles.y,handles.fs);function dials_Callback(hObject, e

9、ventdata, handles)wavplay(handles.q,handles.fs);function paint_Callback(hObject, eventdata, handles)k1=get(handles.fir,value);k2=get(handles.iir,value);if k1k2 handles.l=getappdata(0,l); handles.f=str2num(getappdata(0,f); handles.kp1=str2num(getappdata(0,pk1); handles.kp2=str2num(getappdata(0,pk2);

10、handles.sk1=str2num(getappdata(0,sk1); handles.sk2=str2num(getappdata(0,sk2); handles.apl=str2num(getappdata(0,apl); handles.asl=str2num(getappdata(0,asl); choose=getappdata(0,choose); if choose=3 wp=handles.kp1 handles.kp2; wr=handles.sk1 handles.sk2; else wp=handles.kp1; wr=handles.sk1; end ch=get

11、(handles.listbox1,value); switch ch case 1 N,Wn=buttord(wp/pi,wr/pi,handles.apl,handles.asl); b,a= butter(N,Wn,handles.l); case 2 N,Wn=cheb1ord(wp/pi,wr/pi,handles.apl,handles.asl); b,a= cheby1(N,handles.apl,Wn,handles.l); case 3 end f,w=freqz(b,a); f_mag=abs(f); plot(w/pi,abs(f_mag); handles.a=a; h

12、andles.b=b;else handles.h=getappdata(0,h); handles.n=str2num(getappdata(0,n); handles.wc=str2num(getappdata(0,wc); cl=get(handles.listbox2,value); if handles.h=1 chr=low; else chr=high; end switch cl case 1 b=fir1(handles.n,handles.wc,chr,boxcar(handles.n+1); case 2 b=fir1(handles.n,handles.wc,chr,h

13、amming(handles.n+1); case 3 end stem(b); handles.b=b;endguidata(hObject,handles);function listbox2_Callback(hObject, eventdata, handles)function listbox2_CreateFcn(hObject, eventdata, handles)if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,Backgroun

14、dColor,white);endfunction effect_Callback(hObject, eventdata, handles)q=filtfilt(handles.b,handles.a,handles.nos);figure;subplot(121);plot(q);ts=(handles.len-1)/handles.fs;t=0:1/handles.fs:ts;fy=fft(q);f=(0:length(t)-1)*handles.fs/length(t)-handles.fs/2;subplot(122);plot(f,fftshift(abs(fy);handles.q=q;guidata(hObject,handles);function nois_Callback(hObject, eventdata, handles)function gaosi_Callback(hObject, eventdata, handles)function bai_Callback(hObject, eventdata, handles)nos=0.001*rand(handles.len,1);handles.nos=handles.y+nos;subplot(121);plot(