Advertisement
Guest User

Untitled

a guest
May 22nd, 2018
101
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 0.94 KB | None | 0 0
  1. 101 x=[0,1,2,3,4,5,6,7,8,9,10];
  2. 102 Fy=[10,5,15,10,20,15,25,20,30,25,35];
  3. 103
  4. 104 if length(x)==length(Fy)
  5. 105 n=length(x);
  6. 106 stopien=6;
  7. 107 if stopien<n-1
  8. 108 A=polyfit(x,Fy,stopien);
  9. 109 Wy=zeros(1,n);
  10. 110
  11. 111 for i=1:n
  12. 112 for j=1:stopien+1
  13. 113 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
  14. 114 end
  15. 115 end
  16. 116 plot(x,Fy,'rx');
  17. 117 grid on;
  18. 118 hold on;
  19. 119 plot(x, Wy);
  20. 120
  21. 121 #*******************
  22. 122
  23. 123 disp("Sprawdzenie");
  24. 124 disp("Blad wzgledny");
  25. 125 bw=sum(abs(Fy-Wy))
  26. 126
  27. 127 disp("Blad kwadratowy");
  28. 128 bk=sum((Fy-Wy).^2)
  29. 129
  30. 130 disp("Sredni blad procentowy");
  31. 131 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
  32. 132
  33. 133 disp("Srednie odchylenie od krzywej");
  34. 134 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
  35. 135
  36. 136 opis=sprintf("Bw=%f, Bk=%f, Sbp=%f, Sook=%f",bw, bk, sbp, sook);
  37. 137 title(opis);
  38. 138
  39. 139 disp("********************");
  40. 140
  41. 141 else
  42. 142 disp("Za duzy stopien wielomianu aproksymujacego");
  43. 143 end
  44. 144 end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement