神经网络辨识系统(二)

clc;close all;clear;

y(1)= 0;
k=1;
u(k)= 2*rand-1;
f(k)= 0.6*sin(pi*u(k))+0.3*sin(3*pi*u(k))+0.1*sin(5*pi*u(k));
y(k+1)= 0.3*y(k)+f(k);

for k=1:100
    u(k)= 2*rand-1;
    f(k)= 0.6*sin(pi*u(k))+0.3*sin(3*pi*u(k))+0.1*sin(5*pi*u(k));
    y(k+1)= 0.3*y(k)+0.6*(k-1)+f(k); 
end

uP= u(1:100);
uT= 0.6*sin(pi*uP)+0.3*sin(3*pi*uP)+0.1*sin(5*pi*uP);
netg= newff(minmax(uP),[20,20,1],{'tansig','tansig','purelin'},'trainlm','learngdm','msereg');
[netg,tr] = train(netg,uP,uT);

uQ= -1:0.1:1;
simu= sim(netg,uQ);
figure;
plot(uQ,0.6*sin(pi*uQ)+0.3*sin(3*pi*uQ)+0.1*sin(5*pi*uQ),'b',uQ,simu,'r');
legend('系统输出','系统辨识');

k=1;
u(k)= 2*rand-1;
f(k)= 0.6*sin(pi*u(k))+0.3*sin(3*pi*u(k))+0.1*sin(5*pi*u(k));
y(k+1)= 0.3*y(k)+f(k);

for k=2:100