Delta波：0.5-4Hz;

Theta波：4-8Hz;

Alpha波：8-13Hz;

Beta波：13-32Hz;

Gamma波：>32Hz；

## 代码案例

### Delta波

[header,data] = edfread('Affaf Ikram 20121020 1839.L1.edf');    fp_data=data(1,:);    fp_data=resample(fp_data,100,500); %降采样率为100hz    fc = 600;          %采样频率600Hz    N = 30000;           %采样点数    n = 0:N-1;              f = n*fc/N;         %频率序列    Wn = [0.5*2 4*2]/fc;%设置通带为0.5-4Hz Delta波    [k,l] = butter(2,Wn);%4阶IIR滤波器    result = filtfilt(k,l,fp_data);    tmp = fft(result);    fp1_Delta= ifft(tmp);    figure,    subplot(211),plot(f(1:N/2),abs(tmp(1:N/2)*2/N));    title('滤波后频谱结果，频率:0.5-4Hz');    axis([0,100,-inf,inf])    subplot(212),plot(fp1_Delta);    title('脑波滤波后Delta波结果图像');

### Theta波

[header,data] = edfread('Affaf Ikram 20121020 1839.L1.edf');    fp_data=data(1,:);    fp_data=resample(fp_data,100,500); %降采样率为100hz    fc = 600;          %采样频率600Hz    N = 30000;           %采样点数    n = 0:N-1;              f = n*fc/N;         %频率序列    Wn = [4*2 8*2]/fc;%设置通带为4-8Hz Theta波    [k,l] = butter(2,Wn);%4阶IIR滤波器    result = filtfilt(k,l,fp_data);    tmp = fft(result);    fp1_Delta= ifft(tmp);    figure,    subplot(211),plot(f(1:N/2),abs(tmp(1:N/2)*2/N));    title('滤波后频谱结果，频率:4-8Hz');    axis([0,100,-inf,inf])    subplot(212),plot(fp1_Delta);    title('脑波滤波后Theta波结果图像');

### Alpha波

[header,data] = edfread('Affaf Ikram 20121020 1839.L1.edf');    fp_data=data(1,:);    fp_data=resample(fp_data,100,500); %降采样率为100hz    fc = 600;          %采样频率600Hz    N = 30000;           %采样点数    n = 0:N-1;              f = n*fc/N;         %频率序列    Wn = [8*2 13*2]/fc;%设置通带为8-13Hz Alpha波    [k,l] = butter(2,Wn);%4阶IIR滤波器    result = filtfilt(k,l,fp_data);    tmp = fft(result);    fp1_Delta= ifft(tmp);    figure,    subplot(211),plot(f(1:N/2),abs(tmp(1:N/2)*2/N));    title('滤波后频谱结果，频率:8-13Hz');    axis([0,100,-inf,inf])    subplot(212),plot(fp1_Delta);    title('脑波滤波后Alpha波结果图像');

### Beta波

[header,data] = edfread('Affaf Ikram 20121020 1839.L1.edf');    fp_data=data(1,:);    fp_data=resample(fp_data,100,500); %降采样率为100hz    fc = 600;          %采样频率600Hz    N = 30000;           %采样点数    n = 0:N-1;              f = n*fc/N;         %频率序列    Wn = [13*2 32*2]/fc;%设置通带为13-32Hz Beta波    [k,l] = butter(2,Wn);%4阶IIR滤波器    result = filtfilt(k,l,fp_data);    tmp = fft(result);    fp1_Delta= ifft(tmp);    figure,    subplot(211),plot(f(1:N/2),abs(tmp(1:N/2)*2/N));    title('滤波后频谱结果，频率:13-32Hz');    axis([0,100,-inf,inf])    subplot(212),plot(fp1_Delta);    title('脑波滤波后 Beta波结果图像');

### Gamma波

[header,data] = edfread('Affaf Ikram 20121020 1839.L1.edf');    fp_data=data(1,:);    fp_data=resample(fp_data,100,500); %降采样率为100hz    fc = 600;          %采样频率600Hz    N = 30000;           %采样点数    n = 0:N-1;              f = n*fc/N;         %频率序列    Wn = [32*2 40*2]/fc;%设置通带为32-50Hz Gamma波    [k,l] = butter(2,Wn);%4阶IIR滤波器    result = filtfilt(k,l,fp_data);    tmp = fft(result);    fp1_Delta= ifft(tmp);    figure,    subplot(211),plot(f(1:N/2),abs(tmp(1:N/2)*2/N));    title('滤波后频谱结果，频率:32-50Hz');    axis([0,100,-inf,inf])    subplot(212),plot(fp1_Delta);    title('脑波滤波后Gamma波结果图像');