x = [0, 1, 3, 5, 7, 9, 11,2,4,6,8];
2 Fy = [-2, 3, 1, 6, 8, 9, 10,12,5,4,-1];
3
4 if length(x)==length(Fy)
5 n=length(x);
6 stopien=4;
7 if stopien<n-1
8 A=polyfit(x,Fy,stopien);
9 Wy=zeros(1,n);
10
11
12
13 for i=1:n
14 for j=1:stopien+1
15 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
16 end
17 end
18 plot(x,Fy,'rx');
19 grid on;
20 hold on;
21 plot(x,Wy);
22
23
24 disp("Sprawdzenie");
25 disp("Blad wzgledny");
26 bw=sum(abs(Fy-Wy))
27 disp("Blad kwadratowy");
28 bk=sum((Fy-Wy).^2)
29 disp("Sredni blad procentowy");
30 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
31 disp("Srednie odchylenie od krzywej");
32 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
33 opis=sprintf("Bw=%f, Bk=%f, Sbp=%f, Sook=%f",bw, bk, sbp, sook);
34 title(opis);
35 else
36 disp("Za duzy stopien wielomianu aproksymujacego");
37 end
38 end