www.pudn.com > dspproject.rar > jiangzao.m, change:2016-10-23,size:25037b
function varargout = jiangzao(varargin)
% JIANGZAO MATLAB code for jiangzao.fig
% JIANGZAO, by itself, creates a new JIANGZAO or raises the existing
% singleton*.
%
% H = JIANGZAO returns the handle to a new JIANGZAO or the handle to
% the existing singleton*.
%
% JIANGZAO('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in JIANGZAO.M with the given input arguments.
%
% JIANGZAO('Property','Value',...) creates a new JIANGZAO or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before jiangzao_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to jiangzao_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 jiangzao
% Last Modified by GUIDE v2.5 22-Oct-2016 18:56:42
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @jiangzao_OpeningFcn, ...
'gui_OutputFcn', @jiangzao_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 jiangzao is made visible.
function jiangzao_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 jiangzao (see VARARGIN)
% Choose default command line output for jiangzao
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes jiangzao wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = jiangzao_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 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)
global fpath snr
[y,fs]=wavread(fpath);
fq=fs;%话音回放
n = length (y);
a=str2num(get(handles.edit16,'String'));
Ft=str2num(get(handles.edit6,'String'));;
Fp=str2num(get(handles.edit7,'String'));;
Fs=str2num(get(handles.edit8,'String'));;
wp=2*pi*Fp/Ft;
ws=2*pi*Fs/Ft;
fp=2*Ft*tan(wp/2);
fs=2*Fs*tan(wp/2);
[n11,wn11]=buttord(wp,ws,1,50,'s'); %求低通滤波器的阶数和截止频率
[b11,a11]=butter(n11,wn11,'s'); %求S域的频率响应的参数
[num11,den11]=bilinear(b11,a11,0.5); %利用双线性变换实现频率响应S域到Z域的变换
[h,w]=freqz(num11,den11); %根据参数求出频率响应
z11=filter(num11,den11,y);
sound(z11*a,fq);
n=length(z11);
Y=fft(z11,n); %傅里叶变换
N=length(z11)+1;
Y=fft(z11);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes4,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes2,z11);
snr=SNR_singlech(Y,z11); % 计算维纳滤波后的信噪比
set(handles.text22,'string',snr);%在静态文本框中显示信噪比
% --- 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)
global fpath
[y,fs]=wavread(fpath);
fq=fs;%话音回放
n = length (y)
a=str2num(get(handles.edit16,'String'));
fpl=str2num(get(handles.edit2,'String'));
fpu=str2num(get(handles.edit4,'String'));
fsl=str2num(get(handles.edit3,'String'));
fsu=str2num(get(handles.edit5,'String'));
wp3=[2*fpl/fs,2*fpu/fs];
ws3=[2*fsl/fs,2*fsu/fs];
rp=1;
as=100;
[N3,wp3]=ellipord(wp3,ws3,rp,as); %计算椭圆带通模拟滤波器的阶数和通带边界频率
[B3,A3]=ellip(N3,rp,as,wp3); %计算带通滤波器模拟滤波器系统函数系数
y3=filter(B3,A3,y);
sound(y3*a,fq);
n=length(y3);
Y=fft(y3,n); %傅里叶变换
N=length(y3)+1;
Y=fft(y3);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes4,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes2,y3);
% --- 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)
global fpath
[y,fs]=wavread(fpath);
fq=fs;%话音回放
n = length (y)
a=str2num(get(handles.edit16,'String'));
Ft=str2num(get(handles.edit13,'String'));;
Fp=str2num(get(handles.edit14,'String'));;
Fs=str2num(get(handles.edit15,'String'));;
wp1=tan(pi*Fp/Ft);%高通到低通滤波器参数转换?
ws1=tan(pi*Fs/Ft);
wp=1;
ws=wp1*wp/ws1;
[n13,wn13]=cheb1ord(wp,ws,1,100,'s');%求模拟的低通滤波器阶数和截止频率?
[b13,a13]=cheby1(n13,1,wn13,'s');%求S域的频率响应的参数?
[num,den]=lp2hp(b13,a13,wn13);%将S域低通参数转为高通的?
[num13,den13]=bilinear(num,den,0.5);%利用双线性变换实现频率响应S域到Z域转换?
[h,w]=freqz(num13,den13);
z11=filter(num13,den13,y);
sound(z11*a,fq);
n=length(z11);
Y=fft(z11,n); %傅里叶变换
N=length(z11)+1;
Y=fft(z11);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes4,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes2,z11);
snr=SNR_singlech(Y,z11); % 计算维纳滤波后的信噪比
set(handles.text22,'string',snr);%在静态文本框中显示信噪比
% --- 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)
ButtonName=questdlg('是否开始录音？','MATLAB quest','Yes','No','Yes');
switch ButtonName
case 'Yes'
prompt={'输入录音时间：','输入采样率'};%设置提示字符串
name='Enter trapeaia Data';%设置标题
numlines=1;%指定输入数据的行数
defAns={'6','44100'};%设定默认值
Resize='on';%设定对话框尺寸可调节
answer=inputdlg(prompt,name,numlines,defAns,'on');%创建输入对话框
duration = str2num(answer{1});
fs = str2num(answer{2});
h = waitbar(0,'录音中，请稍后...');
data=audiorecorder(fs,16,2);
recordblocking(data,duration);
y =getaudiodata(data,'int16');
wavwrite(y,fs,'d:\luyinzao.wav');
waitbar(1,h,'录音完毕');
end
% --- 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)
global fpath3
[filename, pathname] = uigetfile('*.wav','选择声音文件');
if isequal(filename,0)||isequal(pathname,0) %判断是否打开音频
return;
end;
fpath3=[pathname filename];
[y,fs]=wavread(fpath3);
audioplayer(y,fs);
n=length(y);
Y=fft(y,n); %傅里叶变换
N=length(y)+1;
Y=fft(y);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes3,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes1,y);
% --- 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)
main;
% --- Executes on button press in pushbutton7.
function pushbutton7_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton7 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global fpath3 %小波降噪
[s,fs]=wavread(fpath3);
N=length(s);
t=s';
[thr,sorh,keepapp]=ddencmp('den','wv',t);
[c,l]=wavedec(t,5,'sym6');
xd=wdencmp('gbl',c,l,'sym6',5,thr,sorh,keepapp);
sound(xd,2*fs);
%snr=SNR_singlech(s,xd); % 计算维纳滤波后的信噪比
%set(handles.text23,'string',snr);%在静态文本框中显示信噪比
n=length(xd);
Y=fft(xd,n); %傅里叶变换
N=length(xd)+1;
Y=fft(xd);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes4,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes2,xd);
function edit6_Callback(hObject, eventdata, handles)
% hObject handle to edit6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit6 as text
% str2double(get(hObject,'String')) returns contents of edit6 as a double
% --- Executes during object creation, after setting all properties.
function edit6_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit7_Callback(hObject, eventdata, handles)
% hObject handle to edit7 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit7 as text
% str2double(get(hObject,'String')) returns contents of edit7 as a double
% --- Executes during object creation, after setting all properties.
function edit7_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit7 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit8_Callback(hObject, eventdata, handles)
% hObject handle to edit8 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit8 as text
% str2double(get(hObject,'String')) returns contents of edit8 as a double
% --- Executes during object creation, after setting all properties.
function edit8_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit8 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit9_Callback(hObject, eventdata, handles)
% hObject handle to edit9 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit9 as text
% str2double(get(hObject,'String')) returns contents of edit9 as a double
% --- Executes during object creation, after setting all properties.
function edit9_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit9 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit2_Callback(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit2 as text
% str2double(get(hObject,'String')) returns contents of edit2 as a double
% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit3_Callback(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit3 as text
% str2double(get(hObject,'String')) returns contents of edit3 as a double
% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit4_Callback(hObject, eventdata, handles)
% hObject handle to edit4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit4 as text
% str2double(get(hObject,'String')) returns contents of edit4 as a double
% --- Executes during object creation, after setting all properties.
function edit4_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit5_Callback(hObject, eventdata, handles)
% hObject handle to edit5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit5 as text
% str2double(get(hObject,'String')) returns contents of edit5 as a double
% --- Executes during object creation, after setting all properties.
function edit5_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit13_Callback(hObject, eventdata, handles)
% hObject handle to edit13 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit13 as text
% str2double(get(hObject,'String')) returns contents of edit13 as a double
% --- Executes during object creation, after setting all properties.
function edit13_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit13 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit14_Callback(hObject, eventdata, handles)
% hObject handle to edit14 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit14 as text
% str2double(get(hObject,'String')) returns contents of edit14 as a double
% --- Executes during object creation, after setting all properties.
function edit14_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit14 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit15_Callback(hObject, eventdata, handles)
% hObject handle to edit15 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit15 as text
% str2double(get(hObject,'String')) returns contents of edit15 as a double
% --- Executes during object creation, after setting all properties.
function edit15_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit15 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in pushbutton8.
function pushbutton8_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton8 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
function edit16_Callback(hObject, eventdata, handles)
% hObject handle to edit16 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit16 as text
% str2double(get(hObject,'String')) returns contents of edit16 as a double
% --- Executes during object creation, after setting all properties.
function edit16_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit16 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in pushbutton9.
function pushbutton9_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton9 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global fpath3 %维纳滤波
[xx,fs]=audioread(fpath3);% 读入数据文件
xx=xx(:,1);
xx=xx-mean(xx); % mean:求矩阵各列均值 消除直流分量
x=xx/max(abs(xx)); % 幅值归一化
N=length(x); % 信号长度
time=(0:N-1)/fs; % 设置时间
Y=fft(x,N); %傅里叶变换
n=length(x)+1;
Y=fft(x);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
IS=.15; % 设置IS
% 调用WienerScalart96m_2函数做维纳滤波
output=WienerScalart96m_2(xx,fs,IS,0.12);
ol=length(output); % 把output补到与x等长
if ol<N
output=[output; zeros(N-ol,1)];
end
snr=SNR_singlech(x,output); % 计算维纳滤波后的信噪比
fp1=1200;
fs1=2000;
wp1=2*fp1/fs;%设定低通滤波器通带截止频率和阻带截止频率
ws1=2*fs1/fs;
rp=1;
as=100;
[N1,wp1]=ellipord(wp1,ws1,rp,as);%计算椭圆低通模拟滤波器的阶数和通带边界频率
[B,A]=ellip(N1,rp,as,wp1);%计算低通滤波器模拟滤波器系统函数系数
output1=filter(B,A,output);
M=length(output); % 信号长度
time1=(0:M-1)/fs; % 设置时间
Y2=fft(output,M); %傅里叶变换
m=length(M)+1;
Y2=fft(output);%做FFT变换，以下为绘图函数
Y2=abs(Y2);
Y2=Y2(1:M/2);
f2=(0:M-1)*fs/M;
f2=f2(1:M/2);
%subplot(2,2,1)
%plot(time,x,'k'); grid; axis tight;
%title('噪音语音波形'); ylabel('幅值'); xlabel('时间/s');
%subplot(2,2,2)
%plot(time,output,'r'); grid; axis tight;
%title('去噪后语音波形'); ylabel('幅值'); xlabel('时间/s');
sound(output1,fs);
plot(handles.axes2,output);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes4,Y2);
set(handles.text22,'string',snr);
function edit17_Callback(hObject, eventdata, handles)
% hObject handle to edit17 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit17 as text
% str2double(get(hObject,'String')) returns contents of edit17 as a double
% --- Executes during object creation, after setting all properties.
function edit17_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit17 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes during object creation, after setting all properties.
function text22_CreateFcn(hObject, eventdata, handles)
% hObject handle to text22 (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 pushbutton11.
function pushbutton11_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton11 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global fpath3 %完美去风声
[Original,fs]=audioread(fpath3);%修改文件路径
ts=1/fs;
N=length(Original)-1;
t=0:1/fs:N/fs;
Nfft=N;
df=fs/Nfft;
fk=(-Nfft/2:Nfft/2-1)*df;
a1=1,a2=-1,b1=1,b2=-1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
SoundLeft=Original(:,1);
SoundRight=Original(:,2);
SoundLeft_f=ts*fftshift(fft(SoundLeft,N));
SoundRight_f=ts*fftshift(fft(SoundRight,N));
% Sound=SoundLeft-SoundRight;
NewLeft=a1*SoundLeft+a2*SoundRight;
NewRight=b1*SoundLeft+b2*SoundRight;
Sound(:,1)=NewLeft;
Sound(:,2)=NewRight;
sound(Sound,fs);
n=length(Sound);
Y=fft(Sound,n); %傅里叶变换
N=length(Sound)+1;
Y=fft(Sound);%做FFT变换，以下为绘图函数
Y=abs(Y);
Y=Y(1:N/2);
f=(0:N-1)*fs/N;
f=f(1:N/2);
plot(handles.axes4,Y);
xlabel('Hz');
ylabel('幅值');
plot(handles.axes2,Sound);