MATLAB 带通滤波器功能实现回调函数

function bandpass_pushbutton10_Callback(~, eventdata, handles) % hObject handle to bandpass_pushbutton10 (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 = [2000,3000]; % fs = [1800,3300]; % Wp = 2fp/Fs; % Ws = 2fs/Fs; % if(Wp >= 1) % Wp = 0.99; % end % if(Ws >= 1) % Ws = 0.99;
% end fp = get(handles.edit3,'string');
fp = str2double(fp)2; fs = get(handles.edit4,'string');
fs = str2double(fs)2; if get(handles.radiobutton1,'value') [n, Wn]=buttord(Wp,Ws, 2, 15); [b, a]=butter(n, Wn,'bandpass'); axes(handles.axes3); [h,w]=freqz(b,a); plot(w/piFs/2,abs(h)); x1=filter(b,a,x1); elseif get(handles.radiobutton4,'value') b2=fir1(30,[fp/Fs fs/Fs],boxcar(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton5,'value') b2=fir1(30,[fp/Fs fs/Fs], triang(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton6,'value') b2=fir1(30,[fp/Fs fs/Fs],hamming(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton7,'value') b2=fir1(30,[fp/Fs fs/Fs],hanning(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton8,'value') b2=fir1(30,[fp/Fs fs/Fs],blackman(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(abs(h))); x1=fftfilt(b2,x1); elseif get(handles.radiobutton9,'value') b2=fir1(30,[fp/Fs fs/Fs],kaiser(31)); axes(handles.axes3); [h,w]=freqz(b2, 1,512); plot(w/piFs/2,20log(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;

% 该函数是用于实现带通滤波器功能的回调函数，即当用户点击GUI界面上的“带通滤波”按钮时，会执行该函数。

% 第1-5行是自动生成的注释，不需要解释。

% 第7-10行是定义全局变量，用于在函数内外传递参数。

% 第12-13行是将当前数据保存到x1中，以便进行滤波操作。

% 第15-23行是根据用户在GUI界面上选择的滤波器类型和截止频率，设计滤波器，同时绘制滤波器的幅频响应曲线。如果用户选择的是Butterworth滤波器，则使用butter函数设计滤波器；如果用户选择的是FIR滤波器，则使用fir1函数设计滤波器。绘制幅频响应曲线使用freqz函数，同时将结果绘制在GUI界面的axes3中。

% 第25-42行是根据选择的滤波器类型，对输入信号进行滤波操作。使用filter函数对信号进行IIR滤波，使用fftfilt函数对信号进行FIR滤波。滤波后得到的结果保存在x1中。

% 第44-48行是绘制滤波后的信号波形图，将结果绘制在GUI界面的axes5中。

% 第50-60行是绘制滤波后的信号频谱图，将结果绘制在GUI界面的axes6中。使用fft函数计算FFT，使用fftshift函数将结果移动到中心，然后绘制幅值谱曲线。

% 第62行是将滤波后的结果保存到全局变量x2中，以便在其他函数中使用。