rfnni_2.m

利用遞迴式模糊類神經網路(recerrent neural networks system identification, RFNN)進行系統辨識, 自己寫的請多包含

%Adaptive PID control based on RBF Identification
clear all;
close all;
clc;
format('long')


xite=0.1;

x=[0; 0];

ci=[-1.5, -.8, 0, .8, 1.5;...
    -1.5, -.8, 0, .8, 1.5];
%ci=ones(2,5);
bi=[ones(1,5); ones(1,5)];
w=0.1*ones(25,1); 
wi=0.1*ones(2,5);

% initial feedback output
O12_1 = zeros(2,5)

error_1 = 0;
error_2 = 0;

tf=2;
ts=0.001;
n=floor(tf/ts);

for k=1:n
   time(k)=k*ts;
   rin(k)=1.0*sin(sin(2*pi*k*ts));
   
  % Input layer
        for i1 = 1:2
           for j1=1:5
              O12(i1,j1) = O12_1(i1,j1)*wi(i1,j1) + x(i1,end);
           end
        end
   % Fuzzification
  for i = 1:2 
   for j=1:5
      f2(i,j)=exp(-norm(O12(i,j)-ci(i,j))^2/(bi(i,j)^2))
   end
  end
   % Rules
 for j = 1:5
    m1(j)=f2(1,j);
    m2(j)=f2(2,j);
 end
  
   for i = 1:5
      for j = 1:5
         ff3(i,j) = m2(i)*m1(j); 
      end
   end
  
   f3 = [ff3(1,:), ff3(2,:), ff3(3,:), ff3(4,:), ff3(5,:)];
   % Defuzzification
   ymout(k)=(w'*f3');         
   
   wof3=w'*f3';
   e(k) = rin(k)-ymout(k);
% Learning 
   d_w=0*w;
   for j=1:25
      d_w(j)=xite*e(k)*f3(1,j);
   end
   w=w+d_w;
   
   d_bi=0*bi;
  for i = 1:2 
   for j=1:5
      d_bi(i,j)=xite*e(k)*wof3*(bi(i,j)^-3)*norm(O12(i,j)-ci(i,j))^2;
   end
  end
   bi=bi+ d_bi;
   %pause
   
   
     for i=1:2
         for j=1:5
      d_ci(i,j)=xite*e(k)*wof3*(O12(i,j)-ci(i,j))*(bi(i,j)^-2);
         end
     end

   ci=ci+d_ci;
 
   d_wi=0*wi;
  for i = 1:2 
   for j=1:5
      d_wi(j)=xite*e(k)*wof3*(O12(i,j)-ci(i,j))*(bi(i,j)^-2)*O12_1(i,j);
   end
  end
   wi=wi+d_wi;
   
   % Return params
   O12_1 = O12;
   x(1,k+1) = e(k);
   x(2,k+1) = e(k)-2*error_1-error_2;

error_2 = error_1;   
error_1 = e(k);
   k
  % pause
end
figure(1);
plot(time,rin,'r',time,ymout,'b');
xlabel('time(s)');ylabel('rin,ymout');
legend('reference', 'NNout');