《基于matlab的肌电信号处理程序:.doc》由会员分享,可在线阅读,更多相关《基于matlab的肌电信号处理程序:.doc(14页珍藏版)》请在金锄头文库上搜索。
1、基于matlab的肌电信号处理程序:function varargout = untitled(varargin)% UNTITLED M-file for untitled.fig% UNTITLED, by itself, creates a new UNTITLED or raises the existing% singleton*.% H = UNTITLED returns the handle to a new UNTITLED or the handle to% the existing singleton*.% UNTITLED(CALLBACK,hObject,eventD
2、ata,handles,.) calls the local% function named CALLBACK in UNTITLED.M with the given input arguments.% UNTITLED(Property,Value,.) creates a new UNTITLED or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before untitled_OpeningFcn gets called. An
3、% unrecognized property name or invalid value makes property application% stop. All inputs are passed to untitled_OpeningFcn via varargin.% *See GUI Options on GUIDEs Tools menu. Choose GUI allows only one% instance to run (singleton).% See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to m
4、odify the response to help untitled % Last Modified by GUIDE v2.5 08-Nov-2014 12:45:18 % Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, untitled_OpeningFcn, . gui_OutputFcn, untitled_OutputFcn, . gui_
5、LayoutFcn, , . gui_Callback, );if nargin & ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);end if nargout varargout1:nargout = gui_mainfcn(gui_State, varargin:);else gui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT % - Executes just before untitled is made
6、visible.function untitled_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
7、 arguments to untitled (see VARARGIN) % Choose default command line output for untitledhandles.output = hObject; % Update handles structureguidata(hObject, handles); % UIWAIT makes untitled wait for user response (see UIRESUME)% uiwait(handles.figure1); % - Outputs from this function are returned to
8、 the command line.function varargout = untitled_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 GUID
9、ATA) % Get default command line output from handles structurevarargout1 = 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
10、% handles structure with handles and user data (see GUIDATA)global s;global a;global t;global count;% count=1;filename,filepath=uigetfile(*.txt,);% str=filepath filename; s=load(str); a=s(:,7); t=s(:,1); axes(handles.axes1);plot(t,a);xlabel();ylabel();title(); if count=2 global s2;global a2;global t
11、2; filename,filepath=uigetfile(*.txt,); str=filepath filename; s2=load(str); a2=s2(:,7); t2=s2(:,1); plot(handles.axes2,t2,a2); end count=1; % - Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata res
12、erved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global a;global t;global s;global y; % global IR;% global SORH;y=fft(a,20000);M =length (a) ;N=length(y);p=size(a);% s1=a(1:20000);C L=wavedec(a,3,db5);cA3=appcoef(C,L,db5,3);cD1=detcoef(C,
13、L,1);cD2=detcoef(C,L,2);cD3=detcoef(C,L,3);thr1=thselect(cD1,rigrsure);thr2=thselect(cD2,rigrsure);thr3=thselect(cD3,rigrsure);TR=thr1,thr2,thr3;SORH=s;XC,CXC,LXC,PERFO,PERF2=wdencmp(lvd,a,. db5,3,TR,SORH);L=p(2);x=a;h=XC;F=0;M=0;for ii=1:L m(ii)=(x(ii)-y(ii)2; t(ii)=y(ii)2; f(ii)=t(ii)/m(ii); F=F+f
14、(ii); M=M+m(ii);end;SNR=10*log10(F);MSE=M/N;SM=SNR/MSE;axes(handles.axes2);plot(XC(1:10000);xlabel();ylabel();title();% plot(handles.axes2,XC(1:20000); % - Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata