基于matlab的肌电信号处理程序Word文档下载推荐.docx
- 文档编号:20536934
- 上传时间:2023-01-23
- 格式:DOCX
- 页数:20
- 大小:17.25KB
基于matlab的肌电信号处理程序Word文档下载推荐.docx
《基于matlab的肌电信号处理程序Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《基于matlab的肌电信号处理程序Word文档下载推荐.docx(20页珍藏版)》请在冰豆网上搜索。
%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)
%a=EMG3(1:
points,2);
%¶
½
ø
ì
Ë
Ù
Fourier±
»
fs=2000;
mag=abs(y);
%Ç
µ
Ã
Õ
·
ù
f=(0:
N-1)/N*fs;
%Æ
Â
ò
axes(handles.axes1);
plot(f,mag);
%plot(handles.axes1,f,mag);
%»
æ
Ö
Æ
Í
xlabel('
¼
ç
N=20000'
ifcount==2
globaly2;
y2=fft(a2,20000);
N=length(y2);
mag=abs(y2);
end
count=1;
%---Executesonbuttonpressinpushbutton4.
functionpushbutton4_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton4(seeGCBO)
y1=fft(XC,20000);
N=length(y1);
mag1=abs(y1);
plot(f,mag1);
%---Executesonbuttonpressinpushbutton5.
functionpushbutton5_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton5(seeGCBO)
%Æ
plot(f,angle(y));
Ï
à
%plot(handles.axes1,f,angle(y));
Ì
Ø
plot(f,angle(y2));
%---Executesonbuttonpressinpushbutton6.
functionpushbutton6_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton6(seeGCBO)
plot(f,angle(y1));
%---Executesonbuttonpressinpushbutton7.
functionpushbutton7_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton7(seeGCBO)
try
delete(allchild(handles.axes1));
delete(allchild(handles.axes2));
set(handles.edit1,'
string'
set(handles.edit3,'
set(handles.edit4,'
set(handles.edit5,'
set(handles.edit6,'
%---Executesonbuttonpressinpushbutton8.
functionpushbutton8_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton8(seeGCBO)
power1=(mag.^2)/20000;
plot(f,power1);
¹
¦
×
globals2;
%---Executesonbuttonpressinpushbutton9.
functionpushbutton9_Callback(hObject,eventdata,handles)
%hObjecthandletopushbutton9(seeGCBO)
s=fft(a,2000);
mag=abs(y1);
power1=(mag.^2)/2000
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 matlab 电信号 处理 程序