a1=load('EEG01.txt');
[c,r]=size(a1);
z=10;%等於幾，繪圖起點從幾開始
s=256*z;%繪圖起點;還有，這里的256是采樣率
d=floor(c/256);
cn=d*256;%d等於多少就是到多少秒
cal=2; %通道數
a=a1(1:cn,1); %取整秒數的時間段
t=(1:c)/256;
noise = 5*sin(2*pi*50*t);
y= conv(noise,a);
%e14保存theta波
e14=[];
n=0;%保存結果矩陣行數
m=0;
for j=1:256:cn
aa=a(j:j+255,:);%第j/256秒的數據塊
N=256;
b=zeros(4,cal);
%去除幅值絕對值超過100uV的數據段
[rr,rc]=find(-100>aa|aa>100);
nc=~isempty(rc);nr=~isempty(rr);
if ~(nc&&nr) %該數據塊幅值絕對值無超過100uV
m=m+1;
[delta,theta,alpha,beta]=Fftfil(aa);% 提取特征波段，這里需要調用另外的FFT函數
%此處不貼了：可以發eegdiy@aliyun.com交流
%計算第1通道的theta波能量
for i=1:N
b(2,1)=b(2,1)+theta(i,1)*theta(i,1);
end
n=n+1;
e14(n,1)=b(2,1)/N;
else
continue
end
b=zeros(4,1);
end
%最后畫出需要的值即可