x=[0,1,2,3,4,5,6,7,8,9,10];
53 Fy=[10,5,15,10,20,15,25,20,30,25,35];
54
55 if length(x)==length(Fy)
56 n=length(x);
57 stopien=5;
58 if stopien<n-1
59 A=polyfit(x,Fy,stopien);
60 Wy=zeros(1,n);
61
62 for i=1:n
63 for j=1:stopien+1
64 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
65 end
66 end
67 plot(x,Fy,'rx');
68 grid on;
69 hold on;
70 plot(x, Wy);
71
72 #*******************
73
74 disp("Sprawdzenie");
75 disp("Blad wzgledny");
76 bw=sum(abs(Fy-Wy))
77
78 disp("Blad kwadratowy");
79 bk=sum((Fy-Wy).^2)
80
81 disp("Sredni blad procentowy");
82 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
83
84 disp("Srednie odchylenie od krzywej");
85 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
86
87 opis=sprintf("Bw=%f, Bk=%f, Sbp=%f, Sook=%f",bw, bk, sbp, sook);
88 title(opis);
89
90 disp("********************");
91
92 else
93 disp("Za duzy stopien wielomianu aproksymujacego");
94 end
95 end