" Н а у к а м о л о д ы х " , 3 0 - 3 1 м а р т а 2 0 1 7 г . , А р з а м а с
П о с в я щ а е т с я 1 0 0 - л е т и ю Р о с т и с л а в а Е в г е н ь е в и ч а А л е к с е е в а
490
Net=newrb(P,T,GOAL,SPREAD),
где P – матрица входных данных; Т – вектор выходных данных, GOAL -
среднеквадратичная ошибка (в нашей модели принята равной 0,0001), SPREAD
– параметр влияния радиально-базисной функции (в нашей модели принят
равным 1). Образец программного кода (Matlab) имеет следующий вид:
exp = length(BD);
Br0n = zeros(1, exp);
e0n = zeros(1, exp);
for
e = 1:exp
B = cat(1,B,BD(e,:));
NRBA;
Br0n(e) = Br1(n);
e0n(e) = (Br1(n) - BD(e, 1))/BD(e, 1);
B(n,:) = [];
A(n,:) = [];
n=n-1;
C(:,n) = [];
T(:,n) = [];
P(:,n) = [];
G(:,n) = [];
Br1(:,n) = [];
es1(n,:) = [];
q1=q1-1;
bq1=bq1-1;
end
ep1 = e0n';
Br9 = Br0n';
[mean(abs(es1)) mean(abs(ep1))]
A = zscore(B);
n = length(B(:,1));
C = A';
T = C(1,:);
C(1,:) = [];
P = C;
net = newrb(P,T,0.00009,1);
net.layers{1}.size;
G = sim(net,P);
Br1 = G*std(B(:,1)) + mean(B(:,1));
q1 = length(Br1);
bq1 = length(B(:,1));
