chengxu.m

采用matlab来产生神经逆模块

its3(1,:)=its1(2,:); ims3(1,:)=ims1(2,:);p3(1,:)=p1(2,:);ws3(1,:)=ws1(2,:);umr3(1,:)=umr1(2,:);utr3(1,:)=utr1(2,:);%dtdf (qudiaoshijianhangde) I/O(hangyiwei)
%stdf (qudiaoshijianhangde) I/O (hangyiwei)

for i=1:1:1997     
dits5(1,i)=(-25*its3(1,i)+48*its3(1,i+1)-368*its3(1,i+2)+16*its3(1,i+3)-3*its3(1,i+4))/0.12
dims5(1,i)=(-25*ims3(1,i)+48*ims3(1,i+1)-368*ims3(1,i+2)+16*ims3(1,i+3)-3*ims3(1,i+4))/0.12

end                                       %qiudaogongshi  shujugeshu-4 
for i=1:1:1997
its5(1,i)=its3(1,i);  ims5(1,i)=ims3(1,i);p5(1,i)=p3(1,i);
ws5(1,i)=ws3(1,i);umr5(1,i)=umr3(1,i);utr5(1,i)=utr3(1,i);

end                                                                    % I/O lietongwei
x=[ims5;dims5;its5;dits5;p5;ws5];


t=[umr5;utr5];                       % I/O hecheng
    m(1,:)=max(abs(x(1,:)));
    m(2,:)=max(abs(x(2,:)));
    m(3,:)=max(abs(x(3,:)));
    m(4,:)=max(abs(x(4,:)));
    m(5,:)=max(abs(x(5,:)));
    m(6,:)=max(abs(x(6,:)));
    m(1,:)
    m(2,:)
    m(3,:)
    m(4,:)
    m(5,:)
    m(6,:)
    P(1,:)=x(1,:)/m(1,:);;P(2,:)=x(2,:)/m(2,:); P(3,:)=x(3,:)/m(3,:); 
    P(4,:)=x(4,:)/m(4,:); P(5,:)=x(5,:)/m(5,:); P(6,:)=x(6,:)/m(6,:); 
    P                                                                          % IN guiyihua
    max(abs(t(1,:)));
    max(abs(t(2,:))); 
    n(1,:)=max(abs(t(1,:)));
    n(2,:)=max(abs(t(2,:)));
    n(1,:)
    n(2,:)
    T(1,:)=t(1,:)/n(1,:);  T(2,:)=t(2,:)/n(2,:) ; 
    T                        % OUT guiyihua   
net=newff(minmax(P),[13 2],{'tansig','purelin'},'trainlm');net=train(net,P,T);gensim(net,-1)         %shengcheng xunlian net