apha_waveleve1.m
上传用户:xhrjvip
上传日期:2020-11-07
资源大小:2k
文件大小:1k
源码类别:

生物技术

开发平台:

Matlab

  1. clc
  2. clear
  3.    filenm='02.eeg';
  4.    numdl=33;
  5.    smpf=210; 
  6.    %i表示第i通道
  7.    i=1;
  8.    %输入时间(分),输出信号段
  9.    %刺激前
  10.    t1=2.20;t2=3.20;
  11.    A_02_bef_sti=Extract_data(filenm,t1,t2,i);
  12.    %刺激前和刺激1的条形图数据分析
  13.    for i=1:512
  14.        x(i)=A_02_bef_sti(1,i);
  15.    end
  16.    for i=1:512
  17.        x1(i)=A_02_sti_sti1(1,i);
  18.    end
  19. [x_delta,x_theta,x_apha,x_beta]=wavelete_energy_decompose(x);
  20. [x1_delta,x1_theta,x1_apha,x1_beta]=wavelete_energy_decompose(x1);
  21. E=[x_delta,x_theta,x_apha,x_beta];
  22. E1=[x1_delta,x1_theta,x1_apha,x1_beta];
  23. data=[E;E1];
  24. b = bar(data);
  25. %ch = get(b,'children');
  26. %set(ch{1},'FaceVertexCData',[1;1;1;1;2;2;2;2;3;3;3;3;4;4;4;4])
  27. %set(ch{2},'FaceVertexCData',[1;1;1;1;2;2;2;2;3;3;3;3;4;4;4;4])
  28. %set(ch{3},'FaceVertexCData',[1;1;1;1;2;2;2;2;3;3;3;3;4;4;4;4])
  29. %set(ch{4},'FaceVertexCData',[1;1;1;1;2;2;2;2;3;3;3;3;4;4;4;4])
  30. %text(1,1,'原始图象 概率分布图象 能量数值图像');
  31. legend('δ','θ','α','β');