fData=textread('MDFWORK.txt');
x=fData(:,2);
y=fData(:,12);
lev=4;
wpt=wpdec(y,lev,'db4','shannon');
for i=1:2^lev
s(i,:)=wprcoef(wpt,[4,i-1]);
end
E1(1:2^lev)=0;
for k=1:2^lev
for j=1:1:length(s(i,:))
E1(k)=E1(k)+s(k,j)*s(k,j)*j;
end
end
S=0;
for k=1:2^lev
S=S+E1(k);
end
for k=1:2^lev
E(k)=E1(k)/S;
end
系统提示第7行有错,不解,求答案呀,谢谢啦