【正文】 47。181。196。249。198。181。204。216。208。212。199。207。223。 grid on。title(39。176。205。204。216。206。214。203。185。181。205。205。168。202。253。214。194。203。178。168。198。247。39。)。xlabel(39。Frequency/Hz39。)。ylabel(39。Magnitude39。)。 else if t==2T=str2double(get(,39。String39。))。 %201。232。214。195。178。201。209。249。214。220。198。218。fs=1/T。 %178。201。209。249。198。181。194。202。206。170。214。220。198。218。181。185。202。253。Wp= str2double(get(,39。String39。))。 Ws= str2double(get(,39。String39。))。 %201。232。214。195。185。233。210。187。187。175。205。168。180。248。186。205。232。180。248。189。216。214。185。198。181。194。202。Ap= str2double(get(,39。String39。))。As= str2double(get(,39。String39。))。 %201。232。214。195。205。168。180。248。238。180。243。186。205。238。208。161。203。165。188。245。 [N,Wc]=buttord(Wp,Ws,Ap,As,39。s39。)。 %181。247。211。195。butter186。175。202。253。200。182。168。176。205。204。216。206。214。203。185。194。203。178。168。198。247。189。202。253。[B,A]=butter(N,Wc, 39。high39。,39。s39。)。 %181。247。211。195。butter186。175。202。253。201。232。188。198。176。205。204。216。206。214。203。185。194。203。178。168。198。247。W=linspace(0,pi,400*pi)。 %214。184。182。168。210。187。182。206。198。181。194。202。214。181。 [D,C]=impinvar(B,A,fs)。 %181。247。211。195。194。246。179。229。207。236。211。166。178。187。177。228。168。 Hz=freqz(D,C,W)。 %181。187。216。198。181。194。202。207。236。211。166。 plot(W/pi,abs(Hz)/abs(Hz(1)))。 %187。230。179。246。176。205。204。216。206。214。203。185。202。253。214。181。205。205。168。194。203。178。168。198。247。181。196。249。198。181。204。216。208。212。199。207。223。 grid on。title(39。176。205。204。216。206。214。203。185。184。223。205。168。202。253。214。194。203。178。168。198。247。39。)。xlabel(39。Frequency/Hz39。)。ylabel(39。Magnitude39。)。else if t==3T=str2double(get(,39。String39。))。 %201。232。214。195。178。201。209。249。214。220。198。218。 c1=str2double(get(,39。String39。))。b1=str2double(get(,39。String39。))。c2=str2double(get(,39。String39。))。b2=str2double(get(,39。String39。))。fs=1/T。 %178。201。209。249。198。181。194。202。206。170。214。220。198。218。181。185。202。253。Wp=[c1 c2]。Ws=[b1 b2]。 %201。232。214。195。185。233。210。187。187。175。205。168。180。248。186。205。232。180。248。189。216。214。185。198。181。194。202。Ap= str2double(get(,39。String39。))。As= str2double(get(,39。String39。))。 %201。232。214。195。205。168。180。248。238。180。243。186。205。238。208。161。203。165。188。245。 [N,Wc]=buttord(Wp,Ws,Ap,As,39。s39。)。 %181。247。211。195。butter186。175。202。253。200。182。168。176。205。204。216。206。214。203。185。194。203。178。168。198。247。189。202。253。[B,A]=butter(N,Wc,39。bandpass39。,39。s39。)。 %181。247。211。195。butter186。175。202。253。201。232。188。198。176。205。204。216。206。214。203。185。194。203。178。168。198。247。W=linspace(0,pi,400*pi)。 %214。184。182。168。210。187。182。206。198。181。194。202。214。181。 [D,C]=impinvar(B,A,fs)。 %181。247。211。195。194。246。179。229。207。236。211。166。178。187。177。228。168。 Hz=freqz(D,C,W)。 %181。187。216。198。181。194。202。207。236。211。166。 plot(W/pi,abs(Hz))。 %187。230。179。246。176。205。204。216。206。214。203。185。202。253。214。184。223。205。168。194。203。178。168。198。247。181。196。249。198。181。204。216。208。212。199。207。223。grid on。title(39。176。205。204。216。206。214。203。185。180。248。205。168。202。253。214。194。203。178。168。198。247。39。)。xlabel(39。Frequency/Hz39。)。ylabel(39。Magnitude39。)。else if t==4 T=str2double(get(,39。String39。))。 fs=1/T。c1=str2double(get(,39。String39。))。b1=str2double(get(,39。String39。))。c2=str2double(get(,39。String39。))。b2=str2double(get(,39。String39。))。Wp=[c1 c2]。Ws=[b1 b2]。 %201。232。214。195。185。233。210。187。187。175。205。168。180。248。186。205。232。180。248。189。216。214。185。198。181。194。202。Ap= str2double(get(,39。String39。))。As= str2double(get(,39。String39。))。 %201。232。214。195。205。168。180。248。238。180。243。186。205。238。208。161。203。165。188。245。 [N,Wc]=buttord(Wp,Ws,Ap,As,39。s39。)。 %181。247。211。195。butter186。175。202。253。200。182。168。176。205。204。216。206。214。203。185。194。203。178。168。198。247。189。202。253。[B,A]=butter(N,Wc,39。stop39。,39。s39。)。 %181。247。211。195。butter186。175。202。253。201。232。188。198。176。205。204。216。206。214。203。185。194。203。178。168。198。247。W=linspace(0,pi,400*pi)。 %214。184。182。168。210。187。182。206。198。181。194。202。214。181。 [D,C]=impinvar(B,A,fs)。 %181。247。211。195。194。246。179。229。207。236。211。166。178。187。177。228。168。 Hz=freqz(D,C,W)。 %181。187。216。198。181。194。202。207。236。211。166。 plot(W/pi,abs(Hz))。 %187。230。179。246。176。205。204。216。206。214。203。185。202。253。214。184。223。205。168。194。203。178。168。198。247。181。196。249。198。181。204。216。208。212。199。207。223。grid on。title(39。176。205。204。216。206。214。203。185。?161。232。202。253。214。194。203。178。168。198。247。39。)。xlabel(39。Frequency/Hz39。)。ylabel(39。Magnitude39。)。 end end endend % Executes on button press in radiobutton1.function radiobutton1_Callback(hObject, eventdata, handles)% hObject handle to radiobutton1 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,39。Value39。) returns toggle state of radiobutton1 % Executes on selection change in kj1. % Executes during object creation, after setting all properties.function kj1_CreateF(hObject, eventdata, handles)% hObject handle to kj1 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles empty handles not created until after all CreateFs called % Hint: popupmenu controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc amp。amp。 isequal(get(hObject,39。BackgroundColor39。), get(0,39。defaultUicontrolBackgroundColor39。)) set(hObject,39。BackgroundColor39。,39。white39。)。end % 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) t=0::9。m=str2double(get(,39。String39。))。n=str2double(get(,39。String39。))。y1=sin(m.*t)。y2=sin(n.*t)。y=y1+y2。 % plot(f,abs(fx(n+1))*2/N)。 grid。 plot(t,y)。 grid。 % Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved to be defined