flch8eg5.m

刚刚接触神经网络的朋友们看了这本书会有一定得帮助的！！

clear
close all
%------------------------------------------------------------------------------------------------------得到样本
[t,y] = ode45(@LorenzDifEqn2,[0:0.01:100],[0.1 0.1 10]);%产生混沌序列
g = y(5001:6027,2); %取其中的1027个 
std_g=std(g);%标准差
SNR=20*log(0.38/std_g);%
%--------------------------------------------------------------------------------
g=g;
%figure;hold;
%plot(g);
s= [g-mean(g)]/[max(g)-min(g)]; % 混沌序列归一化到均值为 0，振幅为 1
%figure;hold;
%plot(s);
%------------------------------------------------------------------------------------------------------重构数据
%First construct the input samples,namely:RPS
num=length(s);%total data length
K=1;%time delay
m=3;%embedding dimension

n=num-(m-1)*K;%Number of phase space state vectors
Sample=zeros(n,m);
for i=1:n
    for j=1:m
        Sample(i,j)=s(i+(j-1)*K);%reconstruction with time delay
    end
end
%---------------------------------------------------------------------------------------------------------
%Then construct target set
target=Sample(:,end);
%train SVM

%train RBF
net = newrb(Sample(1:250,:)',target(2:251)',0.03);%SSE = 0.0308786
%--------------------------------------------------------------------------------------------------------
%president
for k=1:1024
    XT=Sample(k,:)';%
    rbf(k) = sim(net,XT);
end

nnt=1:1024;
figure;hold on;
plot(nnt,target(2:1025),'-',nnt,rbf(1:1024)','.');%图1无信号
%get MSE
y_true=target(2:1025);
y_pred=rbf(1:1024)';
error1=y_true-y_pred;
error(1)=(error1)'*error1;
error(1)=error(1)/1024;
error(1)=sqrt(error(1))