NS L7

clear;
x = linspace(-10, 10, 50);
Q = length(x);
y = x.^3 - 8*x.^2;
r = 10;
t1=y+(rand(1,Q).*y*r/50-y*r/100);
plot(x,y,x,t1,'r-');
legend('первоначальная функция', 'с погрешностями', 2);
title('Исходные функции');

net1 = newff(minmax(x), [3 1], {'tansig' 'purelin'});
net1.trainParam.epochs;
net1 = train(net1, x, t1);
yf = sim(net, x);
figure;
plot(1:Q, yf,'k', 1:Q, t1,'r', 1:Q, y,'g');
legend('сглаженная функция', 'с погрешностями','первоначальная функция',2);
title('Сглаживание функции с помощью HC newff');