gusucode.com > 《MATLAB神经网络超级学习手册》随书光盘源码程序 > code/15/N15_2.m

    clear all
clc
x=0:0.02:5;
T=sin(x*3*pi);
Q=length(T);
P=zeros(4,Q);       %P中存储信号T的前4次值，作为网络输入。
P(1,2:Q)=T(1,1:(Q-1));
P(2,3:Q)=T(1,1:(Q-2));
P(3,4:Q)=T(1,1:(Q-3));
P(4,5:Q)=T(1,1:(Q-4));
plot(x,T)           %绘制信号T曲线
xlabel('时间');
ylabel('目标信号');
title('待预测信号');
net=newlind(P,T);   %根据输入和期望输出直接生成线性网络
a=sim(net,P);       %网络测试
figure(2)
plot(x,a,x,T,'+')
xlabel('时间');
ylabel('仿真和目标');
title('仿真信号和目标信号');
e=T-a
figure(3)
plot(x,e)
hold on
plot([min(x) max(x)],[0 0],'r:')
hold off
xlabel('时间');
ylabel('误差');
title('误差信号');