bijin.m

一个关于BP网络的设计

%函数逼近例子
%bijin.m
N=1;%频率参数
p=[-1:0.05:1];
t=sin(N*pi*p);
plot(p,t,'r*')
title('要逼近的非线性函数');
xlabel('时间');
%建立网络
net=newff(minmax(p),[8,1],{'tansig','purelin'},'trainlm');
y1=sim(net,p);
figure;
plot(p,t,'r*',p,y1,'-')
title('未训练网络输出结果');
xlabel('时间');
ylabel('仿真输出- 原函数输出*');
legend('要逼近的非线性函数','未训练网络输出结果')
%网络训练
net.trainparam.show=10;
net.trainparam.lr=0.05;
net.trainparam.epochs=100;
net.trainparam.goal=0.0001;
net1=train(net,p,t)
%网络测试
y2=sim(net1,p);
figure;
plot(p,t,'r*',p,y1,'b-',p,y2,'ko')
title('训练后网络输出结果');
xlabel('时间');
ylabel('仿真输出');
legend('要逼近的非线性函数','未训练网络输出结果','训练后网络输出结果')