|
a |
|
b/Feature Extraction/extractAR.m |
|
|
1 |
function X=extractAR(x_train,AROrder,startS,endS,wStep,wRange) |
|
|
2 |
|
|
|
3 |
FS=128; |
|
|
4 |
|
|
|
5 |
N=size(x_train,3); |
|
|
6 |
sz=floor((endS-(startS+wRange))/wStep)+1; |
|
|
7 |
X=zeros(sz*FS,2); |
|
|
8 |
cn=0; |
|
|
9 |
for i=1:N |
|
|
10 |
|
|
|
11 |
for sig=startS:wStep:endS-wRange |
|
|
12 |
|
|
|
13 |
sW=sig*FS+1; |
|
|
14 |
eW=(sig+wRange)*FS; |
|
|
15 |
|
|
|
16 |
C3Sig=x_train(sW:eW,1,i); |
|
|
17 |
C4Sig=x_train(sW:eW,3,i); |
|
|
18 |
|
|
|
19 |
c3= arburg(C3Sig, AROrder); |
|
|
20 |
c4= arburg(C4Sig, AROrder); |
|
|
21 |
cn=cn+1; |
|
|
22 |
X(cn,1)=sum(c3.^2)/numel(c3); |
|
|
23 |
X(cn,2)=sum(c4.^2)/numel(c4); |
|
|
24 |
|
|
|
25 |
end |
|
|
26 |
end |
|
|
27 |
|
|
|
28 |
end |
|
|
29 |
|