p3.51.txt

来自「Matlab_Network_and_Application是国防工业出版社董长」· 文本 代码 · 共 37 行

%应用线形网络进行预测，假设信号T持续时间为5秒，采样频率为40次每秒
time=0:0.025:5
T=sin(time*4*pi);
Q=length(T);
%信号T的前5次值，作为网络的输入（构造全0初始矩阵）
P=zeros(5,Q);
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));
P(5,6:Q)=T(1,1:(Q-5));
%显示信号T的曲线
plot(time,T)
xlabel('时间');
ylabel('目标信号');
title('待预测信号');
%设计网络
net=newlind(P,T);
%验证
a=sim(net,P);
%显示信号a与T比较的曲线
plot(time,a,time,T,'+')
xlabel('时间');
ylabel('输出-目标+');
title('输出信号和目标信号');
%计算误差
e=T-a;
%显示误差e的曲线
plot(time,e)
hold on
%:r虚线
plot([min(time) max(time)],[0 0],':r')
hold off
xlabel('时间');
ylabel('误差');
title('误差信号')