flch5eg1.m

采用CMAC网络对对频率为0.1572HZ的谐波进行辨识

%CMAC group program,include four sub-CMAC 
%CMAC  program 1
clear
w1=0.1;w2=0.0;w3=-0.1; w4=0.0
N=144;  N1=20; E=0.0;
b=0.062; k=0; i=0;
for k=1:N   
 
yp=2.5*cos(0.99*k);
for j=1:30
 wa=w1+w2+w3+w4;
 a=b*(1-k/N+1);
 w5=w4+a*[yp-wa];
w1=w2;
 w2=w3;
 w3=w4;
 w4=w5;
end
e=[(yp-wa).^2]/2;
sub1e(k)=e;
 yp1(k)=yp;
 sub1w(k)=wa;
 %if e<=E
 %    break
 % else
 % end
end
sub1wz=sub1w;

 
 %CMAC  program 2
%clear
w1=0.1;w2=0.0;w3=-0.1; w4=0.0;
 N=144; N1=20; E=0.0;
b=0.0612;%b=0.12;
k=1; 
for k=1:N  
   
ypp=10/(pi*3)*cos(1.97*k);
for j=1:30
 wa=w1+w2+w3+w4;
 a=b*(1-k/N+1);
 w5=w4+a*[ypp-wa];
w1=w2;
 w2=w3;
 w3=w4;
 w4=w5;
end
e=[(ypp-wa).^2]/2;
sub2e(k)=e;
 yp2(k)=ypp;
 sub2w(k)=wa;
 %if e<=E
 %    break
 % else
 % end
end
sub2wz=sub2w;



%CMAC  program 3
w1=0.0;w2=0;w3=-0.1; w4=0.1;
 N=144; N1=20; E=0.0; wa0=0;wa1=0;
b=0.062; k=0; i=0;
for k=1:N  
   
yppp=-10/(pi*15)*cos(3.95*k);
for j=1:30
 wa=w1+w2+w3+w4;
 a=b*(1-k/N+1);
 w5=w4+a*[yppp-wa];
w1=w2;
 w2=w3;
 w3=w4;
 w4=w5;
end
%e=yppp-wa;
e=[(yppp-wa).^2]/2;
sub3w(k)=wa;
 yp3(k)=yppp;
 %if e<=E
 %    break
 % else
 % end
sub3e(k)=e;
end
sub3wz=sub3w;


%CMAC  program 4
%clear
w1=0.1;w2=0.0;w3=-0.1; w4=0.0;
 N=144; N1=20; E=0.0;
b=0.0612;%b=0.12;
k=1; 
for k=1:N  
   
ypp=10/(pi*35)*cos(5.93*k);
for j=1:30
 wa=w1+w2+w3+w4;
 a=b*(1-k/N+1);
 w5=w4+a*[ypp-wa];
w1=w2;
 w2=w3;
 w3=w4;
 w4=w5;
end
e=[(ypp-wa).^2]/2;
sub4e(k)=e;
 yp4(k)=ypp;
 sub4w(k)=wa;
 %if e<=E
 %    break
 % else
 % end
end
sub4wz=sub4w;

%CMAC  program zong
%clear
w1=0.1;w2=0.0;w3=-0.1; w4=0.0;
 N=144; N1=20; E=0.0;
b=0.0612;%b=0.12;
k=1; 
for k=1:N  
   
ypp=2.5*cos(0.99*k)+10/(pi*3)*cos(1.97*k)-10/(pi*15)*cos(3.95*k)+10/(pi*35)*cos(5.93*k);
for j=1:30
 wa=w1+w2+w3+w4;
 a=b*(1-k/N+1);
 w5=w4+a*[ypp-wa];
w1=w2;
 w2=w3;
 w3=w4;
 w4=w5;
end
e=[(ypp-wa).^2]/2;
subze(k)=e;
 ypz(k)=ypp;
 subzw(k)=wa;
 %if e<=E
 %    break
 % else
 % end
end
subzwz=subzw;

 %grapher
 i=1:N
 subplot(6,1,1)
 plot(i,yp1,i,sub1w,'rx')
 ylabel('w1')
 legend('y1 ','w1')
 subplot(6,1,2)
 plot(i,yp2,i,sub2w,'rx')
 ylabel('w2')
 legend('y2 ','w2')
 subplot(6,1,3)
 plot(i,yp3,i,sub3w,'rx') 
 ylabel('w3')
legend('y3 ','w3')
 subplot(6,1,4)
 plot(i,yp4,i,sub4w,'rx') 
 ylabel('w4')
legend('y4  ','w4 ')
subplot(6,1,5)
 plot(i,ypz,i,subzw,'rx') 
 ylabel('wz')
legend('yz ','wz')
subplot(6,1,6)
 plot(i,sub1e,'o',i,sub2e,'x',i,sub3e,':',i,sub4e,'-',i,subze,'*') 
 %plot(i,sub1e,'x',i,sub2e,'o')
 xlabel('k')
 ylabel('e')
 legend(' e1 ',' e2',' e3',' e4',' ez')