基于matlab的肌电信号处理程序Word下载.docx

基于matlab的肌电信号处理程序Word下载.docx

《基于matlab的肌电信号处理程序Word下载.docx》由会员分享，可在线阅读，更多相关《基于matlab的肌电信号处理程序Word下载.docx（20页珍藏版）》请在冰豆网上搜索。

Value'

...）createsanewUNTITLEDorraisesthe

%existingsingleton*.Startingfromtheleft,propertyvaluepairsare

%appliedtotheGUIbeforeuntitled_OpeningFcngetscalled.An

%unrecognizedpropertynameorinvalidvaluemakespropertyapplication

%stop.Allinputsarepassedtountitled_OpeningFcnviavarargin.

%*SeeGUIOptionsonGUIDE'

sToolsmenu.Choose"

GUIallowsonlyone

%instancetorun（singleton）"

.

%Seealso:

GUIDE,GUIDATA,GUIHANDLES

%Edittheabovetexttomodifytheresponsetohelpuntitled

%LastModifiedbyGUIDEv2.508-Nov-201412:

45:

18

%Begininitializationcode-DONOTEDIT

gui_Singleton=1;

gui_State=struct（'

gui_Name'

mfilename,...

'

gui_Singleton'

gui_Singleton,...

gui_OpeningFcn'

@untitled_OpeningFcn,...

gui_OutputFcn'

@untitled_OutputFcn,...

gui_LayoutFcn'

[],...

gui_Callback'

[]）;

ifnargin&

&

ischar（varargin{1}）

gui_State.gui_Callback=str2func（varargin{1}）;

end

ifnargout

[varargout{1:

nargout}]=gui_mainfcn（gui_State,varargin{:

}）;

else

gui_mainfcn（gui_State,varargin{:

%Endinitializationcode-DONOTEDIT

%---Executesjustbeforeuntitledismadevisible.

functionuntitled_OpeningFcn（hObject,eventdata,handles,varargin）

%Thisfunctionhasnooutputargs,seeOutputFcn.

%hObjecthandletofigure

%eventdatareserved-tobedefinedinafutureversionofMATLAB

%handlesstructurewithhandlesanduserdata（seeGUIDATA）

%varargincommandlineargumentstountitled（seeVARARGIN）

%Choosedefaultcommandlineoutputforuntitled

handles.output=hObject;

%Updatehandlesstructure

guidata（hObject,handles）;

%UIWAITmakesuntitledwaitforuserresponse（seeUIRESUME）

%uiwait（handles.figure1）;

%---Outputsfromthisfunctionarereturnedtothecommandline.

functionvarargout=untitled_OutputFcn（hObject,eventdata,handles）

%varargoutcellarrayforreturningoutputargs（seeVARARGOUT）;

%Getdefaultcommandlineoutputfromhandlesstructure

varargout{1}=handles.output;

%---Executesonbuttonpressinpushbutton1.

functionpushbutton1_Callback（hObject,eventdata,handles）

%hObjecthandletopushbutton1（seeGCBO）

globals;

globala;

globalt;

globalcount;

%count=1;

[filename,filepath]=uigetfile（'

*.txt'

Ñ

¡

Ô

ñ

Î

Ä

¼

þ

'

）;

%Ñ

Ê

ý

¾

Ý

str=[filepathfilename];

s=load（str）;

a=s（:

7）;

t=s（:

1）;

axes（handles.axes1）;

plot（t,a）;

xlabel（'

±

ä

ylabel（'

Á

¿

title（'

Ê

Ð

Å

º

ifcount==2

globals2;

globala2;

globalt2;

[filename,filepath]=uigetfile（'

s2=load（str）;

a2=s2（:

t2=s2（:

plot（handles.axes2,t2,a2）;

end

count=1;

%---Executesonbuttonpressinpushbutton2.

functionpushbutton2_Callback（hObject,eventdata,handles）

%hObjecthandletopushbutton2（seeGCBO）

globaly;

%globalIR;

%globalSORH;

y=fft（a,20000）;

M=length（a）;

N=length（y）;

p=size（a）;

%s1=a（1:

20000）;

[CL]=wavedec（a,3,'

db5'

cA3=appcoef（C,L,'

3）;

cD1=detcoef（C,L,1）;

cD2=detcoef（C,L,2）;

cD3=detcoef（C,L,3）;

thr1=thselect（cD1,'

rigrsure'

thr2=thselect（cD2,'

thr3=thselect（cD3,'

TR=[thr1,thr2,thr3];

SORH='

s'

;

[XC,CXC,LXC,PERFO,PERF2]=wdencmp（'

lvd'

a,...

3,TR,SORH）;

L=p

（2）;

x=a;

h=XC;

F=0;

M=0;

forii=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）;

plot（XC（1:

10000））;

È

¥

ë

ó

%plot（handles.axes2,XC（1:

20000））;

%---Executesonbuttonpressinpushbutton3.

functionpushbutton3_Callback（hObject,eventdata,handles）

%hObjecthandletopushbutton3（seeGCBO）

global