完整版肌电信号的识别

1、燕山大学数字信号处理课程设计燕山大学课程设计说明书题目：肌电信号分析及动作识别学院（系）：电气工程学院年级专业：10级仪表三班学号：学生姓名：指导教师：教师职称：教授讲师燕山大学数字信号处理课程设计电气工程学院课程设计任务书基层教学单位：自动化仪表系指导教师：谢平张淑清学号学生姓名（专业）班级10级仪表三班设计题目肌电信号分析及动作识别设计技术参数1、独立完成设计任务。2、编程，上机调试。3、连接硬件实验线路，实现所要求的功能。4、完成设计，提交课程设计报告。设计要求1、利用肌电传感器采集人体肌电信号，利用Matlab软件对肌电信号进行时域和频域（包括FFT,功率谱，倒谱）分析，计算均值、标准

2、差、方差、积分肌电值IEMG、均方根RMS等时域指标参数和平均功率频率MPF、中值频率MF等频域指标。比较不同动作下肌电信号特征，根据设定阈值给出动作状态，设计相应的信号分析及显示界面。2、基于肌电信号分析结果，通过串口发送命令给单片机系统，根据肌电信号动作状态控制相应的数码管或LED发光显示。3、扩展：也可通过动作模式驱动电机转动或其他控制输出模块。工作量软件编程与硬件调试相结合，绘制设计流程图，编制相应软件界面，实现单片机控制与信号处理任务的综合应用参考资料1）微型计算机控制系统赖寿宏，机械工业出版社（教材）2）单片机及应用李大友，高等教育出版社（教材）3）信号处理原理及应用谢平等机械工业

3、出版社（教材）4）«Matlab程序设计及其在信号处理中的应用聂祥飞等西南交通大学出版社周次A周第二周应完成内容熟悉伟福单片机编程环境，调试单片机各基本功能模块；熟悉matlab信号处理工具箱，信号处理系统基本功能模块学习和调试单片机系统与信号处理系统综合进行硬件调试，撰写课程设计报告指导教师签字基层教学单位主任签字2燕山大学数字信号处理课程设计第一章摘要2第二章系统总体设计方案3第三章肌电信号的时域参数处理及其分析4第四章肌电信号的频域处理方法及其分析73.1 FFT分析73.2 功率谱分析83.3 倒谱分析93.4 平均功率频率MPFffi中值频率10第五章Matlab程序及GU

4、I11第六章系统整体调试及结果说明24第七章学习心得24参考文献25燕山大学数字信号处理课程设计第一章摘要肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单,病人易接受，有着广泛的应用前景。主要应用领域有：一，仿生学。提出肌肉

5、生理模型来判别肌肉的动作以来，电子假肢的研究进入了新的发展时期，过去电子假肢的控制靠使用者人为开关和选择运动模式来完成，现在则可通过检测人体残肢表面肌电信号提取出肢体的动作特征，来自动控制假肢运动，利用残肢表面肌电信号的肌电假肢研制在国内外都取得较大进展。二，康复工程。如利用表面肌电信号提取出的特征作为功能性电刺激的控制信号，帮助瘫痪的肢体恢复运动功能。通过检测表面肌电信号，并将其作为反馈信号提供给病人和医生，便于进行合理的治疗和训练。三，运动医学。表面肌电信号在运动医学中也可发挥重要作用，通过检测运动员运动时的表面肌电信号，及时反映出肌肉的疲劳和兴奋状态，有助于建立科学的训练方法。本次课程设

6、计的主要任务就是对微弱的肌电信号进行时域和频域的处理及分析，运用数字处理及matlab的知识进行“屈”和“伸”动作识别。然后通过串口将数据发送到单片机下行微机进行显示。燕山大学数字信号处理课程设计第二章系统总体设计方案根据课程设计要求在上微机利用matble分析肌电信号并处理，基于肌电信号分析结果，通过串口发送命令给单片机系统，根据肌电信号动作状态控制相应的数码管显示。并增加了扩展模块，通过动作模式驱动电机转动或其他控制输出模块。通过分析上位机matlab中对信号处理的结果，我们可以得到一系列的信号特征值，其中我们选取了具有代表意义积分肌电值来进行处理，并给出对于屈伸动作的阈值。通过对阈值的判

7、断，使数码管显示积分肌电值，使点阵模块显示相应的“屈”和“伸”字样，使电机根据动作进行正反转，蜂鸣器在“屈”动作是发出鸣响。最后，将上微机的处理数据通过串口通信发送到下微机显示，得到动作识别的要求。燕山大学数字信号处理课程设计第三章肌电信号的时域参数处理及其分析(1)均值：对于一组随机变量来说，均值是一个很重要的数值特征，用来描述一组变量的平均水平。具严格的数学定义非常简单，就是一个随机变量关于概率测度的积分。因此，在此处，均值表示肌电信号的平均水平。公式如下：N(2)方差：方差是各个数据与平均数之差的平方的平均数。在数理统计和概率论中，方差用来度量随机变量和其数学期望(即均值)之间的偏离程度

8、。在此处，研究信号的随机变量和均值之间的偏离程度有着很重要的意义。其求解公式如下：X(y产(3)标准差：标准差也称均方差，是各数据偏离平均数的距离的平均数，它是离均差平方和平均后的方根，用u表示。标准差是方差的算术平方根。标准差能反映一个数据集的离散程度。所以能很好的反映肌电信号的离散程度。求解公式如下：1强=万£%八不二i(4)积分肌电值IEMG:积分肌电值就是对所有信号取绝对值后尽心均值的求解，由于对肌电信号直接求均值，均值近似为零，无法表征信号间的差异。若对肌电信号取绝对值后再进行均值运算后，均值恒大于零，因而可用于提取肌电信号的特征。公式如下：1V-1=讨6叫Ni=Q燕山大学

9、数字信号处理课程设计(5)均方根RMS均方根就是一组数据的平方和除以数据的个数再开方，均方根是最理想的平方滤波方式的典型，让滤波更平滑，更大限度的滤掉噪声。因此，对肌电信号求均方根，可以滤除信号中的噪声，使滤波后的信号更平滑、更明显。公式如下：RMS.士ICTEMG2(f)dt(6)原始信号的时域及上述参数值的Matlab程序allclear;closea=load('qu.txt');/('shen.txt');N=10000;b=0:N-1;axis(0,10000,-1,1);plot(b,a);xlabel('时间(s)');ylabel

10、('被测变量y');gridon;fprintf('n数据基本信息：n')printf('均值=%7.5fn',mean(a);fprintf('标准差=%7.5fn',sqrt(var(a);fprintf('方差=%7.5fn',var(a);fprintf('积分肌电值IEMG=%7.5fn',mean(abs(a);fprintf('均方根有效值RMS=%7.5fn',sqrt(mean(a.A2);屈信号数据基本信息：均值=0.03502标准差=0.05775方差=0.00

11、334积分肌电值IEMG=0.05437均方根RMS=0.17246燕山大学数字信号处理课程设计屈J原始信方10(10200C300040005Q00600070008000900010000时间闾JO.-F“酬翩夏樟伸信号数据基本信息：均值=0.00337标准差=0.24421方差=0.05964积分肌电值IEMG=0.12826均方根RMS=0.19993c伸)原始信号时间15)燕山大学数字信号处理课程设计第四章肌电信号的频域处理方法及其分析3.1 FFT分析：FFT是离散傅立叶变换的快速算法，可以将一个信号变换到频域。有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易

12、看出特征了。这就是很多信号分析采用FFT变换的原因。另外，FFT可以将一个信号的频谱提取出来。采样得到的数字信号FFT变换。N个采样点，经过FFT之后，就可以得到N个点的FFT结果。为了方便进行FFT运算，通常N取2的整数次方。假设采样频率为Fs,信号频率F,采样点数为No那么FFT之后结果就是一个为N点的复数。每一个点就对应着一个频率点。这个点的模值，就是该频率值下的幅度特性。而每个点的相位就是在该频率下的信号的相位。如果要要提高频率分辨率，就需要增加采样点数，也即采样时间。频率分辨率和采样时问是倒数关系。Matlab实现:clear;closea=load(y=fft(a,1024);fs

13、=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;figure;plot(f,mag);xlabel(ylabel(title(all;'qu.txt'）;/'shen.txt'%1024点傅立叶变换娴频谱频率）；幅值）；'（屈/伸）肌电幅频'）;figure;plot(f,angle(y);xlabel(ylabel(title(%相频谱频率）；相位）；'（屈/伸）肌电相频'）;gridon;燕山大学数字信号处理课程设计3.2 功率谱分析：功率谱是信号或噪声的自相关函数的傅里叶变换。如果一确

14、定信号平均功率为有限的，则该信号的自相关函数存在，如随机信号或随机噪声是由二阶随机平稳函数表示的，则其自相关函数存在。为了更好得描述能量信号、功率信号，我们引入能量谱密度和功率谱密度概念。能量谱密度、功率谱密度函数表示信号的能量、功率密度随频率变化的情况。通过研究功率谱密度，可以帮助了解信号的功率分布情况，确定信号的频带等。信号的功率谱即上述FFT分析后，幅频值mag的平方再除以2得到。功率密度谱虽然描述了随机信号的功率在各个不同频率上的分布，但因为它仅与幅度频谱有关，没有相位信息，所以从已知功率谱还难以完整地恢复原来的功率信号。燕山大学数字信号处理课程设计Matlab实现：a=loadfqu

15、.txt'y=fft(a,1024);fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power1=(mag.A2)/2;plot(f,power1););/('shen.txt')%做1024点傅里叶变换%周期图法求功率谱xlabel(ylabel(title('频谱);'功率谱);,(屈/伸)肌电信号功率谱);500频谱3000期百2UOQ与1000(届)肌电信号功率谱4000(1$)肌电信号功率谱10地助日5oa频谱100C10003.3 倒谱分析：倒谱是信号的傅里叶变换谱经对数运算后再进行的傅里叶反变换

16、或者功率谱的对数值的逆傅氏变换称为倒谱。倒频谱函数C(q)(powercepstrum)其数学表达式为：C(q)=|IF(log(s(f)-2。其中，s(f)是信号s(t)的傅里叶变换，log()为取对数，IF为逆傅里叶变换。Matlab实现：前100点的倒谱变化比较明显，所以取前100点波形图。a=load('qu.txty=fft(a,1024);fs=1000;N=length(y););/('shen.txt')%ft1024点傅里叶变换燕山大学数字信号处理课程设计%求倒谱mag=abs(y);f=(0:N-1)/N*fs;w=rceps(a);n=1:100;

17、plot(n,w(1:100);xlabel('时间);ylabel('倒谱);title('(屈/伸)肌电信号倒谱);gridon;（屈）肌电信骂倒谱y=fft(a,1024);fs=1000;）;%10000点傅里叶变换103.4平均功率频率MPF和中值频率MF:平均功率频率是总功率除以总时间中值频率是各个时间段的功率的平均值。公式如下：Matlab实现：a=load('qu.txt燕山大学数字信号处理课程设计N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/2;ss=sum(power);M2=0.5

18、*ss;df=fs/N;M1=0.5*df*(sum(power(1:N-1)+sum(power(2:N);MPF=M1/M2;MF=M2/2;第五章Matlab程序及GUIMATLA是一种面向工程和科学计算的交互式计算软件，它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个简单易用的交互式工作环境中。同时由于MATLA是一个数据分析和处理功能十分强大的工程实用软件，它的信号处理与分析工具箱为肌电信号分析提供了十分丰富的功能函数，利用这些功能函数可以快捷而又方便地完成肌电信号的处理和分析，再通过GUI界面设计，显示信号处理后的结果，方便简明。打开GUI设计界面，根据要求添加相应的按钮，

19、并对按钮进行定义，调整界面整体布局，将界面运行，会产生相应的程序。在程序中，寻找相应的按钮定义，在定义下添加该按钮要实现的功能函数。然后对整个程序进行调试验证。得到可视化GUI界面及功能的实现。11燕山大学数字信号处理课程设计functionvarargout=untitled20(varargin)%UNTITLED20M-fileforuntitled20.fig%UNTITLED20,byitself,createsanewUNTITLED20orraisestheexisting%singleton*.%H=UNTITLED20returnsthehandletoanewUNTITLE

20、D20orthehandleto%theexistingsingleton*.%UNTITLED20('CALLBACK',hObject,eventData,handles,.)callsthelocal%functionnamedCALLBACKinUNTITLED20.Mwiththegiveninputarguments.%UNTITLED20('Property','Value',.)createsanewUNTITLED20orraisesthe%existingsingleton*.Startingfromtheleft,prope

21、rtyvaluepairsare%appliedtotheGUIbeforeuntitled20_OpeningFunctiongetscalled.An%unrecognizedpropertynameorinvalidvaluemakespropertyapplication%stop.Allinputsarepassedtountitled20_OpeningFcnviavarargin.%*SeeGUIOptionsonGUIDE'sToolsmenu.Choose"GUIallowsonlyone%instancetorun(singleton)".%Se

22、ealso:GUIDE,GUIDATA,GUIHANDLES%Edittheabovetexttomodifytheresponsetohelpuntitled20%LastModifiedbyGUIDEv2.501-Jul-201316:55:28%Begininitializationcode-DONOTEDIT12gui_Singleton=1;gui_State=struct('gui_Name''gui_Singleton''gui_OpeningFcn''gui_OutputFcn''gui_LayoutFcn

23、''gui_Callback'燕山大学数字信号处理课程设计mfilename,.,gui_Singleton,.,untitled20_OpeningFcn,untitled20_OutputFcn,，口,.,);ifnargin&&ischar(varargin1)gui_State.gui_Callback=str2func(varargin1);endvarargin:);ifnargoutvarargout1:nargout=gui_mainfcn(gui_State,elsegui_mainfcn(gui_State,varargin:);en

24、d%Endinitializationcode-DONOTEDIT%-Executesjustbeforeuntitled20ismadevisible.functionuntitled20_OpeningFcn(hObject,eventdata,handles,varargin)%Thisfunctionhasnooutputargs,seeOutputFcn.%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(

25、seeGUIDATA)%varargincommandlineargumentstountitled20(seeVARARGIN)%Choosedefaultcommandlineoutputforuntitled20handles.output=hObject;%Updatehandlesstructureguidata(hObject,handles);%UIWAITmakesuntitled20waitforuserresponse(seeUIRESUME)%uiwait(handles.figure1);%-Outputsfromthisfunctionarereturnedtothe

26、commandline.13燕山大学数字信号处理课程设计functionvarargout=untitled20_OutputFcn(hObject,handles)%varargoutcellarrayforreturningoutputargs(see%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversion%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Getdefaultcommandlineoutputfromhandlesstructurevarar

27、gout1=handles.output;%-Executesonbuttonpressinqu.functionqu_Callback(hObject,eventdata,handles)%hObjecthandletoqu(seeGCBO)%eventdatareserved-tobedefinedinafutureversion%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('qu.txt');axes(handles.axes1);plot(a);xlabel('时间(s)');ylab

28、el('被测变量y');title('(屈)原始信号');gridon;A=mean(a);%平均值B=std(a);%标准差C=Var(a);%方差D=mean(abs(a);%积分机电值E=sqrt(mean(a.A2);%匀方根有效值RMSy=fft(a,10000);%10000点傅里叶变换fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/10000;%率谱ss=sum(power);M2=0.5*ss;df=fs/N;eventdata,VARARGOUT);ofMATLABofM

29、ATLABM1=0.5*df*(sum(power(1:N-1)+sum(power(2:N);14燕山大学数字信号处理课程设计F=M1/M2;G=M2/2;%A,B,C,D,E,F=index_calculation(a);%设置显示A1=A*100000;A2=round(A1);%-IfEnable%-Otherwise,='on',executesonmousepressin5pixelexecutesonmousepressin5pixelborderborder.oroverA3=A2/100000;%round函数是对数据取整，以此限制数据的有效位set(hand

30、les.pingjunzhi,'string',num2str(A3);B1=B*100000;B2=round(B1);B3=B2/100000;set(handles.biaozhuncha,'string',num2str(B3);C1=C*100000;C2=round(C1);C3=C2/100000;set(handles.fangcha,'string',num2str(C3);D1=D*100000;D2=round(D1);D3=D2/100000;set(handles.jifenjidianzhi,'string&#

31、39;,num2str(D3);E1=E*100000;E2=round(E1);E3=E2/100000;set(handles.junfanggen,'string',num2str(E3);F1=F*100000;F2=round(F1);F3=F2/100000;set(handles.pingjungonglv,'string',num2str(F3);G1=G*100000;G2=round(G1);G3=G2/100000;set(handles.zhongweipinglv,'string',num2str(G3);pingjun

32、zhi.15燕山大学数字信号处理课程设计functionpingjunzhi_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletopingjunzhi(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',executesonmousepressin5pixelborder.%-Otherwise,executeson

33、mousepressin5pixelborderoroverbiaozhuncha.functionbiaozhuncha_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletobiaozhuncha(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',executesonmousepressin5pixelborde

34、r.%-Otherwise,executesonmousepressin5pixelborderoroverjifenjidianzhi.functionjifenjidianzhi_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletojifenjidianzhi(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',

35、executesonmousepressin5pixelborder.%-Otherwise,executesonmousepressin5pixelborderoroverpingjungonglv.functionpingjungonglv_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletopingjungonglv(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%-IfEnable%-Otherwise,fangcha.%handlesstructu

36、rewithhandlesanduserdata(seeGUIDATA)='on',executesonmousepressin5pixelborder.executesonmousepressin5pixelborderorover16燕山大学数字信号处理课程设计functionfangcha_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletofangcha(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructure

37、withhandlesanduserdata(seeGUIDATA)%-Executesonbuttonpressinshen.functionshen_Callback(hObject,eventdata,handles)%hObjecthandletoshen(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('shen.txt');axes(handles.axes2);plot(a)

38、;xlabel('时间(s)');ylabel('被测变量y');title('(伸)原始信号);gridon;A=mean(a);%平均值B=std(a);%标准差C=Var(a);%方差D=mean(abs(a);%积分机电值E=sqrt(mean(a.A2);%匀方根有效值RMSy=fft(a,10000);%10000点傅里叶变换fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/10000;%周期图法功率谱ss=sum(power);M2=0.5*ss;df=fs/N;M1=

39、0.5*df*(sum(power(1:N-1)+sum(power(2:N);F=M1/M2;G=M2/2;%A,B,C,D,E,F=index_calculation(a);%设置显示A1=A*100000;17燕山大学数字信号处理课程设计A2=round(A1);A3=A2/100000;set(handles.pingjunzhi,'string',num2str(A3);B1=B*100000;B2=round(B1);B3=B2/100000;set(handles.biaozhuncha,'string',num2str(B3);C1=C*1000

40、00;C2=round(C1);C3=C2/100000;set(handles.fangcha,'string',num2str(C3);D1=D*100000;D2=round(D1);D3=D2/100000;set(handles.jifenjidianzhi,'string',num2str(D3);E1=E*100000;E2=round(E1);E3=E2/100000;set(handles.junfanggen,'string',num2str(E3);F1=F*100000;F2=round(F1);F3=F2/100000;

41、set(handles.pingjungonglv,'string',num2str(F3);G1=G*100000;G2=round(G1);G3=G2/100000;set(handles.zhongweipinglv,'string',num2str(G3);%-Executesonbuttonpressinqu1.functionqu1_Callback(hObject,eventdata,handles)MATLAB%hObjecthandletoqu1(seeGCBO)%eventdatareserved-tobedefinedinafutureve

42、rsionof%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('qu.txt');y=fft(a,1024);%故1024点傅里叶变换fs=1000;N=length(y);mag=abs(y);18燕山大学数字信号处理课程设计f=(0:N-1)/N*fs;axes(handles.axes3);plot(f(1:512),mag(1:512);%故幅频谱xlabel('频率');ylabel('幅值');title('(屈)肌电幅频谱)；axes(handles.axe

43、s4);plot(f,angle(y);%做相频谱xlabel('频率');ylabel('幅值');title('(屈)肌电相频谱)；gridon;w=rceps(a);%t倒谱n=1:100;axes(handles.axes5);plot(n,w(1:100);xlabel('时间);ylabel('倒谱');title('(屈)肌电信号倒谱)；gridon;power1=(mag.A2)/2;%期图法求功率谱axes(handles.axes6);plot(f,power1);xlabel('频谱'

44、);ylabel('功率谱');title('(屈)肌电信号功率谱)；gridon;%-Executesonbuttonpressinshen1.functionshen1_Callback(hObject,eventdata,handles)%hObjecthandletoshen1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('shen.txt');y=ff

45、t(a,1024);%故1024点傅里叶变换fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;19燕山大学数字信号处理课程设计axes(handles.axes3);plot(f(1:512),mag(1:512);%故幅频谱xlabel('频率');ylabel('幅值');title('(伸)肌电幅频谱)；axes(handles.axes4);plot(f,angle(y);%做相频谱xlabel('频率');ylabel('幅值');title('(伸)肌电相频谱)

46、；gridonw=rceps(a);%t倒谱n=1:100;axes(handles.axes5);plot(n,w(1:100);xlabel('时间);ylabel('倒谱');title('(伸)肌电信号倒谱)；gridon;power1=(mag.A2)/10000;%周期图法求功率谱axes(handles.axes6);plot(f,power1);xlabel('频谱');ylabel('功率谱');title('(伸)肌电信号功率谱)；gridon;%-Executesonmousepressoveraxe

47、sbackground.functionaxes5_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletoaxes5(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)20燕山大学数字信号处理课程设计%-Executesonbuttonpressinqufa.functionqufa_Callback(hObject,eventdata,handles)%hObj

48、ecthandletoqufa(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalss%a=load('qu.txt');%D=mean(abs(a);%积分机电值%D1=D*100000;%D2=round(D1);%D3=dec2hex(D2);dtmf_series(01);dtmf_series(01);%-Executesonbuttonpressinshenfa.functions

49、henfa_Callback(hObject,eventdata,handles)%hObjecthandletoshenfa(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalss%a=load('shen.txt');%D=mean(abs(a);%积分机电值%D1=D*100000;%D2=round(D1);%D3=dec2hex(D2);dtmf_series(00);dtmf_series(00);functionUntitled_2_Callback(hObject,eventdata,handles)ofMATLAB%hObjecthandletoqu1(seeGCBO)%eventdatareserved-tobedefinedinafutureversion21燕山大学数字信号处理课程设计%handlesstructurewith