脑电beta波的权重计算MATLAB,matlab计算相对功率(示例代码)

1、对脑电数据进行db4四层分解，因为脑电频率是在0-64HZ，分层后如图所示，

细节分量[d1 d2 d3 d4]

近似分量[a4]

重建细节分量和近似分量，然后计算对应频段得相对功率谱，重建出来得四个频段(αβθδ)都有14个通道，所以要计算4频段14通道共56个相对功率

2、代码

function wavelet(signal)

A4Array = zeros(14,5000);

D4Array = zeros(14,5000);

D3Array = zeros(14,5000);

D2Array = zeros(14,5000);

for i=1:14

[C,L] = wavedec(signal(i,1:5000),4,‘db4‘);%函数返回 3 层分解的各组分系数C(连接在一个向量里) ，向量 L 里返回的是各组分的长度。

% [cD1,cD2,cD3,cD4] = detcoef(C,L,[1,2,3,4]);%抽取1234层细节系数

% cA4 = appcoef(C,L,‘d4‘,4);%抽取近似系数

A4 = wrcoef(‘a‘,C,L,‘db4‘,4);%重建4层近似，deta波

A4Array(i,:) = A4;

D4 = wrcoef(‘d‘,C,L,‘db4‘,4);%重建4层细节，si