肌电信号的时域和频域分析

1、 肌电信号的时域和频域分析摘要：肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。 其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应用前景。 本次设计基于matlab用小波变换对肌电信号进行消噪处理，分别选用20N的肌电信号数据和50N的肌电数据进

2、行对比，最后在GUI界面上完成相应的功能处理。 关键字：肌电信号 Matlab 小波去噪 GUI第1章 绪论 肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位 在 时 间 和 空 间 上 的 叠 加，反 映 了 神 经，肌 肉 的 功 能 状 态，在基础医学研究、临床诊断和康复工程中有广泛的应用。其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应

3、用前景。肌电信号本身是一种较微弱的电信号。检测和记录表面肌电信号，需要考虑的主要问题是尽量消除噪声和干扰的影响, 提高信号的保真度1。第2章 肌电信号的时域分析2.1 肌电信号时域图的显示及比较 肌电信号采用两个不同的数据进行比较，通过比较时域图及其特性来进行分析2。其图像如下所示：如上图所示：肌电数据分别是同一个体在20N的力和50N的力所反映的图像。可以看出在不同作用力时，其图像的差别很大。2.2 时域参数 2.2.1 均值对于一个随机变量来说，均值是一个很重要的数值特征。粗略的说，就是来描述一个群体的平均水平。其严格的数学定义非常的简单，就是一个随机变量关于概率测度的积分。这样的积分在测

4、度轮或者实分析里是没有什么直观的解释的。而在概率论里却成为了一个群体的主要指标。在此处，均值表示肌电信号的平均水平。2.2.2 标准差标准差（Standard Deviation），也称均方差（mean square error），是各数据偏离平均数的距离的平均数，它是离均差平方和平均后的方根，用表示。标准差是方差的算术平方根。标准差能反映一个数据集的离散程度。平均数相同的，标准差未必相同。2.2.3 方差 方差是各个数据与平均数之差的平方的平均数。在概率论和数理统计中，方（英文Variance）用来度量随机变量和其数学期望（即均值）之间的偏离程度。在许多实际问题中，研究随机变量和均值之间的偏

5、离程度有着很重要的意义。 选取一个信号，其执行结果如下所示： 其部分程序代码如下所示：clear;clc;s=load(E:肌电信号数据EMGEMG3.txt);%s=load(E:肌电信号数据EMGbs1.txt);a=s(:,7);t=s(:,1);fprintf(n数据基本信息：n);fprintf( 均值=%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(

6、mean(a.2);第3章 肌电信号的时域分析3.1 傅里叶变换傅里叶是离散傅立叶变换的快速算法，可以将一个信号变换到频域。有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后就很容易看出特征了。另外，FFT可以将一个信号的频谱提取出来，这在频谱分析方面也是经常用的。一个模拟信号，经过ADC采样之后，就变成了数字信号。采样得到的数字信号，就可以做FFT变换了。N个采样点，经过FFT之后，就可以得到N个点的FFT结果。本次设计选取N=20000来进行肌电信号的频域分析3。通过傅里叶变换后，两个不同信号的幅频图如下所示： 肌电相频图如下所示：3.2 功率谱分析 能量谱密度、功率谱密度函数表

7、示信号的能量、功率密度随频率变化的情况。通过研究功率谱密度，可以帮助了解信号的功率分布情况，确定信号的频带等。 功率密度谱虽然描述了随机信号的功率在各个不同频率上的分布，但因为它仅与幅度频谱有关，没有相位信息，所以从已知功率谱还难以完整地恢复原来的功率信号。 通过执行相应程序后，其功率谱的显示图如下所示：其部分程序代码如下所示：global 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(

8、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);mag=abs(s);f=(0:N-1)/N*fs;power1=(mag.2)/2000;plot(handles.axes1,power1)第4章 小波去噪分析小波分析方法是一种窗

9、口大小（即窗口面积）固定但其形状可改变，时间窗和频率都可改变的时频局域化分析方法，即在低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，所以被称为数学显微镜。正是这种特性，是小波变换具有对信号的自适应性4。4.1 常用的小波分析方法4.1.1 小波消噪一维的小波消噪过程主要分为以下三个步骤：一维信号的分解：获得各尺度上的细节分量和近似分量；细节分量的阈值化：对1到N尺度的上的每层细节分量选取合适的阈值，进行阈值化处理，保留有用信息，去除噪声成分；一维信号的重构：根据小波分解的第N层近似分量和经阈值化处理后的第1层到第N层细节分量，进行一维信号的小

10、波重构。本文还利用小波的多分辨率特性将心电信号分解为多层，然后进行有用信号的重构，将一些噪声信号去除。 4.1.2 小波函数的选取同一信号，选取不同的小波函数进行处理，将得到不同的效果，所以小波函数的选取显得尤为重要。对于心电信号滤波来说，选择支集长度较短的小波可提高处理的实时性；选取对称性的小波可满足相移为基本线性，使心电信号不失真；选取正则性的小波可使重构以后的信号比较平滑。经过对比降噪效果，在本文中用db5小波基进行消噪处理。4.1.2 阈值选取 阈值的确定是小波收缩消噪最关键的一步，阈值过小，则方差偏大，数据欠平滑；阈值过大，会使数据过平滑，信号的奇异性可能丧失。对小波系数进行阈值操作

11、过程中，有两种方式，其一对每一个小波系数进行阈值操作，其二是成块习俗进行阈值操作。由信号的奇异性理论，心电信号里的噪声具有负的奇异性，其幅度和稠密度随尺度的增大而减小，而信号则相反，因此阈值的选取不能单一。阙值函数体现了对小波系数的不同处理策略，主要分为软阈值函数和硬阈值函数它们的基本思想都是去除小幅值的系数，对大幅值的系数进行保留。它们的去噪效果各有特点(可以从后面的具体例子中看出)：硬阈值法可以很好的保留信号的边缘等局部特征，但去噪的结果具有较大的方差，会出现伪Gibbs现象等视觉失真；软闭值法去噪结果相对平滑，但有较大的偏差，可能会造成边缘模糊等失真现象。4.2 小波去噪的具体分析 为了

12、确定和比较本次设计中小波去噪的效果，先对小波进行一层分解，二层分解，三层分解，四层分解，让它们单独对同一信号进行消噪处理。其结果图如下所示：通过比较上图可以看出，第四幅图像的去噪效果比较明显，因次在本次设计中选择小波四层分解对信号进行去噪处理。分别对20N 的肌电信号和50N的肌电信号进行小波去噪，去噪前后的效果图如下所示：左边是20N小波去噪前后的图，右边是50N作用前后的效果图。通过前后图像的比较可以看出，当选用四层小波分解进行滤波时，已经能去掉大部分噪声的干扰，保存有用的信号。小波去噪部分程序代码如下:global a;global t;global s;global y;% globa

13、l IR;% global SORH;M =length (a) ;N=length(y);p=size(a);s=a(1:20000);C L=wavedec(a,4,db5);cA3=appcoef(C,L,db5,4);cD1=detcoef(C,L,1);cD2=detcoef(C,L,2);cD3=detcoef(C,L,3);cD4=detcoef(C,L,4);thr1=thselect(cD1,rigrsure);thr2=thselect(cD2,rigrsure);thr3=thselect(cD3,rigrsure);thr4=thselect(cD4,rigrsure)

14、;TR=thr1,thr2,thr3,thr4;SORH=s;XC,CXC,LXC,PERFO,PERF2=wdencmp(lvd,a,. db5,4,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;% K=length(d);% t1=(0:K-1)/2000;plot(handles.axes5,XC(1:20000);小波去

15、噪前后的幅频和相频图如下：表面肌电信号一般只有微伏级电压,信号中往往夹带着低频(接近直流)和高频的干扰信号,真正有用的肌电信号大致在10Hz500Hz之间。由上面的对比图可知，消除了噪声信号，有用信号得以保存。其部分程序代码如下所示：global a;global t;%global h;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

16、);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;plot(handles.axes5,f,mag1);去噪前后的功率谱图如下所示： 部分程序代码如下：global a;global t;global s;s=fft(a,2000);C L=wavedec(a,3,db5);cA3=appcoef(

17、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,2000);fs=2000;N=length(y1); mag=abs(y1);f=(0:N-1)/N*fs;powe

18、r2=(mag.2)/2000;plot(handles.axes5,power2);第5章 GUI界面的设计GUI是由窗口、光标、按键、菜单、文字说明等对象（Objects）构成的一个用户界面。用户通过一定的方法（如鼠标或键盘）选择、激活这些图形对象，使计算机产生某种动作或变化，比如实现计算、绘图等。一个好的GUI 能够使程序更加容易的使用。它提供用户一个常见的界面，还提供一些控件，例如，按钮，列表框，滑块，菜单等。用户图形界面应当是易理解且操作是可以预告的，所以当用户进行某一项操作，它知道如何去做。例如，当鼠标在一个按钮上发生了单击事件，用户图形界面初始化它的操作，并在按钮的标签上对这个操

19、作进行描述。本次设计的GUI界面整体图如下所示： 程序运行之后产生的界面图如下所示： 本次设计模块涵盖多个功能，都是通过按钮来进行相应的信号处理。只要M文件在当前目录或在Matlab搜索路径上，在Matlab命令窗口输入对应的M文件，就能打如上图所示的图形用户界面。在此界面上可以进行相关的操作。例如点击open按扭就会打开如下界面：Open按钮的回调函数如下填写：filename,filepath=uigetfile(*.txt,选择文件);%选择数据文件 str=filepath filename;其他的按钮都按照此方法去激活，以实现相应的功能。学习心得 通过这次设计，我对MATLAB有了更

20、深入的理解，学会了在GUI界面通过按钮来实现相应的功能，对以前未接触的小波变换也有了大概了解。在设计的过程中，我也认识到了自己所学知识的不足。这也让我再次认识到知识是无尽的，只有不断的充实自己、完善自己的知识理论体系，才能够更好的胜任自己以后的工作。设计过程中知识的不足也让我更加坚定了终身学习的决心。参考文献1.谢平、王娜、林洪斌等主编，信号处理原理及应用。北京：机械工业出版社， 2008.10 2. 聂祥飞、王海宝、谭泽富主编，Matlab程序设计及其在信号处理中的应用。 成都：西南交通大学出版社，2005 3.吴大正、高西全等主编，Matlab及在电子信息课程中的应用。北京：电子工业出版社

21、，2006.3 4.李培芳、孙晖、李江主编，信号与系统分析基础。北京：清华大学出版社，2006.12 程序附录function varargout = a1fig(varargin)gui_Singleton = 1;gui_State = struct(gui_Name, mfilename, . gui_Singleton, gui_Singleton, . gui_OpeningFcn, a1fig_OpeningFcn, . gui_OutputFcn, a1fig_OutputFcn, . gui_LayoutFcn, , . gui_Callback, );if nargin &

22、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 a1fig is made visible.function a1fig_OpeningFcn(hObject, eve

23、ntdata, 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 arguments to a1fig (see VARARGIN) % Choose defau

24、lt command line output for a1fighandles.output = hObject; % Update handles structureguidata(hObject, handles); % UIWAIT makes a1fig wait for user response (see UIRESUME)% uiwait(handles.figure1); % - Outputs from this function are returned to the command line.function varargout = a1fig_OutputFcn(hOb

25、ject, 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 GUIDATA) % Get default command line output from handles structure

26、varargout1 = 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% handles structure with handles and user data (see GUIDATAgl

27、obal s;global a;global t;filename,filepath=uigetfile(*.txt,选择文件);%选择数据文件 str=filepath filename; s=load(str); a=s(:,7); t=s(:,1); % - Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to b

28、e defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)try delete(allchild(handles.axes1);end % - Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserv

29、ed - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global a;global t;plot(handles.axes1,t,a);title(信号1); % - Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see

30、 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;%a=EMG3(1:points,2);y=fft(a,20000); %对信号进行快速Fourier变换fs=2000;N=length(y);mag=abs(y);%求得Fourier变换后的振幅f=(0:N-1)/N*fs; %频率序列plot(handles.axes1,f

31、,mag);%绘制图形 % - Executes on button press 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;

32、global y;fs=2000;N=length(y);f=(0:N-1)/N*fs; plot(handles.axes1,f,angle(y); % - 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 stru

33、cture with handles and user data (see GUIDATA)%小波去噪global a;global t;global s;global y;% global IR;% global SORH;M =length (a) ;N=length(y);p=size(a);s=a(1:20000);C L=wavedec(a,4,db5);cA3=appcoef(C,L,db5,4);cD1=detcoef(C,L,1);cD2=detcoef(C,L,2);cD3=detcoef(C,L,3);cD4=detcoef(C,L,4);thr1=thselect(cD1

34、,rigrsure);thr2=thselect(cD2,rigrsure);thr3=thselect(cD3,rigrsure);thr4=thselect(cD4,rigrsure);TR=thr1,thr2,thr3,thr4;SORH=s;XC,CXC,LXC,PERFO,PERF2=wdencmp(lvd,a,. db5,4,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*

35、log10(F);MSE=M/N;SM=SNR/MSE;% K=length(d);% t1=(0:K-1)/2000;plot(handles.axes5,XC(1:20000); % - 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 MATL

36、AB% handles structure with handles and user data (see GUIDATA)%去噪后幅频global a;global t;%global h; 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;

37、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;plot(handles.axes5,f,mag1); % - Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)% hObject handle to pushbutton8 (see GCBO)%

38、 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(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,rigr

39、sure);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;plot(handles.axes5,f,angle(y1); % - Executes on button press in pushbutton9.function pushbutton9_Callback(hObject, eventd

40、ata, 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)%功率谱global 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

41、(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); mag=abs(s);f=(0:N-1)/N*fs;power1=(mag.2)/2000;plot(handles.axes1,power1

42、); % - Executes on button press in pushbutton10.function pushbutton10_Callback(hObject, eventdata, handles)% hObject handle to pushbutton10 (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

43、 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);y1=fft(XC,2000);fs=

44、2000;N=length(y1); mag=abs(y1);f=(0:N-1)/N*fs;power2=(mag.2)/2000;plot(handles.axes5,power2); % - 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 future version of

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

46、s 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_CreateFcn(hObject, eventdata, handles)% hObject handle t

47、o edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows.% See ISPC and COMPUTER.if ispc & isequal(get(hObject,BackgroundColor), get(0,d

48、efaultUicontrolBackgroundColor) 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 struct

49、ure with handles and user data (see GUIDATA)global a;x=sqrt(var(a);set(handles.edit2,string,x); function edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (se

50、e GUIDATA) % Hints: get(hObject,String) returns contents of edit2 as text% str2double(get(hObject,String) returns contents of edit2 as a double % - Executes during object creation, after setting all properties.function edit2_CreateFcn(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)%

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

52、groundColor) set(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 GUIDATA) global a;x=(var(a);set(handles.edit3,string,x); function edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (