Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 101 x=[0,1,2,3,4,5,6,7,8,9,10];
- 102 Fy=[10,5,15,10,20,15,25,20,30,25,35];
- 103
- 104 if length(x)==length(Fy)
- 105 n=length(x);
- 106 stopien=6;
- 107 if stopien<n-1
- 108 A=polyfit(x,Fy,stopien);
- 109 Wy=zeros(1,n);
- 110
- 111 for i=1:n
- 112 for j=1:stopien+1
- 113 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
- 114 end
- 115 end
- 116 plot(x,Fy,'rx');
- 117 grid on;
- 118 hold on;
- 119 plot(x, Wy);
- 120
- 121 #*******************
- 122
- 123 disp("Sprawdzenie");
- 124 disp("Blad wzgledny");
- 125 bw=sum(abs(Fy-Wy))
- 126
- 127 disp("Blad kwadratowy");
- 128 bk=sum((Fy-Wy).^2)
- 129
- 130 disp("Sredni blad procentowy");
- 131 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
- 132
- 133 disp("Srednie odchylenie od krzywej");
- 134 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
- 135
- 136 opis=sprintf("Bw=%f, Bk=%f, Sbp=%f, Sook=%f",bw, bk, sbp, sook);
- 137 title(opis);
- 138
- 139 disp("********************");
- 140
- 141 else
- 142 disp("Za duzy stopien wielomianu aproksymujacego");
- 143 end
- 144 end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement