Нейрон с квадратичным сумматором, решение задачи аппроксимации функции

function Lab_1_11
clear
clc
X=[0 1 2 3 4 5 6 7 8 9];
Y=[0 1 4 9 16 25 36 49 64 81];
 
%X1=[0 0.4 2.1 2.8 3.7 4.6 6 7.3 8 9];
%Y1=[0 0.16 4.41 7.84 13.69 21.16 36 53.29 64 81];
 
X2=[0 0.2 1.7 3 3.5 4.7 6.4 7.2 7.8 9];
 
X=X/max(X);
Y=Y/max(Y);
 
X1=X1/max(X1);
Y1=Y1/max(Y1);
 
X2=X2/max(X2);
 
q(1:10, 1:10, 1:10)=0.5;
p(1:10, 1:10)=0.5;
r(1:10)=0.5;
y(1:10)=0;
 
alpha=0.1;
for n=1:1000
    for i=1:10
        U=quadSumm(X,q(:,:,i),p(:,i),r(i));
        y(i)=U;
        
        delta=Y(i)-y(i);
 
        for t=1:10
            for f=1:10
                q(t,f,i)=q(t,f,i)+2*alpha*delta*X(t)*X(f);
            end
            p(t,i)=p(t,i)+alpha*delta*X(t);
        end
        r(i)=r(i)+alpha*delta;
    end
end
for i=1:10
    U=quadSumm(X2,q(:,:,i),p(:,i),r(i));
    y(i)=U;
    display(y);
    %display(U);
end
%display(q);display(p);display(r);
%plot(X,Y,X1,Y1,X2,y);
plot(X2,y);
end
 
 
 
function V=quadSumm(X,q,p,r)
V=0;
for i=1:10
    for j=1:10
 
            V=V+q(i,j)*X(i)*X(j);
 
    end
    V=V+p(i)*X(i);
end
V=V+r;
end