基于某matlab的肌电信号处理程序.doc

《基于某matlab的肌电信号处理程序.doc》由会员分享，可在线阅读，更多相关《基于某matlab的肌电信号处理程序.doc（14页珍藏版）》请在三一办公上搜索。

1、word基于matlab的肌电信号处理程序：function varargout = untitled(varargin)% 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,eventData,handles,.) calls the local

2、% 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_OpeningF gets called. An% unrecognized property name or

3、invalid value makes property application% stop. All inputs are passed to untitled_OpeningF 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 modify the response to help untitled

4、% 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_OpeningF, untitled_OpeningF, . gui_OutputF, untitled_OutputF, . gui_LayoutF, , . gui_Callback, );if nargin & isch

5、ar(varargin1) gui_State.gui_Callback = str2func(varargin1);endif nargout varargout1:nargout = gui_mainf(gui_State, varargin:);else gui_mainf(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before untitled is made visible.function untitled_OpeningF(hObject, eventdata

6、, handles, varargin)% This function has no output args, see OutputF.% 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 mand line arguments to untitled (see VARARGIN)% Choose default mand l

7、ine 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 the mand line.function varargout = untitled_OutputF(hObject, event

8、data, 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 mand line output from handles structurevarargout1 = ha

9、ndles.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 s;global a

10、;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 t2; filename,filepath=uigetfile(*.txt,); str=filepath filename; s2=load(str);

11、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 reserved - to be defined in a future version of MATLAB% handles structure with h

12、andles 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,L,1);cD2=detcoef(C,L,2);cD3=detcoef(C,L,3);thr1=thselect(cD1,rigrsure);thr2=th

13、select(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(ii); M=M+m(ii);end;SNR=10*log10(F);MSE=M/N;SM=SNR/MSE;axes(handles.axes2);plo

14、t(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 reserved - to be defined in a future version of MATLAB% handles structu

15、re with handles and user data (see GUIDATA)global a;global t;global y;global count;%a=EMG3(1:points,2);y=fft(a,20000); %Fourierfs=2000;N=length(y);mag=abs(y);%Fourierf=(0:N-1)/N*fs; %axes(handles.axes1);plot(f,mag);% plot(handles.axes1,f,mag);% xlabel();ylabel();title(N=20000);if count=2global a2;gl

16、obal t2;global y2;%a=EMG3(1:points,2);y2=fft(a2,20000); %Fourierfs=2000;N=length(y2);mag=abs(y2);%Fourierf=(0:N-1)/N*fs; %axes(handles.axes2);plot(f,mag); endcount=1;% - Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4

17、(see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global a;global t;% global IR;% global SORH;C L=wavedec(a,3,db5);cA3=appcoef(C,L,db5,3);cD1=detcoef(C,L,1);cD2=detcoef(C,L,2);cD3=detcoef(C,L,3);thr1=thselect(cD1,r

18、igrsure);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);y1=fft(XC,20000);fs=2000;N=length(y1);mag1=abs(y1);f=(0:N-1)/N*fs;axes(handles.axes2);plot(f,mag1);xlabel();ylabel();title(N=20000);% - Executes on button p

19、ress 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 a;global t;global y;global count;fs=2000;N=leng

20、th(y);f=(0:N-1)/N*fs;%axes(handles.axes1);plot(f,angle(y);xlabel();ylabel();title();% plot(handles.axes1,f,angle(y);%if count=2global a2;global t2;global y2;fs=2000;N=length(y2);f=(0:N-1)/N*fs;%axes(handles.axes2);plot(f,angle(y2);xlabel();ylabel();title();endcount=1;% plot(handles.axes1,f,angle(y);

21、% - 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)global a;global t;C L=wavedec(

22、a,3,db5);cA3=appcoef(C,L,db5,3);cD1=detcoef(C,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);y1=fft(XC,20000);fs=2000;N=length(y1);f=(0:N-1)/

23、N*fs;axes(handles.axes2);plot(f,angle(y1);xlabel();ylabel();title();% - 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 wi

24、th handles and user data (see GUIDATA)try delete(allchild(handles.axes1); delete(allchild(handles.axes2); set(handles.edit1,string,); set(handles.edit1,string,); set(handles.edit3,string,); set(handles.edit4,string,); set(handles.edit5,string,); set(handles.edit6,string,);end% - Executes on button p

25、ress 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)global a;global t;global s;y=fft(a,20000);fs=2000;N=le

26、ngth(y);mag=abs(y);f=(0:N-1)/N*fs;power1=(mag.2)/20000;axes(handles.axes1);plot(f,power1);xlabel();ylabel();title();global count;global a2;global t2;global s2;if count=2y2=fft(a2,20000);fs=2000;N=length(y2);mag=abs(y2);f=(0:N-1)/N*fs;power1=(mag.2)/20000;axes(handles.axes2);plot(f,power1);xlabel();y

27、label();title();endcount=1;% - 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)glo

28、bal a;global t;global s;s=fft(a,2000);C L=wavedec(a,3,db5);cA3=appcoef(C,L,db5,3);cD1=detcoef(C,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

29、);y1=fft(XC,20000);fs=2000;N=length(y1);mag=abs(y1);f=(0:N-1)/N*fs;power1=(mag.2)/20000;%axes(handles.axes2);plot(f,power1);xlabel();ylabel();title();% - Executes on button press in pushbutton10.function pushbutton10_Callback(hObject, eventdata, handles)% hObject handle to pushbutton10 (see GCBO)% e

30、ventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)%global a;set(handles.edit1,string,mean(a);function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a futu

31、re version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,String) returns contents of edit1 as text% str2double(get(hObject,String) returns contents of edit1 as a double% - Executes during object creation, after setting all properties.function edit1_CreateF

32、(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFs called% Hint: edit controls usually have a white background on Windows.% See ISPC and PUTER.if ispc & iseq

33、ual(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);end% - 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

34、future version of MATLAB% handles structure with handles and user data (see GUIDATA)global a;set(handles.edit3,string,sqrt(var(a);function edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles struct

35、ure 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_CreateF(hObject, eventdata, handles)% hObje

36、ct handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFs called% Hint: edit controls usually have a white background on Windows.% See ISPC and PUTER.if ispc & isequal(get(hObject,BackgroundColor), ge

37、t(0,defaultUicontrolBackgroundColor) set(hObject,BackgroundColor,white);end% - Executes on button press in pushbutton12.function pushbutton12_Callback(hObject, eventdata, handles)% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles st

38、ructure with handles and user data (see GUIDATA)global a;set(handles.edit4,string,var(a);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 GUIDA

39、TA)% 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_CreateF(hObject, eventdata, handles)% hObject handle to edit4 (see GCBO)% eventdata

40、reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFs called% Hint: edit controls usually have a white background on Windows.% See ISPC and PUTER.if ispc & isequal(get(hObject,BackgroundColor), get(0,defaultUicontrolBackgroundColor) set(

41、hObject,BackgroundColor,white);end% - Executes on button press in pushbutton13.function pushbutton13_Callback(hObject, eventdata, handles)% hObject handle to pushbutton13 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see G

42、UIDATA)%IEMGglobal a;set(handles.edit5,string,mean(abs(a);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 s