matlab小波网络,小波神经网络MATLAB程序（真正共享）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
load data
M1=20;
epo=15;
A=4;
B=18;
B2=B/2+1
N=500;
M=(A+1)*(B+1);
for a0=1:A+1;
for b0=1:B+1;
i=(B+1)*(a0-1)+b0;
b_init(i)=((b0-B2)/10)/(2^(-A)); a_init(i)=1/(2^(-A));
c_init(i)=(20-A)/2;
end
end
w0=ones(1,M);
for i=1:N
for j=1:M
t=x(i);
t= a_init(j)*t-b_init(j);
%P0(i,j)= (cos(1.75*t)*exp(-t*t/2))/2^c_init(j);
P0(i,j)= ((1-t*t)*exp(-t*t/2))/2^c_init(j);
end
end
%calculation of output of network
for i=1:N
u=0;
for j=1:M
u=u+w0(j)*P0(i,j);%w0?aè¨?μ
end
y0(i)= u;% y(p)= u=??W(j)*phi(p,j)= ??W(j)* |μj(t)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for k=1:M
W(:,k)=P0(:,k);
end
for k=2:M
u=0;
for i=1:k-1
aik(i)=(P0(:,k)'*W(:,i))/(W(:,i)'*W(:,i));
u=u+aik(i) *W(:,i);
end
W(:,k)=P0(:,k)-u;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:M
g(i)= (W(:,i)'*d')/( W(:,i)'* W(:,i));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
u=0;
for i=1:M
u=u+g(i)*(W(:,i)'*W(:,i));
end
DD=u;
for i=1:M
Erro(i)=(g(i)^2)*(W(:,i)'*W(:,i))/DD;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k=1;
while(k<=M1)
u=1;
for i=2:M
if abs(Erro(u))
u=i
else u=u
end
end
I(k)=u;
Erro(u)=0
k=k+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for k=1:M1
u=I(k);
a(k)=a_init(u);
b(k)=b_init(u);
c(k)=c_init(u);
w1(k)=w0(u);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
epoch=1;
error=0.1;
err=0.0001;
lin=0.5;
while (error>=err & epoch<=epo)
for i=1:N
for j=1:M1
t=x(i);
t= a(j)*t-b(j);
%P1(i,j)= (cos(1.75*t)*exp(-t*t/2))/2^c(j);
P1(i,j)= ((1-t*t)*exp(-t*t/2))/2^c(j);
end
end
%calculation of output of network
for i=1:N
u=0;
for j=1:M1
u=u+w1(j)*P1(i,j);
end
y1(i)= u;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
u=0;
for i=1:N
u=u+(d(i)-y1(i))^2;
end
u=u/2;%u=1/2??(d-p)^2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if u>error
lin=lin*0.8;
else
lin=lin*1.2;
end
error=u; %error=u=1/2??(d-p)^2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j=1:M1
u=0;
for i=1:N
u=u+(d(i)-y1(i))*P1(i,j);
end
EW(j)=-u;
end
if epoch==1
SW=-EW;
w1_=w1;
else
SW=-EW+((EW*EW')*SW_)/(EW_*EW_');
end
EW_=EW;
SW_=SW;
w1=w1_+SW*lin;
w1_=w1;
%number of epoch increase by 1
epoch=epoch+1;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(2,1,1);plot(x,d);
subplot(2,1,2);plot(x,y1);